Scatter

Scatter#

#include <Scatter.h>

How to cite this class in this doc#

RST

:ref:`Scatter <class-scatter>`

Markdown

[Scatter](#class-scatter)

Detailed description#

Inherits from#

public : Interprete

Inherited by#

public : ScatterMED

List of Public Methods#

domaine
duplique
get_info
interpreter
lire_domaine
Scatter
taille_memoire
~Scatter

List of Protected Methods#

printOn
read_domain_no_comm
readOn

Complete Member Function Documentation#

domaine

Domaine & Scatter::domaine()

Renvoi le domaine associe.

duplique

int Scatter::duplique() const override

Reimplements:

duplique

get_info

const Type_info * Scatter::get_info() const override

Donne des informations sur le type de l’ Objet_U .

Reimplements:

get_info

interpreter

Entree & Scatter::interpreter(Entree &) override

Lit et complete un domaine parallele selon les motcles lus dans le jeu de donnees.

Format: Scatter [debug] file_name domain_name On lit les sommets, les elements et les sommets et faces de joint, On construit les espaces distants et virtuels en fonction de l’epaisseur de joint.

Reimplements:

interpreter

lire_domaine

void Scatter::lire_domaine(Nom &, Noms &liste_bords_periodiques)

Lit le domaine dans le fichier de nom “nomentree”, de type LecFicDistribueBin ou LecFicDistribue .

Format attendu : Domaine::ReadOn

Scatter

Scatter::Scatter()

taille_memoire

unsigned Scatter::taille_memoire() const override

Reimplements:

taille_memoire

~Scatter

Scatter::~Scatter()

ajouter_joints

void Scatter::ajouter_joints(Domaine &domaine, ArrOfInt &pe_voisins)

Ajoute des joints avec tous les pe de pe_voisins.

Pour que l’ensemble des joints soit symetrique, on en cree aussi un joint sur le processeur destination: Si A ajoute un joint avec B, alors B ajoute un joint avec A. On trie les joints par ordre croissant du numero de PE. ATTENTION: les joints sont donc reordonnes ! On met dans pe_voisins la liste des joints effectivement crees.

calculer_espace_distant

void Scatter::calculer_espace_distant(Domaine &domaine, const int nb_items_reels, const ArrsOfInt &items_to_send, const JOINT_ITEM type_item)

Determination des items distants en fonction d’une liste d’items a envoyer et de listes d’items communs.

exemple: calculer_espace_distant_sommets calculer_espace_distant_faces Pour les sommets: les “items_to_send” sont les sommets des elements distants, Si le processeur A veut que le processeur B connaisse le sommet i, il faut que le processeur qui possede le sommet l’envoie a B. Le processeur qui “possede” le sommet est le plus petit parmi les PEs qui partagent ce sommet (item commun) (requis pour pouvoir faire echange_item_commun et echange_espace_virtuel en une seule passe). De plus, si plusieurs processeurs demandent a envoyer le meme sommet au meme processeur, il ne faut l’inserer qu’une seule fois dans l’espace distant.

calculer_espace_distant_aretes

void Scatter::calculer_espace_distant_aretes(Domaine &domaine, const int nb_aretes_reelles, const IntTab &elem_aretes)

Idem que Scatter::calculer_espace_distant_sommets pour les aretes.

calculer_espace_distant_elements

void Scatter::calculer_espace_distant_elements(Domaine &dom)

Remplissage du tableau “espace_distant()” des elements dans les joints.

C’est ici qu’on determine les elements de joint en fonction de l’epaisseur de joint. Le tableau espace_distant contient les indices locaux des elements distants (a envoyer aux processeurs voisins) Pour un joint d’epaisseur 1, ce sont tous les elements voisins d’un sommet de joint (sommet sur un face de joint ou sommet isole). Pour un joint d’epaisseur n>1, ce sont tous les elements voisins d’un sommet d’un element du joint d’epaisseur n - 1. Le voisinage s’entend sur le domaine global (toutes domaines confondues) Historique: premiere version B.Mathieu le 16/01/2007. Il existe une methode qui determine les elements distants au moment du decoupage (DomaineCutter::construire_elements_distants_ssdom). La methode ci - dessous a ete validee par comparaison avec la methode du decoupeur. Les sorties ont ete verifiees pour des epaisseurs jusqu’a 5 sur des maillages tetra. La difficulte de l’algorithme est d’obtenir les elements virtuels d’epaisseur> 1 qui se trouvent sur des sous - domaines qui ne sont pas en contact direct avec le sous - domaine local. Difficulte resolue par l’algorithme ci - dessous.

calculer_espace_distant_faces

void Scatter::calculer_espace_distant_faces(Domaine &domaine, const int nb_faces_reelles, const IntTab &elem_faces)

Idem que Scatter::calculer_espace_distant_sommets pour les faces.

calculer_espace_distant_sommets

void Scatter::calculer_espace_distant_sommets(Domaine &dom, const Noms &liste_bords_periodiques)

En fonction de l’espace distant des elements, calcule l’espace distant des sommets.

Pour chaque joint, on envoie au processeur voisin l’ensemble des sommets de tous les elements du joint. C’est le processeur proprietaire du sommet (plus petit pe qui le possede) qui le met dans son espace distant. Attention, on cree de nouveaux joints. On remplit les tableaux dom.faces_joint(i).joint_item(JOINT_ITEM::SOMMET).items_distants();

calculer_nb_items_virtuels

void Scatter::calculer_nb_items_virtuels(Joints &joints, const JOINT_ITEM type_item)

Pour un item geometrique “type_item”, remplit le champ nb_items_virtuels_ des joints en fonction du nombre d’items distants :

Le nombre d’items virtuels sur un joint i du processeur j est le nombre d’items distants du joint j sur le processeur i.

calculer_renum_items_communs

void Scatter::calculer_renum_items_communs(Joints &joints, const JOINT_ITEM type_item)

On suppose que chaque joint[i].

joint_item(type_item).items_communs() contient les indices locaux des items de joint communs dans le meme ordre sur les deux processeurs (local et voisin) On remplit renum_items_communs : colonne 0 = contenu du tableau items_communs sur le PE voisin colonne 1 = contenu du tableau items_communs sur le PE local

check_consistancy_remote_items

void Scatter::check_consistancy_remote_items(Domaine &dom, const ArrOfInt &mergedDomaines)

Merged domaines receive joints information from their neighbours to ensure that their common items (vertices) appear in the same order.

If it’s not the case, the merged domaine reorders its common items so that it matches the neighbour’s order When 2 neighbouring domaines have each been merged, only the processor with the lowest rank proceeds to reordering

Chercher_Correspondance

int Scatter::Chercher_Correspondance(const DoubleTab &sommets1, const DoubleTab &sommets2, ArrOfInt &correspondance, const double epsilon)

Construit le tableau “correspondance” tel que Pour 0 < = i <sommets2.

size_array(), Si sommet2(i) existe dans le tableau sommet1, alors sommets2(i, …) = = sommets1(correspondance[i], …) Sinon correspondance[i] = - 1 L’egalite est verifiee a epsilon pres en absolu (soit abs(x1 - x2)<epsilon) L’algorithme est generalement en n1*log(n1) + n2*log(n1) (recherche basee sur un quicksort). En cas d’echec du tri, on utilise un algorithme en n1*n2. Les tableaux sommets1 et sommets2 doivent etre de dimension 2 Le tableau correspondance doit etre de taille sommets2.size_array(). Valeur de retour: nombre de sommets de sommets2 non trouves dans le tableau sommets1.

construire_correspondance_aretes_par_coordonnees

void Scatter::construire_correspondance_aretes_par_coordonnees(Domaine_VF &zvf)

Construction des tableaux joint_item(JOINT_ITEM::ARETE).

items_communs de tous les joints du domaine

construire_correspondance_items_par_coordonnees

void Scatter::construire_correspondance_items_par_coordonnees(Joints &joints, const JOINT_ITEM type_item, const DoubleTab &coord_items)

Methode obsolete (utilisee avec l’ancien decoupeur).

construire_correspondance_sommets_par_coordonnees

void Scatter::construire_correspondance_sommets_par_coordonnees(Domaine &dom)

Construction des tableaux joint_item(JOINT_ITEM::SOMMET).

items_communs de tous les joints du domaine(0) du domaine dom

construire_espace_virtuel_traduction

void Scatter::construire_espace_virtuel_traduction(const MD_Vector &md_indice, const MD_Vector &md_valeur, IntTab &tableau, const int error_is_fatal=1)

Construit la structure items_communs + espaces virtuels d’un tableau contenant des indices d’items geometriques, indexe par un autre type d’item geometrique.

Exemple: tableau indexe par md_indice, contenant des indices md_valeur: type_indice type_valeur exemple de tableau: element sommet domaine.les_elems() face sommet faces_sommets element face elem_faces face element faces_voisins element element ? element arete elem_aretes Nb_valeurs_max est le nombre d’items reels de type “type_valeur”.

construire_md_vector

void Scatter::construire_md_vector(const Domaine &, int nb_items_reels, const JOINT_ITEM, MD_Vector &)

construction d’un MD_Vector_std a partir des informations de joint du domaine pour le type d’item demande.

construire_structures_paralleles

void Scatter::construire_structures_paralleles(Domaine &dom, const Noms &liste_bords_perio)

Construction des structures paralleles du domaine et du domaine (determination des elements distants en fonction de l’epaisseur de joint,.

determination des sommets distants, creation des sommets et des elements virtuels)

corriger_espace_distant_elements_perio

void Scatter::corriger_espace_distant_elements_perio(Domaine &dom, const Noms &liste_bords_periodiques)

Les algorithmes actuels pour le periodique (assembleur P1B, OpDivElem P1B) ont besoin que pour chaque face virtuelle periodique, la face opposee soit.

aussi virtuelle. Ceci n’est pas assure a la sortie de la methode calculer_elements_distants. Cette methode ajoute aux elements distants les elements manquants pour assurer cette condition: Si un element est distant pour un PE donne est voisin d’une face periodique, on ajoute a l’espace distant l’element adjacent a la face opposee.

cree_instance

Objet_U * Scatter::cree_instance()

info

const Type_info * Scatter::info()

init_sequential_domain

template void Scatter::init_sequential_domain(Domaine_32_64< _SIZE_ > &dom)

Create parallel descriptors for the vertex and element arrays of the domain (necessary because Scatter is never invoked in sequential).

In 64bit the corresponding number of items might be big. This is here the main justification for the need of the class MD_Vector_seq which unique useful argument is the total number of items (with type trustIdType). Alternative would have been to make all members of MD_Vector_std compatible with trustIdType …

reordonner_faces_de_joint

void Scatter::reordonner_faces_de_joint(Domaine &dom)

Reordonne les faces de joint de sorte qu’elles apparaissent dans le meme ordre sur chaque couple de processeur voisin.

En pratique, pour un couple pe1 <pe2, pe1 envoie ses faces de joint a pe2 et pe2 les traduit en indices de sommets locaux. Les faces de joint du PE2 ne sont donc pas utilisees.

self_cast

const Scatter & Scatter::self_cast(const Objet_U &)

self_cast

Scatter & Scatter::self_cast(Objet_U &)

trier_les_joints

void Scatter::trier_les_joints(Joints &joints)

Tri des joints dans l’ordre croissant des processeurs.

uninit_sequential_domain

template void Scatter::uninit_sequential_domain(Domaine_32_64< _SIZE_ > &dom)

methode utilisee par les interpretes qui modifient le domaine (sequentiel), detruit les descripteurs des sommets et elements pour permettre la modification de ces tableaux.

printOn

Sortie & Scatter::printOn(Sortie &x) const override

Simple appel a: Interprete::printOn(Sortie&)

Reimplements:

printOn

read_domain_no_comm

void Scatter::read_domain_no_comm(Entree &fic)

Does the exact same thing as the readOn of the class Domaine but without collective communication.

Necessary when the processors don’t have the same numbers of file to read

readOn

Entree & Scatter::readOn(Entree &) override

Simple appel a: Interprete::readOn(Entree&)

Reimplements:

readOn

Attributes Documentation#

info_obj (public)

const Type_info Scatter::info_obj

le_domaine (protected)

TRUST_Ref<Domaine *> Scatter::le_domaine

Scatter

Contents

Scatter#

How to cite this class in this doc#

Detailed description#

Inherits from#

Inherited by#

List of Public Methods#

List of Public Static Methods#

List of Protected Methods#

Complete Member Function Documentation#

Attributes Documentation#