1 // Created on: 1997-11-06
2 // Created by: Roman BORISOV
3 // Copyright (c) 1997-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
6 // This file is part of Open CASCADE Technology software library.
8 // This library is free software; you can redistribute it and/or modify it under
9 // the terms of the GNU Lesser General Public License version 2.1 as published
10 // by the Free Software Foundation, with special exception defined in the file
11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
12 // distribution for complete text of the license and disclaimer of any warranty.
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
17 #include <ProjLib_PrjResolve.ixx>
18 #include <ProjLib_PrjFunc.hxx>
19 #include <math_FunctionSetRoot.hxx>
20 #include <math_NewtonFunctionSetRoot.hxx>
22 ProjLib_PrjResolve::ProjLib_PrjResolve(const Adaptor3d_Curve& C,const Adaptor3d_Surface& S,const Standard_Integer Fix) : myFix(Fix)
24 if (myFix > 3 || myFix < 1) Standard_ConstructionError::Raise();
25 mySolution = gp_Pnt2d(0.,0.);
26 myCurve = (Adaptor3d_CurvePtr)&C;
27 mySurface = (Adaptor3d_SurfacePtr)&S;
30 // void ProjLib_PrjResolve::Perform(const Standard_Real t, const Standard_Real U, const Standard_Real V, const gp_Pnt2d& Tol2d, const gp_Pnt2d& Inf, const gp_Pnt2d& Sup, const Standard_Real FuncTol, const Standard_Boolean StrictInside)
31 void ProjLib_PrjResolve::Perform(const Standard_Real t, const Standard_Real U, const Standard_Real V, const gp_Pnt2d& Tol2d, const gp_Pnt2d& Inf, const gp_Pnt2d& Sup, const Standard_Real FuncTol, const Standard_Boolean )
34 myDone = Standard_False;
35 Standard_Real FixVal = 0.;
36 gp_Pnt2d ExtInf(0.,0.), ExtSup(0.,0.);
37 Standard_Real ExtU = 10*Tol2d.X(), ExtV = 10*Tol2d.Y();
38 math_Vector Tol(1, 2), Start(1, 2), BInf(1, 2), BSup(1, 2);
40 ExtInf.SetCoord(Inf.X() - ExtU, Inf.Y() - ExtV);
41 ExtSup.SetCoord(Sup.X() + ExtU, Sup.Y() + ExtV);
66 ProjLib_PrjFunc F(myCurve, FixVal, mySurface, myFix);
69 // Standard_Integer option = 1;//2;
71 // math_FunctionSetRoot S1 (F, Start,Tol, BInf, BSup);
72 // if (!S1.IsDone()) { return; }
75 math_NewtonFunctionSetRoot SR (F, Tol, 1.e-10);
76 SR.Perform(F, Start, BInf, BSup);
77 // if (!SR.IsDone()) { return; }
79 math_FunctionSetRoot S1 (F, Start,Tol, BInf, BSup);
80 if (!S1.IsDone()) { return; }
85 mySolution.SetXY(F.Solution().XY());
87 // computation of myDone
88 myDone = Standard_True;
90 Standard_Real ExtraU , ExtraV;
91 // if(!StrictInside) {
95 if (mySolution.X() > Inf.X() - Tol2d.X() && mySolution.X() < Inf.X()) mySolution.SetX(Inf.X());
96 if (mySolution.X() > Sup.X() && mySolution.X() < Sup.X() + Tol2d.X()) mySolution.SetX(Sup.X());
97 if (mySolution.Y() > Inf.Y() - Tol2d.Y() && mySolution.Y() < Inf.Y()) mySolution.SetY(Inf.Y());
98 if (mySolution.Y() > Sup.Y() && mySolution.Y() < Sup.Y() + Tol2d.Y()) mySolution.SetY(Sup.Y());
99 if (mySolution.X() < Inf.X() - ExtraU ||
100 mySolution.X() > Sup.X() + ExtraU ||
101 mySolution.Y() < Inf.Y() - ExtraV ||
102 mySolution.Y() > Sup.Y() + ExtraV) myDone = Standard_False;
103 else if (FuncTol > 0) {
104 math_Vector X(1,2,0.), FVal(1,2,0.);
105 X(1) = mySolution.X();
106 X(2) = mySolution.Y();
111 if ((FVal(1)*FVal(1) + FVal(2)*FVal(2)) > FuncTol) myDone = Standard_False;
117 Standard_Boolean ProjLib_PrjResolve::IsDone() const
122 gp_Pnt2d ProjLib_PrjResolve::Solution() const
124 if (!IsDone()) StdFail_NotDone::Raise();