From c4ad6b01cad7d58eae44fbbc21ab03b3c496ff3e Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Thu, 27 Oct 2016 15:21:48 +0200 Subject: [PATCH] Introduce template for preconditioner type in Chebyshev smoother. --- include/deal.II/lac/precondition.h | 420 ++++++++++++------ ...r_08.with_cxx11=on.with_umfpack=off.output | 2 - ...or_08.with_cxx11=on.with_umfpack=on.output | 2 - ...ition_chebyshev_01.with_lapack=true.output | 8 +- ...h_mpi=true.with_p4est=true.mpirun=3.output | 186 ++------ ...h_mpi=true.with_p4est=true.mpirun=7.output | 186 ++------ .../parallel_multigrid_adaptive_01.cc | 3 + ...h_mpi=true.with_p4est=true.mpirun=3.output | 282 +----------- ...h_mpi=true.with_p4est=true.mpirun=7.output | 282 +----------- .../parallel_multigrid_adaptive_02.cc | 3 + ...h_mpi=true.with_p4est=true.mpirun=3.output | 282 +----------- ...h_mpi=true.with_p4est=true.mpirun=7.output | 282 +----------- .../parallel_multigrid_adaptive_03.cc | 4 +- .../parallel_multigrid_adaptive_04.cc | 3 + ...h_mpi=true.with_p4est=true.mpirun=3.output | 286 ------------ tests/matrix_free/parallel_multigrid_mf.cc | 1 - .../step-37.with_lapack=true.output | 42 -- tests/multigrid/step-16-04.output | 97 ---- .../step-16-04.with_lapack=on.output | 43 ++ ...utput => step-16-05.with_lapack=on.output} | 28 -- ...utput => step-16-06.with_lapack=on.output} | 28 -- 21 files changed, 452 insertions(+), 2018 deletions(-) delete mode 100644 tests/multigrid/step-16-04.output create mode 100644 tests/multigrid/step-16-04.with_lapack=on.output rename tests/multigrid/{step-16-05.output => step-16-05.with_lapack=on.output} (55%) rename tests/multigrid/{step-16-06.output => step-16-06.with_lapack=on.output} (55%) diff --git a/include/deal.II/lac/precondition.h b/include/deal.II/lac/precondition.h index dface084ed..7f71ebe525 100644 --- a/include/deal.II/lac/precondition.h +++ b/include/deal.II/lac/precondition.h @@ -23,7 +23,8 @@ #include #include #include -#include +#include +#include #include #include @@ -811,19 +812,73 @@ private: /** * 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 vmult 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. + * + *

Algorithm execution

+ * + * 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 + * [max_eigenvalue/smoothing_range, max_eigenvalue] 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. + * + *

Requirements on the templated classes

+ * + * 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 el(i,i) 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 , class VectorType=Vector > +template , typename VectorType=Vector, typename PreconditionerType=DiagonalMatrix > class PreconditionChebyshev : public Subscriptor { public: @@ -905,8 +960,15 @@ 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 preconditioner; }; PreconditionChebyshev (); @@ -973,7 +1035,7 @@ private: /** * A pointer to the underlying matrix. */ - SmartPointer > matrix_ptr; + SmartPointer > matrix_ptr; /** * Internal vector used for the vmult operation. @@ -986,7 +1048,13 @@ private: mutable VectorType update2; /** - * Stores the additional data provided to the initialize function. + * Internal vector used for the vmult operation. + */ + mutable VectorType update3; + + /** + * Stores the additional data passed to the initialize function, obtained + * through a copy operation. */ AdditionalData data; @@ -1002,10 +1070,17 @@ private: 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. @@ -1021,6 +1096,14 @@ private: */ 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; }; @@ -1604,28 +1687,30 @@ namespace internal // size of the vector // generic part for non-deal.II vectors - template + template 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 @@ -1738,61 +1823,89 @@ namespace internal const VectorUpdater &updater; }; - // selection for deal.II vector + // selection for diagonal matrix around deal.II vector template inline void vector_updates (const ::dealii::Vector &src, - const ::dealii::Vector &matrix_diagonal_inverse, + const DiagonalMatrix<::dealii::Vector > &jacobi, const bool start_zero, const double factor1, const double factor2, ::dealii::Vector &update1, ::dealii::Vector &update2, + ::dealii::Vector &, ::dealii::Vector &dst) { - VectorUpdater upd(src.begin(), matrix_diagonal_inverse.begin(), + VectorUpdater upd(src.begin(), jacobi.get_vector().begin(), start_zero, factor1, factor2, update1.begin(), update2.begin(), dst.begin()); VectorUpdatesRange(upd, src.size()); } - // selection for parallel deal.II vector + // selection for diagonal matrix around parallel deal.II vector template inline void vector_updates (const LinearAlgebra::distributed::Vector &src, - const LinearAlgebra::distributed::Vector &matrix_diagonal_inverse, + const DiagonalMatrix > &jacobi, const bool start_zero, const double factor1, const double factor2, LinearAlgebra::distributed::Vector &update1, LinearAlgebra::distributed::Vector &update2, + LinearAlgebra::distributed::Vector &, LinearAlgebra::distributed::Vector &dst) { - VectorUpdater upd(src.begin(), matrix_diagonal_inverse.begin(), + VectorUpdater upd(src.begin(), jacobi.get_vector().begin(), start_zero, factor1, factor2, update1.begin(), update2.begin(), dst.begin()); VectorUpdatesRange(upd, src.local_size()); } - template - struct DiagonalPreconditioner + template + inline + void + initialize_preconditioner(const MatrixType &matrix, + std_cxx11::shared_ptr &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 + inline + void + initialize_preconditioner(const MatrixType &matrix, + std_cxx11::shared_ptr > &preconditioner, + VectorType &diagonal_inverse) + { + if (preconditioner.get() == NULL || + preconditioner->m() != matrix.m()) + { + if (preconditioner.get() == NULL) + preconditioner.reset(new DiagonalMatrix()); - 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; iget_vector()(i) = 1./matrix.el(i,i); + } + } + } struct EigenvalueTracker { @@ -1809,9 +1922,9 @@ namespace internal -template +template inline -PreconditionChebyshev::AdditionalData:: +PreconditionChebyshev::AdditionalData:: AdditionalData (const unsigned int degree, const double smoothing_range, const bool nonzero_starting, @@ -1829,9 +1942,9 @@ AdditionalData (const unsigned int degree, -template +template inline -PreconditionChebyshev::PreconditionChebyshev () +PreconditionChebyshev::PreconditionChebyshev () : is_initialized (false) { @@ -1844,24 +1957,34 @@ PreconditionChebyshev::PreconditionChebyshev () -template +template inline void -PreconditionChebyshev::initialize +PreconditionChebyshev::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 +inline +void +PreconditionChebyshev::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 @@ -1869,57 +1992,48 @@ PreconditionChebyshev::initialize 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 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(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::epsilon()), + 1e-10, false, false); internal::PreconditionChebyshev::EigenvalueTracker eigenvalue_tracker; - SolverCG solver (control, memory); + SolverCG solver (control); solver.connect_eigenvalues_slot(std_cxx11::bind(&internal::PreconditionChebyshev::EigenvalueTracker::slot, &eigenvalue_tracker, std_cxx11::_1)); - internal::PreconditionChebyshev::DiagonalPreconditioner - 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 { @@ -1930,23 +2044,35 @@ PreconditionChebyshev::initialize 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 *>(this)->delta = (max_eigenvalue-alpha)*0.5; + const_cast *>(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 >::value == false || + (types_are_equal >::value == false + && + types_are_equal >::value == false + )) + update3.reinit (src, true); + + const_cast *>(this)->is_initialized = true; } -template +template inline void -PreconditionChebyshev::do_chebyshev_loop(VectorType &dst, - const VectorType &src) const +PreconditionChebyshev +::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::do_chebyshev_loop(VectorType 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 +template inline void -PreconditionChebyshev::do_transpose_chebyshev_loop(VectorType &dst, - const VectorType &src) const +PreconditionChebyshev +::do_transpose_chebyshev_loop(VectorType &dst, + const VectorType &src) const { double rhok = delta / theta, sigma = theta / delta; for (unsigned int k=0; k::do_transpose_chebyshev_loop(Vector 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 +template inline void -PreconditionChebyshev::vmult (VectorType &dst, - const VectorType &src) const +PreconditionChebyshev +::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 +template inline void -PreconditionChebyshev::Tvmult (VectorType &dst, - const VectorType &src) const +PreconditionChebyshev +::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 +template inline void -PreconditionChebyshev::step (VectorType &dst, - const VectorType &src) const +PreconditionChebyshev +::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 +template inline void -PreconditionChebyshev::Tstep (VectorType &dst, - const VectorType &src) const +PreconditionChebyshev +::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 +template inline -void PreconditionChebyshev::clear () +void +PreconditionChebyshev::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 +template inline -typename PreconditionChebyshev::size_type -PreconditionChebyshev::m () const +typename PreconditionChebyshev::size_type +PreconditionChebyshev::m () const { Assert (matrix_ptr!=0, ExcNotInitialized()); return matrix_ptr->m(); } -template +template inline -typename PreconditionChebyshev::size_type -PreconditionChebyshev::n () const +typename PreconditionChebyshev::size_type +PreconditionChebyshev::n () const { Assert (matrix_ptr!=0, ExcNotInitialized()); return matrix_ptr->n(); diff --git a/tests/lac/linear_operator_08.with_cxx11=on.with_umfpack=off.output b/tests/lac/linear_operator_08.with_cxx11=on.with_umfpack=off.output index 84bdb2e317..55468cb785 100644 --- a/tests/lac/linear_operator_08.with_cxx11=on.with_umfpack=off.output +++ b/tests/lac/linear_operator_08.with_cxx11=on.with_umfpack=off.output @@ -1,8 +1,6 @@ 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 diff --git a/tests/lac/linear_operator_08.with_cxx11=on.with_umfpack=on.output b/tests/lac/linear_operator_08.with_cxx11=on.with_umfpack=on.output index 459bc5bca8..d47fd5d5ed 100644 --- a/tests/lac/linear_operator_08.with_cxx11=on.with_umfpack=on.output +++ b/tests/lac/linear_operator_08.with_cxx11=on.with_umfpack=on.output @@ -1,8 +1,6 @@ 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 diff --git a/tests/lac/precondition_chebyshev_01.with_lapack=true.output b/tests/lac/precondition_chebyshev_01.with_lapack=true.output index f09bc650ad..2d5be039fc 100644 --- a/tests/lac/precondition_chebyshev_01.with_lapack=true.output +++ b/tests/lac/precondition_chebyshev_01.with_lapack=true.output @@ -1,10 +1,6 @@ -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 diff --git a/tests/matrix_free/parallel_multigrid.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=3.output b/tests/matrix_free/parallel_multigrid.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=3.output index 9d5dea07fe..c1311be955 100644 --- a/tests/matrix_free/parallel_multigrid.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=3.output +++ b/tests/matrix_free/parallel_multigrid.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=3.output @@ -1,216 +1,118 @@ 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 @@ -219,21 +121,11 @@ 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 diff --git a/tests/matrix_free/parallel_multigrid.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=7.output b/tests/matrix_free/parallel_multigrid.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=7.output index 9d5dea07fe..c1311be955 100644 --- a/tests/matrix_free/parallel_multigrid.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=7.output +++ b/tests/matrix_free/parallel_multigrid.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=7.output @@ -1,216 +1,118 @@ 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 @@ -219,21 +121,11 @@ 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 diff --git a/tests/matrix_free/parallel_multigrid_adaptive_01.cc b/tests/matrix_free/parallel_multigrid_adaptive_01.cc index 1ac481fb76..1cd256990e 100644 --- a/tests/matrix_free/parallel_multigrid_adaptive_01.cc +++ b/tests/matrix_free/parallel_multigrid_adaptive_01.cc @@ -530,6 +530,9 @@ void do_test (const DoFHandler &dof) preconditioner(dof, mg, mg_transfer); { + // avoid output from inner (coarse-level) solver + deallog.depth_file(3); + ReductionControl control(30, 1e-20, 1e-7); SolverCG > solver(control); solver.solve(fine_matrix, sol, in, preconditioner); diff --git a/tests/matrix_free/parallel_multigrid_adaptive_01.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=3.output b/tests/matrix_free/parallel_multigrid_adaptive_01.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=3.output index e9349c579a..16d865c56a 100644 --- a/tests/matrix_free/parallel_multigrid_adaptive_01.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=3.output +++ b/tests/matrix_free/parallel_multigrid_adaptive_01.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=3.output @@ -1,311 +1,49 @@ 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 diff --git a/tests/matrix_free/parallel_multigrid_adaptive_01.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=7.output b/tests/matrix_free/parallel_multigrid_adaptive_01.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=7.output index e9349c579a..16d865c56a 100644 --- a/tests/matrix_free/parallel_multigrid_adaptive_01.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=7.output +++ b/tests/matrix_free/parallel_multigrid_adaptive_01.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=7.output @@ -1,311 +1,49 @@ 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 diff --git a/tests/matrix_free/parallel_multigrid_adaptive_02.cc b/tests/matrix_free/parallel_multigrid_adaptive_02.cc index 03c6d6e6f1..bf08a47b3b 100644 --- a/tests/matrix_free/parallel_multigrid_adaptive_02.cc +++ b/tests/matrix_free/parallel_multigrid_adaptive_02.cc @@ -385,6 +385,9 @@ void do_test (const DoFHandler &dof) preconditioner(dof, mg, mg_transfer); { + // avoid output from inner (coarse-level) solver + deallog.depth_file(3); + ReductionControl control(30, 1e-20, 1e-7); SolverCG > solver(control); solver.solve(fine_matrix, sol, in, preconditioner); diff --git a/tests/matrix_free/parallel_multigrid_adaptive_02.with_mpi=true.with_p4est=true.mpirun=3.output b/tests/matrix_free/parallel_multigrid_adaptive_02.with_mpi=true.with_p4est=true.mpirun=3.output index e9349c579a..16d865c56a 100644 --- a/tests/matrix_free/parallel_multigrid_adaptive_02.with_mpi=true.with_p4est=true.mpirun=3.output +++ b/tests/matrix_free/parallel_multigrid_adaptive_02.with_mpi=true.with_p4est=true.mpirun=3.output @@ -1,311 +1,49 @@ 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 diff --git a/tests/matrix_free/parallel_multigrid_adaptive_02.with_mpi=true.with_p4est=true.mpirun=7.output b/tests/matrix_free/parallel_multigrid_adaptive_02.with_mpi=true.with_p4est=true.mpirun=7.output index e9349c579a..16d865c56a 100644 --- a/tests/matrix_free/parallel_multigrid_adaptive_02.with_mpi=true.with_p4est=true.mpirun=7.output +++ b/tests/matrix_free/parallel_multigrid_adaptive_02.with_mpi=true.with_p4est=true.mpirun=7.output @@ -1,311 +1,49 @@ 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 diff --git a/tests/matrix_free/parallel_multigrid_adaptive_03.cc b/tests/matrix_free/parallel_multigrid_adaptive_03.cc index 133f617b0d..ea83846e7c 100644 --- a/tests/matrix_free/parallel_multigrid_adaptive_03.cc +++ b/tests/matrix_free/parallel_multigrid_adaptive_03.cc @@ -507,9 +507,6 @@ void do_test (const DoFHandler &dof, const bool threaded) 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 > @@ -531,6 +528,7 @@ void do_test (const DoFHandler &dof, const bool threaded) { // avoid output from inner (coarse-level) solver deallog.depth_file(3); + ReductionControl control(30, 1e-20, 1e-7); SolverCG > solver(control); solver.solve(fine_matrix, sol, in, preconditioner); diff --git a/tests/matrix_free/parallel_multigrid_adaptive_04.cc b/tests/matrix_free/parallel_multigrid_adaptive_04.cc index b189fa73d7..6bf48b8858 100644 --- a/tests/matrix_free/parallel_multigrid_adaptive_04.cc +++ b/tests/matrix_free/parallel_multigrid_adaptive_04.cc @@ -385,6 +385,9 @@ void do_test (const DoFHandler &dof) preconditioner(dof, mg, mg_transfer); { + // avoid output from inner (coarse-level) solver + deallog.depth_file(2); + ReductionControl control(30, 1e-20, 1e-7); SolverCG > solver(control); solver.solve(fine_matrix, sol, in, preconditioner); diff --git a/tests/matrix_free/parallel_multigrid_adaptive_04.with_mpi=true.with_p4est=true.mpirun=3.output b/tests/matrix_free/parallel_multigrid_adaptive_04.with_mpi=true.with_p4est=true.mpirun=3.output index e9349c579a..fa84d2c132 100644 --- a/tests/matrix_free/parallel_multigrid_adaptive_04.with_mpi=true.with_p4est=true.mpirun=3.output +++ b/tests/matrix_free/parallel_multigrid_adaptive_04.with_mpi=true.with_p4est=true.mpirun=3.output @@ -1,311 +1,25 @@ 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 diff --git a/tests/matrix_free/parallel_multigrid_mf.cc b/tests/matrix_free/parallel_multigrid_mf.cc index 3a81c512f4..c7dc7aaef7 100644 --- a/tests/matrix_free/parallel_multigrid_mf.cc +++ b/tests/matrix_free/parallel_multigrid_mf.cc @@ -364,7 +364,6 @@ void do_test (const DoFHandler &dof) // 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 > diff --git a/tests/matrix_free/step-37.with_lapack=true.output b/tests/matrix_free/step-37.with_lapack=true.output index 76f9274e64..f67b52a7ea 100644 --- a/tests/matrix_free/step-37.with_lapack=true.output +++ b/tests/matrix_free/step-37.with_lapack=true.output @@ -1,73 +1,31 @@ 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:: diff --git a/tests/multigrid/step-16-04.output b/tests/multigrid/step-16-04.output deleted file mode 100644 index 7fb681c689..0000000000 --- a/tests/multigrid/step-16-04.output +++ /dev/null @@ -1,97 +0,0 @@ - -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. diff --git a/tests/multigrid/step-16-04.with_lapack=on.output b/tests/multigrid/step-16-04.with_lapack=on.output new file mode 100644 index 0000000000..29325f556a --- /dev/null +++ b/tests/multigrid/step-16-04.with_lapack=on.output @@ -0,0 +1,43 @@ + +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. diff --git a/tests/multigrid/step-16-05.output b/tests/multigrid/step-16-05.with_lapack=on.output similarity index 55% rename from tests/multigrid/step-16-05.output rename to tests/multigrid/step-16-05.with_lapack=on.output index b7f18bfa52..a6986c2c48 100644 --- a/tests/multigrid/step-16-05.output +++ b/tests/multigrid/step-16-05.with_lapack=on.output @@ -3,10 +3,6 @@ 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. @@ -14,12 +10,6 @@ DEAL::Cycle 1: 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. @@ -27,14 +17,6 @@ DEAL::Cycle 2: 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. @@ -42,16 +24,6 @@ DEAL::Cycle 3: 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. diff --git a/tests/multigrid/step-16-06.output b/tests/multigrid/step-16-06.with_lapack=on.output similarity index 55% rename from tests/multigrid/step-16-06.output rename to tests/multigrid/step-16-06.with_lapack=on.output index b7f18bfa52..a6986c2c48 100644 --- a/tests/multigrid/step-16-06.output +++ b/tests/multigrid/step-16-06.with_lapack=on.output @@ -3,10 +3,6 @@ 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. @@ -14,12 +10,6 @@ DEAL::Cycle 1: 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. @@ -27,14 +17,6 @@ DEAL::Cycle 2: 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. @@ -42,16 +24,6 @@ DEAL::Cycle 3: 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. -- 2.39.5