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

Manipulation d'une S.D. B issue d'un problème évolutif

La résolution d'un problème parabolique par une méthode à pas variables conduit à calculer les solutions en des temps a priori inconnus. Le but du module INTERQ est donc d'interpoler ces solutions pour obtenir une structure B contenant des solutions en des temps choisis.

      SUBROUTINE INTERQ(M,XM,NQ,NFB,NIB,PAS,NPAS,T0,TSAV,NFBS,NIBS)
C  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C   SOUS PROGRAMME INTERQ  . SIMPLE PRECISION
C   BUT  : INTERPOLER A L ORDRE NQ LES SOLUTIONS ISSUES DE EVO*GE(RK)
C  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

où

• M et XM représentent le super-tableau,
• NQ est l'ordre de l'interpolation, il code également le choix du mode de sauvegarder des solutions:
  • 1 NQ 10 : interpolation à l'ordre NQ et sauvegarde de la solution au seul temps TSAV,
  • 11 NQ 20 : interpolation à l'ordre NQ - 10 et sauvegarde des solutions à tous les temps,
• NFB, NIB sont le numéro de fichier et le niveau de la structure B issue du module évolutif,
• PAS est le pas de temps,
• NPAS est le nombre de pas de temps (solution initiale comprise),
• T0 est le temps initial,
• TSAV est inutile si NQ > 10 ou est le temps où l'on désire sauver la solution. Ce temps est nécessairement l'une des valeurs atteintes par la suite définie par T0, PAS et NPAS:

  TSAV = T0 + j * PAS pour j=1,NPAS

• NFBS, NIBS sont le numéro de fichier et le niveau de la structure B de sortie.

La résolution d'un problème parabolique par une méthode quelconque peut, dans certains cas, conduire à calculer les solutions dans deux intervalles de temps , et , se chevauchant. Le but du module UNIONB est donc de former la solution dans l'intervalle , en supprimant l'éventuelle intersection.

Plusieurs options sont offertes. Pour les préciser donnons l'exemple suivant:

• la première structure B contient les solutions aux temps = 0, = 10, ... , = 100 et a un tableau associé contenant les valeurs des
• la seconde structure contient les solutions aux temps = 80, = 90, ... , = 200 et a un tableau associé contenant les valeurs des

• la structure concaténée contiendra les solutions aux temps: = 0, = 10, ... , = 80, = 90, ... = 200 si on choisit l'option de troncature dans l'ordre croissant. Son tableau associé contiendra les valeurs de tous ces temps.
• la structure concaténée contiendra les solutions aux temps: = 0, = 10, ... , = 70, = 80, ... = 200 si on choisit l'option de troncature dans l'ordre décroissant. Son tableau associé contiendra les valeurs de tous ces temps.
• la structure concaténée contiendra les solutions mises bout à bout si on choisit cette option.

Le préprocesseur UNIOXX appelle le module UNIONB qui permet cette concaténation.

      SUBROUTINE UNIONB(M,XM,NFB1,NIB1,NFB2,NIB2,JOPTIO,NFBS,NIBS)
C  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C   SOUS PROGRAMME UNIONB  . SIMPLE PRECISION
C   BUT  : RECOLLER 2  S.D.E. B EN UNE SEULE S.D.S. B
C  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

où

• M et XM représentent le super-tableau,
• NFB1, NIB1 sont le numéro de fichier et le niveau de la première structure B,
• NFB2, NIB2 sont le numéro de fichier et le niveau de la seconde structure B,
• JOPTIO est l'option de troncature

  .

  0 : mise bout à bout

  .

  1 : ordre croissant selon le dernier tableau associé,

  .

  2 : ordre décroissant selon le dernier tableau associé,

• NFBS, NIBS sont le numéro de fichier et le niveau de la structure B de sortie.