Changeset d6a338 for src/mg.cpp

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/mg.cpp (modified) (2 diffs)

src/mg.cpp

@@ -133 +133 @@
   if (GetFactory().TestObject("COMM") && GetFactory().TestObject("INTERFACE")) {
 
-    Multigrid* sol = new Multigrid(GetComm(), GetInterface());
+    GetComm()->ComputeDomainDecomposition(*GetInterface());
+
+    Multigrid* sol = new Multigrid(*GetComm(), *GetInterface());
     sol->Register("SOL");
 
-    Multigrid* rhs = new Multigrid(GetComm(), GetInterface());
+    Multigrid* rhs = new Multigrid(*GetComm(), *GetInterface());
     rhs->Register("RHS");
 
@@ -149 +151 @@
       GetCycle()->Generate();
 
-    if (GetFactory().TestObject("COMM"))
-      GetComm()->PostInit(*GetSol(), *GetRhs());
+    GetComm()->PostInit(*GetSol(), *GetRhs());
 
   }