fix mpi comm

Author: asartori86

include/lac/lac_initializer.h

@@ -3,6 +3,7 @@
 
 // This includes all types we know of.
 #include "lac/lac_type.h"
+#include <deal.II/base/mpi.h>
 #include <mpi.h>
 
 /**
@@ -15,20 +16,20 @@ class ScopedLACInitializer
   ScopedLACInitializer(const std::vector<types::global_dof_index> &dofs_per_block,
                        const std::vector<IndexSet> &owned,
                        const std::vector<IndexSet> &relevant,
-                       const MPI_Comm &comm = MPI_COMM_WORLD):
+                       const MPI_Comm extcomm= MPI_COMM_WORLD):
     dofs_per_block(dofs_per_block),
     owned(owned),
     relevant(relevant),
-    comm(comm)
-  {};
+    comm(Utilities::MPI::duplicate_communicator(extcomm))
+  {}
 
   /**
    * Initialize a non ghosted TrilinosWrappers::MPI::BlockVector.
    */
   void operator() (TrilinosWrappers::MPI::BlockVector &v, bool fast=false)
   {
     v.reinit(owned, comm, fast);
-  };
+  }
 
 
   /**
@@ -37,15 +38,15 @@ class ScopedLACInitializer
   void ghosted(TrilinosWrappers::MPI::BlockVector &v, bool fast=false)
   {
     v.reinit(owned, relevant, comm, fast);
-  };
+  }
 
   /**
    * Initialize a serial BlockVector<double>.
    */
   void operator() (BlockVector<double> &v, bool fast=false)
   {
     v.reinit(dofs_per_block, fast);
-  };
+  }
 
 
   /**
@@ -56,7 +57,7 @@ class ScopedLACInitializer
   {
     Assert(false, ExcInternalError("You tried to create a ghosted vector in a serial run."));
     (void)fast;
-  };
+  }
 
   /**
    * Initialize a Trilinos Sparsity Pattern.
@@ -112,7 +113,7 @@ class ScopedLACInitializer
   /**
    * MPI Communicator.
    */
-  const MPI_Comm &comm;
+  MPI_Comm comm;
 };
 
 

include/pidomus.h

@@ -66,7 +66,7 @@ class piDoMUS : public ParameterAcceptor, public SundialsInterface<typename LAC:
 
   piDoMUS (const std::string &name,
            const BaseInterface<dim, spacedim, LAC> &energy,
-           const MPI_Comm &comm = MPI_COMM_WORLD);
+           const MPI_Comm comm = MPI_COMM_WORLD);
 
   virtual void declare_parameters(ParameterHandler &prm);
   virtual void parse_parameters_call_back();
@@ -301,7 +301,7 @@ class piDoMUS : public ParameterAcceptor, public SundialsInterface<typename LAC:
 
   void set_constrained_dofs_to_zero(typename LAC::VectorType &v) const;
 
-  const MPI_Comm &comm;
+  MPI_Comm comm;
   const BaseInterface<dim, spacedim, LAC>    &interface;
 
   unsigned int n_cycles;

source/pidomus.cc

@@ -44,11 +44,11 @@ using namespace deal2lkit;
 template <int dim, int spacedim, typename LAC>
 piDoMUS<dim, spacedim, LAC>::piDoMUS (const std::string &name,
                                       const BaseInterface<dim, spacedim, LAC> &interface,
-                                      const MPI_Comm &communicator)
+                                      const MPI_Comm communicator)
   :
   ParameterAcceptor(name),
   SundialsInterface<typename LAC::VectorType>(communicator),
-  comm(communicator),
+  comm(Utilities::MPI::duplicate_communicator(communicator)),
   interface(interface),
   pcout (std::cout,
          (Utilities::MPI::this_mpi_process(comm)

tests/checkpoint_restart_02.cc

@@ -1,9 +1,7 @@
 #include "pidomus.h"
-#include "interfaces/poisson_problem_signals.h"
+#include "interfaces/poisson_problem.h"
 #include "tests.h"
 
-
-// test that we can write snapshots
 using namespace dealii;
 int main (int argc, char *argv[])
 {
@@ -12,7 +10,6 @@ int main (int argc, char *argv[])
                                                       numbers::invalid_unsigned_int);
 
   MPILogInitAll log;
-  deallog.depth_file(1);
 
   deallog.depth_file(1);
   const int dim = 2;
@@ -29,6 +26,5 @@ int main (int argc, char *argv[])
 
   deallog << "OK" << std::endl;
 
-
   return 0;
 }

tests/parameters/checkpoint_restart_02.prm

@@ -164,7 +164,7 @@ subsection pidomus
   set Print some useful informations about processes = true
   set Refine mesh during transient                   = false
   set Resume computation from snapshot               = false
-  set Save snapshots during simulation               = false
+  set Save snapshots during simulation               = true
   set Snapshot prefix                                = 
   set Threshold for solver's restart                 = 1e-2
   set Time stepper                                   = imex