From b74bb274efbe50ec6de0e72ec5050a777540e8f6 Mon Sep 17 00:00:00 2001 From: guido Date: Thu, 12 May 2005 09:16:31 +0000 Subject: [PATCH] handle quadrature_points_fastest correctly git-svn-id: https://svn.dealii.org/trunk@10666 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/fe/fe_values.cc | 20 ++++++++++---------- 1 file changed, 10 insertions(+), 10 deletions(-) diff --git a/deal.II/deal.II/source/fe/fe_values.cc b/deal.II/deal.II/source/fe/fe_values.cc index 477e360bc3..1eda80854d 100644 --- a/deal.II/deal.II/source/fe/fe_values.cc +++ b/deal.II/deal.II/source/fe/fe_values.cc @@ -849,26 +849,26 @@ void FEValuesBase::get_function_grads ( for (unsigned int point=0; pointis_primitive(shape_func)) - values[point][fe->system_to_component_index(shape_func).first - +mc * n_components] + values[fe->system_to_component_index(shape_func).first + +mc * n_components][point] += fe_function(indices[shape_func+mc*dofs_per_cell]) * shape_grad(shape_func, point); else for (unsigned int c=0; cis_primitive(shape_func)) - values[fe->system_to_component_index(shape_func).first - +mc * n_components][point] + values[point][fe->system_to_component_index(shape_func).first + +mc * n_components] += fe_function(indices[shape_func+mc*dofs_per_cell]) * shape_grad(shape_func, point); else for (unsigned int c=0; c tmp(shape_2nd_derivative(shape_func,point)); tmp *= dof_values(shape_func); - second_derivs[point][fe->system_to_component_index(shape_func).first] + second_derivs[fe->system_to_component_index(shape_func).first][point] += tmp; } else @@ -960,7 +960,7 @@ get_function_2nd_derivatives (const InputVector &fe_func { Tensor<2,dim> tmp = this->shape_2nd_derivative_component(shape_func,point,c); tmp *= dof_values(shape_func); - second_derivs[point][c] += tmp; + second_derivs[c][point] += tmp; } else for (unsigned int point=0; point tmp(shape_2nd_derivative(shape_func,point)); tmp *= dof_values(shape_func); - second_derivs[fe->system_to_component_index(shape_func).first][point] + second_derivs[point][fe->system_to_component_index(shape_func).first] += tmp; } else @@ -977,7 +977,7 @@ get_function_2nd_derivatives (const InputVector &fe_func { Tensor<2,dim> tmp = this->shape_2nd_derivative_component(shape_func,point,c); tmp *= dof_values(shape_func); - second_derivs[c][point] += tmp; + second_derivs[point][c] += tmp; } } -- 2.39.5