ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-12-26 13:50:33 +05:00
Code Issues Packages Projects Releases Wiki Activity

search tree fix

Browse Source
This commit is contained in:
Joachim Schoeberl 2012-10-18 14:18:36 +00:00
parent 1149c0f836
commit 0166fca8d4
1 changed files with 10 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
libsrc/meshing/meshing2.cpp
View File

@ -302,7 +302,8 @@ namespace netgen
box.Set ( mesh[sel[0]] ); box.Set ( mesh[sel[0]] );
box.Add ( mesh[sel[1]] ); box.Add ( mesh[sel[1]] );
box.Add ( mesh[sel[2]] ); box.Add ( mesh[sel[2]] );
surfeltree.Insert (box, i); // surfeltree.Insert (box, i);
surfeltree.Insert (box, seia[i]);
} }
@ -1078,16 +1079,13 @@ namespace netgen
{ {
const Point3d & p = critpoints.Get(i); const Point3d & p = critpoints.Get(i);
for (int jj = 0; jj < intersecttrias.Size(); jj++)
{
// int j = intersecttrias.Get(jj);
// const Element2d & el = mesh.SurfaceElement(j);
/* SurfaceElementIndex j = intersecttrias[jj];
for (j = 1; j <= mesh.GetNSE(); j++) const Element2d & el = mesh[j];
{
*/
int jj;
for (jj = 1; jj <= intersecttrias.Size(); jj++)
{
int j = intersecttrias.Get(jj);
const Element2d & el = mesh.SurfaceElement(j);
int ntrig = (el.GetNP() == 3) ? 1 : 2; int ntrig = (el.GetNP() == 3) ? 1 : 2;
@ -1306,7 +1304,7 @@ namespace netgen
box.Set (mesh[mtri[0]]); box.Set (mesh[mtri[0]]);
box.Add (mesh[mtri[1]]); box.Add (mesh[mtri[1]]);
box.Add (mesh[mtri[2]]); box.Add (mesh[mtri[2]]);
surfeltree.Insert (box, mesh.GetNSE()); surfeltree.Insert (box, mesh.GetNSE()-1);
const Point3d & sep1 = mesh.Point (mtri.PNum(1)); const Point3d & sep1 = mesh.Point (mtri.PNum(1));
const Point3d & sep2 = mesh.Point (mtri.PNum(2)); const Point3d & sep2 = mesh.Point (mtri.PNum(2));