Feature/dynamicpolymorphism

This is the first MR in a series to implement #26 (closed).

The solver interfaces use dynamic polymorphism, but currently it is not possible to fully benefit from this feature. We want to allow the user to take full benefit from a dynamic interface for solver/precoditioner/scalarproduct/linearoperator.

The main change is that the category isn't an enum anymore, but virtual method catergory().

This MR contains the following changes:

• new method category() to check for the same SolverCategory of different components at run time. This is a slight drawback compared to the current compile time checks, but it adds a lot of flexibility as we will be able to write factory classes for solvers and compose/exchange solvers at runtime.
• add SolverCategory::category(object) to be be backwards compatible
• add a cmake switch to enable/disable the backwards compatibility
• update solvers to the new interface
• update preconditioners to the new interface
• update scalar products to the new interface
• update operators to the new interface
• update AMG implementation to the new interface
  • remove static dependency of PropertiesGraphCreator on the category

This MR partially incorporates commits from !5 (closed) and !1 (closed).

dune/istl/schwarz.hh

@@ -98,11 +98,6 @@ namespace Dune {
     //! implementing the same interface.
     typedef C communication_type;
 
-    enum {
-      //! \brief The solver category.
-      category=SolverCategory::overlapping
-    };
-
     /**
      * @brief constructor: just store a reference to a matrix.
      *
@@ -137,6 +132,12 @@ namespace Dune {
       return _A_;
     }
 
+    //! Category of the linear operator (see SolverCategory::Category)
+    virtual SolverCategory::Category category() const
+    {
+      return SolverCategory::overlapping;
+    }
+
   private:
     const matrix_type& _A_;
     const communication_type& communication;
@@ -177,9 +178,6 @@ namespace Dune {
     //! implementing the same interface.
     typedef C communication_type;
 
-    //! define the category
-    enum {category=SolverCategory::overlapping};
-
     /*! \brief Constructor needs to know the grid
      * \param com The communication object for syncing overlap and copy
      * data points. (E.~g. OwnerOverlapCopyCommunication )
@@ -207,6 +205,12 @@ namespace Dune {
       return communication.norm(x);
     }
 
+    //! Category of the scalar product (see SolverCategory::Category)
+    virtual SolverCategory::Category category() const
+    {
+      return SolverCategory::overlapping;
+    }
+
   private:
     const communication_type& communication;
   };
@@ -263,12 +267,6 @@ namespace Dune {
     //! \brief The type of the communication object.
     typedef C communication_type;
 
-    // define the category
-    enum {
-      //! \brief The category the precondtioner is part of.
-      category=SolverCategory::overlapping
-    };
-
     /*! \brief Constructor.
 
        constructor gets all parameters to operate the prec.
@@ -313,6 +311,12 @@ namespace Dune {
      */
     virtual void post (X& x) {}
 
+    //! Category of the preconditioner (see SolverCategory::Category)
+    virtual SolverCategory::Category category() const
+    {
+      return SolverCategory::overlapping;
+    }
+
   private:
     //! \brief The matrix we operate on.
     const matrix_type& _A_;
@@ -374,12 +378,6 @@ namespace Dune {
     //! implementing the same interface.
     typedef C communication_type;
 
-    // define the category
-    enum {
-      //! \brief The category the precondtioner is part of.
-      category=SolverCategory::overlapping
-    };
-
     /*! \brief Constructor.
 
        constructor gets all parameters to operate the prec.
@@ -430,6 +428,12 @@ namespace Dune {
       preconditioner.post(x);
     }
 
+    //! Category of the preconditioner (see SolverCategory::Category)
+    virtual SolverCategory::Category category() const
+    {
+      return SolverCategory::overlapping;
+    }
+
   private:
     //! \brief a sequential preconditioner
     T& preconditioner;