added: broadcast wrapper to ProcessAdm

src/LinAlg/ProcessAdm.C

@@ -53,6 +53,14 @@ ProcessAdm::ProcessAdm(const ProcessAdm& adm) : cout(adm.cout)
 }
 
 
+void ProcessAdm::broadcast (std::vector<double>& vec, int root) const
+{
+#ifdef HAVE_MPI
+  MPI_Bcast(vec.data(), vec.size(), MPI_DOUBLE, root, comm);
+#endif
+}
+
+
 #ifdef HAVE_MPI
 ProcessAdm::ProcessAdm(bool) : cout(std::cout)
 {

src/LinAlg/ProcessAdm.h

@@ -69,6 +69,11 @@ public:
   //! \brief Return if parallel
   int isParallel() const { return parallel; }
 
+  //! \brief Broadcast for a double vector.
+  //! \param vec Double array to broadcast
+  //! \param root Node to broadcast from
+  void broadcast(std::vector<double>& vec, int root) const;
+
 #if defined(HAS_PETSC) || defined(HAVE_MPI)
   //! \brief Return MPI communicator
   MPI_Comm* getCommunicator() { return &comm; }