Suiv.: Symétrisation et recollement d'une S.D. Sup.: 2 Manipulation de solutions Préc.: 2 Manipulation de solutions Index Table des matières

2.1 Construction manuelle d'une S.D. B

Le préprocesseur COSDXX permet de construire une S.D. B à partir :

• d'une valeur donnée : tous les degrés de liberté seront affectés à cette valeur.
• d'une fonction dépendant des coordonnées des noeuds avec:
  • utilisation d'une S.D. NOPO,
  • utilisation des S.D. MAIL et COOR, en particulier dans le cas où le nombre de degrés de liberté par noeud, ND, est non constant.

Dans le cas où on affecte une valeur donnée (via un vecteur), le module COSDB est appelé. Pour affecter une valeur dépendant des coordonnées du noeud, on utilise le module COSNOB qui, pour ND constant, exploite la S.D. NOPO ou on utilise le module COSMCB qui, pour ND quelconque, exploite les S.D. MAIL et COOR.

Le module COSDB correspond au cas où l'on donne un vecteur (U) ou deux vecteurs (U et V) pour construire la S.D. B telle que:

B4(I) = U(I)

ou, pour créer une structure à deux tableaux B4 (deux vecteurs sont alors utilisés):

B4(I) = U(I) et B4(I) = V(I)

      SUBROUTINE COSDB(M,NVEC,U,V,LUV,ND,NFB,NIB,NTAB)
C  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C  UTILITAIRE COSDB : CREER UNE SD B A PARTIR DE UN OU DEUX VECTEURS
C  ----------------
C  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

où

• M est le super-tableau,
• NVEC est le nombre de vecteurs (1 ou 2) indiquant que l'on veut un ou deux tableaux B4. Si NVEC > 0, la structure sera en simple précision, pour NVEC négatif, elle sera en double précision,
• U (V) sont les vecteurs dûment remplis (en fonction du type) contenant les valeurs à affecter,
• LUV est le nombre de mots de chaque vecteur,
• NFB, NIB et NTAB sont le numéro de fichier, le niveau et le nombre de tableaux à associer de la S.D.S. B.

Si NTAB est non nul, on doit donner :

Boucle I de 1 à NTAB

• nom du tableau I sur 4 caractères,
• type de ce tableau (1 = entier, 2 = réel simple précision, 4 = character, 5 = réel double précision, ...)
• nombre de mots de ce tableau
• les variables de ce tableau
• commentaire sur le contenu du tableau I sur 18 mots de 4 caractères.

Le module COSNOB correspond au cas où l'on calcule un vecteur (U) ou deux vecteurs (U et V) pour construire la S.D. B telle que:

B4(I) = VALCL*(I,X,Y,Z)

ou, pour créer une structure à deux tableaux B4 (deux vecteurs sont alors utilisés):

B4(I) = VALCL*(I,X,Y,Z) et B4(I) = VALCL*(I,X,Y,Z)

Les fonctions VALCL* sont en fait VALCLR en simple précision et VALCLD en double précision. Leurs arguments sont I (de 1 à ND, supposé constant, cf. plus bas) le numéro du degré de liberté, X, Y et Z les coordonnées du point traité lues dans le tableau NOP4 de la structure NOPO. Ces fonctions seront données:

• sous la forme d'une FUNCTION dans le cas de l'appel classique du module (FONINT = .FALSE.),
• sous forme interprétée ou sous la forme d'une FUNCTION (au choix) dans le cas d'un appel via COSDXX.

      SUBROUTINE COSNOB(M,XM,NFNOPO,NINOPO,NVEC,U,V,DU,DV,ND,
     +                  NFB,NIB,NTAB,FONINT,VALCLR,VALCLD)
C  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C  UTILITAIRE COSNOB : CREER UNE SD B A PARTIR DE NOPO ET D'UNE FONCTION
C  ----------------   ( MEME BUT QUE COSDB )
C  U(I) ( V(I) ) = FONCTION(X(I),Y(I),Z(I)) AU POINT I DE LA SD NOPO
C  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

où

• M et XM représentent le super-tableau,
• NFNOPO, NINOPO sont le numéro de fichier et le niveau de la structure NOPO,
• NVEC est le nombre de vecteurs (1 ou 2) indiquant que l'on veut un ou deux tableaux B4. Si NVEC > 0, la structure sera en simple précision, pour NVEC négatif, elle sera en double précision,
• U ou DU (V ou DV) sont les vecteurs qui seront utilisés (selon le type) avec les valeurs à affecter,
• ND est le nombre (supposé constant) de degrés de liberté en chaque noeud,
• NFB, NIB et NTAB sont le numéro de fichier, le niveau et le nombre de tableaux à associer de la S.D.S. B,
• FONINT est un logique valant .TRUE. si on utilise les fonctions VALCLR ou VALCLD sous forme interprétée, .FALSE. si on les utilise normalement.

Si NTAB est non nul, on doit donner :

Boucle I de 1 à NTAB

• nom du tableau I sur 4 caractères,
• type de ce tableau (1 = entier, 2 = réel simple précision, 4 = character, 5 = réel double précision, ...)
• nombre de mots de ce tableau
• les variables de ce tableau
• commentaire sur le contenu du tableau I sur 18 mots de 4 caractères.

Selon la valeur de FONINT, on doit écrire les fonctions VALCL* ou non. VALCLR et VALCLD ont la forme suivante:

      FUNCTION VALCLR(I,X,Y,Z)

      DOUBLE PRECISION FUNCTION VALCLD(I,X,Y,Z)

Notons par ailleurs que le nombre de variables des vecteurs est calculé comme ND x NP où NP est le nombre de points du maillage.

Le module COSMCB correspond au cas où l'on calcule un vecteur (U) ou deux vecteurs (U et V) pour construire la S.D. B telle que:

B4(I) = VALCL*(I,X,Y,Z)

ou, pour créer une structure à deux tableaux B4 (deux vecteurs sont alors utilisés):

B4(I) = VALCL*(I,X,Y,Z) et B4(I) = VALCL*(I,X,Y,Z)

Les fonctions VALCL* sont en fait VALCLR en simple précision et VALCLD en double précision. Leurs arguments sont I (de 1 à ND, cf. plus bas) le numéro du degré de liberté, X, Y et Z les coordonnées du point traité lues dans le tableau COO4 de la structure COOR associée à la structure MAIL. Ces fonctions seront données:

• sous la forme d'une FUNCTION dans le cas de l'appel classique du module (FONINT = .FALSE.),
• sous forme interprétée ou sous la forme d'une FUNCTION dans le cas d'un appel via COSDXX.

      SUBROUTINE COSMCB (M,XM,DM,NFMAIL,NIMAIL,NFCOOR,NICOOR,NVEC,
     +   NFB,NIB,NTAB,FONINT,VALCLR,VALCLD,U,V,DU,DV)
C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C  UTILITAIRE COSMCB : CREER UNE SD B A PARTIR DE MAIL ET COOR ET
C  ----------------    D'UNE FONCTION ( MEME BUT QUE COSDB )
C  U(I) ( V(I) ) = FONCTION(X(I),Y(I),Z(I)) AU POINT I DE LA SD COOR
C  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

où

• M, XM et DM représentent le super-tableau,
• NFMAIL, NIMAIL sont le numéro de fichier et le niveau de la structure MAIL,
• NFCOOR, NICOOR sont le numéro de fichier et le niveau de la structure COOR,
• NVEC est le nombre de vecteurs (1 ou 2) indiquant que l'on veut un ou deux tableaux B4. Si NVEC > 0, la structure sera en simple précision, pour NVEC négatif, elle sera en double précision,
• U ou DU (V ou DV) sont les vecteurs qui seront utilisés (selon le type) avec les valeurs à affecter,
• NFB, NIB et NTAB sont le numéro de fichier, le niveau et le nombre de tableaux à associer de la S.D.S. B,
• FONINT est un logique valant .TRUE. si on utilise les fonctions VALCLR ou VALCLD sous forme interprétée, .FALSE. si on les utilise normalement.

Si NTAB est non nul, on doit donner :

Boucle I de 1 à NTAB

• nom du tableau I sur 4 caractères,
• type de ce tableau (1 = entier, 2 = réel simple précision, 4 = character, 5 = réel double précision, ...)
• nombre de mots de ce tableau
• les variables de ce tableau
• commentaire sur le contenu du tableau I sur 18 mots de 4 caractères.

Selon la valeur de FONINT, on doit écrire les fonctions VALCL* ou non. Rappelons que VALCLR et VALCLD ont la forme suivante:

      FUNCTION VALCLR(I,X,Y,Z)

      DOUBLE PRECISION FUNCTION VALCLD(I,X,Y,Z)