Télécharger la plaquette Cast3M


fr en

Télécharger gpocrq.eso

Retour à la liste

Numérotation des lignes : 

gpocrq

  1. C GPOCRQ    SOURCE    GOUNAND   21/06/02    21:16:17     11022          
  2.       SUBROUTINE GPOCRQ(NDIM,
  3.      $     MYBPOL,
  4.      $     IMPR,IRET)
  5.       IMPLICIT REAL*8 (A-H,O-Z)
  6.       IMPLICIT INTEGER (I-N)
  7. C***********************************************************************
  8. C NOM         : GPOCRQ
  9. C DESCRIPTION : Génère les polynômes pour les éléments de
  10. C               Crouzeix-Raviart P1 non conforme de dimension
  11. C               NDIM et l'ajoute à la liste de polynômes.
  12. C
  13. C LANGAGE     : ESOPE
  14. C AUTEUR      : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
  15. C               mél : gounand@semt2.smts.cea.fr
  16. C***********************************************************************
  17. C APPELES          : ISET
  18. C APPELE PAR       : INELQU, INELCU
  19. C***********************************************************************
  20. C ENTREES            : NDIM
  21. C ENTREES/SORTIES    : MYBPOL
  22. C SORTIES            : -
  23. C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  24. C***********************************************************************
  25. C VERSION    : v1, 18/10/2006, version initiale
  26. C HISTORIQUE : v1, 18/10/2006, création
  27. C HISTORIQUE :
  28. C HISTORIQUE :
  29. C***********************************************************************
  30. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  31. C en cas de modification de ce sous-programme afin de faciliter
  32. C la maintenance !
  33. C***********************************************************************
  34.  
  35. -INC PPARAM
  36. -INC CCOPTIO
  37. -INC TNLIN      
  38. *-INC SPOLYNO
  39.       INTEGER NBMON,NDIML
  40.       POINTEUR MYBPOL.POLYNS
  41.       POINTEUR POCRQ.POLYNO
  42. *
  43.       INTEGER NDIM
  44.       INTEGER IMPR,IRET
  45. *
  46.       REAL*8 UN
  47.       PARAMETER(UN=1.D0)
  48.       INTEGER IMON
  49. *
  50. * Executable statements
  51. *
  52.       IF (IMPR.GT.2) WRITE(IOIMP,*) 'Entrée dans gpocrq.eso'
  53. * On génère un polynôme du type :
  54. * x_i^2 - x_{i+1}^2 avec i variant de 1 à NDIM-1
  55.       SEGACT MYBPOL*MOD
  56.       DO IBPOLY=1,NDIM-1
  57.          NDIML=NDIM
  58.          NBMON=2
  59.          SEGINI POCRQ
  60.          POCRQ.COEMON(1)=UN
  61.          POCRQ.EXPMON(IBPOLY,1)=2
  62.          POCRQ.COEMON(2)=-UN
  63.          POCRQ.EXPMON(IBPOLY+1,2)=2
  64.          SEGDES POCRQ
  65.          MYBPOL.LIPOLY(**)=POCRQ
  66.       ENDDO
  67.       SEGDES MYBPOL
  68. *
  69. * Normal termination
  70. *
  71.       IRET=0
  72.       RETURN
  73. *
  74. * Format handling
  75. *
  76. *
  77. * Error handling
  78. *
  79.  9999 CONTINUE
  80.       IRET=1
  81.       WRITE(IOIMP,*) 'An error was detected in subroutine gpocrq'
  82.       RETURN
  83. *
  84. * End of subroutine GPOCRQ
  85. *
  86.       END
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  

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