Commit 78e75ece authored by Michael Turner's avatar Michael Turner
Browse files

dave comment

parent c83d04c9
......@@ -104,28 +104,15 @@ NekDouble CADCurve::loct(Array<OneD, NekDouble> xyz)
Array<OneD, NekDouble> b = Bounds();
gp_Pnt loc(xyz[0]*1000.0, xyz[1]*1000.0, xyz[2]*1000.0);
/*GeomAPI_ProjectPointOnCurve projection(
loc,m_c,b[0],b[1]);
if (projection.NbPoints() == 0)
{
ASSERTL0(false,"failed");
}
else
{
t = projection.Parameter(1);
if(projection.Distance(1) > 1e-6)
{
cout << "large curve projection: " << projection.Distance(1) << endl;
}
}*/
ShapeAnalysis_Curve sac;
gp_Pnt p;
NekDouble d = sac.Project(m_c,loc,1e-7,p,t);
ASSERTL0(p.Distance(loc) < 1e-6, "large loct distance sac");
NekDouble d = sac.Project(m_c,loc,1e-8 ,p,t);
if(p.Distance(loc) > 1e-5)
{
cerr << "large loct distance" << endl;
}
return t;
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment