Point<dim> trial_point;
for (unsigned int i=0; i<dim; ++i)
trial_point[i] = 1.5;
-
+
+ std::vector<Point<dim> > vs(2);
+
for (unsigned int e=0; e<GeometryInfo<dim>::lines_per_cell; ++e)
{
deallog << " Line " << e << ", projected point=";
- if (dim > 1)
- deallog << boundary.project_to_surface (cell->line(e), trial_point);
- else
- deallog << boundary.project_to_surface (cell, trial_point);
+ vs.resize(2);
+ if (dim > 1) {
+ vs[0] = cell->line(e)->vertex(0);
+ vs[1] = cell->line(e)->vertex(1);
+ deallog << boundary.project_to_manifold (vs, trial_point);
+ }
+ else {
+ vs[0] = cell->vertex(0);
+ vs[1] = cell->vertex(1);
+ deallog << boundary.project_to_manifold (vs, trial_point);
+ }
deallog << " (line is from ";
if (dim > 1)
Point<dim> trial_point;
for (unsigned int i=0; i<dim; ++i)
trial_point[i] = 1.5;
-
+
+ std::vector<Point<dim> > vs(GeometryInfo<dim>::vertices_per_cell);
+
for (unsigned int e=0; e<GeometryInfo<dim>::quads_per_cell; ++e)
{
const typename Triangulation<dim>::quad_iterator
*reinterpret_cast<const typename Triangulation<dim>::quad_iterator *>(&cell));
deallog << " Quad " << e << ", projected point=";
+ for(unsigned int i=0; i<GeometryInfo<dim>::vertices_per_cell; ++i)
+ vs[i] = quad->vertex(i);
- const Point<dim> p = boundary.project_to_surface (quad, trial_point);
+ const Point<dim> p = boundary.project_to_manifold (vs, trial_point);
deallog << p;
deallog << " (quad is from ";
deallog << quad->vertex(0);
deallog << "Trial point = " << trial_point << std::endl;
+ std::vector<Point<dim> > vs(GeometryInfo<dim>::vertices_per_cell);
+
for (unsigned int e=0; e<GeometryInfo<dim>::quads_per_cell; ++e)
{
const typename Triangulation<dim>::quad_iterator
*reinterpret_cast<const typename Triangulation<dim>::quad_iterator *>(&cell));
deallog << " Quad " << e << ", projected point=";
+
+ for(unsigned int i=0; i<GeometryInfo<dim>::vertices_per_cell; ++i)
+ vs[i] = quad->vertex(i);
- const Point<dim> p = boundary.project_to_surface (quad, trial_point);
+ const Point<dim> p = boundary.project_to_surface (vs, trial_point);
deallog << p;
deallog << " (quad is from ";
deallog << quad->vertex(0);