mixed dim overlapping epsilon issue
We created the following modified mixeddimoverlapping test: mixeddimoverlappingtest.cc
A 2x2x2 unit cube grid is intersected with a 3-element one-dimensional grid on the z-axis. If I'm correct this should result in 4 intersections (and each intersection has four outside and four inside neighbors/embeddings because they are located on the edge). We get four intersections for the unit cube. Now, scaling the unit cube (and the one-dimensional grid) with epsilon produces the following results:
- eps = 1e-3: (only 3 intersections are found!?)
This is Codim0Extractor on a <3,3> grid!
This is Codim0Extractor on a <1,3> grid!
GridGlue: Constructor succeeded!
>>>> rank 0 coords: 27 and 4
>>>> rank 0 entities: 64 and 6
>>>> rank 0 types: 8 and 3
0 GridGlue::mergePatches : rank 0 / 0
StandardMerge building merged grid...
setup took 0.000171649 seconds.
intersection construction took 0.00623633 seconds.
0 GridGlue::mergePatches : The number of remote intersections is 3
Gluing successful, 3 remote intersections found!
- eps = 1e-8: (following error message)
This is Codim0Extractor on a <3,3> grid!
This is Codim0Extractor on a <1,3> grid!
GridGlue: Constructor succeeded!
>>>> rank 0 coords: 27 and 4
>>>> rank 0 entities: 64 and 6
>>>> rank 0 types: 8 and 3
0 GridGlue::mergePatches : rank 0 / 0
StandardMerge building merged grid...
setup took 0.000171437 seconds.
mixeddimoverlappingtest: /home-local/timok/monolithic/dune-grid-glue/dune/grid-glue/merging/overlappingmerge.cc:114: virtual void Dune::GridGlue::OverlappingMerge<3, 1, 3, double>::computeIntersections(const Dune::GeometryType &, const std::vector<Dune::FieldVector<T, dimworld> > &, std::bitset<(1 << dim1)> &, unsigned int, const Dune::GeometryType &, const std::vector<Dune::FieldVector<T, dimworld> > &, std::bitset<(1 << dim2)> &, unsigned int, std::vector<RemoteSimplicialIntersection> &) [dim1 = 3, dim2 = 1, dimworld = 3, T = double]: Assertion `dimis != 1' failed.
Aborted (core dumped)
- eps = 1e-10: (0 intersections found)
This is Codim0Extractor on a <3,3> grid!
This is Codim0Extractor on a <1,3> grid!
GridGlue: Constructor succeeded!
>>>> rank 0 coords: 27 and 4
>>>> rank 0 entities: 64 and 6
>>>> rank 0 types: 8 and 3
0 GridGlue::mergePatches : rank 0 / 0
StandardMerge building merged grid...
setup took 0.000168485 seconds.
intersection construction took 0.0311377 seconds.
0 GridGlue::mergePatches : The number of remote intersections is 0
Gluing successful, 0 remote intersections found!