for (unsigned int point=0; point<n_q_points; ++point)
phi[point] = psi[point]*normal_vectors[point];
+ // if a coefficient was given: use that
+ // to scale the jump in the gradient
+ if (coefficient != 0)
+ {
+ vector<double> coefficient_values (n_q_points);
+ coefficient->value_list (fe_face_values_cell.get_quadrature_points(),
+ coefficient_values);
+ for (unsigned int point=0; point<n_q_points; ++point)
+ phi[point] *= coefficient_values[point];
+ };
+
if (face->at_boundary() == true)
// neumann boundary face. compute
phi[point] -= g[point];
};
- // if a coefficient was given: use that
- // to scale the jump in the gradient
- if (coefficient != 0)
- {
- vector<double> coefficient_values (n_q_points);
- coefficient->value_list (fe_face_values_cell.get_quadrature_points(),
- coefficient_values);
- for (unsigned int point=0; point<n_q_points; ++point)
- phi[point] *= coefficient_values[point];
- };
-
// now phi contains the following:
// - for an internal face, phi=[a du/dn]