ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-03-18 12:37:56 +05:00
Code Issues Packages Projects Releases Wiki Activity

remove debug output

Browse Source
This commit is contained in:
Christopher Lackner 2025-03-13 18:41:38 +01:00
parent 7aae5369c4
commit b8d722d6a8
1 changed files with 2 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
libsrc/meshing/meshclass.cpp
View File

@ -83,7 +83,6 @@ namespace netgen
bool allowindex) bool allowindex)
{ {
double vlam[3]; double vlam[3];
cout << "find2delement " << p << endl;
ElementIndex velement = ElementIndex::INVALID; ElementIndex velement = ElementIndex::INVALID;
if(mesh.GetNE()) if(mesh.GetNE())
@ -91,7 +90,6 @@ namespace netgen
if(searchtree) if(searchtree)
const_cast<Mesh&>(mesh).BuildElementSearchTree(3); const_cast<Mesh&>(mesh).BuildElementSearchTree(3);
velement = Find3dElement(mesh, p,vlam, nullopt,searchtree ? mesh.GetElementSearchTree() : nullptr,allowindex); velement = Find3dElement(mesh, p,vlam, nullopt,searchtree ? mesh.GetElementSearchTree() : nullptr,allowindex);
cout << "found volume element = " << velement << endl;
} }
//(*testout) << "p " << p << endl; //(*testout) << "p " << p << endl;
@ -106,8 +104,6 @@ namespace netgen
for(auto face : fnrs) for(auto face : fnrs)
faces.Append(topology.GetFace2SurfaceElement(face)); faces.Append(topology.GetFace2SurfaceElement(face));
cout << "faces = " << faces << endl;
for(int i=0; i<faces.Size(); i++) for(int i=0; i<faces.Size(); i++)
{ {
if(!faces[i].IsValid()) if(!faces[i].IsValid())
@ -117,13 +113,10 @@ namespace netgen
continue; continue;
auto & el = mesh[velement]; auto & el = mesh[velement];
cout << "eltype = " << el.GetType() << endl;
if (el.GetType() == TET) if (el.GetType() == TET)
{ {
double lam4[4] = { vlam[0], vlam[1], vlam[2], 1.0-vlam[0]-vlam[1]-vlam[2] }; double lam4[4] = { vlam[0], vlam[1], vlam[2], 1.0-vlam[0]-vlam[1]-vlam[2] };
cout << "lam4 = " << lam4[0] << ", " << lam4[1] << ", " << lam4[2] << ", " << lam4[3] << endl;
double face_lam = lam4[i]; double face_lam = lam4[i];
cout << "face lam = " << face_lam << endl;
if(face_lam < 1e-5) if(face_lam < 1e-5)
{ {
// found volume point very close to a face -> use barycentric coordinates directly // found volume point very close to a face -> use barycentric coordinates directly
@ -132,16 +125,12 @@ namespace netgen
for(auto k : Range(4)) for(auto k : Range(4))
if(sel[j] == el[k]) if(sel[j] == el[k])
lami[j-1] = lam4[k]/(1.0-face_lam); lami[j-1] = lam4[k]/(1.0-face_lam);
cout << "found close tet face = " << faces[i] << ", sel = " << sel << endl;
return SurfaceElementIndex(faces[i]); return SurfaceElementIndex(faces[i]);
} }
} }
if(mesh.PointContainedIn2DElement(p,lami,faces[i],true)) if(mesh.PointContainedIn2DElement(p,lami,faces[i],true))
{ return faces[i];
cout << "apoint contained in 2d el = " << faces[i] << ", sel = " << sel << endl;
return faces[i];
}
} }
} }
@ -172,11 +161,7 @@ namespace netgen
if((allowindex && !contained) || (!allowindex && contained)) continue; if((allowindex && !contained) || (!allowindex && contained)) continue;
} }
if(mesh.PointContainedIn2DElement(p,lami,ii)) if(mesh.PointContainedIn2DElement(p,lami,ii))
{ return ii;
cout << "point contained in 2d el = " << ii << ", sel = " << mesh[ii] << endl;
return ii;
}
} }
return 0; return 0;
} }