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