00001 /* */ 00002 /* Matrix Multiplication */ 00003 /* */ 00004 /* Ver.1.0, Jan.18,1988. */ 00005 /* */ 00006 00007 #include "arith.h" 00008 00009 mat_mult(n,m,l,amat,bmat,pmat) 00010 int n,m,l; 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.0e0; 00018 for(j=0; j<m; j++){ 00019 pmat[i][k] += amat[i][j] * bmat[j][k]; 00020 } 00021 } 00022 } 00023 }