/*
- reduced failing test. Original exception in transform_real_to_unit_cell:
-
-#0 0x00007fffeec931a0 in __cxa_throw () from /usr/lib64/libstdc++.so.6
-#1 0x00007ffff5215d60 in dealii::deal_II_exceptions::internals::issue_error_throw<dealii::Mapping<2, 2>::ExcTransformationFailed> (
- file=0x7ffff6b47ce8 "/w/kwang/kwang/deal.II/source/fe/mapping_q1.cc", line=1803,
- function=0x7ffff6b54e80 "dealii::Point<dim, double> dealii::MappingQ1<dim, spacedim>::transform_real_to_unit_cell_internal(const typename dealii::Triangulation<dim, spacedim>::cell_iterator&, const dealii::Point<spacedim>&, c"..., cond=0x7ffff6b47b3d "false",
- exc_name=0x7ffff6b47fd8 "(typename Mapping<dim,spacedim>::ExcTransformationFailed())", e=...)
- at /w/kwang/kwang/deal.II/include/deal.II/base/exceptions.h:304
-#2 0x00007ffff5201d5f in dealii::MappingQ1<2, 2>::transform_real_to_unit_cell_internal (this=0x7ffff7dd8160, cell={
- triangulation = 0x7fffffffd2b0,
- level = 1,
- index = 4
-}, p={-0.27999999999999992, -0.62999999999999989}, initial_p_unit={1.4871340295728577, -0.21318360089489402}, mdata=...)
- at /w/kwang/kwang/deal.II/source/fe/mapping_q1.cc:1803
-#3 0x00007ffff51ff24a in dealii::MappingQ1<2, 2>::transform_real_to_unit_cell (this=0x7ffff7dd8160, cell={
- triangulation = 0x7fffffffd2b0,
- level = 1,
- index = 4
-}, p={-0.27999999999999992, -0.62999999999999989}) at /w/kwang/kwang/deal.II/source/fe/mapping_q1.cc:1683
-#4 0x00007ffff4037ea2 in dealii::Functions::FEFieldFunction<2, dealii::DoFHandler<2, 2>, dealii::Vector<double> >::compute_point_locations (this=0x7fffffffcc80, points=..., cells=..., qpoints=..., maps=...)
- at /w/kwang/kwang/deal.II/include/deal.II/numerics/fe_field_function.templates.h:443
-#5 0x00007ffff4035808 in dealii::Functions::FEFieldFunction<2, dealii::DoFHandler<2, 2>, dealii::Vector<double> >::vector_value_list (this=0x7fffffffcc80, points=..., values=...) at /w/kwang/kwang/deal.II/include/deal.II/numerics/fe_field_function.templates.h:208
-#6 0x00007ffff40355e5 in dealii::Functions::FEFieldFunction<2, dealii::DoFHandler<2, 2>, dealii::Vector<double> >::value_list (
- this=0x7fffffffcc80, points=..., values=..., component=0)
- at /w/kwang/kwang/deal.II/include/deal.II/numerics/fe_field_function.templates.h:246
-#7 0x0000000000450fd5 in Step6<2>::compute_measure (this=0x7fffffffd2b0) at step-6.cc:627
-#8 0x000000000045874c in Step6<2>::run (this=0x7fffffffd2b0) at step-6.cc:705
-#9 0x000000000045127f in main () at step-6.cc:728
-
-
- Note that test() only looks at the one cell that causes the failure.
-
-test2() is the original code using FEFieldFunction. Note that it depends on the contents of the points in the value_list (see comments).
-
- */
+ FEFieldFunction could not deal with ExcTransformationFailed from transform_real_to_unit_cell()
+ */
#include "../tests.h"
#include <deal.II/numerics/fe_field_function.h>
-template<int dim>
-void test()
-{
- const HyperBallBoundary<dim> boundary_description;
-
- Triangulation<dim> triangulation;
- GridGenerator::hyper_ball (triangulation);
- triangulation.set_boundary (0, boundary_description);
- triangulation.refine_global (1);
- MappingQ1< dim > mapping;
-
- Point<dim> p(-0.27999999999999992, -0.62999999999999989);
- typename Triangulation<dim>::active_cell_iterator cell;
-
- for (cell = triangulation.begin_active();
- cell != triangulation.end(); cell++)
- {
- if (cell->level()==1 && cell->index()==4)
- mapping.transform_real_to_unit_cell(cell, p);
-
-
- }
-
- triangulation.set_boundary (0);
-
- deallog << "OK" << std::endl;
-}
template<int dim>
deallog.depth_console(0);
deallog.threshold_double(1.e-10);
- test<2>();
test2<2>();
return 0;
-DEAL::dim=1, spacedim=1
-DEAL::Vertex displacement: 0.00000
-DEAL::Vertex displacement: 0.0800000
-DEAL::Vertex displacement: 0.160000
-DEAL::Vertex displacement: 0.240000
-DEAL::Vertex displacement: 0.320000
-DEAL::OK
-DEAL::dim=2, spacedim=2
-DEAL::Vertex displacement: 0.00000 0.00000
-DEAL::Vertex displacement: 0.0800000 0.0800000
-DEAL::Vertex displacement: 0.160000 0.160000
-DEAL::Vertex displacement: 0.240000 0.240000
-DEAL::Vertex displacement: 0.320000 0.320000
-DEAL::OK
-DEAL::dim=3, spacedim=3
-DEAL::Vertex displacement: 0.00000 0.00000 0.00000
-DEAL::Vertex displacement: 0.0800000 0.0800000 0.0800000
-DEAL::Vertex displacement: 0.160000 0.160000 0.160000
-DEAL::Vertex displacement: 0.240000 0.240000 0.240000
-DEAL::Vertex displacement: 0.320000 0.320000 0.320000
-DEAL::OK
-DEAL::dim=1, spacedim=2
-DEAL::Vertex displacement: 0.00000 0.00000
-DEAL::Vertex displacement: 0.0800000 0.0800000
-DEAL::Vertex displacement: 0.160000 0.160000
-DEAL::Vertex displacement: 0.240000 0.240000
-DEAL::Vertex displacement: 0.320000 0.320000
-DEAL::OK
-DEAL::dim=1, spacedim=3
-DEAL::Vertex displacement: 0.00000 0.00000 0.00000
-DEAL::Vertex displacement: 0.0800000 0.0800000 0.0800000
-DEAL::Vertex displacement: 0.160000 0.160000 0.160000
-DEAL::Vertex displacement: 0.240000 0.240000 0.240000
-DEAL::Vertex displacement: 0.320000 0.320000 0.320000
-DEAL::OK
-DEAL::dim=2, spacedim=3
-DEAL::Vertex displacement: 0.00000 0.00000 0.00000
-DEAL::Vertex displacement: 0.0800000 0.0800000 0.0800000
-DEAL::Vertex displacement: 0.160000 0.160000 0.160000
-DEAL::Vertex displacement: 0.240000 0.240000 0.240000
-DEAL::Vertex displacement: 0.320000 0.320000 0.320000
DEAL::OK