Télécharger rempmch.eso

Retour à la liste

Numérotation des lignes :

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

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