1 : C TRIANG SOURCE CHAT 09/10/09 21:25:38 6519 2 : SUBROUTINE TRIANG(KRIGI,PREC) 3 : C======================================================================= 4 : C ASSEMBLE LES PETITES MATRICES RIGIDITE ET LES MET SOUS LA FORME 5 : C Lt.D.L. 6 : C IL LE POINTEUR DE LA MATRICE RESULTANTE DANS ICHOLE 7 : C======================================================================= 8 : C 9 : IMPLICIT INTEGER(I-N) 10 : IMPLICIT REAL*8 (A-H,O-Z) 11 : -INC SMRIGID 12 : -INC CCOPTIO 13 : IF(IIMPI.EQ.1)THEN 13 : CALL GIBTEM(XKT) 14 : INTERR(1)=XKT 14 : CALL ERREUR(-259) 15 : ENDIF 16 : IF(IIMPI.EQ.1)WRITE(IOIMP,10) 17 : 10 FORMAT(' L''IMPRESSION PRECDENTE EST AVANT ASSEM1 ') 17 : CALL ASSEM1(KRIGI,MMATRX,INUINX,ITOPOX,IMINIX,IPOX,INCTRX,IITOPX) 18 : IF(IERR.NE.0) RETURN 19 : IF(IIMPI.EQ.1)THEN 19 : CALL GIBTEM(XKT) 20 : INTERR(1)=XKT 20 : CALL ERREUR(-259) 21 : ENDIF 22 : IF(IIMPI.EQ.1)WRITE(IOIMP,11) 23 : 11 FORMAT(' L''IMPRESSION PRECEDENTE EST AVANT ASSEM2') 23 : CALL ASSEM2(KRIGI,ITOPOX,INUINX,IMINIX,MMATRX,IPOX,INCTRX,IITOPX) 24 : IF(IERR.NE.0) RETURN 25 : IF(IIMPI.EQ.1)THEN 25 : CALL GIBTEM(XKT) 26 : INTERR(1)=XKT 26 : CALL ERREUR(-259) 27 : ENDIF 28 : IF(IIMPI.EQ.1)WRITE(IOIMP,12) 29 : 12 FORMAT(' L''IMPRESSION PRECEDENTE EST AVANT CHOLE ') 30 : IF(IERR.NE.0) GO TO 5000 30 : CALL CHOLE(MMATRX,PREC) 31 : IF(IERR.NE.0) RETURN 32 : IF(IIMPI.EQ.1)THEN 32 : CALL GIBTEM(XKT) 33 : INTERR(1)=XKT 33 : CALL ERREUR(-259) 34 : ENDIF 35 : IF(IERR.NE.0) GO TO 5000 36 : IF(IIMPI.EQ.1) WRITE(IOIMP,13) 37 : 13 FORMAT(' L''IMPRESSION PRECEDENTE EST APRES CHOLE') 38 : MRIGID=KRIGI 39 : SEGACT MRIGID*MOD 40 : ICHOLE=MMATRX 41 : SEGDES MRIGID 42 : 5000 CONTINUE 43 : RETURN 44 : END 45 : 46 : 47 : 48 : 49 : 50 : 51 : 52 :
© Cast3M 2003 - All rights reserved.
Disclaimer