share/colnew/fortran/dmzsol.f

   1       SUBROUTINE DMZSOL (KD, MSTAR, N, V, Z, DMZ)
   2 C
   3 C**********************************************************************
   4 C
   5 C   purpose
   6 C          compute dmz in a blockwise manner
   7 C          dmz(i) = dmz(i)  +  v(i) * z(i), i = 1,...,n
   8 C
   9 C**********************************************************************
  10 C
  11       IMPLICIT REAL*8 (A-H,O-Z)
  12       DIMENSION V(KD,1), DMZ(KD,1), Z(1)
  13 C
  14       JZ = 1
  15       DO 30 I = 1, N
  16          DO 20 J = 1, MSTAR
  17             FACT = Z(JZ)
  18             DO 10 L = 1, KD
  19                DMZ(L,I) = DMZ(L,I)  +  FACT * V(L,JZ)
  20    10       CONTINUE
  21             JZ = JZ + 1
  22    20    CONTINUE
  23    30 CONTINUE
  24       RETURN
  25       END