Changeset f57182 for src/samples/discretization_poisson_fd.cpp

Timestamp:

Mar 30, 2013, 2:44:52 PM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

8180d8

Parents:

d13e27

git-author:

Julian Iseringhausen <isering@…> (06/11/12 14:02:16)

git-committer:

Julian Iseringhausen <isering@…> (03/30/13 14:44:52)

Message:

Open boundary conditions.

Conflicts:

lib/vmg/src/Makefile.am
lib/vmg/src/base/factory.cpp
lib/vmg/test/unit_test/library/dirichlet_fas_lr_mpi.cpp

File:

• src/samples/discretization_poisson_fd.cpp (modified) (4 diffs)

src/samples/discretization_poisson_fd.cpp

@@ -18 +18 @@
 
 /**
- * @file   discretization_poisson_fd_collatz.cpp
+ * @file   discretization_poisson_fd.cpp
  * @author Julian Iseringhausen <isering@ins.uni-bonn.de>
  * @date   Mon Apr 18 13:03:47 2011
  *
  * @brief  Discretization of the poisson equation
- *         using the Collatz Mehrstellen Ansatz.
- *         Discretization error: O(h^4)
+ *         using the standard 7-point stencil or
+ *         the Collatz Mehrstellen Ansatz.
+ *         Discretization error: O(h^2) or O(h^4)
  *
  */
@@ -38 +39 @@
 using namespace VMG;
 
-DiscretizationPoissonFD::DiscretizationPoissonFD(const int& order) :
-  Discretization(order)
+void DiscretizationPoissonFD::InitDiscretizationPoissonFD()
 {
   switch (order)
@@ -83 +83 @@
   if (order == 4) {
 
-    Grid& rhs = MG::GetRhsMaxLevel();
-
     Stencil stencil(6.0/12.0);
     stencil.push_back(-1,  0,  0, 1.0/12.0);
@@ -93 +91 @@
     stencil.push_back( 0,  0,  1, 1.0/12.0);
 
-    stencil.Apply(rhs);
+    stencil.Apply(MG::GetRhsMaxLevel());
 
   }