Esope source files

   1 : C PROMA1    SOURCE    CHAT      05/01/13    02:32:58     5004           
   2 : c---------------------------------------------------------------------
   3 : c
   4 :       SUBROUTINE PROMA1 (A, B, C, N1, N2, N3)
   5 : c
   6 : c=====================================================================
   7 : c                                                                    =
   8 : c  This routine evaluates a matrix product.                          =
   9 : c                                                                    =
  10 : c           a(i,j)=b(i,k)*c(k,j)    a=b*c                            =
  11 : c                                                                    =
  12 : c=====================================================================
  13 :       IMPLICIT INTEGER(I-N)
  14 :       integer*4 n1, n2, n3
  15 :       real*8    a (n1,n2), b (n1,n3), c (n3,n2)
  16 : c
  17 :       integer*4 i, j, k
  18 : c
  19 :       do i = 1, n1
  20 :          do j = 1, n2
  21 :             a (i,j) = 0.0 d0
  22 :             do k = 1, n3
  23 :                a (i,j) = a (i,j) + b (i,k) * c (k,j)
  24 :             end do
  25 :           end do
  26 :       end do
  27 : c
  28 :       return
  29 :       end
  30 :  
  31 :