Télécharger cdinij.eso

Retour à la liste

Numérotation des lignes :

cdinij
  1. C CDINIJ SOURCE KK2000 14/04/10 21:15:03 8031
  2.  
  3. SUBROUTINE C_DERINVARIIJJ (SIG,I_S,J2_S,J3_S,N)
  4. IMPLICIT INTEGER(I-N)
  5. integer n,i
  6. complex*16 sig(n),I_s(n),J2_s(n),J3_s(n),aux1,J2
  7. call CDSVIA(sig,J2_s,n)
  8. call C_Invari_J2(sig,n,J2)
  9. J3_s(1)=J2_s(1)*J2_s(1)
  10. J3_s(2)=J2_s(2)*J2_s(2)
  11. J3_s(3)=J2_s(3)*J2_s(3)
  12. if (n.ge.4) then
  13. J3_s(1)=J3_s(1)+J2_s(4)*J2_s(4)
  14. J3_s(2)=J3_s(2)+J2_s(4)*J2_s(4)
  15. J3_s(4)=2.D0*(J2_s(1)+J2_s(2))*J2_s(4)
  16. if (n.eq.6) then
  17. J3_s(1)=J3_s(1)+J2_s(5)*J2_s(5)
  18. J3_s(2)=J3_s(2)+J2_s(6)*J2_s(6)
  19. J3_s(3)=J3_s(3)+J2_s(5)*J2_s(5)+J2_s(6)*J2_s(6)
  20. J3_s(4)=J3_s(4)+J2_s(5)*J2_s(6)
  21. J3_s(5)=(J2_s(1)+J2_s(3))*J2_s(5)+J2_s(4)*J2_s(6)
  22. J3_s(6)=(J2_s(2)+J2_s(3))*J2_s(6)+J2_s(4)*J2_s(5)
  23. endif
  24. endif
  25. aux1=-J2*2.D0/3.D0
  26. do i=1,3
  27. J3_s(i)=J3_s(i)+aux1
  28. I_s(i)=1.D0
  29. enddo
  30. if (n.lt.4) return
  31. do i=4,n
  32. I_s(i)=0.D0
  33. J2_s(i)=2.D0*J2_s(i)
  34. enddo
  35. return
  36. end
  37.  
  38.  
  39.  
  40.  
  41.  

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