]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 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 |