int MUL(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T) /*乘法*/
{ int arow,ccol,tp,p,q,blow,t,i; int ctemp[MAXRC]; if(A->nu!=B->mu) return 0; T->mu=A->mu;T->nu=B->nu;T->tu=0; if(A->tu*B->tu!=0) { for(arow=1;arow<=A->mu;++arow) { for(i=1;i<=T->nu;i++) ctemp[i]=0; T->rpos[arow]=T->tu+1; if(arow<A->mu) tp=A->rpos[arow+1]; else{tp=A->tu+1;} for(p=A->rpos[arow];p<tp;++p) { blow=A->data[p].j; if(blow<B->mu) t=B->rpos[blow+1]; else t=B->tu+1; for(q=B->rpos[blow];q<t;++q) { ccol=B->data[q].j; ctemp[ccol] += A->data[p].e*B->data[q].e; } /*for q*/ } /*求得Q中的crow( arow)行的非零元*/ for(ccol=1;ccol<T->nu;++ccol) { if(ctemp[ccol]) { if(++T->tu>MAXSIZE) return 0; T->data[T->tu].i=arow; T->data[T->tu].j=ccol; T->data[T->tu].e=ctemp[ccol]; } } } Print(A,B,T); return 1; } } void Printt(RLSMatrix *A) /*打印矩阵*/
{ int row,col,i=1,k=0; for(row=1;row<=A->mu;row++) { for(col=1;col<=A->nu;col++) { if(A->data[i].i==row&&A->data[i].j==col&&i<=A->tu) printf("%-4d",A->data[i++].e); else printf("%-4d",k); } printf("\n"); } } void Print(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T) /*打印函数*/
{ clrscr(); if(A->mu!=0) {printf("Print Matrix A:\n");Printt(A); printf("\n");} else printf("Matrix A is blank!\n"); if(B->tu!=0) {printf("Print Matrix B:\n");Printt(B); printf("\n");} else printf("Matrix B is blank!\n"); if(T->tu!=0) {printf("Print Matrix T:\n");Printt(T); printf("\n");} else {printf("Matrix T is blank!\n");getch();} getch(); } void Enter(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T) /*输入函数*/
{ int i,k,num[MAXRC]; char cmd; clrscr(); printf("Enter the mu and nu and tu of Matrix A(mu,nu,tu): \n"); /*输入矩阵A*/ scanf("%d,%d,%d",&A->mu,&A->nu,&A->tu); while(A->mu<=0||A->mu>MAXR||A->nu<=0||A->nu>MAXC) { printf("Enter the mu and nu and tu of Matrix A(mu,nu,tu): \n"); scanf("%d,%d,%d",&A->mu,&A->nu,&A->tu); } printf("Enter the Matrix A(i-j:e):\n"); for(i=1;i<=A->tu;i++) scanf("%d-%d:%d",&A->data[i].i,&A->data[i].j,&A->data[i].e); for(i=1;i<=A->mu;i++) num[i]=0; for(i=1;i<=A->tu;i++) ++num[A->data[i].i]; A->rpos[1]=1; for(i=2;i<=A->mu;i++) A->rpos[i]=A->rpos[i-1]+num[i-1]; printf("Enter the mu and nu and tu of Matrix B(mu,nu,tu): \n"); /*输入矩阵B*/
scanf("%d,%d,%d",&B->mu,&B->nu,&B->tu); while(A->mu<=0||A->mu>MAXR||A->nu<=0||A->nu>MAXC) { printf("Enter the mu and nu and tu of Matrix B(mu,nu,tu): \n"); scanf("%d,%d,%d",&B->mu,&B->nu,&B->tu); } printf("Enter the Matrix B(i-j:e):\n"); for(i=1;i<=B->tu;i++) scanf("%d-%d:%d",&B->data[i].i,&B->data[i].j,&B->data[i].e); for(i=1;i<=B->mu;i++) num[i]=0; for(i=1;i<=B->tu;i++) ++num[B->data[i].i]; B->rpos[1]=1; for(i=2;i<=B->mu;i++) B->rpos[i]=B->rpos[i-1]+num[i-1]; do{ printf("Enter the operation you want:"); /*选择需要的运算*/ cmd=getche(); }while(cmd!='+'&&cmd!='-'&&cmd!='*'); getch(); switch(cmd) { case '+' : if(A->mu==B->mu&&A->nu==B->nu) ADD(A,B,T); break; case '-' : if(A->mu==B->mu&&A->nu==B->nu) SUB(A,B,T); break; case '*' : if(A->nu==B->mu) MUL(A,B,T); break; default : printf("No thia operation!"); break; } } void InitializationMatrix(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T) /*初始化矩阵*/ { A->mu=A->nu=A->tu=0; B->mu=B->nu=B->tu=0; T->mu=T->nu=T->tu=0; } /*************************************主函数************************************/ void InitializationMatrix(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T); void Enter(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T); void Print(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T); void ADD(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T); void SUB(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T); int MUL(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T); void Initialization(); void Myinformation(); void Exit(); void main() { int i,k,num[MAXRC];
char cmd; RLSMatrix *A,*B,*T; A=( RLSMatrix *)malloc(sizeof( RLSMatrix)); /*分配空间*/ B=( RLSMatrix *)malloc(sizeof( RLSMatrix)); T=( RLSMatrix *)malloc(sizeof( RLSMatrix)); Myinformation();
InitializationMatrix(A,B,T); while(1) { Initialization(); do{ gotoxy(10,10); printf("Enter E/e ,Q/q ,P/p to continue:"); cmd=getche(); printf("\n"); }while(cmd!='q'&&cmd!='Q'&&cmd!='e'&&cmd!='E'&&cmd!='p'&&cmd!='P'); switch(cmd) { case 'q' : Exit(); break; case 'Q' : Exit(); break; case 'e' : Enter(A,B,T); break; case 'E' : Enter(A,B,T); break; case 'p' : Print(A,B,T); break; case 'P' : Print(A,B,T); break; } } } void Initialization() /*初始化函数*/
{ clrscr(); printf("*********************************************************************************"); gotoxy(7,2); printf("Enter the RLSMatrix--E"); gotoxy(38,2); printf("Print--P"); gotoxy(55,2); printf("Quit--Q"); gotoxy(80,2); printf("*"); printf("********************************************************************************"); gotoxy(1,22); printf("*********************************************************************************"); gotoxy(10,23); printf("Enter a operation code: E/e ,Q/q ,P/p to continue: "); gotoxy(80,23); printf("*"); printf("********************************************************************************"); } void Myinformation() /*我的信息函数*/ { int i; clrscr(); for(i=1;i<=240;i++) printf("\1"); gotoxy(15,8); printf("This is a Matrix arithmetic calculator!"); gotoxy(15,12); printf("Name: liyongjun."); gotoxy(15,13); printf("Num: 3104006893."); gotoxy(15,14); printf("Grade: 2004."); gotoxy(15,15); printf("Class: cumputer science and technology 10"); gotoxy(1,22); for(i=1;i<=240;i++) printf("\1"); gotoxy(1,18); printf("Print any key to continue..."); getch(); } void Exit() /*退出函数*/
{ char cmd; do{ gotoxy(10,11); printf("Do you want to out?Y/N:"); cmd=getche(); }while(cmd!='y'&&cmd!='Y'&&cmd!='n'&&cmd!='N'); if(cmd=='y'||cmd=='Y') exit(0); }
|
Copyright © 2007-2012 www.chuibin.com 六维论文网 版权所有