Changeset 177495 for src/comm/mpi

Timestamp:

Apr 17, 2013, 6:35:40 PM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

4f3474

Parents:

d6a338

Message:

Implemented boundary values derived from continious problem.

Location:

src/comm/mpi

Files:

• settings.cpp (modified) (3 diffs)
• settings.hpp (modified) (3 diffs)

src/comm/mpi/settings.cpp

@@ -44 +44 @@
 #include <string>
 
+#include "base/index.hpp"
 #include "comm/comm.hpp"
 #include "comm/mpi/settings.hpp"
+#include "discretization/boundary_value_setter.hpp"
 #include "grid/multigrid.hpp"
 #include "grid/tempgrid.hpp"
@@ -218 +220 @@
     AddDatatypeGlobal(FinerGrid(sol(i-1)), sol(i), 1);
   }
+
+  InitializeBoundaryValues();
 }
 
@@ -428 +432 @@
 }
 
+void VMG::MPI::Settings::InitializeBoundaryValues()
+{
+  assert(bv_ranks.size() == 0);
+
+  if (MG::GetFactory().TestObject("BOUNDARY_VALUE_SETTER")) {
+
+    Index coord;
+
+    const int level_index = MG::GetRhs()->MaxLevel() - MG::GetRhs()->GlobalMaxLevel();
+    const std::vector<BoundaryValue>& bvs = MG::GetBoundaryValueSetter()->BoundaryValues();
+    const std::map<Index, std::vector<GlobalIndices> >& global = MG::GetComm()->DecomposedGlobal();
+
+    assert(global.find(0)->second[level_index].BoundaryType() == GlobalMax);
+
+    MPI_Comm comm = CommunicatorLocal((*MG::GetRhs())(MG::GetRhs()->GlobalMaxLevel()));
+
+    bv_ranks.reserve(bvs.size());
+
+    for (std::vector<BoundaryValue>::const_iterator iter_b = bvs.begin(); iter_b != bvs.end(); ++iter_b) {
+      for (std::map<Index, std::vector<GlobalIndices> >::const_iterator iter_g = global.begin(); iter_g != global.end(); ++iter_g) {
+        if (iter_b->GetIndex().IsComponentwiseGreaterOrEqual(iter_g->second[level_index].LocalBegin()) &&
+            iter_b->GetIndex().IsComponentwiseLess(iter_g->second[level_index].LocalEnd())) {
+            bv_ranks.push_back(0);
+            coord = iter_g->first;
+            MPI_Cart_rank(comm, coord.vec(), &bv_ranks.back());
+          break;
+        }
+      }
+    }
+  }
+}
+
 std::string VMG::MPI::Settings::ToString() const
 {

src/comm/mpi/settings.hpp

@@ -70 +70 @@
   VMG::MPI::DatatypesLocal& DatatypesLocal(const Grid& grid);
 
+  const std::vector<int>& BoundaryRanks() const {return bv_ranks;}
+
   std::string ToString() const;
 
@@ -79 +81 @@
   void CreateLocalCommunicator(MPI_Comm& comm_global, const Grid& grid);
 
+  void InitializeBoundaryValues();
+
   std::map<int, MPI_Comm> communicators_global;
   std::map<KeyUnsorted, MPI_Comm> communicators_local;
@@ -86 +90 @@
   std::map<KeyUnsorted, VMG::MPI::DatatypesGlobal> datatypes_global;
   std::map<KeyUnsorted, VMG::MPI::DatatypesLocal> datatypes_local;
+
+  std::vector<int> bv_ranks;
+
 };