cpp1.cpp

一个可进行矩阵计算的模版类,可实现大多数常用的矩阵计算功能

#include <time.h> 
#include"matrix.h"


#ifndef _NO_NAMESPACE
using namespace std;
using namespace math;
#define STD std
#else
#define STD
#endif

#ifndef _NO_TEMPLATE
typedef matrix<double> Matrix;
#else
typedef matrix Matrix;
#endif

#ifndef _NO_EXCEPTION
#  define TRYBEGIN()	try {
#  define CATCHERROR()	} catch (const STD::exception& e) { \
						cerr << "Error: " << e.what() << endl; }
#else
#  define TRYBEGIN()
#  define CATCHERROR()
#endif

int main()
{
	Matrix ma1;
	Matrix ma2;

	ma1.SetSize(4,4);
	ma2.SetSize(4,4);

	ma1(0,0)=0;
	ma1(0,1)=-1;
	ma1(0,2)=0;
	ma1(0,3)=10;
	ma1(1,0)=1;
	ma1(1,1)=0;
	ma1(1,2)=0;
	ma1(1,3)=0;
	ma1(2,0)=0;
	ma1(2,1)=0;
	ma1(2,2)=1;
	ma1(2,3)=0;
	ma1(3,0)=0;
	ma1(3,1)=0;
	ma1(3,2)=0;
	ma1(3,3)=1;
	ma2(0,0)=1;
	ma2(0,1)=0;
	ma2(0,2)=0;
	ma2(0,3)=20;
	ma2(1,0)=0;
	ma2(1,1)=0;
	ma2(1,2)=-1;
	ma2(1,3)=10;
	ma2(2,0)=0;
	ma2(2,1)=1;
	ma2(2,2)=0;
	ma2(2,3)=0;
	ma2(3,0)=0;
	ma2(3,1)=0;
	ma2(3,2)=0;
	ma2(3,3)=1;

	cout<<"matrix1=\n"<<ma1<<endl;
	cout<<"matrix2=\n"<<ma2<<endl;

	Matrix ma3=ma1*ma2;

	cout<<"matrix3=\n"<<ma3<<endl;

	return 0;
}