sp.resize(2);
wp.resize(2);
sp[0] = obj->vertex(0); wp[0] = .5;
- sp[1] = obj->vertex(1); wp[0] = .5;
+ sp[1] = obj->vertex(1); wp[1] = .5;
break;
case 2:
sp.resize(8);
// connection between the new
// points on this face and on the
// opposite face
- std::vector<Point<spacedim> > sp(2);
+ std::vector<Point<spacedim> > ps(2);
std::vector<double> ws(2, .5);
- sp[0] = cell->face(boundary_face)
+ ps[0] = cell->face(boundary_face)
->child(0)->vertex(1);
- sp[1] = cell->face(GeometryInfo<dim>
+ ps[1] = cell->face(GeometryInfo<dim>
::opposite_face[boundary_face])
->child(0)->vertex(1);
const Manifold<spacedim> &
Triangulation<dim, spacedim>::get_manifold (const types::manifold_id m_number) const
{
- Assert(m_number < numbers::invalid_manifold_id,
- ExcIndexRange(m_number,0,numbers::invalid_manifold_id));
-
//look, if there is a boundary stored at
//boundary_id number.
typename std::map<types::manifold_id, SmartPointer<const Manifold<spacedim>, Triangulation<dim, spacedim> > >::const_iterator it