Télécharger la plaquette Cast3M


fr en

Télécharger deri2.eso

Retour à la liste

Numérotation des lignes : 

deri2

  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.  
  56.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales