MINICERN/mathlib/gen/f/eltran.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:34  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE ELTRAN(NM,N,LOW,IGH,A,INT,Z)
  11       INTEGER I,J,N,KL,MM,MP,NM,IGH,LOW,MP1
  12       REAL A(NM,IGH),Z(NM,N)
  13       INTEGER INT(IGH)
  14       DO 80 I = 1, N
  15          DO 60 J = 1, N
  16    60    Z(I,J) = 0.0
  17          Z(I,I) = 1.0
  18    80 CONTINUE
  19       KL = IGH - LOW - 1
  20       IF (KL .LT. 1) GO TO 200
  21       DO 140 MM = 1, KL
  22          MP = IGH - MM
  23          MP1 = MP + 1
  24          DO 100 I = MP1, IGH
  25   100    Z(I,MP) = A(I,MP-1)
  26          I = INT(MP)
  27          IF (I .EQ. MP) GO TO 140
  28          DO 130 J = MP, IGH
  29             Z(MP,J) = Z(I,J)
  30             Z(I,J) = 0.0
  31   130    CONTINUE
  32          Z(I,MP) = 1.0
  33   140 CONTINUE
  34   200 RETURN
  35       END