6strassen矩阵乘法.docVIP

6strassen矩阵乘法 procedure STRASSEN(n,A,B,C); begin if n=2 then MATRIX-MULTIPLY(A，B，C) else begin 将矩阵A和B依(1)式分块; STRASSEN(n/2,A11,B12-B22,M1); STRASSEN(n/2,A11+A12,B22,M2); STRASSEN(n/2,A21+A22,B11,M3); STRASSEN(n/2,A22,B21-B11,M4); STRASSEN(n/2,A11+A22,B11+B22,M5); STRASSEN(n/2,A12-A22,B21+B22,M6); STRASSEN(n/2,A11-A21,B11+B12,M7); end; end; #includestdio.h #define X 32 #define Y 32 int judgement(int n) /*判断用户输入的矩阵的阶数是否为2的N次方*/ { int flag,temp=n; while(temp!=1temp%2==0) { if(temp%2==0) temp/=2; else flag=1; } if(temp==1) flag=0; return flag; } void divide(int *d,int *d11,int *d12,int *d21,int *d22,int n) /*分割矩阵*/ { int i,j; for(i=0;i(2*n);i++) for(j=0;j(2*n);j++) { if(in jn) { *(d11+i*n+j)=*(d+i*2*n+j); } else if(in (j=n) j2*n) { *(d12+i*n+j-n)=*(d+i*2*n+j); } else if(i=n i2*n jn) { *(d21+(i-n)*n+j)=*(d+i*2*n+j); } else if((i=n) i2*n (j=n) j2*n) { *(d22+(i-n)*n+j-n)=*(d+i*2*n+j); } } } int *merge(int *d11,int *d12,int *d21,int *d22,int n) /*合并矩阵*/ { int i,j; int *d; d=(int *)malloc(sizeof(int)*(4*n*n)); for(i=0;i2*n;i++) for(j=0;j2*n;j++) { if(in jn) { *(d+i*2*n+j)=*(d11+i*n+j); } else if(in (j=n) j2*n) { *(d+i*2*n+j)=*(d12+i*n+j-n); } else if(i=n i2*n jn) { *(d+i*2*n+j)=*(d21+(i-n)*n+j); } else if((i=n) i2*n (j=n) j2*n) { *(d+i*2*n+j)=*(d22+(i-n)*n+j-n); } } return d; } int *minus(int *a,int *b,int n) /*矩阵相减*/ { int i,j; int *c; c=(int *)malloc(sizeof(int)* (n*n)); for(i=0;i2*n;i++) for(j=0;j2*n;j++) { *(c+i*n+j)=*(a+i*n+j)-*(b+i*n+j); } return c; } int *plus(int *a,int *b,in