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