Esope source files

   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 :