void check (Triangulation<2> &tria)
{
+ Point<2> p (1./3., 1./2.);
+
Triangulation<2>::active_cell_iterator cell
- = GridTools::find_active_cell_around_point (tria,
- Point<2>(1./3., 1./2.));
+ = GridTools::find_active_cell_around_point (tria, p);
deallog << cell << std::endl;
for (unsigned int v=0; v<GeometryInfo<2>::vertices_per_cell; ++v)
deallog << "<" << cell->vertex(v) << "> ";
deallog << std::endl;
+
+ Assert (p.distance (cell->center()) < cell->diameter()/2,
+ ExcInternalError());
}
void check (Triangulation<3> &tria)
{
+ Point<3> p(1./3.,1./2.,1./5.);
+
Triangulation<3>::active_cell_iterator cell
- = GridTools::find_active_cell_around_point (tria,
- Point<3>(1./3.,
- 1./2.,
- 2./3.));
+ = GridTools::find_active_cell_around_point (tria, p);
deallog << cell << std::endl;
for (unsigned int v=0; v<GeometryInfo<3>::vertices_per_cell; ++v)
deallog << "<" << cell->vertex(v) << "> ";
deallog << std::endl;
+
+ Assert (p.distance (cell->center()) < cell->diameter()/2,
+ ExcInternalError());
}
void check (Triangulation<3> &tria)
{
+ Point<3> p(1./3.,1./2.,-1./5.);
+
Triangulation<3>::active_cell_iterator cell
- = GridTools::find_active_cell_around_point (tria,
- Point<3>(1./3.,
- 1./2.,
- 2./3.));
+ = GridTools::find_active_cell_around_point (tria, p);
deallog << cell << std::endl;
for (unsigned int v=0; v<GeometryInfo<3>::vertices_per_cell; ++v)
deallog << "<" << cell->vertex(v) << "> ";
deallog << std::endl;
+
+ Assert (p.distance (cell->center()) < cell->diameter()/2,
+ ExcInternalError());
}
--- /dev/null
+
+DEAL::2.2
+DEAL::<0.250000 0.250000> <0.500000 0.250000> <0.500000 0.500000> <0.250000 0.500000>
+DEAL::2.75
+DEAL::<0.250000 0.518306> <0.198223 0.405600> <0.396447 0.396447> <0.500000 0.500000>
--- /dev/null
+
+DEAL::2.5
+DEAL::<0.250000 0.250000 0.00000> <0.500000 0.250000 0.00000> <0.500000 0.250000 0.250000> <0.250000 0.250000 0.250000> <0.250000 0.500000 0.00000> <0.500000 0.500000 0.00000> <0.500000 0.500000 0.250000> <0.250000 0.500000 0.250000>
+DEAL::2.445
+DEAL::<0.286905 0.658777 -1.73472e-18> <0.209333 0.437629 -3.46945e-18> <0.203251 0.422066 0.203251> <0.269808 0.604598 0.269808> <0.562887 0.562887 0.00000> <0.418667 0.418667 0.00000> <0.406502 0.406502 0.197169> <0.530301 0.530301 0.245590>
--- /dev/null
+
+DEAL::1.0
+DEAL::<0.00000 0.00000 0.00000> <0.500000 0.00000 0.00000> <0.500000 0.500000 0.00000> <0.00000 0.500000 0.00000> <0.00000 0.00000 -0.500000> <0.500000 0.00000 -0.500000> <0.500000 0.500000 -0.500000> <0.00000 0.500000 -0.500000>
+DEAL::1.16
+DEAL::<0.00000 0.00000 0.00000> <0.500000 0.00000 0.00000> <0.500000 0.500000 0.00000> <0.00000 0.500000 0.00000> <0.00000 0.00000 -0.500000> <0.500000 0.00000 -0.500000> <0.500000 0.500000 -0.500000> <0.00000 0.500000 -0.500000>
+DEAL::1.52
+DEAL::<0.00000 0.577350 -0.577350> <0.00000 0.394338 -0.394338> <-6.93889e-18 0.394338 0.00000> <0.00000 0.577350 0.00000> <0.577350 0.577350 -0.577350> <0.394338 0.394338 -0.394338> <0.394338 0.394338 0.00000> <0.577350 0.577350 0.00000>
--- /dev/null
+
+DEAL::2.2
+DEAL::<0.250000 0.250000> <0.500000 0.250000> <0.500000 0.500000> <0.250000 0.500000>
+DEAL::2.75
+DEAL::<0.250000 0.518306> <0.198223 0.405600> <0.396447 0.396447> <0.500000 0.500000>
--- /dev/null
+
+DEAL::2.5
+DEAL::<0.250000 0.250000 0.00000> <0.500000 0.250000 0.00000> <0.500000 0.250000 0.250000> <0.250000 0.250000 0.250000> <0.250000 0.500000 0.00000> <0.500000 0.500000 0.00000> <0.500000 0.500000 0.250000> <0.250000 0.500000 0.250000>
+DEAL::2.445
+DEAL::<0.286905 0.658777 -1.73472e-18> <0.209333 0.437629 -3.46945e-18> <0.203251 0.422066 0.203251> <0.269808 0.604598 0.269808> <0.562887 0.562887 0.00000> <0.418667 0.418667 0.00000> <0.406502 0.406502 0.197169> <0.530301 0.530301 0.245590>
--- /dev/null
+
+DEAL::1.0
+DEAL::<0.00000 0.00000 0.00000> <0.500000 0.00000 0.00000> <0.500000 0.500000 0.00000> <0.00000 0.500000 0.00000> <0.00000 0.00000 -0.500000> <0.500000 0.00000 -0.500000> <0.500000 0.500000 -0.500000> <0.00000 0.500000 -0.500000>
+DEAL::1.16
+DEAL::<0.00000 0.00000 0.00000> <0.500000 0.00000 0.00000> <0.500000 0.500000 0.00000> <0.00000 0.500000 0.00000> <0.00000 0.00000 -0.500000> <0.500000 0.00000 -0.500000> <0.500000 0.500000 -0.500000> <0.00000 0.500000 -0.500000>
+DEAL::1.52
+DEAL::<0.00000 0.577350 -0.577350> <0.00000 0.394338 -0.394338> <-6.93889e-18 0.394338 0.00000> <0.00000 0.577350 0.00000> <0.577350 0.577350 -0.577350> <0.394338 0.394338 -0.394338> <0.394338 0.394338 0.00000> <0.577350 0.577350 0.00000>
--- /dev/null
+
+DEAL::2.2
+DEAL::<0.250000 0.250000> <0.500000 0.250000> <0.500000 0.500000> <0.250000 0.500000>
+DEAL::2.75
+DEAL::<0.250000 0.518306> <0.198223 0.405600> <0.396447 0.396447> <0.500000 0.500000>
--- /dev/null
+
+DEAL::2.5
+DEAL::<0.250000 0.250000 0.00000> <0.500000 0.250000 0.00000> <0.500000 0.250000 0.250000> <0.250000 0.250000 0.250000> <0.250000 0.500000 0.00000> <0.500000 0.500000 0.00000> <0.500000 0.500000 0.250000> <0.250000 0.500000 0.250000>
+DEAL::2.445
+DEAL::<0.286905 0.658777 -1.73472e-18> <0.209333 0.437629 -3.46945e-18> <0.203251 0.422066 0.203251> <0.269808 0.604598 0.269808> <0.562887 0.562887 0.00000> <0.418667 0.418667 0.00000> <0.406502 0.406502 0.197169> <0.530301 0.530301 0.245590>
--- /dev/null
+
+DEAL::1.0
+DEAL::<0.00000 0.00000 0.00000> <0.500000 0.00000 0.00000> <0.500000 0.500000 0.00000> <0.00000 0.500000 0.00000> <0.00000 0.00000 -0.500000> <0.500000 0.00000 -0.500000> <0.500000 0.500000 -0.500000> <0.00000 0.500000 -0.500000>
+DEAL::1.16
+DEAL::<0.00000 0.00000 0.00000> <0.500000 0.00000 0.00000> <0.500000 0.500000 0.00000> <0.00000 0.500000 0.00000> <0.00000 0.00000 -0.500000> <0.500000 0.00000 -0.500000> <0.500000 0.500000 -0.500000> <0.00000 0.500000 -0.500000>
+DEAL::1.52
+DEAL::<0.00000 0.577350 -0.577350> <0.00000 0.394338 -0.394338> <-6.93889e-18 0.394338 0.00000> <0.00000 0.577350 0.00000> <0.577350 0.577350 -0.577350> <0.394338 0.394338 -0.394338> <0.394338 0.394338 0.00000> <0.577350 0.577350 0.00000>