Télécharger la plaquette Cast3M


fr en

Télécharger active.eso

Retour à la liste

Numérotation des lignes : 

active

  1. C ACTIVE    SOURCE    CB215821  20/11/04    21:15:03     10766          
  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.  
  12. -INC PPARAM
  13. -INC CCOPTIO
  14. C
  15.       CHARACTER*(LOCOMP) MOCOMP
  16.       CHARACTER*4 LISMO(2)
  17.       DATA LISMO/'GEOM','SECA'/
  18. C
  19.       CALL LIRMOT(LISMO,2,IDET,1)
  20.       IF(IERR.NE.0) RETURN
  21. C
  22.       IF(IDET.EQ.1) THEN
  23.       FLOT=0.98D0
  24.       CALL LIRREE(FLOT,0,IRETOU)
  25.       ENDIF
  26. C
  27. C  CAS DES CHPOINTS   - ACCELERATIONS GEOMETRIQUE OU SECANTE
  28. C
  29.       CALL LIROBJ('CHPOINT ',MCHPO1,0,IRETOU)
  30.       IF(IERR.NE.0) RETURN
  31.       IF(IRETOU.EQ.1) THEN
  32.         CALL ACTOBJ('CHPOINT ',MCHPO1,1)
  33.         CALL LIROBJ('CHPOINT ',MCHPO2,1,IRETOU)
  34.         CALL ACTOBJ('CHPOINT ',MCHPO2,1)
  35.         IF(IERR.NE.0) RETURN
  36.         CALL LIROBJ('CHPOINT ',MCHPO3,1,IRETOU)
  37.         CALL ACTOBJ('CHPOINT ',MCHPO3,1)
  38.         IF(IERR.NE.0) RETURN
  39. C
  40.         IF(IDET.EQ.2) THEN
  41.           CALL LIROBJ('CHPOINT ',MCHPO4,1,IRETOU)
  42.           CALL ACTOBJ('CHPOINT ',MCHPO4,1)
  43.           IF(IERR.NE.0) RETURN
  44.         ENDIF
  45. C
  46.         CALL ACTIPO(FLOT,IDET,MCHPO1,MCHPO2,MCHPO3,MCHPO4)
  47. C
  48.         RETURN
  49.       ENDIF
  50. C
  51. C  CAS DES MCHAMLS   - ACCELERATION GEOMETRIQUE UNIQUEMENT
  52. C
  53.       IF(IDET.NE.1) THEN
  54.         CALL ERREUR(21)
  55.         RETURN
  56.       ENDIF
  57. C
  58.       CALL LIROBJ('MCHAML  ',IPCHE1,1,IRETOU)
  59.       CALL ACTOBJ('MCHAML  ',IPCHE1,1)
  60.       IF(IERR.NE.0)RETURN
  61.       CALL LIROBJ('MCHAML  ',IPCHE2,1,IRETOU)
  62.       CALL ACTOBJ('MCHAML  ',IPCHE2,1)
  63.       IF(IERR.NE.0) RETURN
  64.       CALL LIROBJ('MCHAML  ',IPCHE3,1,IRETOU)
  65.       CALL ACTOBJ('MCHAML  ',IPCHE3,1)
  66.       IF(IERR.NE.0)RETURN
  67.       MOCOMP=' '
  68.       CALL LIRCHA(MOCOMP,0,IRETOU)
  69.       IF(IERR.NE.0) RETURN
  70.       CALL ACTICH(FLOT,IPCHE1,IPCHE2,IPCHE3,MOCOMP,IPCHE4)
  71. C
  72.       IF (IERR.NE.0) THEN
  73.         CALL ACTOBJ('MCHAML  ',IPCHE4,1)
  74.         CALL ECROBJ('MCHAML  ',IPCHE4)
  75.       ENDIF
  76.  
  77.       END
  78.  
  79.  
  80.  
  81.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales