From 889685c41e75d22000cea97bf3e7c7ac57026908 Mon Sep 17 00:00:00 2001 From: guido Date: Sun, 15 May 2005 00:57:04 +0000 Subject: [PATCH] order of indices fixed git-svn-id: https://svn.dealii.org/trunk@10679 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/fe/fe_values.cc | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/deal.II/deal.II/source/fe/fe_values.cc b/deal.II/deal.II/source/fe/fe_values.cc index 1eda80854d..38b9403078 100644 --- a/deal.II/deal.II/source/fe/fe_values.cc +++ b/deal.II/deal.II/source/fe/fe_values.cc @@ -736,7 +736,7 @@ FEValuesBase:: get_function_grads (const InputVector &fe_function, std::vector > > &gradients) const { - Assert (n_quadrature_points == gradients.size(), + Assert (gradients.size() == n_quadrature_points, ExcDimensionMismatch(gradients.size(), n_quadrature_points)); const unsigned int n_components = fe->n_components(); @@ -769,7 +769,7 @@ get_function_grads (const InputVector &fe_function, { Tensor<1,dim> tmp = this->shape_grad(shape_func,point); tmp *= dof_values(shape_func); - gradients[fe->system_to_component_index(shape_func).first][point] + gradients[point][fe->system_to_component_index(shape_func).first] += tmp; } else @@ -777,7 +777,7 @@ get_function_grads (const InputVector &fe_function, { Tensor<1,dim> tmp = this->shape_grad_component(shape_func,point,c); tmp *= dof_values(shape_func); - gradients[c][point] += tmp; + gradients[point][c] += tmp; } } -- 2.39.5