Numérotation des lignes :

asbglb
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)
14. CALL ZDANUL(BGLOB,144)
15. CALL ZDANUL(BTMP,144)
16. C PREMIERE LIGNE:
17. DO J=1,8
18. C "MULTIPLICATION AVEC X" ---> UNE PARTIE DE EPS11
19. BTMP(1,J)=B(1,J)*P(1,1)*P(1,1)
20. ! +B(2,J)*P(1,1)*P(2,1)
21. ! +B(3,J)*P(3,1)*P(1,1)
22. END DO
23. DO J=9,16
24. C "MULTIPLICATION AVEC Y" ---> UNE PARTIE DE EPS11
25. BTMP(1,J)=P(2,1)*B(1,J-8)*P(1,1)
26. ! +P(2,1)*B(2,J-8)*P(2,1)
27. ! +P(2,1)*B(3,J-8)*P(3,1)
28. END DO
29. DO J=17,24
30. C "MULTIPLICATION AVEC Z" ---> UNE PARTIE DE EPS11
31. BTMP(1,J)=P(3,1)*B(1,J-16)*P(1,1)
32. ! +P(3,1)*B(2,J-16)*P(2,1)
33. ! +P(3,1)*B(3,J-16)*P(3,1)
34. END DO
35. C 2IEME LIGNE:
36. DO J=1,8
37. C "MULTIPLICATION AVEC X" ---> UNE PARTIE DE EPS22
38. BTMP(2,J)=P(1,2)*B(1,J)*P(1,2)
39. ! +P(1,2)*B(2,J)*P(2,2)
40. ! +P(1,2)*B(3,J)*P(3,2)
41. END DO
42. DO J=9,16
43. C "MULTIPLICATION AVEC Y" ---> UNE PARTIE DE EPS22
44. BTMP(2,J)=P(2,2)*B(1,J-8)*P(1,2)
45. ! +P(2,2)*B(2,J-8)*P(2,2)
46. ! +P(2,2)*B(3,J-8)*P(3,2)
47. END DO
48. DO J=17,24
49. C "MULTIPLICATION AVEC Z" ---> UNE PARTIE DE EPS22
50. BTMP(2,J)=P(3,2)*B(1,J-16)*P(1,2)
51. ! +P(3,2)*B(2,J-16)*P(2,2)
52. ! +P(3,2)*B(3,J-16)*P(3,2)
53. END DO
54. C 3IEME LIGNE:
55. DO J=1,8
56. C "MULTIPLICATION AVEC X" ---> UNE PARTIE DE EPS33
57. BTMP(3,J)=P(1,3)*B(1,J)*P(1,3)
58. ! +P(1,3)*B(2,J)*P(2,3)
59. ! +P(1,3)*B(3,J)*P(3,3)
60. END DO
61. DO J=9,16
62. C "MULTIPLICATION AVEC Y" ---> UNE PARTIE DE EPS33
63. BTMP(3,J)=P(2,3)*B(1,J-8)*P(1,3)
64. ! +P(2,3)*B(2,J-8)*P(2,3)
65. ! +P(2,3)*B(3,J-8)*P(3,3)
66. END DO
67. DO J=17,24
68. C "MULTIPLICATION AVEC Z" ---> UNE PARTIE DE EPS33
69. BTMP(3,J)=P(3,3)*B(1,J-16)*P(1,3)
70. ! +P(3,3)*B(2,J-16)*P(2,3)
71. ! +P(3,3)*B(3,J-16)*P(3,3)
72. END DO
73. C 4IEME LIGNE:
74. DO J=1,8
75. C "MULTIPLICATION AVEC X" ---> UNE PARTIE DE EPS12
76. BTMP(4,J)=B(1,J)*P(1,1)*P(1,2)*2.
77. ! +B(2,J)*P(1,2)*P(2,1)
78. ! +B(2,J)*P(1,1)*P(2,2)
79. ! +B(3,J)*P(1,2)*P(3,1)
80. ! +B(3,J)*P(1,1)*P(3,2)
81. END DO
82. DO J=9,16
83. C "MULTIPLICATION AVEC Y" ---> UNE PARTIE DE EPS12
84. BTMP(4,J)=B(1,J-8)*P(1,2)*P(2,1)
85. ! +B(1,J-8)*P(1,1)*P(2,2)
86. ! +B(2,J-8)*P(2,1)*P(2,2)*2.
87. ! +B(3,J-8)*P(2,2)*P(3,1)
88. ! +B(3,J-8)*P(2,1)*P(3,2)
89. END DO
90. DO J=17,24
91. C "MULTIPLICATION AVEC Z" ---> UNE PARTIE DE EPS12
92. BTMP(4,J)=B(1,J-16)*P(1,2)*P(3,1)
93. ! +B(2,J-16)*P(2,2)*P(3,1)
94. ! +B(1,J-16)*P(1,1)*P(3,2)
95. ! +B(2,J-16)*P(2,1)*P(3,2)
96. ! +B(3,J-16)*P(3,1)*P(3,2)*2.
97. END DO
98. C 5IEME LIGNE:
99. DO J=1,8
100. C "MULTIPLICATION AVEC X" ---> UNE PARTIE DE EPS23
101. BTMP(5,J)=B(1,J)*P(1,2)*P(1,3)*2.
102. ! +B(2,J)*P(1,3)*P(2,2)
103. ! +B(2,J)*P(1,2)*P(2,3)
104. ! +B(3,J)*P(1,3)*P(3,2)
105. ! +B(3,J)*P(1,2)*P(3,3)
106. END DO
107. DO J=9,16
108. C "MULTIPLICATION AVEC Y" ---> UNE PARTIE DE EPS23
109. BTMP(5,J)=B(1,J-8)*P(1,3)*P(2,2)
110. ! +B(1,J-8)*P(1,2)*P(2,3)
111. ! +B(2,J-8)*P(2,2)*P(2,3)*2.
112. ! +B(3,J-8)*P(2,3)*P(3,2)
113. ! +B(3,J-8)*P(2,2)*P(3,3)
114. END DO
115. DO J=17,24
116. C "MULTIPLICATION AVEC Z" ---> UNE PARTIE DE EPS23
117. BTMP(5,J)=B(1,J-16)*P(1,3)*P(3,2)
118. ! +B(2,J-16)*P(2,3)*P(3,2)
119. ! +B(1,J-16)*P(1,2)*P(3,3)
120. ! +B(2,J-16)*P(2,2)*P(3,3)
121. ! +B(3,J-16)*P(3,2)*P(3,3)*2.
122. END DO
123. C 6IEME LIGNE:
124. DO J=1,8
125. C "MULTIPLICATION AVEC X" ---> UNE PARTIE DE EPS13
126. BTMP(6,J)=B(1,J)*P(1,1)*P(1,3)*2
127. ! +B(2,J)*P(1,3)*P(2,1)
128. ! +B(2,J)*P(1,1)*P(2,3)
129. ! +B(3,J)*P(1,3)*P(3,1)
130. ! +B(3,J)*P(1,1)*P(3,3)
131. END DO
132. DO J=9,16
133. C "MULTIPLICATION AVEC Y" ---> UNE PARTIE DE EPS13
134. BTMP(6,J)=B(1,J-8)*P(1,3)*P(2,1)
135. ! +B(1,J-8)*P(1,1)*P(2,3)
136. ! +B(2,J-8)*P(2,1)*P(2,3)*2.
137. ! +B(3,J-8)*P(2,3)*P(3,1)
138. ! +B(3,J-8)*P(2,1)*P(3,3)
139. END DO
140. DO J=17,24
141. C "MULTIPLICATION AVEC Z" ---> UNE PARTIE DE EPS13
142. BTMP(6,J)=B(1,J-16)*P(1,3)*P(3,1)
143. ! +B(2,J-16)*P(2,3)*P(3,1)
144. ! +B(1,J-16)*P(1,1)*P(3,3)
145. ! +B(2,J-16)*P(2,1)*P(3,3)
146. ! +B(3,J-16)*P(3,1)*P(3,3)*2.
147. END DO
148. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
149. DO J=1,24
150. DO I=1,6
151. BGLOB(I,J)=BTMP(I,J)
152. END DO
153. END DO
154. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
155. RETURN
156. END
157.
158.
159.

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