Make global numbering compatible with overlap - not used for now.

MmgTools · Sep 27, 2024 · b567cc3 · b567cc3
1 parent a9ca05a
commit b567cc3
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/src/libparmmg.c b/src/libparmmg.c
@@ -1103,9 +1103,13 @@ int PMMG_Compute_verticesGloNum( PMMG_pParMesh parmesh,MPI_Comm comm ){
  }
 
  /* Store owner in the point flag */
+ MMG5_int np_overlap = 0;
  for( ip = 1; ip <= mesh->np; ip++ ) {
  ppt = &mesh->point[ip];
- if (ppt->tag & MG_OVERLAP) continue;
+ if (ppt->tag & MG_OVERLAP) {
+ ++np_overlap;
+ continue;
+ }
  ppt->flag = parmesh->myrank;
  }
 
@@ -1117,7 +1121,7 @@ int PMMG_Compute_verticesGloNum( PMMG_pParMesh parmesh,MPI_Comm comm ){
  }
 
  /* Count owned nodes */
- nowned = mesh->np;
+ nowned = mesh->np - np_overlap;
  for( idx = 0; idx < int_node_comm->nitem; idx++ ) {
  if( intvalues[idx] != parmesh->myrank ) nowned--;
  }