fixed a bug in set_matrix_mode.

STSetMatMode shall be called only after transformation_data->st
is initialised. Thus, we need to store the desired STMatMode
untill set_context is called from the solver object.

diff --git a/include/deal.II/lac/slepc_spectral_transformation.h b/include/deal.II/lac/slepc_spectral_transformation.h
index 848b6481235ab6f2564105d13186c98e45c83de8..6bf9a4c13bed0eb317b00c78c8e7498b23444f27 100644 (file)
--- a/include/deal.II/lac/slepc_spectral_transformation.h
+++ b/include/deal.II/lac/slepc_spectral_transformation.h
@@ -128,6 +128,8 @@ namespace SLEPcWrappers
     };
 
     std_cxx11::shared_ptr<TransformationData> transformation_data;
+
+    std_cxx11::shared_ptr<STMatMode> mat_mode;
   };
 
   /**

diff --git a/source/lac/slepc_spectral_transformation.cc b/source/lac/slepc_spectral_transformation.cc
index e01696234b28d6c2853fa86a08eba93e77bc397b..0cd4d82d4dd3e5e505c7e693fc14117ae337b35b 100644 (file)
--- a/source/lac/slepc_spectral_transformation.cc
+++ b/source/lac/slepc_spectral_transformation.cc
@@ -50,12 +50,19 @@ namespace SLEPcWrappers
     AssertThrow (ierr == 0, SolverBase::ExcSLEPcError(ierr));
 
     set_transformation_type(transformation_data->st);
+
+    // if mat_mode has been set,
+    // pass it to ST object
+    if (mat_mode.get() != 0)
+      {
+        int ierr = STSetMatMode(transformation_data->st,*(mat_mode.get()) );
+        AssertThrow (ierr == 0, SolverBase::ExcSLEPcError(ierr));
+      }
   }
 
   void TransformationBase::set_matrix_mode(const STMatMode mode)
   {
-    int ierr = STSetMatMode(transformation_data->st,mode);
-    AssertThrow (ierr == 0, SolverBase::ExcSLEPcError(ierr));
+    mat_mode.reset (new STMatMode(mode));
   }
 
   /* ------------------- TransformationShift --------------------- */