Télécharger la plaquette Cast3M


fr en

Télécharger gcfs2.eso

Retour à la liste

Numérotation des lignes : 

gcfs2

  1. C GCFS2     SOURCE    GOUNAND   21/06/02    21:16:02     11022          
  2.       SUBROUTINE GCFS2(PGCOUR,NOPG,DIMSRF,XCOR,POIDS,
  3.      $     IMPR,IRET)
  4.       IMPLICIT REAL*8 (A-H,O-Z)
  5.       IMPLICIT INTEGER (I-N)
  6. C***********************************************************************
  7. C NOM         : GCFS2
  8. C PROJET      : Noyau linéaire NLIN
  9. C DESCRIPTION : Rajoute des points à une méthode d'intégration type
  10. C               Gauss (PGCOUR).
  11. C               Domaine de type [ Fully symmetric region ]
  12. C                       i.e. segment, carré et cube pour nous.
  13. C               Générateur de type [ Fully symmetric ].
  14. C               Rule structure of index 2 i.e.
  15. C               XCOR=(a,0,...,0)
  16. C
  17. C               le nombre de points générés est 2*DIMSRF :
  18. C               (a,0,...,0), (-a,0,...,0)
  19. C               (0,a,...,0), (0,-a,...,0) et les autres permutations
  20. C
  21. C LANGAGE     : ESOPE
  22. C AUTEUR      : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
  23. C               mél : gounand@semt2.smts.cea.fr
  24. C***********************************************************************
  25. C APPELES          : -
  26. C APPELE PAR       : INGACU
  27. C***********************************************************************
  28. C ENTREES            : DIMSRF, XCOR, POIDS
  29. C ENTREES/SORTIES    : PGCOUR (actif en *MOD), NOPG
  30. C SORTIES            : MYPGS
  31. C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  32. C***********************************************************************
  33. C VERSION    : v1, 20/10/99, version initiale
  34. C HISTORIQUE : v1, 20/10/99, création
  35. C HISTORIQUE :
  36. C HISTORIQUE :
  37. C***********************************************************************
  38. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  39. C en cas de modification de ce sous-programme afin de faciliter
  40. C la maintenance !
  41. C***********************************************************************
  42.  
  43. -INC PPARAM
  44. -INC CCOPTIO
  45. -INC TNLIN      
  46. *-INC SPOGAU
  47.       POINTEUR PGCOUR.POGAU
  48. *
  49.       INTEGER NOPG,DIMSRF
  50.       REAL*8 XCOR(DIMSRF)
  51.       REAL*8 POIDS
  52.       INTEGER IMPR,IRET
  53. *
  54.       INTEGER IPERM,ISIGN,ICOOP
  55. *
  56. * Executable statements
  57. *
  58.       IF (IMPR.GT.6) WRITE(IOIMP,*) 'Entrée dans gcfs2'
  59. *
  60. * Boucle sur les permutations
  61.       DO 1 IPERM=1,DIMSRF
  62. * Boucle sur le signe
  63.          DO 12 ISIGN=1,2
  64.             NOPG=NOPG+1
  65. * Boucle sur les coordonnées du point
  66.             DO 122 ICOOP=1,DIMSRF
  67.                IF (ICOOP.EQ.IPERM) THEN
  68.                   PGCOUR.XCOPG(ICOOP,NOPG)=
  69.      $                 DBLE((-1)**ISIGN)
  70.      $                 *XCOR(1)
  71.                ELSE
  72.                   PGCOUR.XCOPG(ICOOP,NOPG)=0.D0
  73.                ENDIF
  74.  122        CONTINUE
  75.             PGCOUR.XPOPG(NOPG)=POIDS
  76.  12      CONTINUE
  77.  1    CONTINUE
  78. *
  79. * Normal termination
  80. *
  81.       IRET=0
  82.       RETURN
  83. *
  84. * Format handling
  85. *
  86. *
  87. * Error handling
  88. *
  89.  9999 CONTINUE
  90.       IRET=1
  91.       WRITE(IOIMP,*) 'An error was detected in subroutine gcfs2'
  92.       RETURN
  93. *
  94. * End of subroutine GCFS2
  95. *
  96.       END
  97.  
  98.  
  99.  
  100.  
  101.  

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