矩阵乘法运算test

//学习CG，刚好想到自己动手实现下矩阵变换吧
//just test and test

#include "gl/glut.h"
#include <iostream>
#include <cmath>
using namespace std;

#define MAX 3
class Matrix
{
public:
int row;
int col;
double num[MAX][MAX];
Matrix()
{
for(int i=0;i<MAX;i++)
for(int j=0;j<MAX;j++)
num[i][j]=0;
row=MAX;
col=MAX;
}
};
// A * B
Matrix MultiMatrix(Matrix &A, Matrix &B)
{
Matrix result;
result.row = A.row;
result.col = B.col;
for(int i=0;i<result.row;i++)
for(int j=0;j<result.col;j++)
for(int k=0;k<A.col;k++)
result.num[i][j] += A.num[i][k] * B.num[k][j];

return result;
}
//缩放变换
Matrix scale(Matrix &A,double x,double y)
{
Matrix s;
s.num[0][0]=x;
s.num[1][1]=y;
s.num[2][2]=1;
return MultiMatrix(s,A);
}
//逆时针旋转变换
Matrix rotate(Matrix &A,double angle)
{
Matrix r;
r.num[0][0]= cos(angle);
r.num[0][1]= -sin(angle);
r.num[1][0]= sin(angle);
r.num[1][1]= cos(angle);
r.num[2][2]= 1;
return MultiMatrix(r,A);
}
//平移变换
Matrix translation(Matrix &A,double x,double y)
{
Matrix t;
t.num[0][0]=1;
t.num[1][1]=1;
t.num[2][0]=-x;
t.num[2][1]=-y;
t.num[2][2]=1;
return MultiMatrix(t,A);
}

int main(int argc,char *argv[])
{
//glutInit(&argc,argv);
//glutInitDisplayMode(GLUT_RGB|GLUT_SINGLE);

//TEST FOR MATRIX
// A:3*2    b:2*3
Matrix A,B;
A.row=3;
A.col=2;
B.row=2;
B.col=3;
for(int i=0;i<3;i++)
for(int j=0;j<2;j++)
cin>>A.num[i][j];
for(int i=0;i<2;i++)
for(int j=0;j<3;j++)
cin>>B.num[i][j];
Matrix C = MultiMatrix(A,B);
for(int i=0;i<C.row;i++)
{
for(int j=0;j<C.col;j++)
cout<<C.num[i][j]<<" ";
cout<<endl;
}
char a;
cin>>a;
return 0;
}