Télécharger maxro3.eso

Retour à la liste

Numérotation des lignes :

  1. C MAXRO3 SOURCE CHAT 05/01/13 01:37:53 5004
  2. SUBROUTINE MAXRO3(ICAS,IDEFO,S,A,R,RT,TRA,NDIM,IDIM,IFOUR)
  3. *
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. *
  7. DIMENSION R(NDIM,*),RT(NDIM,*),A(NDIM,*),TRA(NDIM,*)
  8. DIMENSION S(*)
  9. *
  10.  
  11. * ICAS = 1 : RT A R
  12. * ICAS = 2 : R A RT
  13.  
  14.  
  15.  
  16.  
  17. A(1,1) = S(1)
  18. A(2,2) = S(2)
  19. A(1,2) = S(4)
  20. IF (IDEFO.EQ.1) A(1,2)=A(1,2)/2.D0
  21. A(2,1)=A(1,2)
  22. *
  23. IF (IFOUR.LT.1) GOTO 666
  24. *
  25. A(3,3) = S(3)
  26. A(3,1) = S(5)
  27. A(3,2) = S(6)
  28. *
  29. IF (IDEFO.EQ.1) A(3,1)=A(3,1)/2.D0
  30. IF (IDEFO.EQ.1) A(3,2)=A(3,2)/2.D0
  31. A(1,3)=A(3,1)
  32. A(2,3)=A(3,2)
  33. A(3,3) = S(3)
  34. *
  35. 666 CONTINUE
  36. *
  37.  
  38. IF(ICAS.EQ.1) THEN
  39.  
  40. * t
  41. * >>> Rotation du tenseur : A = R A R <<<
  42. *
  43. CALL MULMAT(TRA,A,R,NDIM,NDIM,NDIM)
  44. CALL MULMAT(A,RT,TRA,NDIM,NDIM,NDIM)
  45.  
  46.  
  47. ELSE IF(ICAS.EQ.2) THEN
  48.  
  49. * t
  50. * >>> Rotation du tenseur : A = R A R <<<
  51. *
  52. CALL MULMAT(TRA,A,RT,NDIM,NDIM,NDIM)
  53. CALL MULMAT(A,R,TRA,NDIM,NDIM,NDIM)
  54.  
  55. ENDIF
  56.  
  57. S(1) = A(1,1)
  58. S(2) = A(2,2)
  59. *
  60. IF (IDEFO.EQ.1) A(1,2)=A(1,2)*2.D0
  61. S(4) = A(1,2)
  62. *
  63. IF (IFOUR.GE.1) THEN
  64. *
  65. S(3) = A(3,3)
  66. *
  67. IF (IDEFO.EQ.1) A(3,1)=A(3,1)*2.D0
  68. IF (IDEFO.EQ.1) A(3,2)=A(3,2)*2.D0
  69. S(5) = A(3,1)
  70. S(6) = A(3,2)
  71. *
  72. ENDIF
  73. *
  74. *
  75. RETURN
  76. END
  77.  
  78.  
  79.  

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