14#include "arcane/utils/PlatformUtils.h"
15#include "arcane/utils/ScopedPtr.h"
16#include "arcane/utils/ITraceMng.h"
17#include "arcane/utils/OStringStream.h"
19#include "arcane/mesh/DynamicMesh.h"
20#include "arcane/mesh/EdgeUniqueIdBuilder.h"
21#include "arcane/mesh/GhostLayerBuilder.h"
22#include "arcane/mesh/OneMeshItemAdder.h"
24#include "arcane/core/IParallelExchanger.h"
25#include "arcane/core/IParallelMng.h"
26#include "arcane/core/ISerializeMessage.h"
27#include "arcane/core/ISerializer.h"
28#include "arcane/core/ParallelMngUtils.h"
29#include "arcane/core/IMeshUniqueIdMng.h"
45, m_mesh(mesh_builder->
mesh())
46, m_mesh_builder(mesh_builder)
56void EdgeUniqueIdBuilder::
57computeEdgesUniqueIds()
63 info() <<
"Using version=" << edge_version <<
" to compute edges unique ids"
64 <<
" mesh=" << m_mesh->
name();
66 if (edge_version == 1)
67 _computeEdgesUniqueIdsParallel3();
68 else if (edge_version == 2)
69 _computeEdgesUniqueIdsParallelV2();
70 else if (edge_version == 3)
71 _computeEdgesUniqueIdsParallel64bit();
72 else if (edge_version == 0)
73 info() <<
"No renumbering for edges";
75 ARCANE_FATAL(
"Invalid valid version '{0}'. Valid values are 0, 1, 2 or 3");
78 Real diff = (
Real)(end_time - begin_time);
79 info() <<
"TIME to compute edge unique ids=" << diff;
87 if (m_mesh_builder->isVerbose()) {
88 info() <<
"NEW EDGES_MAP after re-indexing";
111 , m_nb_back_edge(nb_back_edge)
112 , m_nb_true_boundary_edge(nb_true_boundary_edge)
117 : m_unique_id(NULL_ITEM_ID)
119 , m_nb_true_boundary_edge(0)
125 bool operator<(
const T_CellEdgeInfo& ci)
const
127 return m_unique_id < ci.m_unique_id;
133 Int64 m_nb_back_edge;
134 Int64 m_nb_true_boundary_edge;
137template <
typename DataType>
138class ItemInfoMultiList
147 MyInfo(
const DataType& d,
Integer n)
161 : m_last_index(5000, true)
166 void add(
Int64 node_uid,
const DataType& data)
168 Integer current_index = m_values.size();
173 m_values.add(MyInfo(data, d->value()));
174 d->value() = current_index;
186class Parallel3EdgeUniqueIdBuilder
189 using BoundaryInfosMap = std::unordered_map<Int32, SharedArray<Int64>>;
205 const Int32 m_my_rank = A_NULL_RANK;
206 const Int32 m_nb_rank = A_NULL_RANK;
207 BoundaryInfosMap m_boundary_infos_to_send;
209 std::unordered_map<Int64, SharedArray<Int64>> m_nodes_info;
211 bool m_is_verbose =
false;
216 void _addEdgeBoundaryInfo(
Edge edge);
217 void _computeEdgesUniqueId();
218 void _sendInfosToOtherRanks();
224Parallel3EdgeUniqueIdBuilder::
227, m_mesh(mesh_builder->
mesh())
228, m_mesh_builder(mesh_builder)
229, m_parallel_mng(m_mesh->parallelMng())
230, m_my_rank(m_parallel_mng->commRank())
231, m_nb_rank(m_parallel_mng->commSize())
232, m_uid_to_subdomain_converter(max_node_uid, m_nb_rank)
233, m_is_verbose(m_mesh_builder->isVerbose())
246 for (
const auto& [key, value] : m_boundary_infos_to_send) {
257 s->
setMode(ISerializer::ModeReserve);
268 debug() <<
"END EXCHANGE";
274void Parallel3EdgeUniqueIdBuilder::
275_addEdgeBoundaryInfo(
Edge edge)
280 Int32 dest_rank = -1;
281 if (!m_is_boundary_nodes[first_node.
localId()]) {
282 v = m_nodes_info[first_node_uid];
285 dest_rank = m_uid_to_subdomain_converter.uidToRank(first_node_uid);
286 v = m_boundary_infos_to_send[dest_rank];
288 v.
add(first_node_uid);
292 v.
add(NULL_ITEM_UNIQUE_ID);
293 v.
add(NULL_ITEM_UNIQUE_ID);
296 <<
" n0=" << ItemPrinter(edge.
node(0)) <<
" n1=" << ItemPrinter(edge.
node(1)) <<
" dest_rank=" << dest_rank;
297 for (Node edge_node : edge.
nodes())
298 v.
add(edge_node.uniqueId());
322 Integer nb_local_edge = m_mesh_builder->oneMeshItemAdder()->nbEdge();
323 info() <<
"ComputeEdgesUniqueIdsParallel3 nb_edge=" << nb_local_edge;
326 edges_opposite_cell_uid.
fill(NULL_ITEM_ID);
334 edges_new_uid.
fill(NULL_ITEM_UNIQUE_ID);
359 m_is_boundary_nodes.resize(node_family->
maxLocalId(),
false);
366 if (face_nb_cell == 1) {
368 m_is_boundary_nodes[ilid] =
true;
381 const bool is_new_item_map_impl = ItemInternalMap::UseNewImpl;
382 if (is_new_item_map_impl) {
383 info() <<
"Edge: ItemInternalMap is using new implementation";
389 _addEdgeBoundaryInfo(edge);
394 _addEdgeBoundaryInfo(edge);
398 _computeEdgesUniqueId();
399 _sendInfosToOtherRanks();
403 info() <<
"END OF TEST NEW EDGE COMPUTE";
409void Parallel3EdgeUniqueIdBuilder::
410_computeEdgesUniqueId()
417 Integer nb_receiver = exchanger->nbReceiver();
418 debug() <<
"NB RECEIVER=" << nb_receiver;
420 for (
Integer i = 0; i < nb_receiver; ++i) {
427 received_infos.
resize(nb_info);
432 while (z < nb_info) {
433 Int64 node_uid = received_infos[z + 0];
436 Int32 edge_type = (
Int32)received_infos[z + 3];
454 for (
const auto& [key, value] : m_nodes_info) {
461 while (z < nb_info) {
469 ItemTypeInfo* itt = itm->
typeFromId(edge_type);
478 my_max_edge_node =
math::max(node_nb_edge, my_max_edge_node);
481 debug() <<
"GLOBAL MAX EDGE NODE=" << global_max_edge_node;
484 m_boundary_infos_to_send.
clear();
486 for (
const auto& [key, value] : m_nodes_info) {
492 while (z < nb_info) {
493 Int64 node_uid = a[z + 0];
495 Int64 edge_uid = a[z + 2];
499 ItemTypeInfo* itt = itm->
typeFromId(edge_type);
503 Integer edge_index = node_nb_edge;
504 Int32 edge_new_owner = sender_rank;
505 for (
Integer y = 0; y < node_nb_edge; ++y) {
506 if (memcmp(&a[indexes[y] + 6], &a[z + 6],
sizeof(
Int64) * edge_nb_node) == 0) {
508 edge_new_owner = (
Int32)a[indexes[y] + 1];
511 Int64 edge_new_uid = (node_uid * global_max_edge_node) + edge_index;
512 Int64Array& v = m_boundary_infos_to_send[sender_rank];
516 v.add(edge_new_owner);
526 my_max_edge_node =
math::max(node_nb_edge, my_max_edge_node);
533void Parallel3EdgeUniqueIdBuilder::
534_sendInfosToOtherRanks()
536 const bool is_verbose = m_mesh_builder->isVerbose();
537 IParallelMng* pm = m_parallel_mng;
542 ItemInternalMap& edges_map = m_mesh->edgesMap();
543 Integer nb_receiver = exchanger->nbReceiver();
544 debug() <<
"NB RECEIVER=" << nb_receiver;
546 for (
Integer i = 0; i < nb_receiver; ++i) {
547 ISerializeMessage* sm = exchanger->messageToReceive(i);
549 ISerializer* s = sm->serializer();
551 Int64 nb_info = s->getInt64();
553 info() <<
"RECEIVE NB_INFO=" << nb_info <<
" from=" << orig_rank;
554 received_infos.resize(nb_info);
555 s->getSpan(received_infos);
556 if ((nb_info % 3) != 0)
557 ARCANE_FATAL(
"info size can not be divided by 3 x={0}", nb_info);
558 Int64 nb_item = nb_info / 3;
559 for (
Int64 z = 0; z < nb_item; ++z) {
560 Int64 old_uid = received_infos[(z * 3)];
561 Int64 new_uid = received_infos[(z * 3) + 1];
562 Int32 new_owner =
static_cast<Int32>(received_infos[(z * 3) + 2]);
564 impl::MutableItemBase iedge(edges_map.tryFind(old_uid));
567 iedge.setUniqueId(new_uid);
568 iedge.setOwner(new_owner, m_my_rank);
571 info() <<
"SetEdgeOwner uid=" << new_uid <<
" owner=" << new_owner
573 <<
" n0=" << ItemPrinter(edge.
node(0)) <<
" n1=" << ItemPrinter(edge.
node(1));
588 bool is_verbose = m_mesh_builder->isVerbose();
589 Integer nb_cell = m_mesh_builder->oneMeshItemAdder()->nbCell();
598 if (cell_uid>max_uid)
601 info() <<
"Max uid=" << max_uid;
609 Integer nb_true_boundary_edge = 0;
610 for(
Edge edge : cell.edges()){
613 else if (edge.
nbCell()==1){
614 ++nb_true_boundary_edge;
617 cell_nb_num_back_edge[cell_uid] = nb_num_back_edge;
618 cell_true_boundary_edge[cell_uid] = nb_true_boundary_edge;
622 for(
Integer i=0; i<nb_cell; ++i ){
623 cell_first_edge_uid[i] = current_edge_uid;
624 current_edge_uid += cell_nb_num_back_edge[i] + cell_true_boundary_edge[i];
628 for(
Integer i=0; i<nb_cell; ++i ){
629 info() <<
"Recv: Cell EdgeInfo celluid=" << i
630 <<
" firstedgeuid=" << cell_first_edge_uid[i]
631 <<
" nbback=" << cell_nb_num_back_edge[i]
632 <<
" nbbound=" << cell_true_boundary_edge[i];
639 Integer nb_true_boundary_edge = 0;
640 for(
Edge edge : cell.edges() ){
641 Int64 edge_new_uid = NULL_ITEM_UNIQUE_ID;
644 edge_new_uid = cell_first_edge_uid[cell_uid] + nb_num_back_edge;
647 else if (edge.
nbCell()==1){
648 edge_new_uid = cell_first_edge_uid[cell_uid] + cell_nb_num_back_edge[cell_uid] + nb_true_boundary_edge;
649 ++nb_true_boundary_edge;
651 if (edge_new_uid!=NULL_ITEM_UNIQUE_ID){
665 Int64 opposite_cell_uid = NULL_ITEM_UNIQUE_ID;
666 bool true_boundary =
false;
667 bool internal_other =
false;
670 else if (edge.
nbCell()==1){
671 true_boundary =
true;
674 internal_other =
true;
677 ostr() <<
"NEW LOCAL ID FOR CELLEDGE " << cell_uid <<
' '
678 << index <<
' ' << edge.
uniqueId() <<
" (";
680 ostr() <<
' ' << node.uniqueId();
684 ostr() <<
" internal-other";
686 ostr() <<
" true-boundary";
687 if (opposite_cell_uid!=NULL_ITEM_ID)
688 ostr() <<
" opposite " << opposite_cell_uid;
693 info() << ostr.str();
700void EdgeUniqueIdBuilder::
701_computeEdgesUniqueIdsParallelV2()
722 if (total_max_uid>INT32_MAX)
723 ARCANE_FATAL(
"Max uniqueId() for node is too big v={0} max_allowed={1}",total_max_uid,INT32_MAX);
727 Node node1{edge.
node(1)};
728 Int64 new_uid = (node0.uniqueId().asInt64() * total_max_uid) + node1.
uniqueId().asInt64();
736void EdgeUniqueIdBuilder::
737_computeEdgesUniqueIdsParallel3()
739 IParallelMng* pm = m_mesh->parallelMng();
740 ItemInternalMap& nodes_map = m_mesh->nodesMap();
743 Int64 my_max_node_uid = NULL_ITEM_UNIQUE_ID;
745 Int64 node_uid = item.uniqueId();
746 if (node_uid > my_max_node_uid)
747 my_max_node_uid = node_uid;
750 debug() <<
"NODE_UID_INFO: MY_MAX_UID=" << my_max_node_uid
751 <<
" GLOBAL=" << global_max_node_uid;
753 Parallel3EdgeUniqueIdBuilder builder(
traceMng(), m_mesh_builder, global_max_node_uid);
760void EdgeUniqueIdBuilder::
761_computeEdgesUniqueIdsParallel64bit()
765 ItemInternalMap& edges_map = m_mesh->edgesMap();
767 std::hash<Int64> hasher;
770 Node node0{edge.node(0)};
771 Node node1{edge.node(1)};
772 size_t hash0 = hasher(node0.uniqueId().asInt64());
773 size_t hash1 = hasher(node1.uniqueId().asInt64());
774 hash0 ^= hash1 + 0x9e3779b9 + (hash0 << 6) + (hash0 >> 2);
775 Int64 new_uid = hash0 & 0x7fffffff;
776 edge.mutableItemBase().setUniqueId(new_uid);
#define ARCANE_FATAL(...)
Macro envoyant une exception FatalErrorException.
void fill(const DataType &data)
Remplissage du tableau.
void clear()
Supprime les éléments du tableau.
void resize(Int64 s)
Change le nombre d'éléments du tableau à s.
void reserve(Int64 new_capacity)
Réserve le mémoire pour new_capacity éléments.
void add(ConstReferenceType val)
Ajoute l'élément val à la fin du tableau.
EdgeConnectedListViewType edges() const
Liste des arêtes de la maille.
Vue constante d'un tableau de type T.
constexpr Integer size() const noexcept
Nombre d'éléments du tableau.
Int32 nbCell() const
Nombre de mailles connectées à l'arête.
Int32 nbCell() const
Nombre de mailles de la face (1 ou 2)
Table de hachage pour tableaux associatifs.
Data * lookupAdd(KeyTypeConstRef id, const ValueType &value, bool &is_add)
Recherche ou ajoute la valeur correspondant à la clé id.
bool add(KeyTypeConstRef id, const ValueType &value)
Ajoute la valeur value correspondant à la clé id.
void eachValue(const Lambda &lambda)
Applique le fonctor f à tous les éléments de la collection et utilise x->value() (de type ValueType) ...
Interface d'une famille d'entités.
virtual Int32 maxLocalId() const =0
virtual Integer edgeBuilderVersion() const =0
Version de la numérotation des arêtes.
Echange d'informations entre processeurs.
virtual void addSender(Int32 rank)=0
Ajoute un processeur à envoyer.
virtual Integer nbSender() const =0
Nombre de processeurs auquel on envoie.
virtual bool initializeCommunicationsMessages()=0
Calcule les communications.
virtual ISerializeMessage * messageToSend(Integer i)=0
Message destiné au ième processeur.
virtual void processExchange()=0
Effectue l'échange avec les options par défaut de ParallelExchangerOptions.
Interface du gestionnaire de parallélisme pour un sous-domaine.
virtual char reduce(eReduceType rt, char v)=0
Effectue la réduction de type rt sur le réel v et retourne la valeur.
virtual void barrier()=0
Effectue une barière.
Interface d'un sérialiseur.
void reserveInt64(Int64 n)
Réserve pour n Int64.
@ ModePut
Le sérialiseur attend des reserve()
@ ModeGet
Le sérialiseur attend des get()
virtual Int64 getInt64()=0
Récupère une taille.
virtual void allocateBuffer()=0
Alloue la mémoire du sérialiseur.
virtual void putSpan(Span< const Real > values)
Ajoute le tableau values.
virtual void getSpan(Span< Real > values)
Récupère le tableau values.
virtual void reserveSpan(eBasicDataType dt, Int64 n)=0
Réserve de la mémoire pour n valeurs de dt.
virtual void setMode(eMode new_mode)=0
Positionne le fonctionnement actuel.
virtual void putInt64(Int64 value)=0
Ajoute l'entier value.
Interface du gestionnaire de traces.
virtual void flush()=0
Flush tous les flots.
ItemUniqueId uniqueId() const
Numéro unique de l'entité
ItemBase backCell() const
Maille derrière l'entité (nullItem() si aucune)
Infos sur un type d'entité du maillage.
Integer nbLocalNode() const
Nombre de noeuds de l'entité
Gestionnaire des types d'entités de maillage.
ItemTypeInfo * typeFromId(Integer id) const
Type correspondant au numéro id.
Node node(Int32 i) const
i-ème noeud de l'entité
NodeConnectedListViewType nodes() const
Liste des noeuds de l'entité
NodeLocalIdView nodeIds() const
Liste des noeuds de l'entité
Classe de base d'un élément de maillage.
impl::MutableItemBase mutableItemBase() const
Partie interne modifiable de l'entité.
constexpr Int32 localId() const
Identifiant local de l'entité dans le sous-domaine du processeur.
ItemUniqueId uniqueId() const
Identifiant unique sur tous les domaines.
impl::ItemBase itemBase() const
Partie interne de l'entité.
Int16 type() const
Type de l'entité
Interface d'un message de sérialisation entre IMessagePassingMng.
virtual MessageRank destination() const =0
Rang du destinataire (si isSend() est vrai) ou de l'envoyeur.
virtual ISerializer * serializer()=0
Sérialiseur.
Int32 value() const
Valeur du rang.
Flot de sortie lié à une String.
Référence à une instance.
Vecteur 1D de données avec sémantique par référence.
Classe d'accès aux traces.
TraceAccessor(ITraceMng *m)
Construit un accesseur via le gestionnaire de trace m.
TraceMessageDbg debug(Trace::eDebugLevel=Trace::Medium) const
Flot pour un message de debug.
TraceMessage info() const
Flot pour un message d'information.
ITraceMng * traceMng() const
Gestionnaire de trace.
Vecteur 1D de données avec sémantique par valeur (style STL).
Construction d'un maillage de manière incrémentale.
Implémentation d'un maillage.
IParallelMng * parallelMng() override
Gestionnaire de parallèlisme.
String name() const override
Nom du maillage.
ItemTypeMng * itemTypeMng() const override
Gestionnaire de types d'entités associé
IMeshUniqueIdMng * meshUniqueIdMng() const override
Gestionnare de la numérotation des identifiants uniques.
EdgeUniqueIdBuilder(DynamicMeshIncrementalBuilder *mesh_builder)
Construit une instance pour le maillage mesh.
void _computeEdgesUniqueIdsSequential()
Calcul les numéros uniques de chaque edge en séquentiel.
Tableau associatif de ItemInternal.
void notifyUniqueIdsChanged()
Notifie que les numéros uniques des entités ont changés.
void eachItem(const Lambda &lambda)
Fonction template pour itérer sur les entités de l'instance.
void _exchangeData(IParallelExchanger *exchanger)
void compute()
Calcule les numéros uniques de chaque edge en parallèle.
T max(const T &a, const T &b, const T &c)
Retourne le maximum de trois éléments.
@ ReduceMax
Maximum des valeurs.
Ref< IParallelExchanger > createExchangerRef(IParallelMng *pm)
Retourne une interface pour transférer des messages entre rangs.
Array< Int64 > Int64Array
Tableau dynamique à une dimension d'entiers 64 bits.
UniqueArray< Int64 > Int64UniqueArray
Tableau dynamique à une dimension d'entiers 64 bits.
std::int64_t Int64
Type entier signé sur 64 bits.
Int32 Integer
Type représentant un entier.
@ Int64
Donnée de type entier 64 bits.
ConstArrayView< Int64 > Int64ConstArrayView
Equivalent C d'un tableau à une dimension d'entiers 64 bits.
UniqueArray< Int32 > Int32UniqueArray
Tableau dynamique à une dimension d'entiers 32 bits.
double Real
Type représentant un réel.
UniqueArray< Integer > IntegerUniqueArray
Tableau dynamique à une dimension d'entiers.
std::int32_t Int32
Type entier signé sur 32 bits.