1 : C GRACO1 SOURCE PV 13/10/21 21:15:20 7845 2 : SUBROUTINE GRACO1 ( KRIGI) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : C 6 : C COPIE DE TRIANG POUR FAIRE UN CHOLEVSKI INCOMPLET 7 : C 8 : -INC SMRIGID 9 : -INC CCOPTIO 10 : -INC SMMATRI 11 : REAL*8 XKT 12 : IF(IIMPI.EQ.1)THEN 12 : CALL GIBTEM(XKT) 13 : INTERR(1)=XKT 13 : CALL ERREUR(-259) 14 : ENDIF 15 : IF(IIMPI.EQ.1)WRITE(IOIMP,10) 16 : 10 FORMAT(' L''IMPRESSION PRECDENTE EST AVANT ASSEM1 ') 16 : CALL ASSEM1(KRIGI,MMATRX,INUINX,ITOPOX,IMINIX,IPOX,INCTRX,IITOPX) 17 : IF(IERR.NE.0) RETURN 18 : IF(IIMPI.EQ.1)THEN 18 : CALL GIBTEM(XKT) 19 : INTERR(1)=XKT 19 : CALL ERREUR(-259) 20 : ENDIF 21 : IF(IIMPI.EQ.1)WRITE(IOIMP,11) 22 : 11 FORMAT(' L''IMPRESSION PRECEDENTE EST AVANT ASSEM2') 22 : CALL ASSEM2(KRIGI,ITOPOX,INUINX,IMINIX,MMATRX,IPOX,INCTRX,IITOPX) 23 : IF(IERR.NE.0) RETURN 24 : IF(IIMPI.EQ.1)THEN 24 : CALL GIBTEM(XKT) 25 : INTERR(1)=XKT 25 : CALL ERREUR(-259) 26 : ENDIF 26 : CALL GRACO2(MMATRX) 27 : IF(IERR.NE.0) RETURN 28 : IF(IIMPI.EQ.1)THEN 28 : CALL GIBTEM(XKT) 29 : INTERR(1)=XKT 29 : CALL ERREUR(-259) 30 : ENDIF 31 : IF(IERR.NE.0) GO TO 5000 32 : IF(IIMPI.EQ.1) WRITE(IOIMP,13) 33 : 13 FORMAT(' L''IMPRESSION PRECEDENTE EST APRES GRACO2') 34 : MRIGID=KRIGI 35 : SEGACT MRIGID*MOD 36 : ICHOLE=MMATRX 37 : SEGDES MRIGID 38 : 5000 CONTINUE 39 : RETURN 40 : END 41 : 42 : 43 : 44 : 45 : 46 : 47 : 48 : 49 : 50 :
© Cast3M 2003 - All rights reserved.
Disclaimer