ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-01-27 13:20:34 +05:00
Code Issues Packages Projects Releases Wiki Activity

Merge branch 'elofpoint-tol' into 'master'

Browse Source 
Use tolerance in GetElementOfPoint when using a search tree

I'm getting artifacts at domain edges (see attached), because no tolerance is used when calling `elementsearchtree->GetIntersecting()` inside `GetElementOfPoint`.

![searchtree](/uploads/5aa10da74cac7bef9e134756ea897e49/searchtree.png)

See merge request !29
This commit is contained in:
Joachim Schöberl 2017-01-25 21:25:52 +01:00
commit 52517218c9
1 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
libsrc/meshing/meshclass.cpp
View File

@ -4720,6 +4720,9 @@ namespace netgen
bool build_searchtree, bool build_searchtree,
const bool allowindex) const const bool allowindex) const
{ {
const double pointtol = 1e-12;
netgen::Point<3> pmin = p - Vec<3> (pointtol, pointtol, pointtol);
netgen::Point<3> pmax = p + Vec<3> (pointtol, pointtol, pointtol);
if (dimension == 2) if (dimension == 2)
{ {
int ne; int ne;
@ -4733,7 +4736,7 @@ namespace netgen
{ {
// update if necessary: // update if necessary:
const_cast<Mesh&>(*this).BuildElementSearchTree (); const_cast<Mesh&>(*this).BuildElementSearchTree ();
elementsearchtree->GetIntersecting (p, p, locels); elementsearchtree->GetIntersecting (pmin, pmax, locels);
ne = locels.Size(); ne = locels.Size();
} }
else else
@ -4776,7 +4779,7 @@ namespace netgen
{ {
// update if necessary: // update if necessary:
const_cast<Mesh&>(*this).BuildElementSearchTree (); const_cast<Mesh&>(*this).BuildElementSearchTree ();
elementsearchtree->GetIntersecting (p, p, locels); elementsearchtree->GetIntersecting (pmin, pmax, locels);
ne = locels.Size(); ne = locels.Size();
} }
else else