Changeset 177495 for src/comm/mpi/settings.cpp

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.

File:

• src/comm/mpi/settings.cpp (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
 {