From: Wolfgang Bangerth <bangerth@math.tamu.edu> Date: Sun, 8 Nov 2015 01:49:53 +0000 (-0600) Subject: Make the same kind of changes also to FESystem. X-Git-Tag: v8.4.0-rc2~242^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F1845%2Fhead;p=dealii.git Make the same kind of changes also to FESystem. --- diff --git a/source/fe/fe_system.cc b/source/fe/fe_system.cc index 0dd710f662..51697cfda7 100644 --- a/source/fe/fe_system.cc +++ b/source/fe/fe_system.cc @@ -1197,29 +1197,31 @@ compute_fill (const Mapping<dim,spacedim> &mapping, Assert (this->n_nonzero_components(system_index) == base_fe.n_nonzero_components(base_index), ExcInternalError()); - for (unsigned int s=0; s<this->n_nonzero_components(system_index); ++s) - { - if (base_flags & update_values) - for (unsigned int q=0; q<n_q_points; ++q) - output_data.shape_values[out_index+s][q] = - base_data.shape_values(in_index+s,q); - - if (base_flags & update_gradients) - for (unsigned int q=0; q<n_q_points; ++q) - output_data.shape_gradients[out_index+s][q] = - base_data.shape_gradients[in_index+s][q]; - - if (base_flags & update_hessians) - for (unsigned int q=0; q<n_q_points; ++q) - output_data.shape_hessians[out_index+s][q] = - base_data.shape_hessians[in_index+s][q]; - - if (base_flags & update_3rd_derivatives) - for (unsigned int q=0; q<n_q_points; ++q) - output_data.shape_3rd_derivatives[out_index+s][q] = - base_data.shape_3rd_derivatives[in_index+s][q]; - } + if (base_flags & update_values) + for (unsigned int s=0; s<this->n_nonzero_components(system_index); ++s) + for (unsigned int q=0; q<n_q_points; ++q) + output_data.shape_values[out_index+s][q] = + base_data.shape_values(in_index+s,q); + + if (base_flags & update_gradients) + for (unsigned int s=0; s<this->n_nonzero_components(system_index); ++s) + for (unsigned int q=0; q<n_q_points; ++q) + output_data.shape_gradients[out_index+s][q] = + base_data.shape_gradients[in_index+s][q]; + + if (base_flags & update_hessians) + for (unsigned int s=0; s<this->n_nonzero_components(system_index); ++s) + for (unsigned int q=0; q<n_q_points; ++q) + output_data.shape_hessians[out_index+s][q] = + base_data.shape_hessians[in_index+s][q]; + + if (base_flags & update_3rd_derivatives) + for (unsigned int s=0; s<this->n_nonzero_components(system_index); ++s) + for (unsigned int q=0; q<n_q_points; ++q) + output_data.shape_3rd_derivatives[out_index+s][q] = + base_data.shape_3rd_derivatives[in_index+s][q]; + } } }