Télécharger asbglb.eso

Retour à la liste

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