Télécharger la plaquette Cast3M


fr en

Télécharger vectro.eso

Retour à la liste

Numérotation des lignes : 

vectro

  1. C VECTRO    SOURCE    CHAT      07/10/22    21:18:14     5921
  2.       SUBROUTINE VECTRO(RR,X,NN)
  3.       implicit real*8(a-h,o-z)
  4.       implicit integer(i-n)
  5.  
  6.  
  7. *      SUBROUTINE VECTROT(RR,X,NN)
  8. *
  9. *
  10. *      IMPLICIT NONE
  11. *
  12. *
  13. *
  14. *      REAL *8 RR(3,3)
  15. *      INTEGER NN
  16. *
  17. *      REAL *8 X(24)
  18.        dimension RR(3,3),X(24)
  19. *
  20. *
  21. *
  22. * ---   VARIABLES LOCALES
  23. *
  24. *      REAL *8 AA(24),R(3,3)
  25. *      INTEGER I,J
  26.        dimension AA(24),R(3,3)
  27. *
  28. *
  29.       IF(NN.EQ.2) THEN
  30.         DO I=1,3
  31.            DO J=1,3
  32.               R(I,J)= RR(J,I)
  33.            END DO
  34.         END DO
  35.       ELSEIF(NN.EQ.1) THEN
  36.         DO I=1,3
  37.           DO J=1,3
  38.               R(I,J)= RR(I,J)
  39.           END DO
  40.         END DO
  41.       ELSE
  42.         STOP 'VECTROT'
  43.       ENDIF
  44. *
  45.       AA(1) = R(1,1)*X(1) + R(1,2)*X(2) + R(1,3)*X(3)
  46.       AA(2) = R(2,1)*X(1) + R(2,2)*X(2) + R(2,3)*X(3)
  47.       AA(3) = R(3,1)*X(1) + R(3,2)*X(2) + R(3,3)*X(3)
  48. *
  49.       AA(4) = R(1,1)*X(4) + R(1,2)*X(5) + R(1,3)*X(6)
  50.       AA(5) = R(2,1)*X(4) + R(2,2)*X(5) + R(2,3)*X(6)
  51.       AA(6) = R(3,1)*X(4) + R(3,2)*X(5) + R(3,3)*X(6)
  52. *
  53.       AA(7) = R(1,1)*X(7) + R(1,2)*X(8) + R(1,3)*X(9)
  54.       AA(8) = R(2,1)*X(7) + R(2,2)*X(8) + R(2,3)*X(9)
  55.       AA(9) = R(3,1)*X(7) + R(3,2)*X(8) + R(3,3)*X(9)
  56. *
  57.       AA(10) = R(1,1)*X(10) + R(1,2)*X(11) + R(1,3)*X(12)
  58.       AA(11) = R(2,1)*X(10) + R(2,2)*X(11) + R(2,3)*X(12)
  59.       AA(12) = R(3,1)*X(10) + R(3,2)*X(11) + R(3,3)*X(12)
  60. *
  61.       AA(13) = R(1,1)*X(13) + R(1,2)*X(14) + R(1,3)*X(15)
  62.       AA(14) = R(2,1)*X(13) + R(2,2)*X(14) + R(2,3)*X(15)
  63.       AA(15) = R(3,1)*X(13) + R(3,2)*X(14) + R(3,3)*X(15)
  64. *
  65.       AA(16) = R(1,1)*X(16) + R(1,2)*X(17) + R(1,3)*X(18)
  66.       AA(17) = R(2,1)*X(16) + R(2,2)*X(17) + R(2,3)*X(18)
  67.       AA(18) = R(3,1)*X(16) + R(3,2)*X(17) + R(3,3)*X(18)
  68. *
  69.       AA(19) = R(1,1)*X(19) + R(1,2)*X(20) + R(1,3)*X(21)
  70.       AA(20) = R(2,1)*X(19) + R(2,2)*X(20) + R(2,3)*X(21)
  71.       AA(21) = R(3,1)*X(19) + R(3,2)*X(20) + R(3,3)*X(21)
  72. *
  73.       AA(22) = R(1,1)*X(22) + R(1,2)*X(23) + R(1,3)*X(24)
  74.       AA(23) = R(2,1)*X(22) + R(2,2)*X(23) + R(2,3)*X(24)
  75.       AA(24) = R(3,1)*X(22) + R(3,2)*X(23) + R(3,3)*X(24)
  76. *
  77.       DO I=1,24
  78.           X(I)= AA(I)
  79.       END DO
  80. *
  81. *
  82.       END
  83.  
  84.  
  85.  

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