make sure vectors passed to EPSSetInitialSpace() are nonzero

diff --git a/include/deal.II/lac/slepc_solver.h b/include/deal.II/lac/slepc_solver.h
index ecb5cbf3291e3f69f57b4961a8759ef46a4269b8..47c46eaa8ef145e11a4b9109776c6c48663628e2 100644 (file)
--- a/include/deal.II/lac/slepc_solver.h
+++ b/include/deal.II/lac/slepc_solver.h
@@ -791,6 +791,8 @@ namespace SLEPcWrappers
 
     for (unsigned int i = 0; i < this_initial_space.size(); i++)
       {
+        Assert(this_initial_space[i].l2_norm()>0.0,
+               ExcMessage("Initial vectors should be nonzero."));
         vecs[i] = this_initial_space[i];
       }
 

diff --git a/source/lac/slepc_solver.cc b/source/lac/slepc_solver.cc
index 1a1d9ab707d638cee8e8ca42148afe08fdc7d7e5..d6495ac4691feaa23c0d0b89198c0c1b037e6cb8 100644 (file)
--- a/source/lac/slepc_solver.cc
+++ b/source/lac/slepc_solver.cc
@@ -100,6 +100,9 @@ namespace SLEPcWrappers
   void
   SolverBase::set_initial_vector (const PETScWrappers::VectorBase &this_initial_vector)
   {
+    Assert(this_initial_vector.l2_norm()>0.0,
+           ExcMessage("Initial vector should be nonzero."));
+
     int ierr;
     Vec vec = this_initial_vector;
 #if DEAL_II_PETSC_VERSION_LT(3,1,0)