Parcours_interface

Parcours_interface#

#include <Parcours_interface.h>

How to cite this class in this doc#

RST

:ref:`Parcours\_interface <class-parcours-interface>`

Markdown

[Parcours\_interface](#class-parcours-interface)

Detailed description#

Inherits from#

public : Objet_U

List of Public Methods#

associer_connectivite_frontieres
associer_domaine_dis
calculer_face_sortie_element
calculer_normale_face_bord
calculer_sortie_face_bord
distance_sommet_faces
duplique
get_correction_parcours_thomas
get_erreur_geometrique
get_info
get_parcours_sans_tolerance
Parcours_interface
projeter_vecteur_sur_face
set_correction_parcours_thomas
set_parcours_sans_tolerance
taille_memoire
~Parcours_interface

List of Protected Static Methods#

calcul_inverse_matrice33
calcul_produit_matrice33_vecteur

Complete Member Function Documentation#

associer_connectivite_frontieres

void Parcours_interface::associer_connectivite_frontieres(const Connectivite_frontieres &connect)

associer_domaine_dis

void Parcours_interface::associer_domaine_dis(const Domaine_dis_base &domaine_dis)

Remplissage des variables persistantes de la classe (refdomaine_vf_, nb_faces_elem_, nb_elements_reels_, type_element_,.

equations_plans_faces_).

calculer_face_sortie_element

int Parcours_interface::calculer_face_sortie_element(const Domaine_VF &domaine_vf, const int num_element, double x0, double y0, double z0, double x1, double y1, double z1, double &pos_intersection) const

Pour un point P0 (x0, y0, z0) a l’INTERIEUR de l’element num_element et un autre point P1 (x1, y1, z1), calcule l’intersection du segment (P0,P1)

avec les bords de l’element. Si le point P1 est sur un bord de l’element (a epsilon pres), on considere qu’il est a l’interieur et on ne reporte aucune intersection. Si on trouve une intersection I, on met dans pos_intersection la coordonnee barycentrique de l’intersection definie par I = (1 - pos) * P0 + pos * P1 Si on ne trouve pas d’intersection, pos_intersection est inchange. Valeur de retour: Si une intersection a ete trouvee, numero de la face de sortie dans le domaine_vf (peut servir d’index dans face_voisins par exemple). Sinon, renvoie-1.

calculer_normale_face_bord

void Parcours_interface::calculer_normale_face_bord(int num_face, double x, double y, double z, double &nx_, double &ny_, double &nz_) const

calculer_sortie_face_bord

int Parcours_interface::calculer_sortie_face_bord(const int face_0, const int num_element, double x0, double y0, double z0, double x1, double y1, double z1, double &x, double &y, double &z) const

Methode outil de Maillage_FT_Disc::deplacer_un_point dans le cas d’un marqueur de la ligne de contact.

P0 est la position initiale du marqueur en contact (sur la face_bord) P1 est la position finale visee apres le deplacement Pour un point P0(x0,y0,z0) sur la face “face_bord” et un point P1(x1,y1,z1), on determine la projection orthogonale p(P1) de P1 sur le plan contenant la face, et on calcule l’intersection (x,y,z) du segment [P0,p(P1)] avec les bords de la face. S’il n’y a pas d’intersection, (x,y,z) = p(P1) et la valeur de retour est - 1, sinon la valeur de retour est le numero de l’arete de la face qui est coupee par le segment [P0,p(P1)] (aretes telles qu’elles sont definies dans la classe Connectivite_frontieres ).

distance_sommet_faces

double Parcours_interface::distance_sommet_faces(const Domaine_VF &domaine_vf, const int num_element, double x, double y, double z) const

duplique

int Parcours_interface::duplique() const override

Reimplements:

duplique

get_correction_parcours_thomas

int Parcours_interface::get_correction_parcours_thomas() const

get_erreur_geometrique

double Parcours_interface::get_erreur_geometrique() const

Renvoie une estimation de l’erreur geometrique (valeur homogene a une distance).

Par exemple si le domaine a une dimension caracteristique de 1e-3 et si la precision relative des calculs Erreur_relative_maxi_ est fixee a 1E - 14 (valeur fixee dans le code source) alors get_erreur_geometrique renvoie 1E - 17. Cette valeur est identique sur tous les processeurs et est calculee a partir de la taille du domaine.

get_info

const Type_info * Parcours_interface::get_info() const override

Donne des informations sur le type de l’ Objet_U .

Reimplements:

get_info

get_parcours_sans_tolerance

int Parcours_interface::get_parcours_sans_tolerance() const

Parcours_interface

Parcours_interface::Parcours_interface()

projeter_vecteur_sur_face

void Parcours_interface::projeter_vecteur_sur_face(const int num_face, double &x_, double &y_, double &z_) const

Methode outil utilisee pour le traitement des lignes de contact.

Projection du vecteur x_, y_, z_ sur le plan parallele a la face num_face passant par l’origine (permet d’obtenir la direction de deplacement d’un sommet sur une ligne de contact).

set_correction_parcours_thomas

void Parcours_interface::set_correction_parcours_thomas()

set_parcours_sans_tolerance

void Parcours_interface::set_parcours_sans_tolerance()

taille_memoire

unsigned Parcours_interface::taille_memoire() const override

Reimplements:

taille_memoire

~Parcours_interface

Parcours_interface::~Parcours_interface()

cree_instance

Objet_U * Parcours_interface::cree_instance()

info

const Type_info * Parcours_interface::info()

self_cast

const Parcours_interface & Parcours_interface::self_cast(const Objet_U &)

self_cast

Parcours_interface & Parcours_interface::self_cast(Objet_U &)

calcul_eq_plan

double Parcours_interface::calcul_eq_plan(const Domaine_VF &domaine_vf, const int num_element, const int num_face_element, double &a, double &b, double &c, double &d) const

calcul_intersection_facelem_2D

int Parcours_interface::calcul_intersection_facelem_2D(const Domaine_VF &domaine_vf, Maillage_FT_Disc &maillage, int num_facette, int num_element) const

calcul_intersection_facelem_3D

int Parcours_interface::calcul_intersection_facelem_3D(const Domaine_VF &domaine_vf, Maillage_FT_Disc &maillage, int num_facette, int num_element) const

Cette methode permet de calculer l’intersection entre une facette et un element du maillage eulerien.

Precondition: dimension = 3

coupe_face_rectangulaire

CutFace_Properties Parcours_interface::coupe_face_rectangulaire(const Domaine_VF &domaine_vf, int num_element, int num_face, const DoubleTab &poly_reelles, const FTd_vecteur3 &norme, const ArrOfInt &polygone_plan_coupe, double epsilon) const

Calcul de la contribution d’une facette a l’indicatrice surfacique et au barycentre sur une face d’un element.

Precondition: dimension = 3

eloigner_sommets_des_faces

int Parcours_interface::eloigner_sommets_des_faces(Maillage_FT_Disc &maillage) const

Pour chaque sommet, s’il est trop pres d’une face eulerienne, deplace le sommet pour l’en eloigner.

Mise a jour de l’espace virtuel des sommets

matrice_triangle

void Parcours_interface::matrice_triangle(int num_element, FTd_vecteur2 &origine, FTd_matrice22 &matrice, double &surface) const

parcourir

void Parcours_interface::parcourir(Maillage_FT_Disc &maillage) const

parcours_facette

void Parcours_interface::parcours_facette(const Domaine_VF &domaine_vf, Maillage_FT_Disc &maillage, ArrOfInt &echange_facettes_numfacette, ArrOfInt &echange_facettes_numelement, int num_facette, int element_depart) const

printOn

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

Ecriture de l’objet sur un flot de sortie Methode a surcharger.

Reimplements:

printOn

readOn

Entree & Parcours_interface::readOn(Entree &) override

Lecture d’un Objet_U sur un flot d’entree Methode a surcharger.

Reimplements:

readOn

transformation_2d

void Parcours_interface::transformation_2d(const FTd_vecteur2 &origine, const FTd_matrice22 &matrice, double x, double y, double &u, double &v) const

Applique la transformation calculee par matrice_triangle a une coordonnee (x,y).

Le resultat est stocke dans (u,v). La troisieme coordonnee barycentrique est implicitement egale a 1 - u - v.

update_erreur_coords_max

void Parcours_interface::update_erreur_coords_max()

uzawa2

double Parcours_interface::uzawa2(const Domaine_VF &domaine_vf, const int elem, double &x, double &y, double &z) const

Algorithme base sur une version initiale de Thomas (recode par BM) Ramene le point (x,y,z) a l’interieur de l’element elem du domaine_vf a une.

distance> = Erreur_max_coordonnees_ par un algorithme d’Uzawa. Valeur de retour: distance finale du sommet aux faces de l’element (positive si le sommet est a l’interieur)

volume_barycentre_hexaedre

CutCell_Properties Parcours_interface::volume_barycentre_hexaedre(const Domaine_VF &domaine_vf, int num_element, const DoubleTab &poly_reelles, const FTd_vecteur3 &norme, const FTd_vecteur3 &centre_de_gravite, const ArrOfInt &polygone_plan_coupe, double epsilon) const

Calcul de la contribution de volume d’une facette a la valeur de l’indicatrice dans un element.

C’est une fraction du volume de l’element comprise entre epsilon et 1. - epsilon Precondition: dimension = 3

volume_barycentre_rectangle

double Parcours_interface::volume_barycentre_rectangle(const Domaine_VF &domaine_vf, int num_element, double x0, double y0, double x1, double y1, double epsilon, double liquid_barycentre[3]) const

Calcul de la matrice 2x2 de transformation pour passer d’une coordonnee dans le repere (x,y) global a une coordonnee (u,v,1 - u - v) barycentrique.

dans un element fini triangulaire. Le point 0 du triangle aura pour coordonnees (0,0,1) Le point 1 (1,0,0) Le point 2 (0,1,0)

volume_rectangle

double Parcours_interface::volume_rectangle(const Domaine_VF &domaine_vf, int num_element, double x0, double y0, double x1, double y1, double epsilon) const

volume_tetraedre

double Parcours_interface::volume_tetraedre(const Domaine_VF &domaine_vf, int num_element, int num_facette, const Maillage_FT_Disc &maillage, const DoubleTab &poly_reelles, const FTd_vecteur3 &centre_de_gravite, double epsilon) const

Calcul de la contribution de volume d’une facette a la valeur de l’indicatrice dans un element.

C’est une fraction du volume de l’element comprise entre epsilon et 1. - epsilon Precondition: dimension = 3

volume_tetraedre_reference

double Parcours_interface::volume_tetraedre_reference(const DoubleTab &poly_reelles_ref, const FTd_vecteur3 &norme_ref, const FTd_vecteur3 &centre_de_gravite_ref, double epsilon) const

Calcul de la contribution de volume d’une facette a la valeur de l’indicatrice dans le tetraedre de reference.

Dans ce tetraedre, on utilise le plan OXZ comme plan de projection (projection selon OY donc) Precondition: la surface doit etre un triangle (nb_sommets_poly = = 3)

volume_triangle

double Parcours_interface::volume_triangle(const Domaine_VF &domaine_vf, int num_element, double x0, double y0, double x1, double y1, int plan_coupe0, int plan_coupe1) const

Calcul de la contribution de volume d’une facette a la valeur de l’indicatrice dans un element.

(x0,y0) et (x1,y1) sont les coordonnees de l’extremite du segment (obligatoirement a l’interieur du triangle ou sur un bord). plan_coupe* donne le numero de la face du triangle sur laquelle se trouve chacun des deux sommets du segment ou - 1 si le sommet est strictement a l’interieur du triangle. C’est une fraction du volume de l’element comprise entre Erreur_relative_maxi et 1. - Erreur_relative_maxi sommet2 ^ ordonnee = v |. | . | . | . face 0 (opposee au sommet 0) | . | 0 1. | : A :. | : :B. - - - - - - > abscisse = u (coordonnees dans le triangle de reference : u,v) sommet0 sommet1

calcul_inverse_matrice33

void Parcours_interface::calcul_inverse_matrice33(const FTd_matrice33 &matrice, FTd_matrice33 &matrice_inv)

Cette methode (statique) permet d’inverser une matrice 3x3.

Precondition: des denominateurs non nuls

calcul_produit_matrice33_vecteur

void Parcours_interface::calcul_produit_matrice33_vecteur(const FTd_matrice33 &matrice, const FTd_vecteur3 &vect, FTd_vecteur3 &res)

Cette methode (statique) permet de calculer le produit d’une matrice 3x3 avec un vecteur 3.

Attributes Documentation#

info_obj (public)

const Type_info Parcours_interface::info_obj

compteur_erreur_grossiere (protected)

int Parcours_interface::compteur_erreur_grossiere = - 123

correction_parcours_thomas_ (protected)

int Parcours_interface::correction_parcours_thomas_ = 0

domaine_elem_ptr (protected)

const IntTab* Parcours_interface::domaine_elem_ptr = nullptr

domaine_sommets_ptr (protected)

const DoubleTab* Parcours_interface::domaine_sommets_ptr = nullptr

drapeaux_elements_parcourus_ (protected)

ArrOfBit Parcours_interface::drapeaux_elements_parcourus_

equations_plans_faces_ (protected)

DoubleTabFT Parcours_interface::equations_plans_faces_

Erreur_max_coordonnees_ (protected)

double Parcours_interface::Erreur_max_coordonnees_ = 0.

Erreur_relative_maxi_ (protected)

double Parcours_interface::Erreur_relative_maxi_

nb_elements_reels_ (protected)

int Parcours_interface::nb_elements_reels_ = - 123

nb_faces_elem_ (protected)

int Parcours_interface::nb_faces_elem_ = - 123

nb_sommets_par_face_ (protected)

int Parcours_interface::nb_sommets_par_face_ = - 123

parcours_sans_tolerance_ (protected)

int Parcours_interface::parcours_sans_tolerance_ = 0

refconnect_front_ (protected)

TRUST_Ref<Connectivite_frontieres *> Parcours_interface::refconnect_front_

refdomaine_vf_ (protected)

TRUST_Ref<Domaine_VF *> Parcours_interface::refdomaine_vf_

type_element_ (protected)

enum Parcours_interface::0 Parcours_interface::type_element_

Valeur_max_coordonnees_ (protected)

double Parcours_interface::Valeur_max_coordonnees_ = 0.

Friends#

void Maillage_FT_Disc::parcourir_maillage

Enums#