ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-11-11 16:49:16 +05:00
Code Issues Packages Projects Releases Wiki Activity

pointsearch in 2D

Browse Source
This commit is contained in:
Joachim Schoeberl 2011-08-11 16:13:35 +00:00
parent 5733e491fe
commit 1312e1d1da
1 changed files with 47 additions and 23 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

70
libsrc/meshing/meshclass.cpp
View File

@ -4079,40 +4079,62 @@ namespace netgen
PrintMessage (4, "Rebuild element searchtree");
if (elementsearchtree)
delete elementsearchtree;
delete elementsearchtree;
elementsearchtree = NULL;
Box3d box;
int i, j;
int ne = GetNE();
int ne = (dimension == 2) ? GetNSE() : GetNE();
if (!ne)
{
lock.UnLock();
return;
}
box.SetPoint (Point (VolumeElement(1).PNum(1)));
for (i = 1; i <= ne; i++)
if (dimension == 2)
{
const Element & el = VolumeElement(i);
for (j = 1; j <= el.GetNP(); j++)
box.AddPoint (Point (el.PNum(j)));
box.SetPoint (Point (SurfaceElement(1).PNum(1)));
for (int i = 1; i <= ne; i++)
{
const Element2d & el = SurfaceElement(i);
for (int j = 1; j <= el.GetNP(); j++)
box.AddPoint (Point (el.PNum(j)));
}
box.Increase (1.01 * box.CalcDiam());
elementsearchtree = new Box3dTree (box.PMin(), box.PMax());
for (int i = 1; i <= ne; i++)
{
const Element2d & el = SurfaceElement(i);
box.SetPoint (Point (el.PNum(1)));
for (int j = 1; j <= el.GetNP(); j++)
box.AddPoint (Point (el.PNum(j)));
elementsearchtree -> Insert (box.PMin(), box.PMax(), i);
}
}
box.Increase (1.01 * box.CalcDiam());
elementsearchtree = new Box3dTree (box.PMin(), box.PMax());
for (i = 1; i <= ne; i++)
else
{
const Element & el = VolumeElement(i);
box.SetPoint (Point (el.PNum(1)));
for (j = 1; j <= el.GetNP(); j++)
box.AddPoint (Point (el.PNum(j)));
box.SetPoint (Point (VolumeElement(1).PNum(1)));
for (int i = 1; i <= ne; i++)
{
const Element & el = VolumeElement(i);
for (int j = 1; j <= el.GetNP(); j++)
box.AddPoint (Point (el.PNum(j)));
}
box.Increase (1.01 * box.CalcDiam());
elementsearchtree = new Box3dTree (box.PMin(), box.PMax());
for (int i = 1; i <= ne; i++)
{
const Element & el = VolumeElement(i);
box.SetPoint (Point (el.PNum(1)));
for (int j = 1; j <= el.GetNP(); j++)
box.AddPoint (Point (el.PNum(j)));
elementsearchtree -> Insert (box.PMin(), box.PMax(), i);
elementsearchtree -> Insert (box.PMin(), box.PMax(), i);
}
}
elementsearchtreets = GetTimeStamp();
@ -4501,8 +4523,10 @@ namespace netgen
return ps_startelement;
Array<int> locels;
if (0)
if (elementsearchtree || build_searchtree)
{
// update if necessary:
const_cast<Mesh&>(*this).BuildElementSearchTree ();
elementsearchtree->GetIntersecting (p, p, locels);
ne = locels.Size();
}
@ -4513,7 +4537,7 @@ namespace netgen
{
int ii;
if (0)
if (elementsearchtree)
ii = locels.Get(i);
else
ii = i;