1 : C ACTIVE SOURCE FANDEUR 10/12/14 21:15:23 6812 2 : SUBROUTINE ACTIVE 3 : C-------------------------------------------------------------------- 4 : C OPERATEUR ACCELERATION DE CONVERGENCE 5 : C 6 : C SYNTAXE : VOIR NOTICE 7 : C-------------------------------------------------------------------- 8 : IMPLICIT INTEGER(I-N) 9 : IMPLICIT REAL*8(A-H,O-Z) 10 : C 11 : -INC CCOPTIO 12 : C 13 : CHARACTER*4 MOCOMP 14 : CHARACTER*4 LISMO(2) 15 : DATA LISMO/'GEOM','SECA'/ 16 : C 16 : CALL LIRMOT(LISMO,2,IDET,1) 17 : IF(IERR.NE.0) RETURN 18 : C 19 : IF(IDET.EQ.1) THEN 20 : FLOT=0.98D0 20 : CALL LIRREE(FLOT,0,IRETOU) 21 : ENDIF 22 : C 23 : C CAS DES CHPOINTS - ACCELERATIONS GEOMETRIQUE OU SECANTE 24 : C 24 : CALL LIROBJ('CHPOINT',MCHPO1,0,IRETOU) 25 : IF(IERR.NE.0) RETURN 26 : IF(IRETOU.EQ.1) THEN 26 : CALL LIROBJ('CHPOINT',MCHPO2,1,IRETOU) 27 : IF(IERR.NE.0) RETURN 27 : CALL LIROBJ('CHPOINT',MCHPO3,1,IRETOU) 28 : IF(IERR.NE.0) RETURN 29 : C 30 : IF(IDET.EQ.2) THEN 30 : CALL LIROBJ('CHPOINT',MCHPO4,1,IRETOU) 31 : IF(IERR.NE.0) RETURN 32 : ENDIF 33 : C 33 : CALL ACTIPO(FLOT,IDET,MCHPO1,MCHPO2,MCHPO3,MCHPO4) 34 : C 35 : RETURN 36 : ENDIF 37 : C 38 : C CAS DES MCHAMLS - ACCELERATION GEOMETRIQUE UNIQUEMENT 39 : C 40 : IF(IDET.NE.1) THEN 40 : CALL ERREUR(21) 41 : RETURN 42 : ENDIF 43 : C 43 : CALL LIROBJ('MCHAML',IPCHE1,1,IRETOU) 44 : IF(IERR.NE.0)RETURN 44 : CALL LIROBJ('MCHAML',IPCHE2,1,IRETOU) 45 : IF(IERR.NE.0) RETURN 45 : CALL LIROBJ('MCHAML',IPCHE3,1,IRETOU) 46 : IF(IERR.NE.0)RETURN 47 : MOCOMP=' ' 47 : CALL LIRCHA(MOCOMP,0,IRETOU) 48 : IF(IERR.NE.0) RETURN 48 : CALL ACTICH(FLOT,IPCHE1,IPCHE2,IPCHE3,MOCOMP,IPCHE4) 49 : C 49 : IF (IERR.NE.0) CALL ECROBJ('MCHAML',IPCHE4) 50 : C 51 : RETURN 52 : END 53 : 54 :
© Cast3M 2003 - All rights reserved.
Disclaimer