Merge pull request #4458 from nfehn/master

Initialization of Chebyshev smoother

diff --cc include/deal.II/lac/precondition.h
index 4b39bb9ec947157e9f15677dff20c4e0ce14e278,b8d833aee451762b45976f7a830e93657ea8c618..bed1220a3e908fcbcfc57c682b4a996174cbc470
--- 1/include/deal.II/lac/precondition.h
--- 2/include/deal.II/lac/precondition.h
+++ b/include/deal.II/lac/precondition.h
@@@ -2218,12 -2224,18 +2224,12 @@@ PreconditionChebyshev<MatrixType,Vector
          const VectorType &src) const
  {
    Threads::Mutex::ScopedLock lock(mutex);
-   if (is_initialized == false)
+   if (eigenvalues_are_initialized == false)
      estimate_eigenvalues(src);
  
 -  if (!dst.all_zero())
 -    {
 -      matrix_ptr->vmult (update2, dst);
 -      internal::PreconditionChebyshev::vector_updates
 -      (src, *data.preconditioner, false, 0., 1./theta, update1, update2, update3, dst);
 -    }
 -  else
 -    internal::PreconditionChebyshev::vector_updates
 -    (src, *data.preconditioner, true, 0., 1./theta, update1, update2, update3, dst);
 +  matrix_ptr->vmult (update2, dst);
 +  internal::PreconditionChebyshev::vector_updates
 +  (src, *data.preconditioner, false, 0., 1./theta, update1, update2, update3, dst);
  
    do_chebyshev_loop(dst, src);
  }
@@@ -2238,12 -2250,18 +2244,12 @@@ PreconditionChebyshev<MatrixType,Vector
           const VectorType &src) const
  {
    Threads::Mutex::ScopedLock lock(mutex);
-   if (is_initialized == false)
+   if (eigenvalues_are_initialized == false)
      estimate_eigenvalues(src);
  
 -  if (!dst.all_zero())
 -    {
 -      matrix_ptr->Tvmult (update2, dst);
 -      internal::PreconditionChebyshev::vector_updates
 -      (src, *data.preconditioner, false, 0., 1./theta, update1, update2, update3, dst);
 -    }
 -  else
 -    internal::PreconditionChebyshev::vector_updates
 -    (src, *data.preconditioner, true, 0., 1./theta, update1, update2, update3, dst);
 +  matrix_ptr->Tvmult (update2, dst);
 +  internal::PreconditionChebyshev::vector_updates
 +  (src, *data.preconditioner, false, 0., 1./theta, update1, update2, update3, dst);
  
    do_transpose_chebyshev_loop(dst, src);
  }