home *** CD-ROM | disk | FTP | other *** search
/ Usenet 1994 October / usenetsourcesnewsgroupsinfomagicoctober1994disk2.iso / misc / volume34 / newmat06 / part07 / tmt5.cxx < prev    next >
C/C++ Source or Header  |  1992-12-06  |  1KB  |  51 lines
  1.  
  2. //#define WANT_STREAM
  3.  
  4.  
  5. #include "include.h"
  6.  
  7. #include "newmat.h"
  8.  
  9.  
  10. /**************************** test program ******************************/
  11.  
  12. void Print(const Matrix& X);
  13. void Print(const UpperTriangularMatrix& X);
  14. void Print(const DiagonalMatrix& X);
  15. void Print(const SymmetricMatrix& X);
  16. void Print(const LowerTriangularMatrix& X);
  17.  
  18. void trymat5()
  19. {
  20. //   cout << "\nFifth test of Matrix package\n";
  21.    Tracer et("Fifth test of Matrix package");
  22.    Exception::PrintTrace(TRUE);
  23.  
  24.    int i,j;
  25.  
  26.    Matrix A(5,6);
  27.    for (i=1;i<=5;i++) for (j=1;j<=6;j++) A(i,j)=1+i*j+i*i+j*j;
  28.    ColumnVector CV(6);
  29.    for (i=1;i<=6;i++) CV(i)=i*i+3;
  30.    ColumnVector CV2(5); for (i=1;i<=5;i++) CV2(i)=1.0;
  31.    ColumnVector CV1=CV;
  32.  
  33.    { 
  34.       CV=A*CV;
  35.       RowVector RV=CV.t(); // RowVector RV; RV=CV.t();
  36.       RV=RV-1.0;
  37.       CV=(RV*A).t()+A.t()*CV2; CV1=(A.t()*A)*CV1 - CV;
  38.       Print(CV1);
  39.    }
  40.  
  41.    CV1.ReDimension(6);
  42.    CV2.ReDimension(6);
  43.    CV.ReDimension(6);
  44.    for (i=1;i<=6;i++) { CV1(i)=i*3+1; CV2(i)=10-i; CV(i)=11+i*2; }
  45.    ColumnVector CX=CV2-CV; { CX=CX+CV1; Print(CX); }  
  46.    Print(ColumnVector(CV1+CV2-CV));
  47.    RowVector RV=CV.t(); RowVector RV1=CV1.t();   
  48.    RowVector R=RV-RV1; Print(RowVector(R-CV2.t()));
  49. //   cout << "\nEnd of fifth test\n";
  50. }
  51.