Sources Esope | Cast3M

rempmch
C REMPMCH   SOURCE    CB215821  26/08/24    21:18:12     12622          
C REMPMT    SOURCE    CHAT      05/01/13    02:50:12     5004
      SUBROUTINE REMPMCH
     &  (REL1,REL2,REL3,REL4,REL5,REL6,REL7,REL8,
     &  REL9,LRE,RE)
C====================================================================
C
C     TRANSVASE REL (LRE,LRE) MATRICE DE RIGIDITE TRIANGULAIRE
C     INFERIEURE DANS RE(lre,lre) STOCKE LIGNE PAR LIGNE
C
C====================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION REL1(LRE,LRE),REL2(LRE,LRE),REL3(LRE,LRE),REL4(LRE,LRE)
      DIMENSION REL5(LRE,LRE),REL6(LRE,LRE),REL7(LRE,LRE),REL8(LRE,LRE)
      DIMENSION REL9(LRE,LRE)
      DIMENSION RE(3*LRE,3*LRE)
 
*      JA=1
      DO 901 IA=1,LRE
      DO 100 IB=1,LRE
      RE(IA,IB)=REL1(IA,IB)
 100  CONTINUE
  901 CONTINUE
 
        DO 902 IA=LRE+1,2*LRE
      DO 200 IB=1,LRE
      RE(IA,IB)=REL4(IA-LRE,IB)
 200  CONTINUE
  902   CONTINUE
 
        DO 903 IA=(2*LRE)+1,3*LRE
      DO 300 IB=1,LRE
      RE(IA,IB)=REL7(IA-2*LRE,IB)
 300  CONTINUE
  903   CONTINUE
 
 
      DO 904 IA=1,LRE
      DO 400 IB=LRE+1,2*LRE
      RE(IA,IB)=REL2(IA,IB-LRE)
 400  CONTINUE
  904 CONTINUE
 
        DO 905 IA=LRE+1,2*LRE
      DO 500 IB=LRE+1,2*LRE
      RE(IA,IB)=REL5(IA-LRE,IB-LRE)
 500  CONTINUE
  905   CONTINUE
 
        DO 906 IA=(2*LRE)+1,3*LRE
      DO 600 IB=LRE+1,2*LRE
      RE(IA,IB)=REL8(IA-2*LRE,IB-LRE)
 600  CONTINUE
  906   CONTINUE
 
 
      DO 907 IA=1,LRE
      DO 700 IB=2*LRE+1,3*LRE
      RE(IA,IB)=REL3(IA,IB-2*LRE)
 700  CONTINUE
  907 CONTINUE
 
        DO 908 IA=LRE+1,2*LRE
      DO 800 IB=2*LRE+1,3*LRE
      RE(IA,IB)=REL6(IA-LRE,IB-2*LRE)
 800  CONTINUE
  908   CONTINUE
 
        DO 909 IA=(2*LRE)+1,3*LRE
      DO 900 IB=2*LRE+1,3*LRE
      RE(IA,IB)=REL9(IA-2*LRE,IB-2*LRE)
 900  CONTINUE
  909   CONTINUE
      RETURN
      END