omega in PreconditionBlockSOR can be changed

git-svn-id: https://svn.dealii.org/trunk@2707 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/lac/include/lac/precondition_block.h b/deal.II/lac/include/lac/precondition_block.h
index 5c64cf7236feb7917b6a6d4b8c6573d81a72922c..1875e3fd44460ac2d5350203fefff7324929ed2d 100644 (file)
--- a/deal.II/lac/include/lac/precondition_block.h
+++ b/deal.II/lac/include/lac/precondition_block.h
@@ -223,6 +223,11 @@ class PreconditionBlockSOR : public PreconditionBlock<number,inverse_type>
                                      */
     virtual ~PreconditionBlockSOR();
 
+                                    /**
+                                     * Set the relaxation parameter.
+                                     */
+    void set_omega(number omega);
+    
                                     /**
                                      * Execute block SOR preconditioning.
                                      * Make sure that the right
@@ -248,7 +253,7 @@ class PreconditionBlockSOR : public PreconditionBlock<number,inverse_type>
                                     /**
                                      * Damping parameter.
                                      */
-    const number omega;
+    number omega;
 };
 
 

diff --git a/deal.II/lac/include/lac/precondition_block.templates.h b/deal.II/lac/include/lac/precondition_block.templates.h
index 3084c1fb19e934b9ecf3969954fc48915becbc5b..1a299727cf3df9a151cfe90bbcfdf982e8e4855c 100644 (file)
--- a/deal.II/lac/include/lac/precondition_block.templates.h
+++ b/deal.II/lac/include/lac/precondition_block.templates.h
@@ -135,6 +135,17 @@ template<typename number, typename inverse_type>
 PreconditionBlockSOR<number,inverse_type>::~PreconditionBlockSOR(){}
 
 
+
+template<typename number, typename inverse_type>
+void
+PreconditionBlockSOR<number,inverse_type>::set_omega(number om)
+{
+  omega = om;
+}
+
+
+
+
 template <typename number, typename inverse_type>
 template <typename number2>
 void PreconditionBlockSOR<number,inverse_type>::operator() (Vector<number2>       &dst,