for (unsigned int k=0; k<this->dofs_per_cell; ++k)
{
- for (unsigned int i=0;i<quadrature.n_quadrature_points;++i)
- if (flags & update_values)
+ if (flags & update_values)
+ for (unsigned int i=0; i<quadrature.n_quadrature_points; ++i)
data.shape_values(k,i) = fe_data.shape_values[k][i+offset];
if (flags & update_gradients)
// offset determines which data set
// to take (all data sets for all
// sub-faces are stored contiguously)
- const unsigned int offset = (face * GeometryInfo<dim>::subfaces_per_face + subface)
- * quadrature.n_quadrature_points;
+ const unsigned int offset = ((face * GeometryInfo<dim>::subfaces_per_face + subface)
+ * quadrature.n_quadrature_points);
const UpdateFlags flags(fe_data.update_once | fe_data.update_each);
for (unsigned int k=0; k<this->dofs_per_cell; ++k)
{
- for (unsigned int i=0;i<quadrature.n_quadrature_points;++i)
- if (flags & update_values)
+ if (flags & update_values)
+ for (unsigned int i=0; i<quadrature.n_quadrature_points; ++i)
data.shape_values(k,i) = fe_data.shape_values[k][i+offset];
if (flags & update_gradients)