17 using namespace NEWMAT;
52 Tracer et(
"Fifth test of Matrix package");
58 for (i=1;i<=5;i++)
for (j=1;j<=6;j++) A(i,j)=1+i*j+i*i+j*j;
60 for (i=1;i<=6;i++) CV(i)=i*i+3;
68 CV=(RV*A).t()+A.
t()*CV2; CV1=(A.
t()*A)*CV1 - CV;
75 for (i=1;i<=6;i++) { CV1(i)=i*3+1; CV2(i)=10-i; CV(i)=11+i*2; }
84 for (i=1;i<=10;i++) RV(i) = i*i;
86 RV1 << 1 << 4 << 9 << 16 << 25 << 36 << 49 << 64 << 81 << 100;
89 et.
ReName(
"Fifth test of Matrix package - almost at end");
97 X(1,1) -= 11; X(1,2) -= 12; X(1,3) -= 13;
98 X(2,1) -= 21; X(2,2) -= 22; X(2,3) -= 23;
101 Y(1,1) -= 11; Y(2,1) -= 12; Y(3,1) -= 13;
102 Y(1,2) -= 21; Y(2,2) -= 22; Y(3,2) -= 23;
105 et.
ReName(
"Fifth test of Matrix package - at end");
114 for (i=1; i<=31; i++)
for (j=i; j<=31; j++) UT(i,j) = i+j+(i-j)*(i-2*j);
ReturnMatrix Returner6(const GenericMatrix &GM)
ReturnMatrix Returner5(const GenericMatrix &GM)
void Inject(const GeneralMatrix &GM)
ReturnMatrix Returner1(const GenericMatrix &GM)
ReturnMatrix Returner2(const GenericMatrix &GM)
Upper triangular band matrix.
TransposedMatrix t() const
ReturnMatrix Returner0(const GenericMatrix &GM)
ReturnMatrix Returner3(const GenericMatrix &GM)
ReturnMatrix Returner4(const GenericMatrix &GM)
The usual rectangular matrix.
Rectangular matrix for use with Numerical Recipes in C.
Lower triangular band matrix.
void Print(const Matrix &X)
A matrix which can be of any GeneralMatrix type.
void ReName(const char *)
ReturnMatrix Returner7(const GenericMatrix &GM)