matrixbk_test.cc

Go to the documentation of this file.

/*
 * Copyright (c) 1999-2003 Bert Kampes
 * Copyright (c) 1999-2003 Delft University of Technology, The Netherlands
 *
 * This file is part of Doris, the Delft o-o radar interferometric software.
 *
 * Doris program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * Doris is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 * Publications that contain results produced by the Doris software should
 * contain an acknowledgment. (For example: The interferometric processing
 * was performed using the freely available Doris software package developed
 * by the Delft Institute for Earth-Oriented Space Research (DEOS), Delft
 * University of Technology, or include a reference to: Bert Kampes and
 * Stefania Usai. \"Doris: The Delft Object-oriented Radar Interferometric
 * software.\" In: proceedings 2nd ITC ORS symposium, August 1999. (cdrom)).
 *
 */
/****************************************************************
 * $Source: /users/kampes/DEVELOP/DORIS/doris/src/RCS/matrixbk_test.cc,v $      *
 * $Revision: 3.8 $                                             *
 * $Date: 2005/03/07 10:38:40 $                                 *
 * $Author: kampes $                                            *
 *                                                              *
 * Test program for matrix class.                               *
 * Compile with Doris Makefile.                                 *
 #%// BK 17-Aug-2000                                            *
 ****************************************************************/


#include "constants.hh"                 // typedefs
#include "matrixbk.hh"                  // matrix functions
bk_messages TRACE;
bk_messages DEBUG;
bk_messages INFO;
bk_messages PROGRESS;
bk_messages WARNING;
bk_messages ERROR;
bk_messages matDEBUG;
bk_messages matERROR;
// ______ used in matrix_bk.cc to keep track of total allocated ______
uint totalallocated=0;  // [B] matrixdebuging





/****************************************************************
 *    main                                                      *
 #%// BK 17-Aug-2000                                            *
 ****************************************************************/
int32 main(
        int argc,
        char* argv[])
  {
  const int L = 4;      // numlines A
  const int P = 2;      // numcols  A
  int i,j;

  matrix<complr4> A(L,P);
  matrix<complr4> B(P,L+1);
  for (i=0; i<A.lines(); ++i)
    for (j=0; j<A.pixels(); ++j)
      A(i,j) = complr4(float(i)/2.3+j,float(j)/2+i);
  for (i=0; i<B.lines(); ++i)
    for (j=0; j<B.pixels(); ++j)
      B(i,j) = complr4(float(j+i)/2.3+i,float(i)/(2.2+i-j*j));

  // ______ Show matrices for testing ______
  cout << "matrix A(" << A.lines() << "," << A.pixels() << "):\n";
  A.showdata();
  cout << "matrix B(" << B.lines() << "," << B.pixels() << "):\n";
  B.showdata();
  cout << "\n\n---------------------\n\n";

  // ______ Test - ______
  B = complr4(2.) * -B;
  cout << "matrix 2*-B(" << B.lines() << "," << B.pixels() << "):\n";
  B.showdata();
  cout << "\n\n---------------------\n\n";

  matrix<complr4> Bconj = conj(B);
  cout << "conjugated:\n";
  Bconj.showdata();
  cout << "\n\n---------------------\n\n";

  // ______ Test * ______
  matrix<complr4> C = A*B;
  cout << "matrix C(" << C.lines() << "," << C.pixels() << ") = A*B\n";
  C.showdata();
  cout << "\n\n---------------------\n\n";

  // ______ Test 1d fft ______
  fft(A,1);
  cout << "1d fft over columns of matrix A(" << A.lines() << "," << A.pixels() << "):\n";
  A.showdata();
  cout << "\n\n---------------------\n\n";

//    // ______ Test operator = ______
//    for (i=0; i<10; ++i)
//      {
//      cout << i << ": A=B\n";
//      A=B;
//      cout << i << ": B=B\n";
//      B=B;
//      cout << i << ": C=A\n";
//      matrix<complr4> C=A;
//      cout << "\n";
//      }

  // ______ Test find function ______
  // matrix<int32> indexNaN = A.find(2);

  // ______ Test mypow function ______
  matrix<float> R(L,P);
  for (i=0; i<R.lines(); ++i)
    for (j=0; j<R.pixels(); ++j)
      R(i,j) = float(i/2.3+j*i*j);
  
  cout << "\n\n-RRRRRRRRRRRRRRRRRRR-\n\n";
  R.showdata();
  R.mypow(1.5);
  R.showdata();
// ?? but does B*= A^C work?
// ?? but does B*= A.mypow(C) work?

  // ______ Test CR4*R4 function ______
  cout << "\n\n-RRRRRRRRRRRRRRRRRRR-\n\n";
  matrix<complr4> Q(L,P);
  for (i=0; i<Q.lines(); ++i)
    for (j=0; j<Q.pixels(); ++j)
      Q(i,j) = complr4(float(i)/2.3+j,float(j)/2+i);
  Q.showdata();
  R.showdata();
  Q*=R;
  Q.showdata();

  // ______ Test smooth function ______
  // R.smooth(1);
  // ______ Test pow ^ operator ______
  // R = R^2.0;

  // ______ Test wshift function ______
  matrix<real4> AA(1,7);
  for (i=0; i<AA.lines(); ++i)
    for (j=0; j<AA.pixels(); ++j)
      AA(i,j) = i+j;
  AA.showdata();
  cout << "wshift AA, -2\n";
  wshift(AA,-2);
  AA.showdata();

  // ______ Test diagxmat R4*CR4 ______
  matrix<complr4> QQ(7,7);
  for (i=0; i<QQ.lines(); ++i)
    for (j=0; j<QQ.pixels(); ++j)
      QQ(i,j) = complr4(float(i)/2.3+j,float(j)/2+i);
  QQ.showdata();
  matrix<complr4> BB = diagxmat(AA,QQ);
  cout << "BB=diagxmat(AA,QQ)\n";
  BB.showdata();

  // ______ Test cos/sin complr4 operator ______
  cout << "VVVVVVVVVVVVVVVVVVVVVVVVVVVVV\n";
  matrix<real8> r8A(3,5);
  r8A = 5.3;
  cout << "r8A=5.3:\n";
  r8A.showdata();

  for (i=0; i<r8A.lines(); ++i)
    for (j=0; j<r8A.pixels(); ++j)
      r8A(i,j) = (3.14/180.0)*((0.5+i-j)*100);  // rad
  cout << "r8A=radians:\n";
  r8A.showdata();
  matrix<real8> cosr8A = cos(r8A);
  cout << "\ncosr8A:\n";
  cosr8A.showdata();
  //
  matrix<real8> r8B = r8A+0.2;
  matrix<real8> sinr8B = sin(r8B);
  cout << "\nsinr8B:\n";
  sinr8B.showdata();
  //
  matrix<complr4> cr4AB = mat2cr4(r8A,r8B);
  cout << "\ncomplr4(r8A,r8B):\n";
  cr4AB.showdata();



  cout << "\n\nNormal termination.\n\n";
  return 0;
  } // END main

---