From 646114346dba96728647c0d6026a918b8a7c99fb Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Mon, 6 Apr 1998 09:05:16 +0000 Subject: [PATCH] Bug fixes and enhancements. git-svn-id: https://svn.dealii.org/trunk@141 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/numerics/base.cc | 31 ++++++++++++++----------- 1 file changed, 18 insertions(+), 13 deletions(-) diff --git a/deal.II/deal.II/source/numerics/base.cc b/deal.II/deal.II/source/numerics/base.cc index 7b555ba916..bae1b0ee03 100644 --- a/deal.II/deal.II/source/numerics/base.cc +++ b/deal.II/deal.II/source/numerics/base.cc @@ -159,6 +159,7 @@ void ProblemBase::integrate_difference (const Function &exact_sol FEValues::UpdateStruct update_flags; update_flags.q_points = true; + update_flags.jacobians = true; update_flags.JxW_values = true; FEValues fe_values(fe, q, update_flags); @@ -172,7 +173,7 @@ void ProblemBase::integrate_difference (const Function &exact_sol // incrementing is faster than conversion...) DoFHandler::active_cell_iterator cell = dof_handler->begin_active(), endc = dof_handler->end(); - Triangulation::active_cell_iterator tria_cell = dof_handler->get_tria().begin(); + Triangulation::active_cell_iterator tria_cell = dof_handler->get_tria().begin_active(); for (; cell != endc; ++cell, ++tria_cell) { @@ -210,7 +211,7 @@ void ProblemBase::integrate_difference (const Function &exact_sol const dFMatrix & shape_values = fe_values.get_shape_values(); vector dof_values; cell->get_dof_values (solution, dof_values); - + vector psi; // in praxi: first compute @@ -219,10 +220,15 @@ void ProblemBase::integrate_difference (const Function &exact_sol psi); // then subtract finite element // solution - const vector &JxW_values = fe_values.get_JxW_values(); for (unsigned int j=0; j