This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / d / epditr.F
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