mf.apply_inverse (inverse, in);
- SolverControl control(1000, 1e-10);
+ SolverControl control(1000, 1e-12);
+ // do not output iteration numbers to log file because these are too
+ // sensitive
+ std::ostringstream stream;
+ deallog.attach(stream);
SolverCG<Vector<number> > solver(control);
solver.solve (mf, reference, in, PreconditionIdentity());
+ deallog.attach(logfile);
inverse -= reference;
const double diff_norm = inverse.linfty_norm() / reference.linfty_norm();
DEAL:2d::Testing FE_DGQ<2>(1)
-DEAL:2d:cg::Starting value 6.32
-DEAL:2d:cg::Convergence step 48 value 0
+
DEAL:2d::Norm of difference: 0
DEAL:2d::
DEAL:2d:float::Testing FE_DGQ<2>(1)
-DEAL:2d:float:cg::Starting value 6.15
-DEAL:2d:float:cg::Convergence step 57 value 0
+
DEAL:2d:float::Norm of difference: 0
DEAL:2d:float::
DEAL:2d::Testing FE_DGQ<2>(2)
-DEAL:2d:cg::Starting value 8.62
-DEAL:2d:cg::Convergence step 105 value 0
+
DEAL:2d::Norm of difference: 0
DEAL:2d::
DEAL:2d:float::Testing FE_DGQ<2>(2)
-DEAL:2d:float:cg::Starting value 9.05
-DEAL:2d:float:cg::Convergence step 133 value 0
+
DEAL:2d:float::Norm of difference: 0
DEAL:2d:float::
DEAL:2d::Testing FE_DGQ<2>(4)
-DEAL:2d:cg::Starting value 14.6
-DEAL:2d:cg::Convergence step 311 value 0
+
DEAL:2d::Norm of difference: 0
DEAL:2d::
DEAL:2d:float::Testing FE_DGQ<2>(4)
-DEAL:2d:float:cg::Starting value 14.7
-DEAL:2d:float:cg::Convergence step 376 value 0
+
DEAL:2d:float::Norm of difference: 0
DEAL:2d:float::
DEAL:3d::Testing FE_DGQ<3>(1)
-DEAL:3d:cg::Starting value 13.7
-DEAL:3d:cg::Convergence step 124 value 0
+
DEAL:3d::Norm of difference: 0
DEAL:3d::
DEAL:3d::Testing FE_DGQ<3>(2)
-DEAL:3d:cg::Starting value 13.5
-DEAL:3d:cg::Convergence step 241 value 0
+
DEAL:3d::Norm of difference: 0
DEAL:3d::
mf.apply_inverse (inverse, in);
SolverControl control(1000, 1e-12);
+ std::ostringstream stream;
+ deallog.attach(stream);
SolverCG<Vector<number> > solver(control);
solver.solve (mf, reference, in, PreconditionIdentity());
+ deallog.attach(logfile);
inverse -= reference;
const double diff_norm = inverse.linfty_norm() / reference.linfty_norm();
DEAL:2d::Testing FESystem<2>[FE_DGQ<2>(1)^2]
-DEAL:2d:cg::Starting value 7.90
-DEAL:2d:cg::Convergence step 7 value 0
+
DEAL:2d::Norm of difference: 0
DEAL:2d::
DEAL:2d::Testing FESystem<2>[FE_DGQ<2>(2)^2]
-DEAL:2d:cg::Starting value 11.6
-DEAL:2d:cg::Convergence step 15 value 0
+
DEAL:2d::Norm of difference: 0
DEAL:2d::
DEAL:2d::Testing FESystem<2>[FE_DGQ<2>(4)^2]
-DEAL:2d:cg::Starting value 19.1
-DEAL:2d:cg::Convergence step 46 value 0
+
DEAL:2d::Norm of difference: 0
DEAL:2d::
DEAL:3d::Testing FESystem<3>[FE_DGQ<3>(1)^3]
-DEAL:3d:cg::Starting value 25.0
-DEAL:3d:cg::Convergence step 10 value 0
+
DEAL:3d::Norm of difference: 0
DEAL:3d::
DEAL:3d::Testing FESystem<3>[FE_DGQ<3>(2)^3]
-DEAL:3d:cg::Starting value 24.3
-DEAL:3d:cg::Convergence step 32 value 0
+
DEAL:3d::Norm of difference: 0
DEAL:3d::
mf.apply_inverse (inverse, in);
SolverControl control(10000, 1e-12);
+ std::ostringstream stream;
+ deallog.attach(stream);
SolverCG<Vector<number> > solver(control);
solver.solve (mf, reference, in, PreconditionIdentity());
+ deallog.attach(logfile);
inverse -= reference;
const double diff_norm = inverse.linfty_norm() / reference.linfty_norm();
DEAL:2d::Testing FESystem<2>[FE_DGQ<2>(1)^3]
-DEAL:2d:cg::Starting value 10.6
-DEAL:2d:cg::Convergence step 108 value 0
+
DEAL:2d::Norm of difference: 0
DEAL:2d::
DEAL:2d::Testing FESystem<2>[FE_DGQ<2>(2)^3]
-DEAL:2d:cg::Starting value 15.1
-DEAL:2d:cg::Convergence step 252 value 0
+
DEAL:2d::Norm of difference: 0
DEAL:2d::
DEAL:2d::Testing FESystem<2>[FE_DGQ<2>(4)^3]
-DEAL:2d:cg::Starting value 25.4
-DEAL:2d:cg::Convergence step 699 value 0
+
DEAL:2d::Norm of difference: 0
DEAL:2d::
DEAL:3d::Testing FESystem<3>[FE_DGQ<3>(1)^3]
-DEAL:3d:cg::Starting value 23.7
-DEAL:3d:cg::Convergence step 412 value 0
+
DEAL:3d::Norm of difference: 0
DEAL:3d::
DEAL:3d::Testing FESystem<3>[FE_DGQ<3>(2)^3]
-DEAL:3d:cg::Starting value 23.7
-DEAL:3d:cg::Convergence step 1019 value 0
+
DEAL:3d::Norm of difference: 0
DEAL:3d::