From f7509fac41a898679833d65b61746ea01fa8c895 Mon Sep 17 00:00:00 2001 From: bangerth Date: Tue, 28 Oct 2008 19:46:24 +0000 Subject: [PATCH] Take over some changes from step-31. git-svn-id: https://svn.dealii.org/trunk@17377 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-32/step-32.cc | 17 ++++++++++------- 1 file changed, 10 insertions(+), 7 deletions(-) diff --git a/deal.II/examples/step-32/step-32.cc b/deal.II/examples/step-32/step-32.cc index 73415567ba..cbe3838b69 100644 --- a/deal.II/examples/step-32/step-32.cc +++ b/deal.II/examples/step-32/step-32.cc @@ -1212,19 +1212,20 @@ void BoussinesqFlowProblem::solve () pcout << " Solving..." << std::endl; { - LinearSolvers::InverseMatrix + const LinearSolvers::InverseMatrix mp_inverse (stokes_preconditioner_matrix.block(1,1), *Mp_preconditioner); - LinearSolvers::BlockSchurPreconditioner + const LinearSolvers::BlockSchurPreconditioner preconditioner (stokes_matrix, mp_inverse, *Amg_preconditioner); SolverControl solver_control (stokes_matrix.m(), 1e-6*stokes_rhs.l2_norm()); - SolverGMRES gmres(solver_control, - SolverGMRES::AdditionalData(100)); + SolverGMRES + gmres (solver_control, + SolverGMRES::AdditionalData(100)); gmres.solve(stokes_matrix, stokes_solution, stokes_rhs, preconditioner); @@ -1244,6 +1245,9 @@ void BoussinesqFlowProblem::solve () GridTools::minimal_cell_diameter(triangulation) / std::max (get_maximal_velocity(), 1.e-5); + std::cout << " " << "Time step: " << time_step + << std::endl; + temperature_solution = old_temperature_solution; @@ -1465,7 +1469,6 @@ void BoussinesqFlowProblem::run () { pcout << "Timestep " << timestep_number << ": t=" << time - << ", dt=" << time_step << std::endl; assemble_stokes_system (); -- 2.39.5