Hack: Direkter Zugriff zur Matrix eines Operators

[[Imported from SVN: r597]]

fem/feoperator.hh

@@ -86,6 +86,7 @@ namespace Dune {
     //! Stores the destination
     DiscFunctionType * dest_;
 
+  public:
     //! ???
     void assemble ( ) const
     {
@@ -246,7 +247,7 @@ namespace Dune {
     }
 
   public:
-    //! Operator mode
+    //! %Operator mode
     enum OpMode { ON_THE_FLY, ASSEMBLED };
 
     //! ???
@@ -291,6 +292,24 @@ namespace Dune {
       }
     }
 
+    //! ???
+    void assembleMatrix() const
+    {
+      if ( opMode_ == ASSEMBLED )
+      {
+        if ( !matrix_assembled_ )
+        {
+          matrix_ = this->newEmptyMatrix( );
+          assemble();
+        }
+        //matrix_->apply( arg, dest );
+      }
+      else
+      {
+        //multiplyOnTheFly( arg, dest );
+      }
+    }
+
   public:
     //! Store argument and destination
     void prepareGlobal(const DiscFunctionType &Arg, DiscFunctionType & Dest )
@@ -375,7 +394,7 @@ namespace Dune {
     } // end applyLocal
 
 
-    //! ???
+    //! Eliminates the Dirichlet rows and columns
     template <class EntityType>
     void finalizeLocal ( EntityType &en ) const
     {