www.pudn.com > matrix-.rar > BRMUL.C

 
  void brmul(a,b,m,n,k,c) 
  int m,n,k; 
  double a[],b[],c[]; 
  { int i,j,l,u; 
    for (i=0; i<=m-1; i++) 
    for (j=0; j<=k-1; j++) 
      { u=i*k+j; c[u]=0.0; 
        for (l=0; l<=n-1; l++) 
          c[u]=c[u]+a[i*n+l]*b[l*k+j]; 
      } 
    return; 
  }