#include <deal.II/base/config.h>
#include <deal.II/base/subscriptor.h>
#include <deal.II/base/point.h>
-#include <deal.II/grid/tria.h>
DEAL_II_NAMESPACE_OPEN
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.templates.h>
+#include <deal.II/grid/manifold.h>
+#include <deal.II/grid/tria_boundary.h>
#include <deal.II/distributed/tria.h>
#include <cmath>
// otherwise we will run into troubles...
if( (structdim < dim) ||
(dim < spacedim) )
- return this->tria->get_manifold
- (numbers::invalid_manifold_id);
+ return get_boundary();
+// return this->tria->get_manifold
+// (numbers::invalid_manifold_id);
/* *dynamic_cast<const Manifold<spacedim> *> */
/* (&this->tria->get_boundary(mi)); */
else
#include <deal.II/grid/tria_iterator.templates.h>
#include <deal.II/base/geometry_info.h>
#include <deal.II/grid/grid_tools.h>
+#include <deal.II/grid/manifold.h>
+#include <deal.II/grid/tria_boundary.h>
#include <deal.II/fe/mapping_q1.h>
#include <cmath>
return Point<spacedim>();
}
+template <int dim, int spacedim>
+Point<spacedim>
+Boundary<dim, spacedim>::
+get_new_point (const std::vector<Point<spacedim> > &,
+ const std::vector<double> &) const
+{
+ Assert (false, ExcPureFunctionCalled());
+ return Point<spacedim>();
+}
template <int dim, int spacedim>