Télécharger cinvj3.eso

Retour à la liste

Numérotation des lignes :

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

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