From: Martin Kronbichler Date: Fri, 9 Jan 2009 10:42:00 +0000 (+0000) Subject: Simplified a the refine grid function a bit. X-Git-Tag: v8.0.0~8167 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=b03a9c7eac930c8f6981c770a8ddf4a1d0785761;p=dealii.git Simplified a the refine grid function a bit. git-svn-id: https://svn.dealii.org/trunk@18149 0785d39b-7218-0410-832d-ea1e28bc413d --- 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;