<ol>
<li> New: Added two optional parameters to TriaAccessor::center()
- and a new method TriaAccessor::point(). They allow to query for a
+ and a new method TriaAccessor::intermediate_point().
+ They allow to query for a
geometrically coherent center, or ask for arbitrary points on the
underlying Manifold, given the dim coordinates in the reference
element.
/**
* Returns a point belonging to the Manifold<dim,spacedim> where
* this object lives, given its parametric coordinates on the
- * reference #strucdim cell. This function queries the underlying
+ * reference #structdim cell. This function queries the underlying
* manifold object, and can be used to obtain the exact geometrical
* location of arbitrary points on this object.
+ *
+ * Notice that the argument @p coordinates are the coordinates on
+ * the <emph>reference cell</emph>, given in reference
+ * coordinates. In other words, the argument provides a weighting
+ * between the different vertices. For example, for lines, calling
+ * this function with argument Point<1>(.5), is equivalent to asking
+ * the line for its center.
*/
- Point<spacedim> point(const Point<structdim> &coordinates) const;
+ Point<spacedim> intermediate_point(const Point<structdim> &coordinates) const;
/**
* Center of the object. The center of an object is defined to be
// triangulation.vertices[next_unused_vertex] = new_point;
triangulation.vertices[next_unused_vertex] =
cell->get_manifold().get_new_point
- (Manifolds::get_default_quadrature
- <typename Triangulation<dim,spacedim>::active_cell_iterator,
- spacedim>(cell));
+ (Manifolds::get_default_quadrature(cell));
// if the user_flag is set, i.e. if the
// cell is at the boundary, use a
// function
triangulation.vertices[next_unused_vertex] =
quad->get_manifold().get_new_point
- (Manifolds::get_default_quadrature
- <typename Triangulation<dim,spacedim>::quad_iterator,
- spacedim>(quad, true));
+ (Manifolds::get_default_quadrature(quad, true));
}
triangulation.vertices_used[next_unused_vertex] = true;
// now that we created the right point, make up
template <int structdim, int dim, int spacedim>
Point<spacedim>
-TriaAccessor<structdim, dim, spacedim>::point (const Point<structdim> &coordinates) const
+TriaAccessor<structdim, dim, spacedim>::intermediate_point (const Point<structdim> &coordinates) const
{
// We use an FE_Q<structdim>(1) to extract the "weights" of each
// vertex, used to get a point from the manifold.