push_forward_gradient(const Point<spacedim> &chart_point) const override;
/**
- * Return the (normalized) normal vector at the point @p p.
+ * @copydoc Manifold::normal_vector()
*/
virtual Tensor<1, spacedim>
normal_vector(
const Point<spacedim> &x2) const override;
/**
- * Return the (normalized) normal vector at the point @p p.
+ * @copydoc Manifold::normal_vector()
*/
virtual Tensor<1, spacedim>
normal_vector(
// (tangential to the sphere). In this case, the normal vector is
// easy to compute since it is proportional to the vector from the
// center to the point 'p'.
- if (spherical_face_is_horizontal<dim,spacedim>(face, center))
+ if (spherical_face_is_horizontal<dim, spacedim>(face, center))
{
// So, if this is a "horizontal" face, then just compute the normal
// vector as the one from the center to the point 'p', adequately
// base class.
return Manifold<dim, spacedim>::normal_vector(face, p);
- return Tensor<1,spacedim>();
+ return Tensor<1, spacedim>();
}
// (tangential to the sphere). In this case, the normal vector is
// easy to compute since it is proportional to the vector from the
// center to the point 'p'.
- if (spherical_face_is_horizontal<dim,spacedim>(face, center))
+ if (spherical_face_is_horizontal<dim, spacedim>(face, center))
{
// So, if this is a "horizontal" face, then just compute the normal
// vector as the one from the center to the point 'p', adequately
// base class.
return Manifold<dim, spacedim>::normal_vector(face, p);
- return Tensor<1,spacedim>();
+ return Tensor<1, spacedim>();
}
// (tangential to the sphere). In this case, the normal vector is
// easy to compute since it is proportional to the vector from the
// center to the point 'p'.
- if (spherical_face_is_horizontal<dim,spacedim>(face,center))
+ if (spherical_face_is_horizontal<dim, spacedim>(face, center))
{
// So, if this is a "horizontal" face, then just compute the normal
// vector as the one from the center to the point 'p', adequately
// scaled.
- for (unsigned int vertex = 0; vertex < GeometryInfo<spacedim>::vertices_per_face; ++vertex)
+ for (unsigned int vertex = 0;
+ vertex < GeometryInfo<spacedim>::vertices_per_face;
+ ++vertex)
face_vertex_normals[vertex] = face->vertex(vertex) - center;
}
else