switch (dim)
{
case 2:
- {
- for (; result != end; ++result, ++tang1)
- cross_product (*result, *tang1);
- break;
- };
+ {
+ for (; result != end; ++result, ++tang1)
+ cross_product (*result, *tang1);
+ break;
+ };
case 3:
- {
- Assert (face_no+GeometryInfo<dim>::faces_per_cell <
- data.unit_tangentials.n_rows(),
- ExcInternalError());
- Assert (data.aux[1].size() <=
- data.unit_tangentials[face_no+GeometryInfo<dim>::faces_per_cell].size(),
- ExcInternalError());
- transform_contravariant(data.aux[1].begin(),
- data.aux[1].end(),
- data.unit_tangentials[
- face_no+GeometryInfo<dim>::faces_per_cell].begin(),
- data);
- const Tensor<1,dim> *tang2 = data.aux[1].begin();
- for (;result != end; ++result, ++tang1, ++tang2)
- cross_product (*result, *tang1, *tang2);
- break;
- };
+ {
+ Assert (face_no+GeometryInfo<dim>::faces_per_cell <
+ data.unit_tangentials.n_rows(),
+ ExcInternalError());
+ Assert (data.aux[1].size() <=
+ data.unit_tangentials[face_no+GeometryInfo<dim>::faces_per_cell].size(),
+ ExcInternalError());
+ transform_contravariant(data.aux[1].begin(),
+ data.aux[1].end(),
+ data.unit_tangentials[
+ face_no+GeometryInfo<dim>::faces_per_cell].begin(),
+ data);
+ const Tensor<1,dim> *tang2 = data.aux[1].begin();
+ for (;result != end; ++result, ++tang1, ++tang2)
+ cross_product (*result, *tang1, *tang2);
+ break;
+ };
default:
- Assert(false, ExcNotImplemented());
+ Assert(false, ExcNotImplemented());
}
if (update_flags & (update_normal_vectors