From 5a7c626287a1f1e435016142e363d43a6ce0fa32 Mon Sep 17 00:00:00 2001 From: wolf Date: Sat, 5 Jun 1999 20:53:03 +0000 Subject: [PATCH] Hopefully fix a bug in arrays that were accessed in the wrong order. git-svn-id: https://svn.dealii.org/trunk@1370 0785d39b-7218-0410-832d-ea1e28bc413d --- .../source/numerics/error_estimator.cc | 29 ++++++++++--------- 1 file changed, 16 insertions(+), 13 deletions(-) diff --git a/deal.II/deal.II/source/numerics/error_estimator.cc b/deal.II/deal.II/source/numerics/error_estimator.cc index 587ff50813..a1ecfd31b6 100644 --- a/deal.II/deal.II/source/numerics/error_estimator.cc +++ b/deal.II/deal.II/source/numerics/error_estimator.cc @@ -317,11 +317,11 @@ integrate_over_regular_face (const active_cell_iterator &cell, fe_face_values_cell.get_function_grads (solution, psi); else { - vector > > tmp (n_components, - psi); + vector > > tmp (n_q_points, + vector >(n_components)); fe_face_values_cell.get_function_grads (solution, tmp); - - psi.swap (tmp[selected_component]); + for (unsigned int i=0; i > > tmp (n_components, - neighbor_psi); + vector > > tmp (n_q_points, + vector >(n_components)); fe_face_values_neighbor.get_function_grads (solution, tmp); - neighbor_psi.swap (tmp[selected_component]); + for (unsigned int i=0; i > > tmp (n_components, - psi); + vector > > tmp (n_q_points, + vector >(n_components)); fe_subface_values.get_function_grads (solution, tmp); - psi.swap (tmp[selected_component]); + for (unsigned int i=0; i > > tmp (n_components, - neighbor_psi); + vector > > tmp (n_q_points, + vector >(n_components)); fe_face_values.get_function_grads (solution, tmp); - neighbor_psi.swap (tmp[selected_component]); + for (unsigned int i=0; i