#1658 HierarchicSearch suffers from round-off errors
Metadata
Property | Value |
---|---|
Reported by | Anton Schiela (schiela@zib.de) |
Reported at | May 29, 2015 08:29 |
Type | Bug Report |
Version | 2.3 |
Operating System | Unspecified / All |
Description
HierarchicSearch occasionally throws an exception, because a coordinate is found in the father element, but in none of the children:
cf. hierarchicsearch.hh:
DUNE_THROW(Exception,"{" << className(*this) << "} Unexpected "
"internal Error: none of the children of the entity "
"{" << formatEntityInformation(entity) << "} contains "
"coordinate (" << global << "). Children are: "
"[" << children.str() << "].");
This usually occurs, if the coordinate is very close to a vertex of an element, so that at least
one of the local coordinates is very close to zero (approx. +-1e-14), and round-off errors couse a "wrong" result
of checkInside
, i.e. checkInside
is true
for the father, while it is false
for all children.