Changeset 4571da for src/base/discretization.hpp

Timestamp:

Apr 27, 2012, 11:34:57 PM (14 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

1a92cf

Parents:

b2154a3

Message:

vmg: Implement fourth-order discretization of the Poisson equation.

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

File:

• src/base/discretization.hpp (modified) (3 diffs)

src/base/discretization.hpp

@@ -28 +28 @@
 {
 public:
-  Discretization() : stencil(1.0) {}
-  Discretization(const Stencil& stencil_) : stencil(stencil_) {}
+  Discretization() :
+    stencil(1.0),
+    order(2)
+  {}
+
+  Discretization(const int& order) :
+    stencil(1.0),
+    order(order)
+  {}
+
+  Discretization(const Stencil& stencil_, const int& order) :
+    stencil(stencil_),
+    order(order)
+  {}
 
   Discretization(std::string id) :
     Object(id),
-    stencil(1.0)
+    stencil(1.0),
+    order(2)
   {}
 
-  Discretization(std::string id, const Stencil& stencil_) :
+  Discretization(std::string id, const int& order) :
     Object(id),
-    stencil(stencil_)
+    stencil(1.0),
+    order(order)
+  {}
+
+  Discretization(std::string id, const Stencil& stencil_, const int& order) :
+    Object(id),
+    stencil(stencil_),
+    order(order)
   {}
 
@@ -45 +65 @@
   virtual vmg_float OperatorPrefactor(const Grid& grid) const = 0; ///< Returns the prefactor of the operator.
 
-  /** 
+  virtual void ModifyRightHandSide() {}
+
+  /**
    * This function gets called whenever boundary points at inner boundaries are needed.
    * Inner boundaries occur when using adaptive grid refinement.
-   * 
+   *
    * @param sol_fine Solution vector / fine level
    * @param rhs_fine Right handside vector / fine level
@@ -60 +82 @@
 
 private:
-  virtual void SetInnerBoundaryCompute(Grid& sol_fine, Grid& rhs_fine, Grid& sol_coarse) const = 0;
+  virtual void SetInnerBoundaryCompute(Grid& sol_fine, Grid& rhs_fine, Grid& sol_coarse) const {}
 
 protected:
   VMG::Stencil stencil;
+  int order;
 };