Esope source files

   1 : C DERI2     SOURCE    CHAT      05/01/12    22:42:50     5004           
   2 :       SUBROUTINE DERI2(D2FSIG,DPHI,P,S,NSTRS,BETJEF)
   3 : C
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8 (A-H,O-Z)
   6 :       DIMENSION S(4),P(4,4),D2FSIG(4,4)
   7 :       DIMENSION AB(4),AC(4),DP(4,4),P1(4,4)
   8 : *
   9 :       SEGMENT BETJEF
  10 :        REAL*8 AA,BETA,RB,ALPHA,EX,PXY,GFC,GFT,CAR,ETA,TDEF,
  11 :      &          TCON,DPSTF1,DPSTF2,TETA,PDT
  12 :        INTEGER ICT,ICC,IMOD,IVISS,ITER,
  13 :      &                ISIM,IBB1,IGAU1,IZON
  14 :       ENDSEGMENT
  15 : *
  16 : *     COMMON /DBETJEF/AA,BETA,RB,ALPHA,EX,PXY,GFC,GFT,CAR,ETA,TDEF,
  17 : *    &          TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,ITER,
  18 : *    &                ISIM,IBB1,IGAU1,IZON
  19 : C             *************************************************
  20 : C             **  CALCUL DE LA DERIVE SECONDE DU CRITERE     **
  21 : C             *************************************************
  22 : C-------------------------------------------------------------------
  23 :       DO 1 I=1,NSTRS
  24 :         AB(I)=0.D0
  25 :       DO 1 J=1,NSTRS
  26 :         AB(I)=AB(I)+P(I,J)*S(J)
  27 :     1 CONTINUE
  28 : C
  29 :       DO 6 I=1,NSTRS
  30 :         AC(I)=0.D0
  31 :       DO 6 J=1,NSTRS
  32 :         AC(I)=AC(I)+P(J,I)*S(J)
  33 :     6 CONTINUE
  34 : C
  35 :       DO 7 I=1,NSTRS
  36 :       DO 7 J=1,NSTRS
  37 :        DP(I,J)=0.D0
  38 :        DP(I,J)=AB(I)*AC(J)
  39 :     7 CONTINUE
  40 : C
  41 :       DO 9 I=1,NSTRS
  42 :       DO 9 J=1,NSTRS
  43 :        P1(I,J)=P(I,J)/2.D0/DPHI
  44 :        DP(I,J)=DP(I,J)/4.D0/DPHI/DPHI/DPHI
  45 :        D2FSIG(I,J)=P1(I,J)-DP(I,J)
  46 :     9 CONTINUE
  47 : C-------------------------------------------------------------------
  48 :       RETURN
  49 :       END
  50 : 
  51 : 
  52 : 
  53 : 
  54 :  
  55 :