checking MPI_THREAD_MULTIPLE when MPI pre-initialized

fschlimb · fschlimb · commit ccef45d0d262 · 2022-03-11T05:04:54.000-06:00
diff --git a/src/MPITransceiver.cpp b/src/MPITransceiver.cpp
@@ -16,6 +16,9 @@ MPITransceiver::MPITransceiver()
                                      "Please use a thread-safe MPI implementation.");
         }
     } else {
+        MPI_Query_thread(&flag);
+        if(flag != MPI_THREAD_MULTIPLE)
+            throw(std::logic_error("MPI had been initialized incorrectly: not MPI_THREAD_MULTIPLE"));
         std::cerr << "MPI already initialized\n";
     }
     int nranks, rank;

Original file line number	Diff line number	Diff line change
`@@ -16,6 +16,9 @@ MPITransceiver::MPITransceiver()`
`16`	`16`	`"Please use a thread-safe MPI implementation.");`
`17`	`17`	`}`
`18`	`18`	`} else {`
	`19`	`+ MPI_Query_thread(&flag);`
	`20`	`+ if(flag != MPI_THREAD_MULTIPLE)`
	`21`	`+ throw(std::logic_error("MPI had been initialized incorrectly: not MPI_THREAD_MULTIPLE"));`
`19`	`22`	`std::cerr << "MPI already initialized\n";`
`20`	`23`	`}`
`21`	`24`	`int nranks, rank;`