reduce volumetric elements

seamplex · Sep 29, 2023 · 59d9c60 · 59d9c60
1 parent 60aedae
commit 59d9c60
Showing 1 changed file with 5 additions and 2 deletions.
diff --git a/src/pdes/build.c b/src/pdes/build.c
@@ -89,7 +89,10 @@ int feenox_problem_build(void) {
     }
   }
 
-  if (volumetric_elements == 0) {
+  int valid_local = (volumetric_elements > 0);
+  int valid_global = 0;
+  MPI_Allreduce(&valid_local, &valid_global, 1, MPIU_INT, MPIU_SUM, PETSC_COMM_WORLD);
+  if (valid_local == 0) {
     feenox_push_error_message("no volumetric elements found in '%s'", feenox.pde.mesh->file->name);
     return FEENOX_ERROR;
   }
@@ -413,4 +416,4 @@ int feenox_problem_rhs_set(element_t *e, unsigned int q, double *value) {
 #endif
 
   return FEENOX_OK;
-}
+}