if(myU-myUinfium < 2*aDelta)
{
F1=F;
- const Standard_Real U1 = myU, U2 = myU+aDelta, U3 = myU+2*aDelta;
+ //const Standard_Real U1 = myU;
+ const Standard_Real U2 = myU + aDelta;
+ const Standard_Real U3 = myU + aDelta * 2.0;
if(!((Value(U2,F2)) && (Value(U3,F3))))
{
else
{
F3 = F;
- const Standard_Real U1 = myU-2*aDelta, U2 = myU-aDelta, U3 = myU;
-
+ const Standard_Real U1 = myU - aDelta * 2.0;
+ const Standard_Real U2 = myU - aDelta;
+ //const Standard_Real U3 = myU;
+
if(!((Value(U2,F2)) && (Value(U1,F1))))
{
#ifdef DEB
n = TheCurveTool::NbIntervals(*((TheCurve*)myC), GeomAbs_C2);
TColStd_Array1OfReal theInter(1, n+1);
Standard_Boolean isPeriodic = TheCurveTool::IsPeriodic(*((TheCurve*)myC));
- Standard_Real aPeriodicShift = 0.;
TheCurveTool::Intervals(*((TheCurve*)myC), theInter, GeomAbs_C2);
mysample = Max(mysample/n, 17);
TheVector V1;
aCache2->CalculatePoints();
// Calcul des distances
- const Curve1& aCurve1 = *((Curve1*)(myF.CurvePtr (1)));
- const Curve2& aCurve2 = *((Curve1*)(myF.CurvePtr (2)));
const Handle(ArrayOfPnt)& aPntArray1 = aCache1->Points();
const Handle(ArrayOfPnt)& aPntArray2 = aCache2->Points();
Standard_Integer NoU, NoV;
Tol(2) = mytolv;
math_Vector UV(1, 2);
-
theParams.Parameter(UV(1), UV(2));
math_Vector UVinf(1,2), UVsup(1,2);
UVsup(1) = myusup;
UVsup(2) = myvsup;
- math_Vector errors(1,2);
- math_Vector root(1, 2);
-
- Standard_Integer aNbMaxIter = 100;
+ const Standard_Integer aNbMaxIter = 100;
+ math_FunctionSetRoot S (myF, UV, Tol, UVinf, UVsup, aNbMaxIter);
- gp_Pnt PStart = theParams.Value();
-
- math_FunctionSetRoot S (myF,UV,Tol,UVinf,UVsup, aNbMaxIter);
-
myDone = Standard_True;
}