diff --git a/src/compute-disc-error.cc b/src/compute-disc-error.cc
index 9f57695c46c167e66e344fcbdef38d50fb1da62c..821407435221d8f4faf018b09725f71456962009 100644
--- a/src/compute-disc-error.cc
+++ b/src/compute-disc-error.cc
@@ -310,6 +310,8 @@ void measureEOC(const std::shared_ptr<GridType> grid,
 
 int main (int argc, char *argv[]) try
 {
+  MPIHelper::instance(argc, argv);
+
   // Start Python interpreter
   Python::start();
   Python::Reference main = Python::import("__main__");
diff --git a/src/harmonicmaps.cc b/src/harmonicmaps.cc
index b84095f9f51c6f1d9e2186fb758cee3e60092233..fd2e24a6f2651a60b50ada51883bb2963e0e0c26 100644
--- a/src/harmonicmaps.cc
+++ b/src/harmonicmaps.cc
@@ -69,6 +69,8 @@ using namespace Dune;
 
 int main (int argc, char *argv[]) try
 {
+    MPIHelper::instance(argc, argv);
+
     //feenableexcept(FE_INVALID);
     // Start Python interpreter
     Python::start();