Changeset dfed1c for src/commands/com_force_discrete_compatibility.cpp

Timestamp:

Nov 22, 2011, 9:22:10 PM (14 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

facba0

Parents:

66f24d

Message:

Major vmg update.

git-svn-id: ​https://svn.version.fz-juelich.de/scafacos/trunk@1136 5161e1c8-67bf-11de-9fd5-51895aff932f

File:

• src/commands/com_force_discrete_compatibility.cpp (modified) (3 diffs)

src/commands/com_force_discrete_compatibility.cpp

@@ -6 +6 @@
  * @brief  May be used to explicitly force the discrete
  *         compatibility condition \sum_{i,j,k} f_{i,j,k} = 0.
- *         This should not be necessary since this is also
- *         handled in the library for arbitrary compatibility
- *         conditions.
+ *
  */
 
@@ -27 +25 @@
   Request Run(Command::argument_vector arguments)
   {
-    (*MG::GetRhs())().ForceDiscreteCompatibilityCondition();
+    MPE_EVENT_BEGIN()
+
+    Grid& rhs = (*MG::GetRhs())(MG::GetRhs()->MaxLevel());
+    Comm& comm = *MG::GetComm();
+
+    if (comm.BoundaryConditions()[0] == Periodic &&
+        comm.BoundaryConditions()[1] == Periodic &&
+        comm.BoundaryConditions()[2] == Periodic)
+      rhs.ForceDiscreteCompatibilityCondition();
+
+    MPE_EVENT_END()
 
     return Continue;
@@ -36 +44 @@
 };
 
-CREATE_INITIALIZER(VMGCommandForceDiscreteCompatibility);
+CREATE_INITIALIZER(VMGCommandForceDiscreteCompatibility)