Télécharger rempth.eso

Retour à la liste

Numérotation des lignes :

rempth
  1. C REMPTH SOURCE CHAT 08/10/02 21:15:46 6168
  2. SUBROUTINE REMPTH
  3. & (REL1,REL2,REL3,REL4,REL5,REL6,REL7,REL8,
  4. & REL9,RE10,RE11,RE12,LRE,RE)
  5. C====================================================================
  6. C
  7. C TRANSVASE REL (LRE,LRE) MATRICE DE RIGIDITE TRIANGULAIRE
  8. C INFERIEURE DANS RE(lre,lre) STOCKE LIGNE PAR LIGNE
  9. C
  10. C====================================================================
  11. IMPLICIT INTEGER(I-N)
  12. IMPLICIT REAL*8(A-H,O-Z)
  13. DIMENSION REL1(LRE,LRE),REL2(LRE,LRE),REL3(LRE,LRE),REL4(LRE,LRE)
  14. DIMENSION REL5(LRE,LRE),REL6(LRE,LRE),REL7(LRE,LRE),REL8(LRE,LRE)
  15. DIMENSION REL9(LRE,LRE),RE10(LRE,LRE),RE11(LRE,LRE),RE12(LRE,LRE)
  16. DIMENSION RE(3*LRE,3*LRE)
  17. C
  18. DO 100 IA=1,LRE
  19. DO 100 IB=1,LRE
  20. RE(IA,IB)=REL1(IA,IB)
  21. 100 CONTINUE
  22.  
  23. DO 200 IA=LRE+1,2*LRE
  24. DO 200 IB=1,LRE
  25. RE(IA,IB)=REL4(IA-LRE,IB)
  26. 200 CONTINUE
  27.  
  28. DO 300 IA=(2*LRE)+1,3*LRE
  29. DO 300 IB=1,LRE
  30. RE(IA,IB)=REL7(IA-2*LRE,IB) +RE10(IA-2*LRE,IB)
  31. 300 CONTINUE
  32.  
  33. DO 400 IA=1,LRE
  34. DO 400 IB=LRE+1,2*LRE
  35. RE(IA,IB)=REL2(IA,IB-LRE)
  36. 400 CONTINUE
  37.  
  38. DO 500 IA=LRE+1,2*LRE
  39. DO 500 IB=LRE+1,2*LRE
  40. RE(IA,IB)=REL5(IA-LRE,IB-LRE)
  41. 500 CONTINUE
  42.  
  43. DO 600 IA=(2*LRE)+1,3*LRE
  44. DO 600 IB=LRE+1,2*LRE
  45. RE(IA,IB)=REL8(IA-2*LRE,IB-LRE)+RE11(IA-2*LRE,IB-LRE)
  46. 600 CONTINUE
  47.  
  48. DO 700 IA=1,LRE
  49. DO 700 IB=2*LRE+1,3*LRE
  50. RE(IA,IB)=REL3(IA,IB-2*LRE)
  51. 700 CONTINUE
  52.  
  53. DO 800 IA=LRE+1,2*LRE
  54. DO 800 IB=2*LRE+1,3*LRE
  55. RE(IA,IB)=REL6(IA-LRE,IB-2*LRE)
  56. 800 CONTINUE
  57.  
  58. DO 900 IA=(2*LRE)+1,3*LRE
  59. DO 900 IB=2*LRE+1,3*LRE
  60. RE(IA,IB)=REL9(IA-2*LRE,IB-2*LRE)+RE12(IA-2*LRE,IB-2*LRE)
  61. 900 CONTINUE
  62.  
  63. RETURN
  64. END
  65.  
  66.  
  67.  

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