From e33e40b7dfe46ca66adf2bcf36a1143fe6512e83 Mon Sep 17 00:00:00 2001 From: kronbichler Date: Fri, 9 Jan 2009 10:42:00 +0000 Subject: [PATCH] Simplified a the refine grid function a bit. git-svn-id: https://svn.dealii.org/trunk@18149 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-31/step-31.cc | 11 ++++------- 1 file changed, 4 insertions(+), 7 deletions(-) diff --git a/deal.II/examples/step-31/step-31.cc b/deal.II/examples/step-31/step-31.cc index 3944da2da0..7e96df6b6b 100644 --- a/deal.II/examples/step-31/step-31.cc +++ b/deal.II/examples/step-31/step-31.cc @@ -347,7 +347,7 @@ namespace LinearSolvers vmult (VectorType &dst, const VectorType &src) const { - SolverControl solver_control (src.size(), 1e-7*src.l2_norm()); + SolverControl solver_control (src.size(), 5e-8*src.l2_norm()); SolverCG cg (solver_control); dst = 0; @@ -491,7 +491,7 @@ namespace LinearSolvers stokes_matrix (&S), m_inverse (&Mpinv), a_preconditioner (Apreconditioner), - tmp (stokes_matrix->block(1,1).matrix->RowMap()) + tmp (stokes_matrix->block(1,1).m()) {} @@ -2652,12 +2652,9 @@ void BoussinesqFlowProblem::refine_mesh (const unsigned int max_grid_level) // data vectors for refinement (in this // order). std::vector x_temperature (2); - x_temperature[0].reinit (temperature_solution); x_temperature[0] = temperature_solution; - x_temperature[1].reinit (temperature_solution); x_temperature[1] = old_temperature_solution; - TrilinosWrappers::BlockVector x_stokes(2); - x_stokes = stokes_solution; + TrilinosWrappers::BlockVector x_stokes = stokes_solution; SolutionTransfer temperature_trans(temperature_dof_handler); @@ -2702,7 +2699,7 @@ void BoussinesqFlowProblem::refine_mesh (const unsigned int max_grid_level) temperature_solution = tmp[0]; old_temperature_solution = tmp[1]; - + stokes_trans.interpolate (x_stokes, stokes_solution); rebuild_stokes_matrix = true; -- 2.39.5