SUBROUTINE GMPRD ( A, B, R, N, M, L )

Argument Definitions (+ indicates altered content)
DOUBLE PRECISION   A(1),       B(1),      +R(1)
INTEGER            L,          M,          N
Description
.................................................................. 
 
   subroutine gmprd 
 
   purpose 
      multiply two general matrices to form a resultant general 
      matrix 
 
   usage 
      call gmprd(a,b,r,n,m,l) 
 
   description of parameters 
      a - name of first input matrix 
      b - name of second input matrix 
      r - name of output matrix 
      n - number of rows in a 
      m - number of columns in a and rows in b 
      l - number of columns in b 
 
   remarks 
      all matrices must be stored as general matrices 
      matrix r cannot be in the same location as matrix a 
      matrix r cannot be in the same location as matrix b 
      number of columns of matrix a must be equal to number of row 
      of matrix b 
 
   subroutines and function subprograms required 
      none 
 
   method 
      the m by l matrix b is premultiplied by the n by m matrix a 
      and the result is stored in the n by l matrix r. 
 
..................................................................
Source file:gmprd.f
Local Variables (+ indicates altered content)
INTEGER           +I,         +IB,        +IK,        +IR
INTEGER           +J,         +JI,        +K