[updated 3.May.2010]

GEOMETRIE !impression de niveau 10
LECTURE SUR FICHIER 
'/Library/Melina/maillages/t0101/16x16'
%
! impression de niveau 10
%
% Description de l'inconnue
%
 INCONNUE 'u' de type SCALAIRE INTERPOLATION de LAGRANGE
 (AUX ABSCISSES DE GAUSS-LOBATTO) P3
%
% Formule de quadrature adaptee (2k+1 pour Qk)
%
QUADRATURE DE GAUSS DE DEGRE 5
%
% Description des termes
%
CALCUL sur le DOMAINE 'Omega' !impression de niveau 10
!  Matrice de rigidite
  TERME ELEMENT FINI 'Matris' NIVEAU 0
!     impression de niveau 10
      INCONNUE 'u'
      INTEGRAND 'GRAdGRAD'
!
! Deux possibilite de calcul du second membre de l'equation
! voir DONNEE 'Calcul_2membre'
!
! 1) Interpolation de la donnee F de -Delta u = F
!   Matrice de masse pour le calcul du second membre
!   du pb approche par interpolation de la donnee de second membre
!   et produit 'Masse' x '-deltau'
  TERME ELEMENT FINI 'Masse' NIVEAU 0
! impression de niveau 10
       INCONNUE 'u'
       INTEGRAND 'UV'
  TERME VALEURS NODALES 'Fval' NIVEAU 0
       INCONNUE 'u'
       DONNEE '-deltau' FONCTION
! 2) Integration numerique
!    Second membre (de surface) calcule par integration numerique
   TERME ELEMENT FINI 'SecmbF' NIVEAU 0
! impression de niveau 10
       INCONNUE 'u'
       INTEGRAND 'V'
       DONNEE '-deltau' FONCTION
!
! Calcul de la solution exacte
!
  TERME VALEURS NODALES 'SOLEXA' NIVEAU 0
! impression de niveau 10
      INCONNUE 'u'
      DONNEE 'uexacte' FONCTION
!
! Le bord Gamma_i porte soit une condition de Dirichlet
! soit une condition de Neuman (voir DONNEE 'CLi' plus bas)
!
CALCUL sur le DOMAINE 'Gamma1' 
! impression de niveau 10
  CONDITION ESSENTIELLE 'Dirich' NIVEAU 1 de type 'u=g'
      DONNEE 'uexacte'
!
! Deux possibilite de calcul du second membre de Neuman
! voir DONNEE 'Calcul_2membre'
!
! 1) Interpolation de la donnee g de du/dn = g
*!   Matrice de masse de bord pour du calcul du second membre
*!   du pb approche par interpolation de la donnee de Neuman
*!   et produit 'Masse' x 'du/dn'
  TERME ELEMENT FINI 'Masse' NIVEAU 1
! impression de niveau 10
      INCONNUE 'u'
      INTEGRAND 'UV'
  TERME VALEURS NODALES 'dudn' NIVEAU 1 DE TYPE 'N.F'
      INCONNUE 'u'
      DONNEE 'gradu' FONCTION REELLE
! 2) Integration numerique
   TERME ELEMENT FINI 'SecmbN' NIVEAU 1
! impression de niveau 10
      INCONNUE 'u' 
      INTEGRAND 'V' 
      DONNEE 'dudn' FONCTION REELLE
!
CALCUL sur le DOMAINE 'Gamma2' 
!impression niveau 10
  CONDITION ESSENTIELLE 'Dirich' NIVEAU 2 de type 'u=g'
      DONNEE 'uexacte'
  TERME ELEMENT FINI 'SecmbN' NIVEAU 2
      INCONNUE 'u' INTEGRAND 'V' DONNEE 'dudn' FONCTION REELLE
!
CALCUL sur le DOMAINE 'Gamma3'
! impression niveau 10
  CONDITION ESSENTIELLE 'Dirich' NIVEAU 3 de type 'u=g'
      DONNEE 'uexacte'
  TERME ELEMENT FINI 'SecmbN' NIVEAU 3
      INCONNUE 'u' INTEGRAND 'V' DONNEE 'dudn' FONCTION REELLE
!
CALCUL sur le DOMAINE 'Gamma4' 
!impression niveau 10
  CONDITION ESSENTIELLE 'Dirich' NIVEAU 4 de type 'u=g'
      DONNEE 'uexacte'
  TERME ELEMENT FINI 'SecmbN' NIVEAU 4
      INCONNUE 'u' INTEGRAND 'V' DONNEE 'dudn' FONCTION REELLE
FIN
!impression de niveau 10
DONNEE 'Domaines' CONSTANTE ENTIERE 4
!    Condition de 'Dirichlet' ou 'NeumanNH' sur le Bord i
       'CL1' CONSTANTE CARACTERE 'NeumanNH'
       'CL2' CONSTANTE CARACTERE 'Dirichlet'
       'CL3' CONSTANTE CARACTERE 'Dirichlet'
       'CL4' CONSTANTE CARACTERE 'Dirichlet'
*       'Calcul_2membre' CONSTANTE CARACTERE 'Integration'
       'Calcul_2membre' CONSTANTE CARACTERE 'Interpolation'
       'Resolution' CONSTANTE CARACTERE 'Choleski'
       'Solution_exacte' CONSTANTE CARACTERE 'SOLEXA'
       'Imp_ass' CONSTANTE ENTIERE 0
       'Imp_sys' CONSTANTE ENTIERE 0
       'Imp_res' CONSTANTE ENTIERE 0
       'Imp_sol' CONSTANTE ENTIERE 0
FIN