]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.1 1996/04/01 15:02:18 mclareni | |
6 | * Mathlib gen | |
7 | * | |
8 | * | |
9 | #include "gen/pilot.h" | |
10 | SUBROUTINE EPDITR(BETA,DMAX,UCORR,ICORR,GNORM) | |
11 | C | |
12 | C THE EXTRAPOLATED-LIEBMANN ITERATION | |
13 | C | |
14 | DIMENSION DUM1(7025),UMAT(3250),DUM2(20),NUMPT(2900),DIFCO(2900,5) | |
15 | COMMON DUM1 , NC , NR , KODBAS , DX , DY | |
16 | COMMON JOBNUM , UMAT , NEQU , DUM2 , NUMPT , DIFCO | |
17 | DMAX=0 | |
18 | GNORM=0 | |
19 | DO 2 I=1,NEQU | |
20 | J=NUMPT(I) | |
21 | K=J+NC | |
22 | L=J-NC | |
23 | RES=DIFCO(I,5)-UMAT(J) | |
24 | IF(DIFCO(I,1))22,21,22 | |
25 | 22 RES=RES+DIFCO(I,1)*UMAT(J+1) | |
26 | 21 IF(DIFCO(I,2))24,23,24 | |
27 | 24 RES=RES+DIFCO(I,2)*UMAT(K) | |
28 | 23 IF(DIFCO(I,3))26,25,26 | |
29 | 26 RES=RES+DIFCO(I,3)*UMAT(J-1) | |
30 | 25 IF(DIFCO(I,4))28,27,28 | |
31 | 28 RES=RES+DIFCO(I,4)*UMAT(L) | |
32 | 27 CONTINUE | |
33 | DELTA=BETA*RES | |
34 | UMAT(J)=UMAT(J)+DELTA | |
35 | GNORM=GNORM+ABS(DELTA) | |
36 | IF(ABS(DELTA)-DMAX)2,2,1 | |
37 | 1 DMAX=ABS(DELTA) | |
38 | UCORR=UMAT(J) | |
39 | ICORR=I | |
40 | 2 CONTINUE | |
41 | RETURN | |
42 | END |