MINICERN/mathlib/gen/f/bakvec.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:32  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE BAKVEC(NM,N,T,E,M,Z,IERR)
  11       INTEGER I,J,M,N,NM,IERR
  12       REAL T(NM,3),E(N),Z(NM,M)
  13       IERR = 0
  14       E(1) = 1.0
  15       IF (N .EQ. 1) GO TO 1001
  16       DO 100 I = 2, N
  17          IF (E(I) .NE. 0.0) GO TO 80
  18          IF (T(I,1) .NE. 0.0 .OR. T(I-1,3) .NE. 0.0) GO TO 1000
  19          E(I) = 1.0
  20          GO TO 100
  21    80    E(I) = E(I-1) * E(I) / T(I-1,3)
  22   100 CONTINUE
  23       DO 120 J = 1, M
  24          DO 120 I = 2, N
  25          Z(I,J) = Z(I,J) * E(I)
  26   120 CONTINUE
  27       GO TO 1001
  28  1000 IERR = 2 * N + I
  29  1001 RETURN
  30       END