Changeset d6a338 for src/comm/comm.hpp

Timestamp:

Apr 16, 2013, 11:27:31 AM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

177495

Parents:

e58835d

Message:

Save the whole domain decomposition instead of just the part for the current process.

File:

• src/comm/comm.hpp (modified) (4 diffs)

src/comm/comm.hpp

@@ -29 +29 @@
 #define COMM_HPP_
 
+#include <map>
 #include <list>
 
@@ -34 +35 @@
 #include "base/index.hpp"
 #include "base/object.hpp"
+#include "comm/domain_decomposition.hpp"
 #include "thirdparty/pugixml/pugixml.hpp"
 
@@ -40 +42 @@
 
 class DomainDecomposition;
+class GlobalIndices;
 class Grid;
+class Interface;
 class Multigrid;
 class Stencil;
@@ -129 +133 @@
   const Boundary& BoundaryConditions() const {return bound;}
 
-  DomainDecomposition& GetDomainDecomposition() {return *dd;}
+  vmg_float ComputeResidualNorm(Multigrid& sol, Multigrid& rhs);
 
-  vmg_float ComputeResidualNorm(Multigrid& sol, Multigrid& rhs);
+  void ComputeDomainDecomposition(const Interface& interface) {dd->Compute(*this, interface, decomposed_global);}
+  const std::vector<GlobalIndices>& DecomposedGlobalMe() const {return decomposed_global.find(GlobalPos())->second;}
+  const std::map<Index, std::vector<GlobalIndices> >& DecomposedGlobal() const {return decomposed_global;}
+
 
 protected:
   const std::string& OutputPath();
   int OutputCount() {return ++output_count;}
+
+  std::map<Index, std::vector<GlobalIndices> > decomposed_global;
 
   Boundary bound;