bicgstab.solve(M, solution, rhs, PreconditionIdentity());
}
- solution.print(deallog);
+ solution.print(deallog.get_file_stream());
Vector<double> res(4);
M.residual(res, solution, rhs);
SolverBicgstab<> bicgstab(control);
bicgstab.solve(M1, solution, rhs, PreconditionIdentity());
}
- solution.print(deallog);
+ solution.print(deallog.get_file_stream());
M1.residual(res, solution, rhs);
deallog << "residual=" << res.l2_norm() << std::endl;
}
DEAL:Bicgstab::Starting value 2.000
DEAL:Bicgstab::Convergence step 4 value 0
-DEAL::1.000 0.1000 0.09091 0.02381
+1.000 0.1000 0.09091 0.02381
DEAL::residual=0
DEAL:Bicgstab::Starting value 2.000e+10
DEAL:Bicgstab::Convergence step 4 value 0.001939
-DEAL::1.000 0.1000 0.09091 0.02381
+1.000 0.1000 0.09091 0.02381
DEAL::residual=0.001937
VectorType v(N);
v = 1.0;
auto niter = kinsol.solve(v);
- deallog << v << std::endl;
+ v.print(deallog.get_file_stream());
deallog << "Converged in " << niter << " iterations." << std::endl;
}
-DEAL::1.00000 2.00000 3.00000 4.00000 5.00000 6.00000 7.00000 8.00000 9.00000 10.0000
-DEAL::
+1.00000 2.00000 3.00000 4.00000 5.00000 6.00000 7.00000 8.00000 9.00000 10.0000
DEAL::Converged in 23 iterations.