]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Introduce template for preconditioner type in Chebyshev smoother.
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Thu, 27 Oct 2016 13:21:48 +0000 (15:21 +0200)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Sat, 29 Oct 2016 08:41:41 +0000 (10:41 +0200)
21 files changed:
include/deal.II/lac/precondition.h
tests/lac/linear_operator_08.with_cxx11=on.with_umfpack=off.output
tests/lac/linear_operator_08.with_cxx11=on.with_umfpack=on.output
tests/lac/precondition_chebyshev_01.with_lapack=true.output
tests/matrix_free/parallel_multigrid.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=3.output
tests/matrix_free/parallel_multigrid.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=7.output
tests/matrix_free/parallel_multigrid_adaptive_01.cc
tests/matrix_free/parallel_multigrid_adaptive_01.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=3.output
tests/matrix_free/parallel_multigrid_adaptive_01.with_trilinos=true.with_mpi=true.with_p4est=true.mpirun=7.output
tests/matrix_free/parallel_multigrid_adaptive_02.cc
tests/matrix_free/parallel_multigrid_adaptive_02.with_mpi=true.with_p4est=true.mpirun=3.output
tests/matrix_free/parallel_multigrid_adaptive_02.with_mpi=true.with_p4est=true.mpirun=7.output
tests/matrix_free/parallel_multigrid_adaptive_03.cc
tests/matrix_free/parallel_multigrid_adaptive_04.cc
tests/matrix_free/parallel_multigrid_adaptive_04.with_mpi=true.with_p4est=true.mpirun=3.output
tests/matrix_free/parallel_multigrid_mf.cc
tests/matrix_free/step-37.with_lapack=true.output
tests/multigrid/step-16-04.output [deleted file]
tests/multigrid/step-16-04.with_lapack=on.output [new file with mode: 0644]
tests/multigrid/step-16-05.with_lapack=on.output [moved from tests/multigrid/step-16-05.output with 55% similarity]
tests/multigrid/step-16-06.with_lapack=on.output [moved from tests/multigrid/step-16-06.output with 55% similarity]

index dface084ed06951dcbee55f67f98e509a2d37ee5..7f71ebe5258a0f93c4bd39883ddc23b0a1b640d1 100644 (file)
@@ -23,7 +23,8 @@
 #include <deal.II/base/utilities.h>
 #include <deal.II/base/parallel.h>
 #include <deal.II/base/template_constraints.h>
-#include <deal.II/lac/tridiagonal_matrix.h>
+#include <deal.II/base/thread_management.h>
+#include <deal.II/lac/diagonal_matrix.h>
 #include <deal.II/lac/solver_cg.h>
 #include <deal.II/lac/vector_memory.h>
 
@@ -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 <tt>vmult</tt> implementation.
- * However, it does not need to explicitly know the matrix entries.
+ * matrices. This preconditioner is based on an iteration of an inner
+ * preconditioner of type @p PreconditionerType with coefficients that are
+ * adapted to optimally cover an eigenvalue range between the largest
+ * eigenvalue down to a given lower eigenvalue specified by the optional
+ * parameter @p smoothing_range. The typical use case for the preconditioner
+ * is a Jacobi preconditioner specified through DiagonalMatrix, which is also
+ * the default value for the preconditioner. Note that if the degree variable
+ * is set to zero, the Chebyshev iteration corresponds to a Jacobi
+ * preconditioner (or the underlying preconditioner type) with relaxation
+ * parameter according to the specified smoothing range.
  *
- * This class is useful e.g. in multigrid smoother objects, since it is
- * trivially %parallel (assuming that matrix-vector products are %parallel).
+ * Besides the default choice of a pointwise Jacobi preconditioner, this class
+ * also allows for more advanced types of preconditioners, for example
+ * iterating block-Jacobi preconditioners in DG methods.
  *
- * @author Martin Kronbichler, 2009; extension for full compatibility with
+ * Apart from the inner preconditioner object, this iteration does not need
+ * access to matrix entries, which makes it an ideal ingredient for
+ * matrix-free computations. In that context, this class can be used as a
+ * multigrid smoother that is trivially %parallel (assuming that matrix-vector
+ * products are %parallel and the inner preconditioner is %parallel). Its use
+ * is demonstrated in the step-37 tutorial program.
+ *
+ * <h4>Algorithm execution</h4>
+ *
+ * The Chebyshev method relies on an estimate of the eigenvalues of the matrix
+ * which are computed during the first invocation of vmult(). The algorithm
+ * invokes a conjugate gradient solver so symmetry and positive definiteness
+ * of the (preconditioned) matrix system are strong requirements. The
+ * computation of eigenvalues needs to be deferred until the first vmult()
+ * invocation because temporary vectors of the same layout as the source and
+ * destination vectors are necessary for these computations and this
+ * information gets only available through vmult().
+ *
+ * The estimation of eigenvalues can also be bypassed by setting
+ * PreconditionChebyshev::AdditionalData::eig_cg_n_iterations to zero and
+ * providing sensible values for the largest eigenvalues in the field
+ * PreconditionChebyshev::AdditionalData::max_eigenvalue. If the range
+ * <tt>[max_eigenvalue/smoothing_range, max_eigenvalue]</tt> contains all
+ * eigenvalues of the preconditioned matrix system and the degree (i.e.,
+ * number of iterations) is high enough, this class can also be used as a
+ * direct solver. For an error estimation of the Chebyshev iteration that can
+ * be used to determine the number of iteration, see e.g. section 5.1 of
+ * R. S. Varga, Matrix iterative analysis, 2nd ed., Springer, 2009.
+ *
+ * <h4>Requirements on the templated classes</h4>
+ *
+ * The class MatrixType must be derived from Subscriptor because a
+ * SmartPointer to MatrixType is held in the class. In particular, this means
+ * that the matrix object needs to persist during the lifetime of
+ * PreconditionChebyshev. The preconditioner is held in a shared_ptr that is
+ * copied into the AdditionalData member variable of the class, so the
+ * variable used for initialization can safely be discarded after calling
+ * initialize(). Both the matrix and the preconditioner need to provide @p
+ * vmult functions for the matrix-vector product and @p m functions for
+ * accessing the number of rows in the (square) matrix. Furthermore, the
+ * matrix must provide <tt>el(i,i)</tt> methods for accessing the matrix
+ * diagonal in case the preconditioner type is a diagonal matrix. Even though
+ * it is highly recommended to pass the inverse diagonal entries inside a
+ * separate preconditioner object for implementing the Jacobi method (and the
+ * only possibility when computing in %parallel with MPI), there is a backward
+ * compatibility function that can extract the diagonal in case of a serial
+ * computation.
+ *
+ * @author Martin Kronbichler, 2009, 2016; extension for full compatibility with
  * LinearOperator class: Jean-Paul Pelteret, 2015
  */
-template <typename MatrixType=SparseMatrix<double>, class VectorType=Vector<double> >
+template <typename MatrixType=SparseMatrix<double>, typename VectorType=Vector<double>, typename PreconditionerType=DiagonalMatrix<VectorType> >
 class PreconditionChebyshev : public Subscriptor
 {
 public:
@@ -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<PreconditionerType> preconditioner;
   };
 
   PreconditionChebyshev ();
@@ -973,7 +1035,7 @@ private:
   /**
    * A pointer to the underlying matrix.
    */
-  SmartPointer<const MatrixType,PreconditionChebyshev<MatrixType,VectorType> > matrix_ptr;
+  SmartPointer<const MatrixType,PreconditionChebyshev<MatrixType,VectorType,PreconditionerType> > matrix_ptr;
 
   /**
    * Internal vector used for the <tt>vmult</tt> operation.
@@ -986,7 +1048,13 @@ private:
   mutable VectorType update2;
 
   /**
-   * Stores the additional data provided to the initialize function.
+   * Internal vector used for the <tt>vmult</tt> operation.
+   */
+  mutable VectorType update3;
+
+  /**
+   * Stores the additional data passed to the initialize function, obtained
+   * through a copy operation.
    */
   AdditionalData data;
 
@@ -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 <typename VectorType>
+    template <typename VectorType, typename PreconditionerType>
     inline
     void
-    vector_updates (const VectorType &src,
-                    const VectorType &matrix_diagonal_inverse,
-                    const bool       start_zero,
-                    const double     factor1,
-                    const double     factor2,
-                    VectorType       &update1,
-                    VectorType       &update2,
-                    VectorType       &dst)
+    vector_updates (const VectorType         &src,
+                    const PreconditionerType &preconditioner,
+                    const bool                start_zero,
+                    const double              factor1,
+                    const double              factor2,
+                    VectorType               &update1,
+                    VectorType               &update2,
+                    VectorType               &update3,
+                    VectorType               &dst)
     {
       if (start_zero)
         {
-          dst.equ (factor2, src);
-          dst.scale (matrix_diagonal_inverse);
+          update1.equ (factor2, src);
+          preconditioner.vmult(dst, update1);
           update1.equ(-1.,dst);
         }
       else
         {
           update2 -= src;
-          update2.scale (matrix_diagonal_inverse);
+          preconditioner.vmult(update3, update2);
+          update2 = update3;
           if (factor1 == 0.)
             update1.equ(factor2, update2);
           else
@@ -1738,61 +1823,89 @@ namespace internal
       const VectorUpdater<Number> &updater;
     };
 
-    // selection for deal.II vector
+    // selection for diagonal matrix around deal.II vector
     template <typename Number>
     inline
     void
     vector_updates (const ::dealii::Vector<Number> &src,
-                    const ::dealii::Vector<Number> &matrix_diagonal_inverse,
+                    const DiagonalMatrix<::dealii::Vector<Number> > &jacobi,
                     const bool    start_zero,
                     const double  factor1,
                     const double  factor2,
                     ::dealii::Vector<Number> &update1,
                     ::dealii::Vector<Number> &update2,
+                    ::dealii::Vector<Number> &,
                     ::dealii::Vector<Number> &dst)
     {
-      VectorUpdater<Number> upd(src.begin(), matrix_diagonal_inverse.begin(),
+      VectorUpdater<Number> upd(src.begin(), jacobi.get_vector().begin(),
                                 start_zero, factor1, factor2,
                                 update1.begin(), update2.begin(), dst.begin());
       VectorUpdatesRange<Number>(upd, src.size());
     }
 
-    // selection for parallel deal.II vector
+    // selection for diagonal matrix around parallel deal.II vector
     template <typename Number>
     inline
     void
     vector_updates (const LinearAlgebra::distributed::Vector<Number> &src,
-                    const LinearAlgebra::distributed::Vector<Number> &matrix_diagonal_inverse,
+                    const DiagonalMatrix<LinearAlgebra::distributed::Vector<Number> > &jacobi,
                     const bool    start_zero,
                     const double  factor1,
                     const double  factor2,
                     LinearAlgebra::distributed::Vector<Number> &update1,
                     LinearAlgebra::distributed::Vector<Number> &update2,
+                    LinearAlgebra::distributed::Vector<Number> &,
                     LinearAlgebra::distributed::Vector<Number> &dst)
     {
-      VectorUpdater<Number> upd(src.begin(), matrix_diagonal_inverse.begin(),
+      VectorUpdater<Number> upd(src.begin(), jacobi.get_vector().begin(),
                                 start_zero, factor1, factor2,
                                 update1.begin(), update2.begin(), dst.begin());
       VectorUpdatesRange<Number>(upd, src.local_size());
     }
 
-    template <typename VectorType>
-    struct DiagonalPreconditioner
+    template <typename MatrixType, typename VectorType, typename PreconditionerType>
+    inline
+    void
+    initialize_preconditioner(const MatrixType                          &matrix,
+                              std_cxx11::shared_ptr<PreconditionerType> &preconditioner,
+                              VectorType &)
     {
-      DiagonalPreconditioner (const VectorType &vector)
-        :
-        diagonal_vector(vector)
-      {}
+      (void)matrix;
+      (void)preconditioner;
+      AssertThrow(preconditioner.get() != NULL, ExcNotInitialized());
+      AssertDimension(preconditioner->m(), matrix.m());
+    }
 
-      void vmult (VectorType       &dst,
-                  const VectorType &src) const
-      {
-        dst = src;
-        dst.scale(diagonal_vector);
-      }
+    template <typename MatrixType, typename VectorType>
+    inline
+    void
+    initialize_preconditioner(const MatrixType                                   &matrix,
+                              std_cxx11::shared_ptr<DiagonalMatrix<VectorType> > &preconditioner,
+                              VectorType                                         &diagonal_inverse)
+    {
+      if (preconditioner.get() == NULL ||
+          preconditioner->m() != matrix.m())
+        {
+          if (preconditioner.get() == NULL)
+            preconditioner.reset(new DiagonalMatrix<VectorType>());
 
-      const VectorType &diagonal_vector;
-    };
+          Assert(preconditioner->m() == 0,
+                 ExcMessage("Preconditioner appears to be initialized but not sized correctly"));
+
+          // Check if we can initialize from vector that then gets set to zero
+          // as the matrix will own the memory
+          preconditioner->reinit(diagonal_inverse);
+          diagonal_inverse.reinit(0);
+
+          // This part only works in serial
+          if (preconditioner->m() != matrix.m())
+            {
+              preconditioner->get_vector().reinit(matrix.m());
+              for (typename VectorType::size_type i=0; i<matrix.m(); ++i)
+                preconditioner->get_vector()(i) = 1./matrix.el(i,i);
+            }
+        }
+    }
 
     struct EigenvalueTracker
     {
@@ -1809,9 +1922,9 @@ namespace internal
 
 
 
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, class VectorType, typename PreconditionerType>
 inline
-PreconditionChebyshev<MatrixType,VectorType>::AdditionalData::
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::AdditionalData::
 AdditionalData (const unsigned int degree,
                 const double       smoothing_range,
                 const bool         nonzero_starting,
@@ -1829,9 +1942,9 @@ AdditionalData (const unsigned int degree,
 
 
 
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
 inline
-PreconditionChebyshev<MatrixType,VectorType>::PreconditionChebyshev ()
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::PreconditionChebyshev ()
   :
   is_initialized (false)
 {
@@ -1844,24 +1957,34 @@ PreconditionChebyshev<MatrixType,VectorType>::PreconditionChebyshev ()
 
 
 
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
 inline
 void
-PreconditionChebyshev<MatrixType,VectorType>::initialize
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::initialize
 (const MatrixType     &matrix,
  const AdditionalData &additional_data)
 {
   matrix_ptr = &matrix;
   data = additional_data;
-  if (data.matrix_diagonal_inverse.size() != matrix.m())
-    {
-      Assert(data.matrix_diagonal_inverse.size() == 0,
-             ExcMessage("Matrix diagonal vector set but not sized correctly"));
-      data.matrix_diagonal_inverse.reinit(matrix.m());
-      for (unsigned int i=0; i<matrix.m(); ++i)
-        data.matrix_diagonal_inverse(i) = 1./matrix.el(i,i);
-    }
+  internal::PreconditionChebyshev::initialize_preconditioner(matrix,
+                                                             data.preconditioner,
+                                                             data.matrix_diagonal_inverse);
+  is_initialized = false;
+}
+
+
+
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
+inline
+void
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::estimate_eigenvalues
+(const VectorType &src) const
+{
+  Assert(is_initialized == false, ExcInternalError());
+  Assert(data.preconditioner.get() != NULL, ExcNotInitialized());
 
+  update1.reinit(src);
+  update2 = src;
 
   // calculate largest eigenvalue using a hand-tuned CG iteration on the
   // matrix weighted by its diagonal. we start with a vector that consists of
@@ -1869,57 +1992,48 @@ PreconditionChebyshev<MatrixType,VectorType>::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<VectorType> memory;
-      VectorType *rhs = memory.alloc();
-      VectorType *dummy = memory.alloc();
-      rhs->reinit(data.matrix_diagonal_inverse);
-      dummy->reinit(data.matrix_diagonal_inverse);
-
-      // heuristically, a right hand side close to a constant has been shown
-      // to quickly reveal the largest eigenvalue. however, avoid to use the
-      // exact constant because that might be not in the range space of some
-      // matrices (purely Neumann matrices with constant mode filtered out by
-      // orthogonal projection in the matrix-vector product)
-      *rhs = 1./std::sqrt(static_cast<double>(matrix.m()));
-      if (rhs->locally_owned_elements().is_element(0))
-        (*rhs)(0) = 0.;
-      rhs->compress(VectorOperation::insert);
+      ReductionControl control (data.eig_cg_n_iterations,
+                                std::sqrt(std::numeric_limits<typename VectorType::value_type>::epsilon()),
+                                1e-10, false, false);
 
       internal::PreconditionChebyshev::EigenvalueTracker eigenvalue_tracker;
-      SolverCG<VectorType> solver (control, memory);
+      SolverCG<VectorType> solver (control);
       solver.connect_eigenvalues_slot(std_cxx11::bind(&internal::PreconditionChebyshev::EigenvalueTracker::slot,
                                                       &eigenvalue_tracker,
                                                       std_cxx11::_1));
-      internal::PreconditionChebyshev::DiagonalPreconditioner<VectorType>
-      preconditioner(data.matrix_diagonal_inverse);
+
+      // catch case the rhs vector is exactly zero (no eigenvalue estimate
+      // possible) by setting a different initial guess
+      if (update2 * update2 == typename VectorType::value_type())
+        {
+          update2 = 1.;
+          if (update2.locally_owned_elements().is_element(0))
+            update2(0) = 0.;
+        }
+
       try
         {
-          solver.solve(matrix, *dummy, *rhs, preconditioner);
+          solver.solve(*matrix_ptr, update1, update2, *data.preconditioner);
         }
       catch (SolverControl::NoConvergence &)
         {
         }
 
-      memory.free(dummy);
-      memory.free(rhs);
-
       // read the eigenvalues from the attached eigenvalue tracker
       if (eigenvalue_tracker.values.empty())
         min_eigenvalue = max_eigenvalue = 1;
       else
         {
           min_eigenvalue = eigenvalue_tracker.values.front();
-          max_eigenvalue = eigenvalue_tracker.values.back();
-        }
 
-      // include a safety factor since the CG method will in general not be
-      // converged
-      max_eigenvalue *= 1.2;
+          // include a safety factor since the CG method will in general not be
+          // converged
+          max_eigenvalue = 1.2*eigenvalue_tracker.values.back();
+        }
     }
   else
     {
@@ -1930,23 +2044,35 @@ PreconditionChebyshev<MatrixType,VectorType>::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<PreconditionChebyshev<MatrixType,VectorType,PreconditionerType> *>(this)->delta = (max_eigenvalue-alpha)*0.5;
+  const_cast<PreconditionChebyshev<MatrixType,VectorType,PreconditionerType> *>(this)->theta = (max_eigenvalue+alpha)*0.5;
+
+  // We do not need the third auxiliary vector in case we have a
+  // DiagonalMatrix as preconditioner and use deal.II's own vectors
+  if (types_are_equal<PreconditionerType,DiagonalMatrix<VectorType> >::value == false ||
+      (types_are_equal<VectorType,dealii::Vector<typename VectorType::value_type> >::value == false
+       &&
+       types_are_equal<VectorType,LinearAlgebra::distributed::Vector<typename VectorType::value_type> >::value == false
+      ))
+    update3.reinit (src, true);
+
+  const_cast<PreconditionChebyshev<MatrixType,VectorType,PreconditionerType> *>(this)->is_initialized = true;
 }
 
 
 
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
 inline
 void
-PreconditionChebyshev<MatrixType,VectorType>::do_chebyshev_loop(VectorType       &dst,
-    const VectorType &src) const
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>
+::do_chebyshev_loop(VectorType       &dst,
+                    const VectorType &src) const
 {
+  // if delta is zero, we do not need to iterate because the updates will be
+  // zero
+  if (std::abs(delta) < 1e-40)
+    return;
+
   double rhok  = delta / theta,  sigma = theta / delta;
   for (unsigned int k=0; k<data.degree; ++k)
     {
@@ -1955,18 +2081,18 @@ PreconditionChebyshev<MatrixType,VectorType>::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 <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
 inline
 void
-PreconditionChebyshev<MatrixType,VectorType>::do_transpose_chebyshev_loop(VectorType       &dst,
-    const VectorType &src) const
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>
+::do_transpose_chebyshev_loop(VectorType       &dst,
+                              const VectorType &src) const
 {
   double rhok  = delta / theta,  sigma = theta / delta;
   for (unsigned int k=0; k<data.degree; ++k)
@@ -1976,119 +2102,131 @@ PreconditionChebyshev<MatrixType,VectorType>::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 <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
 inline
 void
-PreconditionChebyshev<MatrixType,VectorType>::vmult (VectorType       &dst,
-                                                     const VectorType &src) const
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>
+::vmult (VectorType       &dst,
+         const VectorType &src) const
 {
-  Assert (is_initialized, ExcMessage("Preconditioner not initialized"));
+  Threads::Mutex::ScopedLock lock(mutex);
+  if (is_initialized == false)
+    estimate_eigenvalues(src);
+
   internal::PreconditionChebyshev::vector_updates
-  (src, data.matrix_diagonal_inverse, true, 0., 1./theta, update1,
-   update2, dst);
+  (src, *data.preconditioner, true, 0., 1./theta, update1, update2, update3, dst);
 
   do_chebyshev_loop(dst, src);
 }
 
 
 
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
 inline
 void
-PreconditionChebyshev<MatrixType,VectorType>::Tvmult (VectorType       &dst,
-                                                      const VectorType &src) const
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>
+::Tvmult (VectorType       &dst,
+          const VectorType &src) const
 {
-  Assert (is_initialized, ExcMessage("Preconditioner not initialized"));
+  Threads::Mutex::ScopedLock lock(mutex);
+  if (is_initialized == false)
+    estimate_eigenvalues(src);
+
   internal::PreconditionChebyshev::vector_updates
-  (src, data.matrix_diagonal_inverse, true, 0., 1./theta, update1,
-   update2, dst);
+  (src, *data.preconditioner, true, 0., 1./theta, update1, update2, update3, dst);
 
   do_transpose_chebyshev_loop(dst, src);
 }
 
 
 
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
 inline
 void
-PreconditionChebyshev<MatrixType,VectorType>::step (VectorType       &dst,
-                                                    const VectorType &src) const
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>
+::step (VectorType       &dst,
+        const VectorType &src) const
 {
-  Assert (is_initialized, ExcMessage("Preconditioner not initialized"));
+  Threads::Mutex::ScopedLock lock(mutex);
+  if (is_initialized == false)
+    estimate_eigenvalues(src);
+
   if (!dst.all_zero())
     {
       matrix_ptr->vmult (update2, dst);
       internal::PreconditionChebyshev::vector_updates
-      (src, data.matrix_diagonal_inverse, false, 0., 1./theta, update1,
-       update2, dst);
+      (src, *data.preconditioner, false, 0., 1./theta, update1, update2, update3, dst);
     }
   else
     internal::PreconditionChebyshev::vector_updates
-    (src, data.matrix_diagonal_inverse, true, 0., 1./theta, update1,
-     update2, dst);
+    (src, *data.preconditioner, true, 0., 1./theta, update1, update2, update3, dst);
 
   do_chebyshev_loop(dst, src);
 }
 
 
 
-template <typename MatrixType, class VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
 inline
 void
-PreconditionChebyshev<MatrixType,VectorType>::Tstep (VectorType       &dst,
-                                                     const VectorType &src) const
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>
+::Tstep (VectorType       &dst,
+         const VectorType &src) const
 {
-  Assert (is_initialized, ExcMessage("Preconditioner not initialized"));
+  Threads::Mutex::ScopedLock lock(mutex);
+  if (is_initialized == false)
+    estimate_eigenvalues(src);
+
   if (!dst.all_zero())
     {
       matrix_ptr->Tvmult (update2, dst);
       internal::PreconditionChebyshev::vector_updates
-      (src, data.matrix_diagonal_inverse, false, 0., 1./theta, update1,
-       update2, dst);
+      (src, *data.preconditioner, false, 0., 1./theta, update1, update2, update3, dst);
     }
   else
     internal::PreconditionChebyshev::vector_updates
-    (src, data.matrix_diagonal_inverse, true, 0., 1./theta, update1,
-     update2, dst);
+    (src, *data.preconditioner, true, 0., 1./theta, update1, update2, update3, dst);
 
   do_transpose_chebyshev_loop(dst, src);
 }
 
 
 
-template <typename MatrixType, typename VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
 inline
-void PreconditionChebyshev<MatrixType,VectorType>::clear ()
+void
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::clear ()
 {
   is_initialized = false;
   matrix_ptr = 0;
   data.matrix_diagonal_inverse.reinit(0);
+  data.preconditioner.reset();
   update1.reinit(0);
   update2.reinit(0);
+  update3.reinit(0);
 }
 
 
-template <typename MatrixType, typename VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
 inline
-typename PreconditionChebyshev<MatrixType,VectorType>::size_type
-PreconditionChebyshev<MatrixType,VectorType>::m () const
+typename PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::size_type
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::m () const
 {
   Assert (matrix_ptr!=0, ExcNotInitialized());
   return matrix_ptr->m();
 }
 
 
-template <typename MatrixType, typename VectorType>
+template <typename MatrixType, typename VectorType, typename PreconditionerType>
 inline
-typename PreconditionChebyshev<MatrixType,VectorType>::size_type
-PreconditionChebyshev<MatrixType,VectorType>::n () const
+typename PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::size_type
+PreconditionChebyshev<MatrixType,VectorType,PreconditionerType>::n () const
 {
   Assert (matrix_ptr!=0, ExcNotInitialized());
   return matrix_ptr->n();
index 84bdb2e31720e3f81d988bdc765c267722735d82..55468cb785450465244d8741d005fb7e05f06825 100644 (file)
@@ -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
index 459bc5bca8d87e695c2138450ef91c94c536969a..d47fd5d5edbc298d3275eb05d284c5e24273295b 100644 (file)
@@ -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
index f09bc650ad51ab7749f420acc6aaa643326b38af..2d5be039fca20753513ebd42e5babd67d2407f93 100644 (file)
@@ -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 
index 9d5dea07fe1ec2bd0258abc050c744c68c765e60..c1311be95595a179f2c2e22b8fdae328006c7eb1 100644 (file)
 
 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
index 9d5dea07fe1ec2bd0258abc050c744c68c765e60..c1311be95595a179f2c2e22b8fdae328006c7eb1 100644 (file)
 
 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
index 1ac481fb7673a14b6601ecb4be6cf56cdf7a4dbc..1cd256990e81a6653f3e76440d5466661685210b 100644 (file)
@@ -530,6 +530,9 @@ void do_test (const DoFHandler<dim>  &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<LinearAlgebra::distributed::Vector<double> > solver(control);
     solver.solve(fine_matrix, sol, in, preconditioner);
index e9349c579a69ba876f226c1924c28f5b6e7f0f47..16d865c56ac6bd81a95c708c9f02413c5b1d35d6 100644 (file)
 
 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
index e9349c579a69ba876f226c1924c28f5b6e7f0f47..16d865c56ac6bd81a95c708c9f02413c5b1d35d6 100644 (file)
 
 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
index 03c6d6e6f15ea37e0259ee83636ad53790878788..bf08a47b3bd9e5944cdfbbf41704e30b4dfd5e52 100644 (file)
@@ -385,6 +385,9 @@ void do_test (const DoFHandler<dim>  &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<LinearAlgebra::distributed::Vector<double> > solver(control);
     solver.solve(fine_matrix, sol, in, preconditioner);
index e9349c579a69ba876f226c1924c28f5b6e7f0f47..16d865c56ac6bd81a95c708c9f02413c5b1d35d6 100644 (file)
 
 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
index e9349c579a69ba876f226c1924c28f5b6e7f0f47..16d865c56ac6bd81a95c708c9f02413c5b1d35d6 100644 (file)
 
 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
index 133f617b0d4321e0fb3117452c742a183a22d795..ea83846e7cce0ea4f1b85f56d9d57a8c7247304a 100644 (file)
@@ -507,9 +507,6 @@ void do_test (const DoFHandler<dim>  &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<LinearAlgebra::distributed::Vector<number> >
@@ -531,6 +528,7 @@ void do_test (const DoFHandler<dim>  &dof, const bool threaded)
   {
     // avoid output from inner (coarse-level) solver
     deallog.depth_file(3);
+
     ReductionControl control(30, 1e-20, 1e-7);
     SolverCG<LinearAlgebra::distributed::Vector<double> > solver(control);
     solver.solve(fine_matrix, sol, in, preconditioner);
index b189fa73d7eabf8980b68760e863c6c29f110b30..6bf48b88584f7908033c943113d73d32c564a618 100644 (file)
@@ -385,6 +385,9 @@ void do_test (const DoFHandler<dim>  &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<LinearAlgebra::distributed::Vector<double> > solver(control);
     solver.solve(fine_matrix, sol, in, preconditioner);
index e9349c579a69ba876f226c1924c28f5b6e7f0f47..fa84d2c1324956516a918526d75164429e94b5ae 100644 (file)
 
 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
index 3a81c512f40f0a61a83b8b30b67693518d484ae5..c7dc7aaef7ebcd06705330bf9d87682c0d5745e7 100644 (file)
@@ -364,7 +364,6 @@ void do_test (const DoFHandler<dim>  &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<LinearAlgebra::distributed::Vector<number> >
index 76f9274e64735ef43805fd1198fffd85c18ae8c4..f67b52a7ea0b7f94bd614cfb6ca27532952691b3 100644 (file)
@@ -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 (file)
index 7fb681c..0000000
+++ /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 (file)
index 0000000..29325f5
--- /dev/null
@@ -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.
similarity index 55%
rename from tests/multigrid/step-16-05.output
rename to tests/multigrid/step-16-05.with_lapack=on.output
index b7f18bfa5216362f5406d9bace2b7b8143a6769f..a6986c2c48eb1216bc1a406a74c2b42d818c8bd2 100644 (file)
@@ -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.
similarity index 55%
rename from tests/multigrid/step-16-06.output
rename to tests/multigrid/step-16-06.with_lapack=on.output
index b7f18bfa5216362f5406d9bace2b7b8143a6769f..a6986c2c48eb1216bc1a406a74c2b42d818c8bd2 100644 (file)
@@ -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.

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.