Sources Esope | Cast3M

proma1
C PROMA1    SOURCE    CB215821  17/07/21    21:15:27     9513           
c---------------------------------------------------------------------
c
      SUBROUTINE PROMA1 (A, B, C, N1, N2, N3)
c
c=====================================================================
c                                                                    =
c  This routine evaluates a matrix product.                          =
c                                                                    =
c           a(i,j)=b(i,k)*c(k,j)    a=b*c                            =
c                                                                    =
c=====================================================================
      IMPLICIT INTEGER(I-N)
      integer   n1, n2, n3
      real*8    a (n1,n2), b (n1,n3), c (n3,n2)
c
      integer   i, j, k
c
      do i = 1, n1
         do j = 1, n2
            a (i,j) = 0.d0
            do k = 1, n3
               a (i,j) = a (i,j) + b (i,k) * c (k,j)
            enddo
          enddo
      enddo
c
      return
      end