#include <stdio.h>
#include <math.h>
#define ND  1000
extern double A[ND][ND], B[ND][ND], C[ND][ND] ;
//=================================================================C
void MULT1(int N, int M, int L)
//=================================================================C
//  Real-Dense Matrix Multiplication                               C
//   C = A*B  with JKI-Type                                        C
//-----------------------------------------------------------------C
//    A,B,C     Global define Matrix                               C
//    N         I*4, In,  Matrix Size of A(Column), C(Column)      C
//    M         I*4, In,  Matrix Size of B(Row), C(Row)            C
//    L         I*4, In,  Matrix Size of A(Row), B(Column)         C
//-----------------------------------------------------------------C
//    Written by Yasunori Ushiro ,  2007/05/30                     C
//        ( Tokyo Polytechnic University )                         C
//     後 保範（東京工芸大学）                                     C
//=================================================================C
{ int     i, j, k ;
  double  BW ;
//   Clear C
  for (j=0; j<M; j++) {
    for (i=0; i<N; i++) {
      C[i][j] = 0.0 ; }
   }
//   C=C+A*B
  for (j=0; j<M; j++) {
    for (k=0; k<L; k++) {
      BW = B[k][j] ;
      for (i=0; i<N; i++) {
        C[i][j] = C[i][j] + A[i][k]*BW ; }
   }  }
 }