Télécharger cinvat.eso

Retour à la liste

Numérotation des lignes :

cinvat
  1. C CINVAT SOURCE KK2000 14/04/09 21:15:08 8027
  2.  
  3. SUBROUTINE CINVAT(SIG,N,T)
  4. IMPLICIT INTEGER(I-N)
  5. integer n
  6. complex*16 sig(n),t,J2,J3
  7. call C_Invari_J2(sig,n,J2)
  8. t=(0.D0,0.D0)
  9. if (abs(J2).gt.1.D-17) then
  10. call C_Invari_J3(sig,n,J3)
  11. t=CMPLX(3.D0*sqrt(3.D0))*J3/(CMPLX(2.D0)*J2**CMPLX(1.5D0))
  12. c if (abs(t).lt.-1.D0) t=-1.D0
  13. c if (abs(t).gt.1.D0) t= 1.D0
  14. t=-(0.D0,1.D0)*log(t+sqrt((1.D0,0.D0)-t*t)*(0.D0,1.D0))/
  15. & (3.D0,0.D0)
  16. endif
  17. return
  18. end
  19.  
  20.  
  21.  
  22.  
  23.  

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