Sources Esope | Cast3M

kprior
C KPRIOR    SOURCE    CHAT      05/01/13    01:04:48     5004
      SUBROUTINE  KPRIOR(NES,NSA,NSB,A,B,U,V,IVU)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C
      DIMENSION A(NES,NSA),B(NES,NSB),U(1),V(1)
C
      EPS = 1E-5
C
      LS1 = 0
      DO 1 IS=1,NSB
      F=0.
      DO 11 IES = 1,NES
      F = F+B(IES,IS)*U(IES)
 11   CONTINUE
      F = F + U(NES+1)
      LS2=KSIG(F,EPS)
      IF (IS.GE.2.AND.LS1.NE.LS2.AND.LS2.NE.0.AND.LS1.NE.0) GOTO 4
      IF(LS2.NE.0) LS1=LS2
 1    CONTINUE
C
      IF(LS1.EQ.1) THEN
         LS1 = 0
         DO 2 IS=1,NSA
          F=0.
           DO 21 IES = 1,NES
           F = F+A(IES,IS)*V(IES)
 21        CONTINUE
           F = F + V(NES+1)
           LS2=KSIG(F,EPS)
           IF (IS.GE.2.AND.LS1.NE.LS2.AND.LS2.NE.0.AND.LS1.NE.0) THEN
             IVU = 2
             RETURN
           ENDIF
           IF(LS2.NE.0) LS1=LS2
 2         CONTINUE
C
           IF(LS1.EQ.1)  THEN
             IVU=1
             RETURN
           ELSE
             IVU = 0
             RETURN
           ENDIF
C
      ELSE
        IVU = 0
        RETURN
      ENDIF
C
 4    CONTINUE
           F=0.
           DO 41 IES = 1,NES
           F = F+A(IES,1)*V(IES)
 41        CONTINUE
           F = F + V(NES+1)
           LS2 = KSIG(F,EPS)
           IF (LS2.EQ.-1) THEN
              IVU = 0
              RETURN
           ELSE
              IVU = 3
              RETURN
           ENDIF
      END                                                               00000070