santiy checks

KiT-RT · Sep 13, 2024 · 3156758 · 3156758
1 parent 3fcc2ea
commit 3156758
Show file tree

Hide file tree

Showing 3 changed files with 8 additions and 2 deletions.
diff --git a/src/common/mesh.cpp b/src/common/mesh.cpp
@@ -582,10 +582,10 @@ unsigned Mesh::GetCellOfKoordinate( const double x, const double y ) const {
     unsigned koordinate_cell_id = std::numeric_limits<unsigned>::max();
     bool found                  = false;
 
-    // #pragma omp parallel for shared( found )
+#pragma omp parallel for shared( found )
     for( unsigned idx_cell = 0; idx_cell < _numCells; idx_cell++ ) {
         if( IsPointInsideCell( idx_cell, x, y ) ) {
-            // #pragma omp critical
+            #pragma omp critical
             {
                 if( !found ) {
                     koordinate_cell_id = idx_cell;

diff --git a/src/main.cpp b/src/main.cpp
@@ -32,7 +32,12 @@ int main( int argc, char** argv ) {
 // Py_SetProgramName( program );
 #ifdef BUILD_MPI
     MPI_Init( &argc, &argv );
+    printf( "MPI initialized\n" );
 #endif
+#ifndef BUILD_MPI
+    printf( "MPI not initialized\n" );
+#endif
+
     std::string filename = ParseArguments( argc, argv );
 
     // CD  Load Settings from File

diff --git a/src/solvers/snsolver_hpc.cpp b/src/solvers/snsolver_hpc.cpp
@@ -87,6 +87,7 @@ SNSolverHPC::SNSolverHPC( Config* settings ) {
     auto nodes           = _mesh->GetNodes();
     auto cells           = _mesh->GetCells();
 
+#pragma omp parallel for
     for( unsigned idx_cell = 0; idx_cell < _nCells; idx_cell++ ) {
         _areas[idx_cell] = areas[idx_cell];
     }