#include <deal.II/base/utilities.h>
#include <deal.II/base/parallel.h>
#include <deal.II/base/template_constraints.h>
-#include <deal.II/lac/tridiagonal_matrix.h>
+#include <deal.II/base/thread_management.h>
+#include <deal.II/lac/diagonal_matrix.h>
#include <deal.II/lac/solver_cg.h>
#include <deal.II/lac/vector_memory.h>
/**
* Preconditioning with a Chebyshev polynomial for symmetric positive definite
- * matrices. This preconditioner is similar to a Jacobi preconditioner if the
- * degree variable is set to one, otherwise some higher order polynomial
- * corrections are used. This preconditioner needs access to the diagonal of
- * the matrix it acts on and needs a respective <tt>vmult</tt> implementation.
- * However, it does not need to explicitly know the matrix entries.
+ * matrices. This preconditioner is based on an iteration of an inner
+ * preconditioner of type @p PreconditionerType with coefficients that are
+ * adapted to optimally cover an eigenvalue range between the largest
+ * eigenvalue down to a given lower eigenvalue specified by the optional
+ * parameter @p smoothing_range. The typical use case for the preconditioner
+ * is a Jacobi preconditioner specified through DiagonalMatrix, which is also
+ * the default value for the preconditioner. Note that if the degree variable
+ * is set to zero, the Chebyshev iteration corresponds to a Jacobi
+ * preconditioner (or the underlying preconditioner type) with relaxation
+ * parameter according to the specified smoothing range.
*
- * This class is useful e.g. in multigrid smoother objects, since it is
- * trivially %parallel (assuming that matrix-vector products are %parallel).
+ * Besides the default choice of a pointwise Jacobi preconditioner, this class
+ * also allows for more advanced types of preconditioners, for example
+ * iterating block-Jacobi preconditioners in DG methods.
*
- * @author Martin Kronbichler, 2009; extension for full compatibility with
+ * Apart from the inner preconditioner object, this iteration does not need
+ * access to matrix entries, which makes it an ideal ingredient for
+ * matrix-free computations. In that context, this class can be used as a
+ * multigrid smoother that is trivially %parallel (assuming that matrix-vector
+ * products are %parallel and the inner preconditioner is %parallel). Its use
+ * is demonstrated in the step-37 tutorial program.
+ *
+ * <h4>Algorithm execution</h4>
+ *
+ * The Chebyshev method relies on an estimate of the eigenvalues of the matrix
+ * which are computed during the first invocation of vmult(). The algorithm
+ * invokes a conjugate gradient solver so symmetry and positive definiteness
+ * of the (preconditioned) matrix system are strong requirements. The
+ * computation of eigenvalues needs to be deferred until the first vmult()
+ * invocation because temporary vectors of the same layout as the source and
+ * destination vectors are necessary for these computations and this
+ * information gets only available through vmult().
+ *
+ * The estimation of eigenvalues can also be bypassed by setting
+ * PreconditionChebyshev::AdditionalData::eig_cg_n_iterations to zero and
+ * providing sensible values for the largest eigenvalues in the field
+ * PreconditionChebyshev::AdditionalData::max_eigenvalue. If the range
+ * <tt>[max_eigenvalue/smoothing_range, max_eigenvalue]</tt> contains all
+ * eigenvalues of the preconditioned matrix system and the degree (i.e.,
+ * number of iterations) is high enough, this class can also be used as a
+ * direct solver. For an error estimation of the Chebyshev iteration that can
+ * be used to determine the number of iteration, see e.g. section 5.1 of
+ * R. S. Varga, Matrix iterative analysis, 2nd ed., Springer, 2009.
+ *
+ * <h4>Requirements on the templated classes</h4>
+ *
+ * The class MatrixType must be derived from Subscriptor because a
+ * SmartPointer to MatrixType is held in the class. In particular, this means
+ * that the matrix object needs to persist during the lifetime of
+ * PreconditionChebyshev. The preconditioner is held in a shared_ptr that is
+ * copied into the AdditionalData member variable of the class, so the
+ * variable used for initialization can safely be discarded after calling
+ * initialize(). Both the matrix and the preconditioner need to provide @p
+ * vmult functions for the matrix-vector product and @p m functions for
+ * accessing the number of rows in the (square) matrix. Furthermore, the
+ * matrix must provide <tt>el(i,i)</tt> methods for accessing the matrix
+ * diagonal in case the preconditioner type is a diagonal matrix. Even though
+ * it is highly recommended to pass the inverse diagonal entries inside a
+ * separate preconditioner object for implementing the Jacobi method (and the
+ * only possibility when computing in %parallel with MPI), there is a backward
+ * compatibility function that can extract the diagonal in case of a serial
+ * computation.
+ *
+ * @author Martin Kronbichler, 2009, 2016; extension for full compatibility with
* LinearOperator class: Jean-Paul Pelteret, 2015
*/
-template <typename MatrixType=SparseMatrix<double>, class VectorType=Vector<double> >
+template <typename MatrixType=SparseMatrix<double>, typename VectorType=Vector<double>, typename PreconditionerType=DiagonalMatrix<VectorType> >
class PreconditionChebyshev : public Subscriptor
{
public:
/**
* Stores the inverse of the diagonal of the underlying matrix.
+ *
+ * @deprecated Set the variable @p preconditioner defined below instead.
*/
VectorType matrix_diagonal_inverse;
+
+ /**
+ * Stores the preconditioner object that the Chebyshev is wrapped around.
+ */
+ std_cxx11::shared_ptr<PreconditionerType> preconditioner;
};
PreconditionChebyshev ();
/**
* A pointer to the underlying matrix.
*/
- SmartPointer<const MatrixType,PreconditionChebyshev<MatrixType,VectorType> > matrix_ptr;
+ SmartPointer<const MatrixType,PreconditionChebyshev<MatrixType,VectorType,PreconditionerType> > matrix_ptr;
/**
* Internal vector used for the <tt>vmult</tt> operation.
mutable VectorType update2;
/**
- * Stores the additional data provided to the initialize function.
+ * Internal vector used for the <tt>vmult</tt> operation.
+ */
+ mutable VectorType update3;
+
+ /**
+ * Stores the additional data passed to the initialize function, obtained
+ * through a copy operation.
*/
AdditionalData data;
double delta;
/**
- * Stores whether the preconditioner has been set up.
+ * Stores whether the preconditioner has been set up and eigenvalues have
+ * been computed.
*/
bool is_initialized;
+ /**
+ * A mutex to avoid that multiple vmult() invocations by different threads
+ * overwrite the temporary vectors.
+ */
+ mutable Threads::Mutex mutex;
+
/**
* Runs the inner loop of the Chebyshev preconditioner that is the same for
* vmult() and step() methods.
*/
void do_transpose_chebyshev_loop(VectorType &dst,
const VectorType &src) const;
+
+ /**
+ * Initializes the factors theta and delta based on an eigenvalue
+ * computation. If the user set provided values for the largest eigenvalue
+ * in AdditionalData, no computation is performed and the information given
+ * by the user is used.
+ */
+ void estimate_eigenvalues(const VectorType &src) const;
};
// size of the vector
// generic part for non-deal.II vectors
- template <typename VectorType>
+ template <typename VectorType, typename PreconditionerType>
inline
void
- vector_updates (const VectorType &src,
- const VectorType &matrix_diagonal_inverse,
- const bool start_zero,
- const double factor1,
- const double factor2,
- VectorType &update1,
- VectorType &update2,
- VectorType &dst)
+ vector_updates (const VectorType &src,
+ const PreconditionerType &preconditioner,
+ const bool start_zero,
+ const double factor1,
+ const double factor2,
+ VectorType &update1,
+ VectorType &update2,
+ VectorType &update3,
+ VectorType &dst)
{
if (start_zero)
{
- dst.equ (factor2, src);
- dst.scale (matrix_diagonal_inverse);
+ update1.equ (factor2, src);
+ preconditioner.vmult(dst, update1);
update1.equ(-1.,dst);
}
else
{
update2 -= src;
- update2.scale (matrix_diagonal_inverse);
+ preconditioner.vmult(update3, update2);
+ update2 = update3;
if (factor1 == 0.)
update1.equ(factor2, update2);
else
const VectorUpdater<Number> &updater;
};
- // selection for deal.II vector
+ // selection for diagonal matrix around deal.II vector
template <typename Number>
inline
void
vector_updates (const ::dealii::Vector<Number> &src,
- const ::dealii::Vector<Number> &matrix_diagonal_inverse,
+ const DiagonalMatrix<::dealii::Vector<Number> > &jacobi,
const bool start_zero,
const double factor1,
const double factor2,
::dealii::Vector<Number> &update1,
::dealii::Vector<Number> &update2,
+ ::dealii::Vector<Number> &,
::dealii::Vector<Number> &dst)
{
- VectorUpdater<Number> upd(src.begin(), matrix_diagonal_inverse.begin(),
+ VectorUpdater<Number> upd(src.begin(), jacobi.get_vector().begin(),
start_zero, factor1, factor2,
update1.begin(), update2.begin(), dst.begin());
VectorUpdatesRange<Number>(upd, src.size());
}
- // selection for parallel deal.II vector
+ // selection for diagonal matrix around parallel deal.II vector
template <typename Number>
inline
void
vector_updates (const LinearAlgebra::distributed::Vector<Number> &src,
- const LinearAlgebra::distributed::Vector<Number> &matrix_diagonal_inverse,
+ const DiagonalMatrix<LinearAlgebra::distributed::Vector<Number> > &jacobi,
const bool start_zero,
const double factor1,
const double factor2,
LinearAlgebra::distributed::Vector<Number> &update1,
LinearAlgebra::distributed::Vector<Number> &update2,
+ LinearAlgebra::distributed::Vector<Number> &,
LinearAlgebra::distributed::Vector<Number> &dst)
{
- VectorUpdater<Number> upd(src.begin(), matrix_diagonal_inverse.begin(),
+ VectorUpdater<Number> upd(src.begin(), jacobi.get_vector().begin(),
start_zero, factor1, factor2,
update1.begin(), update2.begin(), dst.begin());
VectorUpdatesRange<Number>(upd, src.local_size());
}
- template <typename VectorType>
- struct DiagonalPreconditioner
+ template <typename MatrixType, typename VectorType, typename PreconditionerType>
+ inline
+ void
+ initialize_preconditioner(const MatrixType &matrix,
+ std_cxx11::shared_ptr<PreconditionerType> &preconditioner,
+ VectorType &)
{
- DiagonalPreconditioner (const VectorType &vector)
- :
- diagonal_vector(vector)
- {}
+ (void)matrix;
+ (void)preconditioner;
+ AssertThrow(preconditioner.get() != NULL, ExcNotInitialized());
+ AssertDimension(preconditioner->m(), matrix.m());
+ }
- void vmult (VectorType &dst,
- const VectorType &src) const
- {
- dst = src;
- dst.scale(diagonal_vector);
- }
+ template <typename MatrixType, typename VectorType>
+ inline
+ void
+ initialize_preconditioner(const MatrixType &matrix,
+ std_cxx11::shared_ptr<DiagonalMatrix<VectorType> > &preconditioner,
+ VectorType &diagonal_inverse)
+ {
+ if (preconditioner.get() == NULL ||
+ preconditioner->m() != matrix.m())
+ {
+ if (preconditioner.get() == NULL)
+ preconditioner.reset(new DiagonalMatrix<VectorType>());
- const VectorType &diagonal_vector;
- };
+ Assert(preconditioner->m() == 0,
+ ExcMessage("Preconditioner appears to be initialized but not sized correctly"));
+
+ // Check if we can initialize from vector that then gets set to zero
+ // as the matrix will own the memory
+ preconditioner->reinit(diagonal_inverse);
+ diagonal_inverse.reinit(0);
+
+ // This part only works in serial
+ if (preconditioner->m() != matrix.m())
+ {
+ preconditioner->get_vector().reinit(matrix.m());
+ for (typename VectorType::size_type i=0; i<matrix.m(); ++i)
+ preconditioner->get_vector()(i) = 1./matrix.el(i,i);
+ }
+ }
+ }
struct EigenvalueTracker
{
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, class VectorType, typename PreconditionerType>
inline
-PreconditionChebyshev<MatrixType,VectorType>::AdditionalData::
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::AdditionalData::
AdditionalData (const unsigned int degree,
const double smoothing_range,
const bool nonzero_starting,
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
inline
-PreconditionChebyshev<MatrixType,VectorType>::PreconditionChebyshev ()
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::PreconditionChebyshev ()
:
is_initialized (false)
{
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
inline
void
-PreconditionChebyshev<MatrixType,VectorType>::initialize
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::initialize
(const MatrixType &matrix,
const AdditionalData &additional_data)
{
matrix_ptr = &matrix;
data = additional_data;
- if (data.matrix_diagonal_inverse.size() != matrix.m())
- {
- Assert(data.matrix_diagonal_inverse.size() == 0,
- ExcMessage("Matrix diagonal vector set but not sized correctly"));
- data.matrix_diagonal_inverse.reinit(matrix.m());
- for (unsigned int i=0; i<matrix.m(); ++i)
- data.matrix_diagonal_inverse(i) = 1./matrix.el(i,i);
- }
+ internal::PreconditionChebyshev::initialize_preconditioner(matrix,
+ data.preconditioner,
+ data.matrix_diagonal_inverse);
+ is_initialized = false;
+}
+
+
+
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
+inline
+void
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::estimate_eigenvalues
+(const VectorType &src) const
+{
+ Assert(is_initialized == false, ExcInternalError());
+ Assert(data.preconditioner.get() != NULL, ExcNotInitialized());
+ update1.reinit(src);
+ update2 = src;
// calculate largest eigenvalue using a hand-tuned CG iteration on the
// matrix weighted by its diagonal. we start with a vector that consists of
double max_eigenvalue, min_eigenvalue;
if (data.eig_cg_n_iterations > 0)
{
- Assert (additional_data.eig_cg_n_iterations > 2,
+ Assert (data.eig_cg_n_iterations > 2,
ExcMessage ("Need to set at least two iterations to find eigenvalues."));
// set a very strict tolerance to force at least two iterations
- ReductionControl control (data.eig_cg_n_iterations, 1e-35, 1e-10);
- GrowingVectorMemory<VectorType> memory;
- VectorType *rhs = memory.alloc();
- VectorType *dummy = memory.alloc();
- rhs->reinit(data.matrix_diagonal_inverse);
- dummy->reinit(data.matrix_diagonal_inverse);
-
- // heuristically, a right hand side close to a constant has been shown
- // to quickly reveal the largest eigenvalue. however, avoid to use the
- // exact constant because that might be not in the range space of some
- // matrices (purely Neumann matrices with constant mode filtered out by
- // orthogonal projection in the matrix-vector product)
- *rhs = 1./std::sqrt(static_cast<double>(matrix.m()));
- if (rhs->locally_owned_elements().is_element(0))
- (*rhs)(0) = 0.;
- rhs->compress(VectorOperation::insert);
+ ReductionControl control (data.eig_cg_n_iterations,
+ std::sqrt(std::numeric_limits<typename VectorType::value_type>::epsilon()),
+ 1e-10, false, false);
internal::PreconditionChebyshev::EigenvalueTracker eigenvalue_tracker;
- SolverCG<VectorType> solver (control, memory);
+ SolverCG<VectorType> solver (control);
solver.connect_eigenvalues_slot(std_cxx11::bind(&internal::PreconditionChebyshev::EigenvalueTracker::slot,
&eigenvalue_tracker,
std_cxx11::_1));
- internal::PreconditionChebyshev::DiagonalPreconditioner<VectorType>
- preconditioner(data.matrix_diagonal_inverse);
+
+ // catch case the rhs vector is exactly zero (no eigenvalue estimate
+ // possible) by setting a different initial guess
+ if (update2 * update2 == typename VectorType::value_type())
+ {
+ update2 = 1.;
+ if (update2.locally_owned_elements().is_element(0))
+ update2(0) = 0.;
+ }
+
try
{
- solver.solve(matrix, *dummy, *rhs, preconditioner);
+ solver.solve(*matrix_ptr, update1, update2, *data.preconditioner);
}
catch (SolverControl::NoConvergence &)
{
}
- memory.free(dummy);
- memory.free(rhs);
-
// read the eigenvalues from the attached eigenvalue tracker
if (eigenvalue_tracker.values.empty())
min_eigenvalue = max_eigenvalue = 1;
else
{
min_eigenvalue = eigenvalue_tracker.values.front();
- max_eigenvalue = eigenvalue_tracker.values.back();
- }
- // include a safety factor since the CG method will in general not be
- // converged
- max_eigenvalue *= 1.2;
+ // include a safety factor since the CG method will in general not be
+ // converged
+ max_eigenvalue = 1.2*eigenvalue_tracker.values.back();
+ }
}
else
{
const double alpha = (data.smoothing_range > 1. ?
max_eigenvalue / data.smoothing_range :
std::min(0.9*max_eigenvalue, min_eigenvalue));
- delta = (max_eigenvalue-alpha)*0.5;
- theta = (max_eigenvalue+alpha)*0.5;
-
- update1.reinit (data.matrix_diagonal_inverse, true);
- update2.reinit (data.matrix_diagonal_inverse, true);
-
- is_initialized = true;
+ const_cast<PreconditionChebyshev<MatrixType,VectorType,PreconditionerType> *>(this)->delta = (max_eigenvalue-alpha)*0.5;
+ const_cast<PreconditionChebyshev<MatrixType,VectorType,PreconditionerType> *>(this)->theta = (max_eigenvalue+alpha)*0.5;
+
+ // We do not need the third auxiliary vector in case we have a
+ // DiagonalMatrix as preconditioner and use deal.II's own vectors
+ if (types_are_equal<PreconditionerType,DiagonalMatrix<VectorType> >::value == false ||
+ (types_are_equal<VectorType,dealii::Vector<typename VectorType::value_type> >::value == false
+ &&
+ types_are_equal<VectorType,LinearAlgebra::distributed::Vector<typename VectorType::value_type> >::value == false
+ ))
+ update3.reinit (src, true);
+
+ const_cast<PreconditionChebyshev<MatrixType,VectorType,PreconditionerType> *>(this)->is_initialized = true;
}
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
inline
void
-PreconditionChebyshev<MatrixType,VectorType>::do_chebyshev_loop(VectorType &dst,
- const VectorType &src) const
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>
+::do_chebyshev_loop(VectorType &dst,
+ const VectorType &src) const
{
+ // if delta is zero, we do not need to iterate because the updates will be
+ // zero
+ if (std::abs(delta) < 1e-40)
+ return;
+
double rhok = delta / theta, sigma = theta / delta;
for (unsigned int k=0; k<data.degree; ++k)
{
const double factor1 = rhokp * rhok, factor2 = 2.*rhokp/delta;
rhok = rhokp;
internal::PreconditionChebyshev::vector_updates
- (src, data.matrix_diagonal_inverse, false, factor1, factor2, update1,
- update2, dst);
+ (src, *data.preconditioner, false, factor1, factor2, update1, update2, update3, dst);
}
}
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
inline
void
-PreconditionChebyshev<MatrixType,VectorType>::do_transpose_chebyshev_loop(VectorType &dst,
- const VectorType &src) const
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>
+::do_transpose_chebyshev_loop(VectorType &dst,
+ const VectorType &src) const
{
double rhok = delta / theta, sigma = theta / delta;
for (unsigned int k=0; k<data.degree; ++k)
const double factor1 = rhokp * rhok, factor2 = 2.*rhokp/delta;
rhok = rhokp;
internal::PreconditionChebyshev::vector_updates
- (src, data.matrix_diagonal_inverse, false, factor1, factor2, update1,
- update2, dst);
+ (src, *data.preconditioner, false, factor1, factor2, update1, update2, update3, dst);
}
}
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
inline
void
-PreconditionChebyshev<MatrixType,VectorType>::vmult (VectorType &dst,
- const VectorType &src) const
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>
+::vmult (VectorType &dst,
+ const VectorType &src) const
{
- Assert (is_initialized, ExcMessage("Preconditioner not initialized"));
+ Threads::Mutex::ScopedLock lock(mutex);
+ if (is_initialized == false)
+ estimate_eigenvalues(src);
+
internal::PreconditionChebyshev::vector_updates
- (src, data.matrix_diagonal_inverse, true, 0., 1./theta, update1,
- update2, dst);
+ (src, *data.preconditioner, true, 0., 1./theta, update1, update2, update3, dst);
do_chebyshev_loop(dst, src);
}
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
inline
void
-PreconditionChebyshev<MatrixType,VectorType>::Tvmult (VectorType &dst,
- const VectorType &src) const
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>
+::Tvmult (VectorType &dst,
+ const VectorType &src) const
{
- Assert (is_initialized, ExcMessage("Preconditioner not initialized"));
+ Threads::Mutex::ScopedLock lock(mutex);
+ if (is_initialized == false)
+ estimate_eigenvalues(src);
+
internal::PreconditionChebyshev::vector_updates
- (src, data.matrix_diagonal_inverse, true, 0., 1./theta, update1,
- update2, dst);
+ (src, *data.preconditioner, true, 0., 1./theta, update1, update2, update3, dst);
do_transpose_chebyshev_loop(dst, src);
}
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
inline
void
-PreconditionChebyshev<MatrixType,VectorType>::step (VectorType &dst,
- const VectorType &src) const
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>
+::step (VectorType &dst,
+ const VectorType &src) const
{
- Assert (is_initialized, ExcMessage("Preconditioner not initialized"));
+ Threads::Mutex::ScopedLock lock(mutex);
+ if (is_initialized == false)
+ estimate_eigenvalues(src);
+
if (!dst.all_zero())
{
matrix_ptr->vmult (update2, dst);
internal::PreconditionChebyshev::vector_updates
- (src, data.matrix_diagonal_inverse, false, 0., 1./theta, update1,
- update2, dst);
+ (src, *data.preconditioner, false, 0., 1./theta, update1, update2, update3, dst);
}
else
internal::PreconditionChebyshev::vector_updates
- (src, data.matrix_diagonal_inverse, true, 0., 1./theta, update1,
- update2, dst);
+ (src, *data.preconditioner, true, 0., 1./theta, update1, update2, update3, dst);
do_chebyshev_loop(dst, src);
}
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
inline
void
-PreconditionChebyshev<MatrixType,VectorType>::Tstep (VectorType &dst,
- const VectorType &src) const
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>
+::Tstep (VectorType &dst,
+ const VectorType &src) const
{
- Assert (is_initialized, ExcMessage("Preconditioner not initialized"));
+ Threads::Mutex::ScopedLock lock(mutex);
+ if (is_initialized == false)
+ estimate_eigenvalues(src);
+
if (!dst.all_zero())
{
matrix_ptr->Tvmult (update2, dst);
internal::PreconditionChebyshev::vector_updates
- (src, data.matrix_diagonal_inverse, false, 0., 1./theta, update1,
- update2, dst);
+ (src, *data.preconditioner, false, 0., 1./theta, update1, update2, update3, dst);
}
else
internal::PreconditionChebyshev::vector_updates
- (src, data.matrix_diagonal_inverse, true, 0., 1./theta, update1,
- update2, dst);
+ (src, *data.preconditioner, true, 0., 1./theta, update1, update2, update3, dst);
do_transpose_chebyshev_loop(dst, src);
}
-template <typename MatrixType, typename VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
inline
-void PreconditionChebyshev<MatrixType,VectorType>::clear ()
+void
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::clear ()
{
is_initialized = false;
matrix_ptr = 0;
data.matrix_diagonal_inverse.reinit(0);
+ data.preconditioner.reset();
update1.reinit(0);
update2.reinit(0);
+ update3.reinit(0);
}
-template <typename MatrixType, typename VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
inline
-typename PreconditionChebyshev<MatrixType,VectorType>::size_type
-PreconditionChebyshev<MatrixType,VectorType>::m () const
+typename PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::size_type
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::m () const
{
Assert (matrix_ptr!=0, ExcNotInitialized());
return matrix_ptr->m();
}
-template <typename MatrixType, typename VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
inline
-typename PreconditionChebyshev<MatrixType,VectorType>::size_type
-PreconditionChebyshev<MatrixType,VectorType>::n () const
+typename PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::size_type
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::n () const
{
Assert (matrix_ptr!=0, ExcNotInitialized());
return matrix_ptr->n();
DEAL::Preconditioners
DEAL:Preconditioners::PreconditionChebyshev
-DEAL:Preconditioners:cg::Starting value 0.9486832981
-DEAL:Preconditioners:cg::Convergence step 1 value 0.000000000
DEAL:Preconditioners:Exact inverse:cg::Starting value 16.88194302
DEAL:Preconditioners:Exact inverse:cg::Convergence step 1 value 2.601335705e-15
DEAL:Preconditioners::PreconditionIdentity
DEAL::Preconditioners
DEAL:Preconditioners::PreconditionChebyshev
-DEAL:Preconditioners:cg::Starting value 0.9486832981
-DEAL:Preconditioners:cg::Convergence step 1 value 0.000000000
DEAL:Preconditioners:Exact inverse:cg::Starting value 16.88194302
DEAL:Preconditioners:Exact inverse:cg::Convergence step 1 value 2.601335705e-15
DEAL:Preconditioners::PreconditionIdentity
-DEAL:cg::Starting value 0.95
-DEAL:cg::Convergence step 1 value 0
DEAL::Exact inverse: 0.84 0.20 0.26 0.20 0.18 0.03 0.05 0.10 0.03 0.06
-DEAL::Check vmult orig: 1.09 0.26 0.34 0.26 0.24 0.04 0.06 0.12 0.04 0.07
-DEAL::Check Tvmult orig: 1.09 0.26 0.34 0.26 0.24 0.04 0.06 0.12 0.04 0.07
-DEAL:cg::Starting value 0.95
-DEAL:cg::Failure step 8 value 0.00
+DEAL::Check vmult orig: 0.57 0.13 0.18 0.14 0.12 0.02 0.03 0.07 0.02 0.04
+DEAL::Check Tvmult orig: 0.57 0.13 0.18 0.14 0.12 0.02 0.03 0.07 0.02 0.04
DEAL::Check vmult diag: 0.82 0.26 0.33 0.21 0.15 0.02 0.03 0.09 0.03 0.07
DEAL::Check Tvmult diag: 0.82 0.26 0.33 0.21 0.15 0.02 0.03 0.09 0.03 0.07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 81
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
DEAL:2d:cg::Starting value 9.000
-DEAL:2d:cg:cg::Starting value 0.004670
+DEAL:2d:cg:cg::Starting value 0.009355
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 3.232e-06
+DEAL:2d:cg:cg::Starting value 6.428e-06
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 3.457e-08
+DEAL:2d:cg:cg::Starting value 6.896e-08
DEAL:2d:cg:cg::Convergence step 1 value 0
DEAL:2d:cg::Convergence step 3 value 4.626e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 289
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
DEAL:2d:cg::Starting value 17.00
-DEAL:2d:cg:cg::Starting value 0.01350
+DEAL:2d:cg:cg::Starting value 0.02710
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 2.583e-06
+DEAL:2d:cg:cg::Starting value 5.173e-06
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 4.246e-10
+DEAL:2d:cg:cg::Starting value 7.863e-10
DEAL:2d:cg:cg::Convergence step 1 value 0
DEAL:2d:cg::Convergence step 3 value 3.600e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 1089
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9995
-DEAL:2d:cg::Failure step 15 value 0.1302
DEAL:2d:cg::Starting value 33.00
-DEAL:2d:cg:cg::Starting value 0.04871
+DEAL:2d:cg:cg::Starting value 0.09750
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 1.396e-06
+DEAL:2d:cg:cg::Starting value 2.792e-06
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 5.548e-08
+DEAL:2d:cg:cg::Starting value 1.110e-07
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 8.781e-11
+DEAL:2d:cg:cg::Starting value 1.758e-10
DEAL:2d:cg:cg::Convergence step 1 value 0
DEAL:2d:cg::Convergence step 4 value 1.484e-08
DEAL:2d::Testing FE_Q<2>(2)
DEAL:2d::Number of degrees of freedom: 289
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 11 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 0.001560
DEAL:2d:cg::Starting value 17.00
-DEAL:2d:cg:cg::Starting value 0.5165
+DEAL:2d:cg:cg::Starting value 0.6907
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 7.798e-05
+DEAL:2d:cg:cg::Starting value 0.0001064
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 4.794e-07
+DEAL:2d:cg:cg::Starting value 6.439e-07
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 1.177e-09
+DEAL:2d:cg:cg::Starting value 1.579e-09
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg::Convergence step 4 value 1.404e-08
+DEAL:2d:cg::Convergence step 4 value 1.405e-08
DEAL:2d::Testing FE_Q<2>(2)
DEAL:2d::Number of degrees of freedom: 1089
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 11 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 0.001560
-DEAL:2d:cg::Starting value 0.9995
-DEAL:2d:cg::Failure step 15 value 0.5714
DEAL:2d:cg::Starting value 33.00
-DEAL:2d:cg:cg::Starting value 1.948
+DEAL:2d:cg:cg::Starting value 2.631
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 7.876e-05
+DEAL:2d:cg:cg::Starting value 0.0001068
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 9.527e-07
+DEAL:2d:cg:cg::Starting value 1.290e-06
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 1.973e-09
+DEAL:2d:cg:cg::Starting value 2.685e-09
DEAL:2d:cg:cg::Convergence step 2 value 0
DEAL:2d:cg::Convergence step 4 value 1.627e-08
DEAL:2d::Testing FE_Q<2>(2)
DEAL:2d::Number of degrees of freedom: 4225
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 11 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 0.001560
-DEAL:2d:cg::Starting value 0.9995
-DEAL:2d:cg::Failure step 15 value 0.5714
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 2.322
DEAL:2d:cg::Starting value 65.00
-DEAL:2d:cg:cg::Starting value 7.684
+DEAL:2d:cg:cg::Starting value 10.37
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 0.0004225
+DEAL:2d:cg:cg::Starting value 0.0005721
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 5.550e-07
+DEAL:2d:cg:cg::Starting value 7.479e-07
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 8.109e-09
+DEAL:2d:cg:cg::Starting value 1.098e-08
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg::Convergence step 4 value 3.009e-08
+DEAL:2d:cg::Convergence step 4 value 3.011e-08
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 125
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
DEAL:3d:cg::Starting value 11.18
-DEAL:3d:cg:cg::Starting value 0.07042
+DEAL:3d:cg:cg::Starting value 0.1410
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 3.438e-05
+DEAL:3d:cg:cg::Starting value 6.992e-05
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 2.189e-09
+DEAL:3d:cg:cg::Starting value 4.858e-09
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg::Convergence step 3 value 8.878e-09
+DEAL:3d:cg::Convergence step 3 value 9.022e-09
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 729
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
DEAL:3d:cg::Starting value 27.00
-DEAL:3d:cg:cg::Starting value 0.01748
+DEAL:3d:cg:cg::Starting value 0.03499
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0001434
+DEAL:3d:cg:cg::Starting value 0.0002870
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 7.360e-08
+DEAL:3d:cg:cg::Starting value 1.473e-07
DEAL:3d:cg:cg::Convergence step 1 value 0
DEAL:3d:cg::Convergence step 3 value 1.327e-06
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 4913
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 2.351e-05
DEAL:3d:cg::Starting value 70.09
-DEAL:3d:cg:cg::Starting value 0.03090
+DEAL:3d:cg:cg::Starting value 0.06180
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 3.917e-05
+DEAL:3d:cg:cg::Starting value 7.869e-05
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 5.305e-07
+DEAL:3d:cg:cg::Starting value 1.067e-06
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 7.196e-10
+DEAL:3d:cg:cg::Starting value 1.447e-09
DEAL:3d:cg:cg::Convergence step 1 value 0
DEAL:3d:cg::Convergence step 4 value 6.374e-08
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 729
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
DEAL:3d:cg::Starting value 27.00
DEAL:3d:cg:cg::Starting value 1.155
DEAL:3d:cg:cg::Convergence step 2 value 0
-DEAL:3d:cg:cg::Starting value 0.0001365
+DEAL:3d:cg:cg::Starting value 0.0001367
DEAL:3d:cg:cg::Convergence step 2 value 0
DEAL:3d:cg:cg::Starting value 1.471e-07
DEAL:3d:cg:cg::Convergence step 2 value 0
DEAL:3d:cg::Convergence step 3 value 5.042e-08
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 4913
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
DEAL:3d:cg::Starting value 70.09
-DEAL:3d:cg:cg::Starting value 6.474
+DEAL:3d:cg:cg::Starting value 6.475
DEAL:3d:cg:cg::Convergence step 2 value 0
DEAL:3d:cg:cg::Starting value 0.0005841
DEAL:3d:cg:cg::Convergence step 2 value 0
DEAL:3d:cg::Convergence step 3 value 2.560e-06
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 35937
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
-DEAL:3d:cg::Starting value 1.000
-DEAL:3d:cg::Failure step 15 value 0.3991
DEAL:3d:cg::Starting value 189.6
DEAL:3d:cg:cg::Starting value 47.62
DEAL:3d:cg:cg::Convergence step 2 value 0
-DEAL:3d:cg:cg::Starting value 0.006356
+DEAL:3d:cg:cg::Starting value 0.006358
DEAL:3d:cg:cg::Convergence step 2 value 0
-DEAL:3d:cg:cg::Starting value 8.271e-06
+DEAL:3d:cg:cg::Starting value 8.274e-06
DEAL:3d:cg:cg::Convergence step 2 value 0
DEAL:3d:cg::Convergence step 3 value 1.007e-05
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 81
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
DEAL:2d:cg::Starting value 9.000
-DEAL:2d:cg:cg::Starting value 0.004670
+DEAL:2d:cg:cg::Starting value 0.009355
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 3.232e-06
+DEAL:2d:cg:cg::Starting value 6.428e-06
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 3.457e-08
+DEAL:2d:cg:cg::Starting value 6.896e-08
DEAL:2d:cg:cg::Convergence step 1 value 0
DEAL:2d:cg::Convergence step 3 value 4.626e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 289
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
DEAL:2d:cg::Starting value 17.00
-DEAL:2d:cg:cg::Starting value 0.01350
+DEAL:2d:cg:cg::Starting value 0.02710
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 2.583e-06
+DEAL:2d:cg:cg::Starting value 5.173e-06
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 4.246e-10
+DEAL:2d:cg:cg::Starting value 7.863e-10
DEAL:2d:cg:cg::Convergence step 1 value 0
DEAL:2d:cg::Convergence step 3 value 3.600e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 1089
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9995
-DEAL:2d:cg::Failure step 15 value 0.1302
DEAL:2d:cg::Starting value 33.00
-DEAL:2d:cg:cg::Starting value 0.04871
+DEAL:2d:cg:cg::Starting value 0.09750
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 1.396e-06
+DEAL:2d:cg:cg::Starting value 2.792e-06
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 5.548e-08
+DEAL:2d:cg:cg::Starting value 1.110e-07
DEAL:2d:cg:cg::Convergence step 1 value 0
-DEAL:2d:cg:cg::Starting value 8.781e-11
+DEAL:2d:cg:cg::Starting value 1.758e-10
DEAL:2d:cg:cg::Convergence step 1 value 0
DEAL:2d:cg::Convergence step 4 value 1.484e-08
DEAL:2d::Testing FE_Q<2>(2)
DEAL:2d::Number of degrees of freedom: 289
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 11 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 0.001560
DEAL:2d:cg::Starting value 17.00
-DEAL:2d:cg:cg::Starting value 0.5165
+DEAL:2d:cg:cg::Starting value 0.6907
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 7.798e-05
+DEAL:2d:cg:cg::Starting value 0.0001064
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 4.794e-07
+DEAL:2d:cg:cg::Starting value 6.439e-07
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 1.177e-09
+DEAL:2d:cg:cg::Starting value 1.579e-09
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg::Convergence step 4 value 1.404e-08
+DEAL:2d:cg::Convergence step 4 value 1.405e-08
DEAL:2d::Testing FE_Q<2>(2)
DEAL:2d::Number of degrees of freedom: 1089
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 11 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 0.001560
-DEAL:2d:cg::Starting value 0.9995
-DEAL:2d:cg::Failure step 15 value 0.5714
DEAL:2d:cg::Starting value 33.00
-DEAL:2d:cg:cg::Starting value 1.948
+DEAL:2d:cg:cg::Starting value 2.631
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 7.876e-05
+DEAL:2d:cg:cg::Starting value 0.0001068
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 9.527e-07
+DEAL:2d:cg:cg::Starting value 1.290e-06
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 1.973e-09
+DEAL:2d:cg:cg::Starting value 2.685e-09
DEAL:2d:cg:cg::Convergence step 2 value 0
DEAL:2d:cg::Convergence step 4 value 1.627e-08
DEAL:2d::Testing FE_Q<2>(2)
DEAL:2d::Number of degrees of freedom: 4225
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 11 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 0.001560
-DEAL:2d:cg::Starting value 0.9995
-DEAL:2d:cg::Failure step 15 value 0.5714
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 2.322
DEAL:2d:cg::Starting value 65.00
-DEAL:2d:cg:cg::Starting value 7.684
+DEAL:2d:cg:cg::Starting value 10.37
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 0.0004225
+DEAL:2d:cg:cg::Starting value 0.0005721
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 5.550e-07
+DEAL:2d:cg:cg::Starting value 7.479e-07
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg:cg::Starting value 8.109e-09
+DEAL:2d:cg:cg::Starting value 1.098e-08
DEAL:2d:cg:cg::Convergence step 2 value 0
-DEAL:2d:cg::Convergence step 4 value 3.009e-08
+DEAL:2d:cg::Convergence step 4 value 3.011e-08
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 125
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
DEAL:3d:cg::Starting value 11.18
-DEAL:3d:cg:cg::Starting value 0.07042
+DEAL:3d:cg:cg::Starting value 0.1410
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 3.438e-05
+DEAL:3d:cg:cg::Starting value 6.992e-05
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 2.189e-09
+DEAL:3d:cg:cg::Starting value 4.858e-09
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg::Convergence step 3 value 8.878e-09
+DEAL:3d:cg::Convergence step 3 value 9.022e-09
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 729
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
DEAL:3d:cg::Starting value 27.00
-DEAL:3d:cg:cg::Starting value 0.01748
+DEAL:3d:cg:cg::Starting value 0.03499
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0001434
+DEAL:3d:cg:cg::Starting value 0.0002870
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 7.360e-08
+DEAL:3d:cg:cg::Starting value 1.473e-07
DEAL:3d:cg:cg::Convergence step 1 value 0
DEAL:3d:cg::Convergence step 3 value 1.327e-06
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 4913
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 2.351e-05
DEAL:3d:cg::Starting value 70.09
-DEAL:3d:cg:cg::Starting value 0.03090
+DEAL:3d:cg:cg::Starting value 0.06180
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 3.917e-05
+DEAL:3d:cg:cg::Starting value 7.869e-05
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 5.305e-07
+DEAL:3d:cg:cg::Starting value 1.067e-06
DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 7.196e-10
+DEAL:3d:cg:cg::Starting value 1.447e-09
DEAL:3d:cg:cg::Convergence step 1 value 0
DEAL:3d:cg::Convergence step 4 value 6.374e-08
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 729
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
DEAL:3d:cg::Starting value 27.00
DEAL:3d:cg:cg::Starting value 1.155
DEAL:3d:cg:cg::Convergence step 2 value 0
-DEAL:3d:cg:cg::Starting value 0.0001365
+DEAL:3d:cg:cg::Starting value 0.0001367
DEAL:3d:cg:cg::Convergence step 2 value 0
DEAL:3d:cg:cg::Starting value 1.471e-07
DEAL:3d:cg:cg::Convergence step 2 value 0
DEAL:3d:cg::Convergence step 3 value 5.042e-08
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 4913
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
DEAL:3d:cg::Starting value 70.09
-DEAL:3d:cg:cg::Starting value 6.474
+DEAL:3d:cg:cg::Starting value 6.475
DEAL:3d:cg:cg::Convergence step 2 value 0
DEAL:3d:cg:cg::Starting value 0.0005841
DEAL:3d:cg:cg::Convergence step 2 value 0
DEAL:3d:cg::Convergence step 3 value 2.560e-06
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 35937
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
-DEAL:3d:cg::Starting value 1.000
-DEAL:3d:cg::Failure step 15 value 0.3991
DEAL:3d:cg::Starting value 189.6
DEAL:3d:cg:cg::Starting value 47.62
DEAL:3d:cg:cg::Convergence step 2 value 0
-DEAL:3d:cg:cg::Starting value 0.006356
+DEAL:3d:cg:cg::Starting value 0.006358
DEAL:3d:cg:cg::Convergence step 2 value 0
-DEAL:3d:cg:cg::Starting value 8.271e-06
+DEAL:3d:cg:cg::Starting value 8.274e-06
DEAL:3d:cg:cg::Convergence step 2 value 0
DEAL:3d:cg::Convergence step 3 value 1.007e-05
preconditioner(dof, mg, mg_transfer);
{
+ // avoid output from inner (coarse-level) solver
+ deallog.depth_file(3);
+
ReductionControl control(30, 1e-20, 1e-7);
SolverCG<LinearAlgebra::distributed::Vector<double> > solver(control);
solver.solve(fine_matrix, sol, in, preconditioner);
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 507
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9984
-DEAL:2d:cg::Failure step 15 value 0.001385
DEAL:2d:cg::Starting value 21.93
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 5 value 7.812e-07
+DEAL:2d:cg::Convergence step 5 value 7.781e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 881
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9988
-DEAL:2d:cg::Failure step 15 value 0.002581
-DEAL:2d:cg::Starting value 0.9989
-DEAL:2d:cg::Failure step 15 value 0.0001487
DEAL:2d:cg::Starting value 27.77
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 9.245e-07
+DEAL:2d:cg::Convergence step 6 value 9.215e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 2147
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9988
-DEAL:2d:cg::Failure step 15 value 0.003288
-DEAL:2d:cg::Starting value 0.9993
-DEAL:2d:cg::Failure step 15 value 0.005495
-DEAL:2d:cg::Starting value 0.9997
-DEAL:2d:cg::Failure step 15 value 0.05100
DEAL:2d:cg::Starting value 43.26
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 1.145e-06
+DEAL:2d:cg::Convergence step 6 value 1.142e-06
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 6517
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9989
-DEAL:2d:cg::Failure step 15 value 0.005519
-DEAL:2d:cg::Starting value 0.9994
-DEAL:2d:cg::Failure step 15 value 0.01098
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 0.09572
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.7527
DEAL:2d:cg::Starting value 76.92
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 1.559e-06
+DEAL:2d:cg::Convergence step 6 value 1.563e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 4259
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.763
DEAL:2d:cg::Starting value 64.26
-DEAL:2d:cg:cg::Starting value 108.9
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.03777
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.0007755
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 7.094e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 2.710e-07
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 1.631e-06
+DEAL:2d:cg::Convergence step 5 value 1.632e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 7457
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.672
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.9896
DEAL:2d:cg::Starting value 83.11
-DEAL:2d:cg:cg::Starting value 165.9
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.07398
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.003378
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 2.891e-05
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 1.842e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 7.691e-06
+DEAL:2d:cg::Convergence step 5 value 7.695e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 18535
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 1.886
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.8395
-DEAL:2d:cg::Starting value 1.000
-DEAL:2d:cg::Failure step 15 value 2.204
DEAL:2d:cg::Starting value 131.0
-DEAL:2d:cg:cg::Starting value 384.0
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.2113
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.002203
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 3.160e-05
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 1.004e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 6.045e-06
+DEAL:2d:cg::Convergence step 5 value 6.059e-06
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 1103
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9989
-DEAL:3d:cg::Convergence step 14 value 6.535e-11
DEAL:3d:cg::Starting value 31.21
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
DEAL:3d:cg::Convergence step 6 value 3.615e-08
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 3694
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9994
-DEAL:3d:cg::Failure step 15 value 1.839e-08
-DEAL:3d:cg::Starting value 0.9998
-DEAL:3d:cg::Failure step 15 value 7.469e-05
DEAL:3d:cg::Starting value 55.00
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
DEAL:3d:cg::Convergence step 7 value 3.941e-06
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 9566
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9995
-DEAL:3d:cg::Failure step 15 value 1.479e-07
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.0004644
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 1.830e-10
DEAL:3d:cg::Starting value 76.18
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg::Convergence step 8 value 5.376e-06
+DEAL:3d:cg::Convergence step 8 value 5.362e-06
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 7494
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
-DEAL:3d:cg::Starting value 0.9998
-DEAL:3d:cg::Failure step 15 value 0.0008241
DEAL:3d:cg::Starting value 82.90
-DEAL:3d:cg:cg::Starting value 7.225
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.01157
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0009262
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0001114
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 3.295e-06
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 2.088e-07
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 2.026e-08
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg::Convergence step 7 value 8.696e-07
+DEAL:3d:cg::Convergence step 7 value 8.682e-07
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 26548
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.008306
-DEAL:3d:cg::Starting value 1.000
-DEAL:3d:cg::Failure step 15 value 0.1026
DEAL:3d:cg::Starting value 152.6
-DEAL:3d:cg:cg::Starting value 17.46
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.008363
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.001190
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0001843
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 8.789e-06
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 7.697e-07
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 6.053e-08
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg::Convergence step 7 value 3.971e-06
+DEAL:3d:cg::Convergence step 7 value 3.969e-06
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 507
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9984
-DEAL:2d:cg::Failure step 15 value 0.001385
DEAL:2d:cg::Starting value 21.93
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 5 value 7.812e-07
+DEAL:2d:cg::Convergence step 5 value 7.781e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 881
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9988
-DEAL:2d:cg::Failure step 15 value 0.002581
-DEAL:2d:cg::Starting value 0.9989
-DEAL:2d:cg::Failure step 15 value 0.0001487
DEAL:2d:cg::Starting value 27.77
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 9.245e-07
+DEAL:2d:cg::Convergence step 6 value 9.215e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 2147
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9988
-DEAL:2d:cg::Failure step 15 value 0.003288
-DEAL:2d:cg::Starting value 0.9993
-DEAL:2d:cg::Failure step 15 value 0.005495
-DEAL:2d:cg::Starting value 0.9997
-DEAL:2d:cg::Failure step 15 value 0.05100
DEAL:2d:cg::Starting value 43.26
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 1.145e-06
+DEAL:2d:cg::Convergence step 6 value 1.142e-06
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 6517
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9989
-DEAL:2d:cg::Failure step 15 value 0.005519
-DEAL:2d:cg::Starting value 0.9994
-DEAL:2d:cg::Failure step 15 value 0.01098
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 0.09572
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.7527
DEAL:2d:cg::Starting value 76.92
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 1.559e-06
+DEAL:2d:cg::Convergence step 6 value 1.563e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 4259
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.763
DEAL:2d:cg::Starting value 64.26
-DEAL:2d:cg:cg::Starting value 108.9
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.03777
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.0007755
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 7.094e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 2.710e-07
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 1.631e-06
+DEAL:2d:cg::Convergence step 5 value 1.632e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 7457
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.672
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.9896
DEAL:2d:cg::Starting value 83.11
-DEAL:2d:cg:cg::Starting value 165.9
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.07398
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.003378
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 2.891e-05
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 1.842e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 7.691e-06
+DEAL:2d:cg::Convergence step 5 value 7.695e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 18535
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 1.886
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.8395
-DEAL:2d:cg::Starting value 1.000
-DEAL:2d:cg::Failure step 15 value 2.204
DEAL:2d:cg::Starting value 131.0
-DEAL:2d:cg:cg::Starting value 384.0
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.2113
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.002203
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 3.160e-05
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 1.004e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 6.045e-06
+DEAL:2d:cg::Convergence step 5 value 6.059e-06
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 1103
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9989
-DEAL:3d:cg::Convergence step 14 value 6.535e-11
DEAL:3d:cg::Starting value 31.21
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
DEAL:3d:cg::Convergence step 6 value 3.615e-08
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 3694
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9994
-DEAL:3d:cg::Failure step 15 value 1.839e-08
-DEAL:3d:cg::Starting value 0.9998
-DEAL:3d:cg::Failure step 15 value 7.469e-05
DEAL:3d:cg::Starting value 55.00
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
DEAL:3d:cg::Convergence step 7 value 3.941e-06
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 9566
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9995
-DEAL:3d:cg::Failure step 15 value 1.479e-07
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.0004644
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 1.830e-10
DEAL:3d:cg::Starting value 76.18
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg::Convergence step 8 value 5.376e-06
+DEAL:3d:cg::Convergence step 8 value 5.362e-06
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 7494
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
-DEAL:3d:cg::Starting value 0.9998
-DEAL:3d:cg::Failure step 15 value 0.0008241
DEAL:3d:cg::Starting value 82.90
-DEAL:3d:cg:cg::Starting value 7.225
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.01157
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0009262
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0001114
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 3.295e-06
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 2.088e-07
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 2.026e-08
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg::Convergence step 7 value 8.696e-07
+DEAL:3d:cg::Convergence step 7 value 8.682e-07
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 26548
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.008306
-DEAL:3d:cg::Starting value 1.000
-DEAL:3d:cg::Failure step 15 value 0.1026
DEAL:3d:cg::Starting value 152.6
-DEAL:3d:cg:cg::Starting value 17.46
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.008363
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.001190
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0001843
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 8.789e-06
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 7.697e-07
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 6.053e-08
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg::Convergence step 7 value 3.971e-06
+DEAL:3d:cg::Convergence step 7 value 3.969e-06
preconditioner(dof, mg, mg_transfer);
{
+ // avoid output from inner (coarse-level) solver
+ deallog.depth_file(3);
+
ReductionControl control(30, 1e-20, 1e-7);
SolverCG<LinearAlgebra::distributed::Vector<double> > solver(control);
solver.solve(fine_matrix, sol, in, preconditioner);
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 507
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9984
-DEAL:2d:cg::Failure step 15 value 0.001385
DEAL:2d:cg::Starting value 21.93
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 5 value 7.812e-07
+DEAL:2d:cg::Convergence step 5 value 7.781e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 881
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9988
-DEAL:2d:cg::Failure step 15 value 0.002581
-DEAL:2d:cg::Starting value 0.9989
-DEAL:2d:cg::Failure step 15 value 0.0001487
DEAL:2d:cg::Starting value 27.77
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 9.245e-07
+DEAL:2d:cg::Convergence step 6 value 9.215e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 2147
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9988
-DEAL:2d:cg::Failure step 15 value 0.003288
-DEAL:2d:cg::Starting value 0.9993
-DEAL:2d:cg::Failure step 15 value 0.005495
-DEAL:2d:cg::Starting value 0.9997
-DEAL:2d:cg::Failure step 15 value 0.05100
DEAL:2d:cg::Starting value 43.26
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 1.145e-06
+DEAL:2d:cg::Convergence step 6 value 1.142e-06
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 6517
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9989
-DEAL:2d:cg::Failure step 15 value 0.005519
-DEAL:2d:cg::Starting value 0.9994
-DEAL:2d:cg::Failure step 15 value 0.01098
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 0.09572
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.7527
DEAL:2d:cg::Starting value 76.92
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 1.559e-06
+DEAL:2d:cg::Convergence step 6 value 1.563e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 4259
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.763
DEAL:2d:cg::Starting value 64.26
-DEAL:2d:cg:cg::Starting value 108.9
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.03777
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.0007755
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 7.094e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 2.710e-07
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 1.631e-06
+DEAL:2d:cg::Convergence step 5 value 1.632e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 7457
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.672
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.9896
DEAL:2d:cg::Starting value 83.11
-DEAL:2d:cg:cg::Starting value 165.9
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.07398
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.003378
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 2.891e-05
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 1.842e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 7.691e-06
+DEAL:2d:cg::Convergence step 5 value 7.695e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 18535
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 1.886
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.8395
-DEAL:2d:cg::Starting value 1.000
-DEAL:2d:cg::Failure step 15 value 2.204
DEAL:2d:cg::Starting value 131.0
-DEAL:2d:cg:cg::Starting value 384.0
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.2113
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.002203
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 3.160e-05
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 1.004e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 6.045e-06
+DEAL:2d:cg::Convergence step 5 value 6.059e-06
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 1103
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9989
-DEAL:3d:cg::Convergence step 14 value 6.535e-11
DEAL:3d:cg::Starting value 31.21
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
DEAL:3d:cg::Convergence step 6 value 3.615e-08
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 3694
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9994
-DEAL:3d:cg::Failure step 15 value 1.839e-08
-DEAL:3d:cg::Starting value 0.9998
-DEAL:3d:cg::Failure step 15 value 7.469e-05
DEAL:3d:cg::Starting value 55.00
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
DEAL:3d:cg::Convergence step 7 value 3.941e-06
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 9566
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9995
-DEAL:3d:cg::Failure step 15 value 1.479e-07
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.0004644
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 1.830e-10
DEAL:3d:cg::Starting value 76.18
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg::Convergence step 8 value 5.376e-06
+DEAL:3d:cg::Convergence step 8 value 5.362e-06
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 7494
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
-DEAL:3d:cg::Starting value 0.9998
-DEAL:3d:cg::Failure step 15 value 0.0008241
DEAL:3d:cg::Starting value 82.90
-DEAL:3d:cg:cg::Starting value 7.225
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.01157
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0009262
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0001114
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 3.295e-06
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 2.088e-07
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 2.026e-08
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg::Convergence step 7 value 8.696e-07
+DEAL:3d:cg::Convergence step 7 value 8.682e-07
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 26548
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.008306
-DEAL:3d:cg::Starting value 1.000
-DEAL:3d:cg::Failure step 15 value 0.1026
DEAL:3d:cg::Starting value 152.6
-DEAL:3d:cg:cg::Starting value 17.46
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.008363
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.001190
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0001843
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 8.789e-06
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 7.697e-07
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 6.053e-08
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg::Convergence step 7 value 3.971e-06
+DEAL:3d:cg::Convergence step 7 value 3.969e-06
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 507
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9984
-DEAL:2d:cg::Failure step 15 value 0.001385
DEAL:2d:cg::Starting value 21.93
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 5 value 7.812e-07
+DEAL:2d:cg::Convergence step 5 value 7.781e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 881
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9988
-DEAL:2d:cg::Failure step 15 value 0.002581
-DEAL:2d:cg::Starting value 0.9989
-DEAL:2d:cg::Failure step 15 value 0.0001487
DEAL:2d:cg::Starting value 27.77
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 9.245e-07
+DEAL:2d:cg::Convergence step 6 value 9.215e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 2147
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9988
-DEAL:2d:cg::Failure step 15 value 0.003288
-DEAL:2d:cg::Starting value 0.9993
-DEAL:2d:cg::Failure step 15 value 0.005495
-DEAL:2d:cg::Starting value 0.9997
-DEAL:2d:cg::Failure step 15 value 0.05100
DEAL:2d:cg::Starting value 43.26
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 1.145e-06
+DEAL:2d:cg::Convergence step 6 value 1.142e-06
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 6517
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9989
-DEAL:2d:cg::Failure step 15 value 0.005519
-DEAL:2d:cg::Starting value 0.9994
-DEAL:2d:cg::Failure step 15 value 0.01098
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 0.09572
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.7527
DEAL:2d:cg::Starting value 76.92
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 1.559e-06
+DEAL:2d:cg::Convergence step 6 value 1.563e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 4259
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.763
DEAL:2d:cg::Starting value 64.26
-DEAL:2d:cg:cg::Starting value 108.9
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.03777
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.0007755
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 7.094e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 2.710e-07
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 1.631e-06
+DEAL:2d:cg::Convergence step 5 value 1.632e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 7457
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.672
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.9896
DEAL:2d:cg::Starting value 83.11
-DEAL:2d:cg:cg::Starting value 165.9
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.07398
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.003378
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 2.891e-05
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 1.842e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 7.691e-06
+DEAL:2d:cg::Convergence step 5 value 7.695e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 18535
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 1.886
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.8395
-DEAL:2d:cg::Starting value 1.000
-DEAL:2d:cg::Failure step 15 value 2.204
DEAL:2d:cg::Starting value 131.0
-DEAL:2d:cg:cg::Starting value 384.0
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.2113
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.002203
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 3.160e-05
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 1.004e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 6.045e-06
+DEAL:2d:cg::Convergence step 5 value 6.059e-06
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 1103
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9989
-DEAL:3d:cg::Convergence step 14 value 6.535e-11
DEAL:3d:cg::Starting value 31.21
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
DEAL:3d:cg::Convergence step 6 value 3.615e-08
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 3694
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9994
-DEAL:3d:cg::Failure step 15 value 1.839e-08
-DEAL:3d:cg::Starting value 0.9998
-DEAL:3d:cg::Failure step 15 value 7.469e-05
DEAL:3d:cg::Starting value 55.00
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
DEAL:3d:cg::Convergence step 7 value 3.941e-06
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 9566
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9995
-DEAL:3d:cg::Failure step 15 value 1.479e-07
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.0004644
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 1.830e-10
DEAL:3d:cg::Starting value 76.18
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg::Convergence step 8 value 5.376e-06
+DEAL:3d:cg::Convergence step 8 value 5.362e-06
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 7494
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
-DEAL:3d:cg::Starting value 0.9998
-DEAL:3d:cg::Failure step 15 value 0.0008241
DEAL:3d:cg::Starting value 82.90
-DEAL:3d:cg:cg::Starting value 7.225
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.01157
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0009262
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0001114
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 3.295e-06
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 2.088e-07
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 2.026e-08
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg::Convergence step 7 value 8.696e-07
+DEAL:3d:cg::Convergence step 7 value 8.682e-07
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 26548
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.008306
-DEAL:3d:cg::Starting value 1.000
-DEAL:3d:cg::Failure step 15 value 0.1026
DEAL:3d:cg::Starting value 152.6
-DEAL:3d:cg:cg::Starting value 17.46
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.008363
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.001190
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0001843
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 8.789e-06
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 7.697e-07
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 6.053e-08
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg::Convergence step 7 value 3.971e-06
+DEAL:3d:cg::Convergence step 7 value 3.969e-06
mg_matrices[level].get_matrix_diagonal_inverse();
}
- // temporarily disable deallog for the setup of the preconditioner that
- // involves a CG solver for eigenvalue estimation
- deallog.depth_file(0);
mg_smoother.initialize(mg_matrices, smoother_data);
mg::Matrix<LinearAlgebra::distributed::Vector<number> >
{
// avoid output from inner (coarse-level) solver
deallog.depth_file(3);
+
ReductionControl control(30, 1e-20, 1e-7);
SolverCG<LinearAlgebra::distributed::Vector<double> > solver(control);
solver.solve(fine_matrix, sol, in, preconditioner);
preconditioner(dof, mg, mg_transfer);
{
+ // avoid output from inner (coarse-level) solver
+ deallog.depth_file(2);
+
ReductionControl control(30, 1e-20, 1e-7);
SolverCG<LinearAlgebra::distributed::Vector<double> > solver(control);
solver.solve(fine_matrix, sol, in, preconditioner);
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 507
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9984
-DEAL:2d:cg::Failure step 15 value 0.001385
-DEAL:2d:cg::Starting value 21.93
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 5 value 7.812e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 881
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9988
-DEAL:2d:cg::Failure step 15 value 0.002581
-DEAL:2d:cg::Starting value 0.9989
-DEAL:2d:cg::Failure step 15 value 0.0001487
-DEAL:2d:cg::Starting value 27.77
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 9.245e-07
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 2147
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9988
-DEAL:2d:cg::Failure step 15 value 0.003288
-DEAL:2d:cg::Starting value 0.9993
-DEAL:2d:cg::Failure step 15 value 0.005495
-DEAL:2d:cg::Starting value 0.9997
-DEAL:2d:cg::Failure step 15 value 0.05100
-DEAL:2d:cg::Starting value 43.26
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 1.145e-06
DEAL:2d::Testing FE_Q<2>(1)
DEAL:2d::Number of degrees of freedom: 6517
-DEAL:2d:cg::Starting value 0.8660
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9428
-DEAL:2d:cg::Convergence step 1 value 0
-DEAL:2d:cg::Starting value 0.9798
-DEAL:2d:cg::Convergence step 4 value 0
-DEAL:2d:cg::Starting value 0.9938
-DEAL:2d:cg::Convergence step 10 value 0
-DEAL:2d:cg::Starting value 0.9983
-DEAL:2d:cg::Failure step 15 value 2.432e-05
-DEAL:2d:cg::Starting value 0.9989
-DEAL:2d:cg::Failure step 15 value 0.005519
-DEAL:2d:cg::Starting value 0.9994
-DEAL:2d:cg::Failure step 15 value 0.01098
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 0.09572
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.7527
-DEAL:2d:cg::Starting value 76.92
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg:cg::Convergence step 0 value 0
-DEAL:2d:cg::Convergence step 6 value 1.559e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 4259
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.763
-DEAL:2d:cg::Starting value 64.26
-DEAL:2d:cg:cg::Starting value 108.9
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.03777
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.0007755
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 7.094e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 2.710e-07
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 1.631e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 7457
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.672
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.9896
-DEAL:2d:cg::Starting value 83.11
-DEAL:2d:cg:cg::Starting value 165.9
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.07398
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.003378
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 2.891e-05
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 1.842e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 7.691e-06
DEAL:2d::Testing FE_Q<2>(3)
DEAL:2d::Number of degrees of freedom: 18535
-DEAL:2d:cg::Starting value 0.9682
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.9897
-DEAL:2d:cg::Convergence step 7 value 0
-DEAL:2d:cg::Starting value 0.9970
-DEAL:2d:cg::Failure step 15 value 2.456e-05
-DEAL:2d:cg::Starting value 0.9992
-DEAL:2d:cg::Failure step 15 value 0.2156
-DEAL:2d:cg::Starting value 0.9998
-DEAL:2d:cg::Failure step 15 value 1.976
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 1.886
-DEAL:2d:cg::Starting value 0.9999
-DEAL:2d:cg::Failure step 15 value 0.8395
-DEAL:2d:cg::Starting value 1.000
-DEAL:2d:cg::Failure step 15 value 2.204
-DEAL:2d:cg::Starting value 131.0
-DEAL:2d:cg:cg::Starting value 384.0
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.2113
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 0.002203
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 3.160e-05
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg:cg::Starting value 1.004e-06
-DEAL:2d:cg:cg::Convergence step 3 value 0
-DEAL:2d:cg::Convergence step 5 value 6.045e-06
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 1103
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9989
-DEAL:3d:cg::Convergence step 14 value 6.535e-11
-DEAL:3d:cg::Starting value 31.21
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg::Convergence step 6 value 3.615e-08
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 3694
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9994
-DEAL:3d:cg::Failure step 15 value 1.839e-08
-DEAL:3d:cg::Starting value 0.9998
-DEAL:3d:cg::Failure step 15 value 7.469e-05
-DEAL:3d:cg::Starting value 55.00
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg::Convergence step 7 value 3.941e-06
DEAL:3d::Testing FE_Q<3>(1)
DEAL:3d::Number of degrees of freedom: 9566
-DEAL:3d:cg::Starting value 0.9354
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 13 value 9.706e-11
-DEAL:3d:cg::Starting value 0.9995
-DEAL:3d:cg::Failure step 15 value 1.479e-07
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.0004644
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 1.830e-10
-DEAL:3d:cg::Starting value 76.18
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg:cg::Convergence step 0 value 0
-DEAL:3d:cg::Convergence step 8 value 5.376e-06
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 7494
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
-DEAL:3d:cg::Starting value 0.9998
-DEAL:3d:cg::Failure step 15 value 0.0008241
-DEAL:3d:cg::Starting value 82.90
-DEAL:3d:cg:cg::Starting value 7.225
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.01157
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0009262
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0001114
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 3.295e-06
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 2.088e-07
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 2.026e-08
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg::Convergence step 7 value 8.696e-07
DEAL:3d::Testing FE_Q<3>(2)
DEAL:3d::Number of degrees of freedom: 26548
-DEAL:3d:cg::Starting value 0.9813
-DEAL:3d:cg::Convergence step 1 value 0
-DEAL:3d:cg::Starting value 0.9960
-DEAL:3d:cg::Convergence step 5 value 0
-DEAL:3d:cg::Starting value 0.9993
-DEAL:3d:cg::Convergence step 15 value 5.276e-11
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.001572
-DEAL:3d:cg::Starting value 0.9999
-DEAL:3d:cg::Failure step 15 value 0.008306
-DEAL:3d:cg::Starting value 1.000
-DEAL:3d:cg::Failure step 15 value 0.1026
-DEAL:3d:cg::Starting value 152.6
-DEAL:3d:cg:cg::Starting value 17.46
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.008363
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.001190
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 0.0001843
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 8.789e-06
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 7.697e-07
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg:cg::Starting value 6.053e-08
-DEAL:3d:cg:cg::Convergence step 1 value 0
-DEAL:3d:cg::Convergence step 7 value 3.971e-06
// temporarily disable deallog for the setup of the preconditioner that
// involves a CG solver for eigenvalue estimation
- deallog.depth_file(0);
mg_smoother.initialize(mg_matrices, smoother_data);
mg::Matrix<LinearAlgebra::distributed::Vector<number> >
DEAL:2d::Cycle 0
DEAL:2d::Number of degrees of freedom: 81
-DEAL:2d:cg::Starting value 0.943
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.980
-DEAL:2d:cg::Convergence step 8 value 0
-DEAL:2d:cg::Starting value 0.994
-DEAL:2d:cg::Failure step 10 value 0.00110
DEAL:2d:cg::Starting value 0.132
DEAL:2d:cg::Convergence step 4 value 0
DEAL:2d::
DEAL:2d::Cycle 1
DEAL:2d::Number of degrees of freedom: 289
-DEAL:2d:cg::Starting value 0.943
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.980
-DEAL:2d:cg::Convergence step 8 value 0
-DEAL:2d:cg::Starting value 0.994
-DEAL:2d:cg::Failure step 10 value 0.00110
-DEAL:2d:cg::Starting value 0.998
-DEAL:2d:cg::Failure step 10 value 0.153
DEAL:2d:cg::Starting value 0.0677
DEAL:2d:cg::Convergence step 5 value 0
DEAL:2d::
DEAL:2d::Cycle 2
DEAL:2d::Number of degrees of freedom: 1089
-DEAL:2d:cg::Starting value 0.943
-DEAL:2d:cg::Convergence step 2 value 0
-DEAL:2d:cg::Starting value 0.980
-DEAL:2d:cg::Convergence step 8 value 0
-DEAL:2d:cg::Starting value 0.994
-DEAL:2d:cg::Failure step 10 value 0.00110
-DEAL:2d:cg::Starting value 0.998
-DEAL:2d:cg::Failure step 10 value 0.153
-DEAL:2d:cg::Starting value 1.00
-DEAL:2d:cg::Failure step 10 value 1.07
DEAL:2d:cg::Starting value 0.0343
DEAL:2d:cg::Convergence step 5 value 0
DEAL:2d::
DEAL:3d::Cycle 0
DEAL:3d::Number of degrees of freedom: 125
-DEAL:3d:cg::Starting value 0.981
-DEAL:3d:cg::Convergence step 2 value 0
-DEAL:3d:cg::Starting value 0.996
-DEAL:3d:cg::Convergence step 10 value 0
DEAL:3d:cg::Starting value 0.125
DEAL:3d:cg::Convergence step 5 value 0
DEAL:3d::
DEAL:3d::Cycle 1
DEAL:3d::Number of degrees of freedom: 729
-DEAL:3d:cg::Starting value 0.981
-DEAL:3d:cg::Convergence step 2 value 0
-DEAL:3d:cg::Starting value 0.996
-DEAL:3d:cg::Convergence step 10 value 0
-DEAL:3d:cg::Starting value 0.999
-DEAL:3d:cg::Failure step 10 value 0.000695
DEAL:3d:cg::Starting value 0.0479
DEAL:3d:cg::Convergence step 4 value 0
DEAL:3d::
DEAL:3d::Cycle 2
DEAL:3d::Number of degrees of freedom: 4913
-DEAL:3d:cg::Starting value 0.981
-DEAL:3d:cg::Convergence step 2 value 0
-DEAL:3d:cg::Starting value 0.996
-DEAL:3d:cg::Convergence step 10 value 0
-DEAL:3d:cg::Starting value 0.999
-DEAL:3d:cg::Failure step 10 value 0.000695
-DEAL:3d:cg::Starting value 1.00
-DEAL:3d:cg::Failure step 10 value 0.0952
DEAL:3d:cg::Starting value 0.0176
DEAL:3d:cg::Convergence step 4 value 0
DEAL:3d::
+++ /dev/null
-
-DEAL::Cycle 0:
-DEAL:: Number of active cells: 64
-DEAL:: Number of degrees of freedom: 81 (by level: 25, 81)
-DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
-DEAL:cg::Starting value 0.1094
-DEAL:cg::Convergence step 10 value 0
-DEAL:: 10 CG iterations needed to obtain convergence.
-DEAL::Cycle 1:
-DEAL:: Number of active cells: 127
-DEAL:: Number of degrees of freedom: 160 (by level: 25, 81, 117)
-DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
-DEAL:cg::Starting value 0.9957
-DEAL:cg::Convergence step 15 value 0
-DEAL:cg::Starting value 0.09101
-DEAL:cg::Convergence step 13 value 0
-DEAL:: 13 CG iterations needed to obtain convergence.
-DEAL::Cycle 2:
-DEAL:: Number of active cells: 256
-DEAL:: Number of degrees of freedom: 301 (by level: 25, 81, 245, 71)
-DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
-DEAL:cg::Starting value 0.9980
-DEAL:cg::Failure step 20 value 7.015e-05
-DEAL:cg::Starting value 0.9929
-DEAL:cg::Convergence step 11 value 0
-DEAL:cg::Starting value 0.06376
-DEAL:cg::Convergence step 12 value 0
-DEAL:: 12 CG iterations needed to obtain convergence.
-DEAL::Cycle 3:
-DEAL:: Number of active cells: 493
-DEAL:: Number of degrees of freedom: 600 (by level: 25, 81, 265, 376, 119)
-DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
-DEAL:cg::Starting value 0.9981
-DEAL:cg::Failure step 20 value 0.0001644
-DEAL:cg::Starting value 0.9987
-DEAL:cg::Failure step 20 value 5.735e-09
-DEAL:cg::Starting value 0.9958
-DEAL:cg::Convergence step 12 value 0
-DEAL:cg::Starting value 0.05706
-DEAL:cg::Convergence step 15 value 0
-DEAL:: 15 CG iterations needed to obtain convergence.
-DEAL::Cycle 4:
-DEAL:: Number of active cells: 979
-DEAL:: Number of degrees of freedom: 1093 (by level: 25, 81, 273, 791, 269, 88)
-DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
-DEAL:cg::Starting value 0.9982
-DEAL:cg::Failure step 20 value 0.0002935
-DEAL:cg::Starting value 0.9994
-DEAL:cg::Failure step 20 value 0.02813
-DEAL:cg::Starting value 0.9981
-DEAL:cg::Failure step 20 value 1.498e-09
-DEAL:cg::Starting value 0.9943
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.04157
-DEAL:cg::Convergence step 15 value 0
-DEAL:: 15 CG iterations needed to obtain convergence.
-DEAL::Cycle 5:
-DEAL:: Number of active cells: 1912
-DEAL:: Number of degrees of freedom: 2184 (by level: 25, 81, 277, 961, 1070, 501, 172)
-DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
-DEAL:cg::Starting value 0.9982
-DEAL:cg::Failure step 20 value 0.0004177
-DEAL:cg::Starting value 0.9995
-DEAL:cg::Failure step 20 value 0.05428
-DEAL:cg::Starting value 0.9995
-DEAL:cg::Failure step 20 value 0.0001063
-DEAL:cg::Starting value 0.9990
-DEAL:cg::Failure step 20 value 1.860e-09
-DEAL:cg::Starting value 0.9971
-DEAL:cg::Convergence step 13 value 0
-DEAL:cg::Starting value 0.03510
-DEAL:cg::Convergence step 15 value 0
-DEAL:: 15 CG iterations needed to obtain convergence.
--- /dev/null
+
+DEAL::Cycle 0:
+DEAL:: Number of active cells: 64
+DEAL:: Number of degrees of freedom: 81 (by level: 25, 81)
+DEAL:: Size of coarse grid matrix: 25
+DEAL:cg::Starting value 0.1094
+DEAL:cg::Convergence step 10 value 0
+DEAL:: 10 CG iterations needed to obtain convergence.
+DEAL::Cycle 1:
+DEAL:: Number of active cells: 127
+DEAL:: Number of degrees of freedom: 160 (by level: 25, 81, 117)
+DEAL:: Size of coarse grid matrix: 25
+DEAL:cg::Starting value 0.09101
+DEAL:cg::Convergence step 13 value 0
+DEAL:: 13 CG iterations needed to obtain convergence.
+DEAL::Cycle 2:
+DEAL:: Number of active cells: 256
+DEAL:: Number of degrees of freedom: 301 (by level: 25, 81, 245, 71)
+DEAL:: Size of coarse grid matrix: 25
+DEAL:cg::Starting value 0.06376
+DEAL:cg::Convergence step 12 value 0
+DEAL:: 12 CG iterations needed to obtain convergence.
+DEAL::Cycle 3:
+DEAL:: Number of active cells: 493
+DEAL:: Number of degrees of freedom: 600 (by level: 25, 81, 265, 376, 119)
+DEAL:: Size of coarse grid matrix: 25
+DEAL:cg::Starting value 0.05706
+DEAL:cg::Convergence step 15 value 0
+DEAL:: 15 CG iterations needed to obtain convergence.
+DEAL::Cycle 4:
+DEAL:: Number of active cells: 979
+DEAL:: Number of degrees of freedom: 1093 (by level: 25, 81, 273, 791, 269, 88)
+DEAL:: Size of coarse grid matrix: 25
+DEAL:cg::Starting value 0.04157
+DEAL:cg::Convergence step 15 value 0
+DEAL:: 15 CG iterations needed to obtain convergence.
+DEAL::Cycle 5:
+DEAL:: Number of active cells: 1912
+DEAL:: Number of degrees of freedom: 2184 (by level: 25, 81, 277, 961, 1070, 501, 172)
+DEAL:: Size of coarse grid matrix: 25
+DEAL:cg::Starting value 0.03510
+DEAL:cg::Convergence step 15 value 0
+DEAL:: 15 CG iterations needed to obtain convergence.
DEAL:: Number of active cells: 64
DEAL:: Number of degrees of freedom: 81 (by level: 25, 81)
DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
DEAL:cg::Starting value 0.1094
DEAL:cg::Convergence step 10 value 0
DEAL:: 10 CG iterations needed to obtain convergence.
DEAL:: Number of active cells: 256
DEAL:: Number of degrees of freedom: 289 (by level: 25, 81, 289)
DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
-DEAL:cg::Starting value 0.9983
-DEAL:cg::Failure step 20 value 0.001167
DEAL:cg::Starting value 0.05859
DEAL:cg::Convergence step 10 value 0
DEAL:: 10 CG iterations needed to obtain convergence.
DEAL:: Number of active cells: 1024
DEAL:: Number of degrees of freedom: 1089 (by level: 25, 81, 289, 1089)
DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
-DEAL:cg::Starting value 0.9983
-DEAL:cg::Failure step 20 value 0.001167
-DEAL:cg::Starting value 0.9995
-DEAL:cg::Failure step 20 value 0.1128
DEAL:cg::Starting value 0.03027
DEAL:cg::Convergence step 11 value 0
DEAL:: 11 CG iterations needed to obtain convergence.
DEAL:: Number of active cells: 4096
DEAL:: Number of degrees of freedom: 4225 (by level: 25, 81, 289, 1089, 4225)
DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
-DEAL:cg::Starting value 0.9983
-DEAL:cg::Failure step 20 value 0.001167
-DEAL:cg::Starting value 0.9995
-DEAL:cg::Failure step 20 value 0.1128
-DEAL:cg::Starting value 0.9999
-DEAL:cg::Failure step 20 value 1.273
DEAL:cg::Starting value 0.01538
DEAL:cg::Convergence step 11 value 0
DEAL:: 11 CG iterations needed to obtain convergence.
DEAL:: Number of active cells: 64
DEAL:: Number of degrees of freedom: 81 (by level: 25, 81)
DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
DEAL:cg::Starting value 0.1094
DEAL:cg::Convergence step 10 value 0
DEAL:: 10 CG iterations needed to obtain convergence.
DEAL:: Number of active cells: 256
DEAL:: Number of degrees of freedom: 289 (by level: 25, 81, 289)
DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
-DEAL:cg::Starting value 0.9983
-DEAL:cg::Failure step 20 value 0.001167
DEAL:cg::Starting value 0.05859
DEAL:cg::Convergence step 10 value 0
DEAL:: 10 CG iterations needed to obtain convergence.
DEAL:: Number of active cells: 1024
DEAL:: Number of degrees of freedom: 1089 (by level: 25, 81, 289, 1089)
DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
-DEAL:cg::Starting value 0.9983
-DEAL:cg::Failure step 20 value 0.001167
-DEAL:cg::Starting value 0.9995
-DEAL:cg::Failure step 20 value 0.1128
DEAL:cg::Starting value 0.03027
DEAL:cg::Convergence step 11 value 0
DEAL:: 11 CG iterations needed to obtain convergence.
DEAL:: Number of active cells: 4096
DEAL:: Number of degrees of freedom: 4225 (by level: 25, 81, 289, 1089, 4225)
DEAL:: Size of coarse grid matrix: 25
-DEAL:cg::Starting value 0.9798
-DEAL:cg::Convergence step 7 value 0
-DEAL:cg::Starting value 0.9938
-DEAL:cg::Convergence step 18 value 0
-DEAL:cg::Starting value 0.9983
-DEAL:cg::Failure step 20 value 0.001167
-DEAL:cg::Starting value 0.9995
-DEAL:cg::Failure step 20 value 0.1128
-DEAL:cg::Starting value 0.9999
-DEAL:cg::Failure step 20 value 1.273
DEAL:cg::Starting value 0.01538
DEAL:cg::Convergence step 11 value 0
DEAL:: 11 CG iterations needed to obtain convergence.
projects / dealii.git / commitdiff