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