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 :
© Cast3M 2003 - All rights reserved.
Disclaimer