From: frohne Date: Thu, 20 Oct 2011 18:54:19 +0000 (+0000) Subject: make step-41 work X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=b30b8f3295d5aabe1d87551684cef008b4a6a27d;p=dealii-svn.git make step-41 work git-svn-id: https://svn.dealii.org/trunk@24647 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-41/step-41.cc b/deal.II/examples/step-41/step-41.cc index 34f366248b..a790098d59 100644 --- a/deal.II/examples/step-41/step-41.cc +++ b/deal.II/examples/step-41/step-41.cc @@ -710,10 +710,9 @@ template void Step4::solve () { ReductionControl reduction_control (100, 1e-12, 1e-2); - SolverCG solver (reduction_control); -// SolverBicgstab<> solver_bicgstab (reduction_control); - TrilinosWrappers::PreconditionSSOR precondition; - precondition.initialize (system_matrix, 1.2); + SolverCG solver (reduction_control); + TrilinosWrappers::PreconditionAMG precondition; + precondition.initialize (system_matrix); solver.solve (system_matrix, solution, system_rhs, precondition); @@ -768,10 +767,10 @@ void Step4::output_results (TrilinosWrappers::Vector vector_to_plot, const (title + ".vtk").c_str ()); data_out.write_vtk (output_vtk); - std::ofstream output_gnuplot (dim == 2 ? - (title + ".gp").c_str () : - (title + ".gp").c_str ()); - data_out.write_gnuplot (output_gnuplot); +// std::ofstream output_gnuplot (dim == 2 ? +// (title + ".gp").c_str () : +// (title + ".gp").c_str ()); +// data_out.write_gnuplot (output_gnuplot); } @@ -827,14 +826,9 @@ void Step4::run () // << solution (k) << ", " // << system_rhs.l2_norm () // <::run () if (resid_vector (k) > 0) resid_vector (k) = 0; - std::ostringstream filename_residuum; - filename_residuum << "residuum_"; - filename_residuum << i; - output_results (resid_vector, filename_residuum.str ()); + std::cout<< "Create Output:" <