00001 /* */ 00002 /* Matrix Multiplication with an entory point */ 00003 /* */ 00004 /* Ver.1.0, Jun,30,1988. */ 00005 /* */ 00006 00007 #include "arith.h" 00008 00009 mat_mult2(n,m,l,amat,ai,aj,bmat,bi,bj,pmat) 00010 int n,m,l,ai,aj,bi,bj; 00011 MATRIX amat,bmat,pmat; 00012 { 00013 int i,j,k; 00014 00015 for(i=0; i<n; i++){ 00016 for(k=0; k<l; k++){ 00017 pmat[i][k]=0; 00018 for(j=0; j<m; j++){ 00019 pmat[i][k] += amat[ai+i][aj+j] * bmat[bi+j][bj+k]; 00020 } 00021 } 00022 } 00023 }