Télécharger ddinij.eso

Retour à la liste

Numérotation des lignes :

  1. C DDINIJ SOURCE KK2000 14/04/09 21:15:14 8027
  2.  
  3. SUBROUTINE DDINIJ (SIG,PIJKL,WIJKL,N)
  4. IMPLICIT INTEGER(I-N)
  5. integer n,i,j
  6. real*8 sig(n),pijkl(6,6),wijkl(6,6),desv(6)
  7. do i=1,6
  8. do j=1,6
  9. wijkl(i,j)=0.D0
  10. pijkl(i,j)=0.D0
  11. enddo
  12. enddo
  13. call Desviador(sig,desv,n)
  14. do i=1,3
  15. do j=1,3
  16. pijkl(i,j)=-1.D0/3.D0
  17. enddo
  18. pijkl(i,i)=2.D0/3.D0
  19. enddo
  20. wijkl(1,1)=desv(1)
  21. wijkl(2,2)=desv(2)
  22. wijkl(3,3)=desv(3)
  23. wijkl(2,1)=desv(3)
  24. wijkl(3,1)=desv(2)
  25. wijkl(3,2)=desv(1)
  26. if (n.ge.4) then
  27. do i=4,n
  28. pijkl(i,i)=2.D0
  29. enddo
  30. wijkl(4,1)=desv(4)
  31. wijkl(4,2)=desv(4)
  32. wijkl(4,3)=-2.D0*desv(4)
  33. wijkl(4,4)=-3.D0*desv(3)
  34. endif
  35. do i=1,n
  36. do j=1,i
  37. wijkl(i,j)=2.D0/3.D0*wijkl(i,j)
  38. enddo
  39. enddo
  40. do i=2,n
  41. do j=1,i-1
  42. wijkl(j,i)=wijkl(i,j)
  43. enddo
  44. enddo
  45. return
  46. end
  47.  
  48.  
  49.  
  50.  
  51.  

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