1 : C ASBGLB SOURCE CHAT 07/10/22 21:15:05 5921 2 : SUBROUTINE ASBGLB(BGLOB,B,P) 3 : IMPLICIT REAL*8(A-H,O-Z) 4 : IMPLICIT INTEGER (I-N) 5 : 6 : * SUBROUTINE ASSE_BGLOB(BGLOB,B,P) 7 : C ANTOINE LEGAY LE 30/01/2001 8 : C TRANSFORME B(6,24) DANS LOCAL EN BGLOB(6,24) DANS GLOBAL 9 : C AVEC P(3,3) MATRICE DE PASSAGE 10 : * REAL*8 BGLOB(6,24),B(3,8),P(3,3),BTMP(6,24) 11 : * INTEGER J 12 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 13 : dimension BGLOB(6,24),B(3,8),P(3,3),BTMP(6,24) 13 : CALL ZDANUL(BGLOB,144) 13 : CALL ZDANUL(BTMP,144) 14 : C PREMIERE LIGNE: 15 : DO J=1,8 16 : C "MULTIPLICATION AVEC X" ---> UNE PARTIE DE EPS11 17 : BTMP(1,J)=B(1,J)*P(1,1)*P(1,1) 18 : ! +B(2,J)*P(1,1)*P(2,1) 19 : ! +B(3,J)*P(3,1)*P(1,1) 20 : END DO 21 : DO J=9,16 22 : C "MULTIPLICATION AVEC Y" ---> UNE PARTIE DE EPS11 23 : BTMP(1,J)=P(2,1)*B(1,J-8)*P(1,1) 24 : ! +P(2,1)*B(2,J-8)*P(2,1) 25 : ! +P(2,1)*B(3,J-8)*P(3,1) 26 : END DO 27 : DO J=17,24 28 : C "MULTIPLICATION AVEC Z" ---> UNE PARTIE DE EPS11 29 : BTMP(1,J)=P(3,1)*B(1,J-16)*P(1,1) 30 : ! +P(3,1)*B(2,J-16)*P(2,1) 31 : ! +P(3,1)*B(3,J-16)*P(3,1) 32 : END DO 33 : C 2IEME LIGNE: 34 : DO J=1,8 35 : C "MULTIPLICATION AVEC X" ---> UNE PARTIE DE EPS22 36 : BTMP(2,J)=P(1,2)*B(1,J)*P(1,2) 37 : ! +P(1,2)*B(2,J)*P(2,2) 38 : ! +P(1,2)*B(3,J)*P(3,2) 39 : END DO 40 : DO J=9,16 41 : C "MULTIPLICATION AVEC Y" ---> UNE PARTIE DE EPS22 42 : BTMP(2,J)=P(2,2)*B(1,J-8)*P(1,2) 43 : ! +P(2,2)*B(2,J-8)*P(2,2) 44 : ! +P(2,2)*B(3,J-8)*P(3,2) 45 : END DO 46 : DO J=17,24 47 : C "MULTIPLICATION AVEC Z" ---> UNE PARTIE DE EPS22 48 : BTMP(2,J)=P(3,2)*B(1,J-16)*P(1,2) 49 : ! +P(3,2)*B(2,J-16)*P(2,2) 50 : ! +P(3,2)*B(3,J-16)*P(3,2) 51 : END DO 52 : C 3IEME LIGNE: 53 : DO J=1,8 54 : C "MULTIPLICATION AVEC X" ---> UNE PARTIE DE EPS33 55 : BTMP(3,J)=P(1,3)*B(1,J)*P(1,3) 56 : ! +P(1,3)*B(2,J)*P(2,3) 57 : ! +P(1,3)*B(3,J)*P(3,3) 58 : END DO 59 : DO J=9,16 60 : C "MULTIPLICATION AVEC Y" ---> UNE PARTIE DE EPS33 61 : BTMP(3,J)=P(2,3)*B(1,J-8)*P(1,3) 62 : ! +P(2,3)*B(2,J-8)*P(2,3) 63 : ! +P(2,3)*B(3,J-8)*P(3,3) 64 : END DO 65 : DO J=17,24 66 : C "MULTIPLICATION AVEC Z" ---> UNE PARTIE DE EPS33 67 : BTMP(3,J)=P(3,3)*B(1,J-16)*P(1,3) 68 : ! +P(3,3)*B(2,J-16)*P(2,3) 69 : ! +P(3,3)*B(3,J-16)*P(3,3) 70 : END DO 71 : C 4IEME LIGNE: 72 : DO J=1,8 73 : C "MULTIPLICATION AVEC X" ---> UNE PARTIE DE EPS12 74 : BTMP(4,J)=B(1,J)*P(1,1)*P(1,2)*2. 75 : ! +B(2,J)*P(1,2)*P(2,1) 76 : ! +B(2,J)*P(1,1)*P(2,2) 77 : ! +B(3,J)*P(1,2)*P(3,1) 78 : ! +B(3,J)*P(1,1)*P(3,2) 79 : END DO 80 : DO J=9,16 81 : C "MULTIPLICATION AVEC Y" ---> UNE PARTIE DE EPS12 82 : BTMP(4,J)=B(1,J-8)*P(1,2)*P(2,1) 83 : ! +B(1,J-8)*P(1,1)*P(2,2) 84 : ! +B(2,J-8)*P(2,1)*P(2,2)*2. 85 : ! +B(3,J-8)*P(2,2)*P(3,1) 86 : ! +B(3,J-8)*P(2,1)*P(3,2) 87 : END DO 88 : DO J=17,24 89 : C "MULTIPLICATION AVEC Z" ---> UNE PARTIE DE EPS12 90 : BTMP(4,J)=B(1,J-16)*P(1,2)*P(3,1) 91 : ! +B(2,J-16)*P(2,2)*P(3,1) 92 : ! +B(1,J-16)*P(1,1)*P(3,2) 93 : ! +B(2,J-16)*P(2,1)*P(3,2) 94 : ! +B(3,J-16)*P(3,1)*P(3,2)*2. 95 : END DO 96 : C 5IEME LIGNE: 97 : DO J=1,8 98 : C "MULTIPLICATION AVEC X" ---> UNE PARTIE DE EPS23 99 : BTMP(5,J)=B(1,J)*P(1,2)*P(1,3)*2. 100 : ! +B(2,J)*P(1,3)*P(2,2) 101 : ! +B(2,J)*P(1,2)*P(2,3) 102 : ! +B(3,J)*P(1,3)*P(3,2) 103 : ! +B(3,J)*P(1,2)*P(3,3) 104 : END DO 105 : DO J=9,16 106 : C "MULTIPLICATION AVEC Y" ---> UNE PARTIE DE EPS23 107 : BTMP(5,J)=B(1,J-8)*P(1,3)*P(2,2) 108 : ! +B(1,J-8)*P(1,2)*P(2,3) 109 : ! +B(2,J-8)*P(2,2)*P(2,3)*2. 110 : ! +B(3,J-8)*P(2,3)*P(3,2) 111 : ! +B(3,J-8)*P(2,2)*P(3,3) 112 : END DO 113 : DO J=17,24 114 : C "MULTIPLICATION AVEC Z" ---> UNE PARTIE DE EPS23 115 : BTMP(5,J)=B(1,J-16)*P(1,3)*P(3,2) 116 : ! +B(2,J-16)*P(2,3)*P(3,2) 117 : ! +B(1,J-16)*P(1,2)*P(3,3) 118 : ! +B(2,J-16)*P(2,2)*P(3,3) 119 : ! +B(3,J-16)*P(3,2)*P(3,3)*2. 120 : END DO 121 : C 6IEME LIGNE: 122 : DO J=1,8 123 : C "MULTIPLICATION AVEC X" ---> UNE PARTIE DE EPS13 124 : BTMP(6,J)=B(1,J)*P(1,1)*P(1,3)*2 125 : ! +B(2,J)*P(1,3)*P(2,1) 126 : ! +B(2,J)*P(1,1)*P(2,3) 127 : ! +B(3,J)*P(1,3)*P(3,1) 128 : ! +B(3,J)*P(1,1)*P(3,3) 129 : END DO 130 : DO J=9,16 131 : C "MULTIPLICATION AVEC Y" ---> UNE PARTIE DE EPS13 132 : BTMP(6,J)=B(1,J-8)*P(1,3)*P(2,1) 133 : ! +B(1,J-8)*P(1,1)*P(2,3) 134 : ! +B(2,J-8)*P(2,1)*P(2,3)*2. 135 : ! +B(3,J-8)*P(2,3)*P(3,1) 136 : ! +B(3,J-8)*P(2,1)*P(3,3) 137 : END DO 138 : DO J=17,24 139 : C "MULTIPLICATION AVEC Z" ---> UNE PARTIE DE EPS13 140 : BTMP(6,J)=B(1,J-16)*P(1,3)*P(3,1) 141 : ! +B(2,J-16)*P(2,3)*P(3,1) 142 : ! +B(1,J-16)*P(1,1)*P(3,3) 143 : ! +B(2,J-16)*P(2,1)*P(3,3) 144 : ! +B(3,J-16)*P(3,1)*P(3,3)*2. 145 : END DO 146 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 147 : DO J=1,24 148 : DO I=1,6 149 : BGLOB(I,J)=BTMP(I,J) 150 : END DO 151 : END DO 152 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 153 : RETURN 154 : END 155 : 156 :
© Cast3M 2003 - All rights reserved.
Disclaimer