// Copyright (c) 1995-1999 Matra Datavision // Copyright (c) 1999-2014 OPEN CASCADE SAS // // This file is part of Open CASCADE Technology software library. // // This library is free software; you can redistribute it and/or modify it under // the terms of the GNU Lesser General Public License version 2.1 as published // by the Free Software Foundation, with special exception defined in the file // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT // distribution for complete text of the license and disclaimer of any warranty. // // Alternatively, this file may be used under the terms of Open CASCADE // commercial license or contractual agreement. #ifndef OCCT_DEBUG #define No_Standard_RangeError #define No_Standard_OutOfRange #endif #include #include #include #include #include #include IntImp_Int2S::IntImp_Int2S(const ThePSurface& surf1, const ThePSurface& surf2, const Standard_Real TolTangency ) : done(Standard_True), empty(Standard_True), tangent(Standard_False), myZerParFunc(surf1,surf2), tol(TolTangency*TolTangency) { ua0 = ThePSurfaceTool::FirstUParameter(surf1); //-- ThePSurfaceTool::UIntervalFirst(surf1); va0 = ThePSurfaceTool::FirstVParameter(surf1); //-- ThePSurfaceTool::VIntervalFirst(surf1); ua1 = ThePSurfaceTool::LastUParameter(surf1); //-- ThePSurfaceTool::UIntervalLast(surf1); va1 = ThePSurfaceTool::LastVParameter(surf1); //-- ThePSurfaceTool::VIntervalLast(surf1); ub0 = ThePSurfaceTool::FirstUParameter(surf2); //-- ThePSurfaceTool::UIntervalFirst(surf2); vb0 = ThePSurfaceTool::FirstVParameter(surf2); //-- ThePSurfaceTool::VIntervalFirst(surf2); ub1 = ThePSurfaceTool::LastUParameter(surf2); //-- ThePSurfaceTool::UIntervalLast(surf2); vb1 = ThePSurfaceTool::LastVParameter(surf2); //-- ThePSurfaceTool::VIntervalLast(surf2); ures1 = ThePSurfaceTool::UResolution(surf1,Precision::Confusion()); vres1 = ThePSurfaceTool::VResolution(surf1,Precision::Confusion()); ures2 = ThePSurfaceTool::UResolution(surf2,Precision::Confusion()); vres2 = ThePSurfaceTool::VResolution(surf2,Precision::Confusion()); } IntImp_Int2S::IntImp_Int2S(const TColStd_Array1OfReal& Param, const ThePSurface& surf1, const ThePSurface& surf2, const Standard_Real TolTangency ) : done(Standard_True), empty(Standard_True), myZerParFunc(surf1,surf2), tol(TolTangency*TolTangency) { math_FunctionSetRoot Rsnld(myZerParFunc,15); //-- Modif lbr 18 MAI ????????????? ua0 = ThePSurfaceTool::FirstUParameter(surf1); //-- ThePSurfaceTool::UIntervalFirst(surf1); va0 = ThePSurfaceTool::FirstVParameter(surf1); //-- ThePSurfaceTool::VIntervalFirst(surf1); ua1 = ThePSurfaceTool::LastUParameter(surf1); //-- ThePSurfaceTool::UIntervalLast(surf1); va1 = ThePSurfaceTool::LastVParameter(surf1); //-- ThePSurfaceTool::VIntervalLast(surf1); ub0 = ThePSurfaceTool::FirstUParameter(surf2); //-- ThePSurfaceTool::UIntervalFirst(surf2); vb0 = ThePSurfaceTool::FirstVParameter(surf2); //-- ThePSurfaceTool::VIntervalFirst(surf2); ub1 = ThePSurfaceTool::LastUParameter(surf2); //-- ThePSurfaceTool::UIntervalLast(surf2); vb1 = ThePSurfaceTool::LastVParameter(surf2); //-- ThePSurfaceTool::VIntervalLast(surf2); ures1 = ThePSurfaceTool::UResolution(surf1,Precision::Confusion()); vres1 = ThePSurfaceTool::VResolution(surf1,Precision::Confusion()); ures2 = ThePSurfaceTool::UResolution(surf2,Precision::Confusion()); vres2 = ThePSurfaceTool::VResolution(surf2,Precision::Confusion()); Perform(Param,Rsnld); } IntImp_ConstIsoparametric IntImp_Int2S:: Perform(const TColStd_Array1OfReal& Param, math_FunctionSetRoot& Rsnld, const IntImp_ConstIsoparametric ChoixIso) { Standard_Real BornInfBuf[3] = {}, BornSupBuf[3] = {}, ToleranceBuf[3] = {}, UVapBuf[3] = {}; Standard_Real UvresBuf[4] = {}; math_Vector BornInf (BornInfBuf, 1, 3), BornSup (BornSupBuf, 1, 3), Tolerance (ToleranceBuf, 1, 3), UVap (UVapBuf, 1, 3); TColStd_Array1OfReal Uvres (UvresBuf[0], 1, 4); IntImp_ConstIsoparametric BestChoix; myZerParFunc.ComputeParameters(ChoixIso,Param,UVap,BornInf,BornSup,Tolerance); Rsnld.SetTolerance(Tolerance); Rsnld.Perform(myZerParFunc,UVap,BornInf,BornSup); BestChoix = ChoixIso; if (Rsnld.IsDone()) { if (Abs(myZerParFunc.Root()) <= tol) { //distance des 2 points // dans la tolerance Rsnld.Root(UVap); empty = Standard_False; tangent = myZerParFunc.IsTangent(UVap,Uvres,BestChoix); pint.SetValue(myZerParFunc.Point(),Uvres(1),Uvres(2),Uvres(3),Uvres(4)); if (!tangent) { d3d = myZerParFunc.Direction(); d2d1 = myZerParFunc.DirectionOnS1(); d2d2 = myZerParFunc.DirectionOnS2(); } } else { empty = Standard_True; } } else { empty = Standard_True; } return BestChoix; } IntImp_ConstIsoparametric IntImp_Int2S:: Perform(const TColStd_Array1OfReal& Param, math_FunctionSetRoot& Rsnld) { gp_Vec DPUV[4]; gp_Pnt P1,P2; Standard_Real Epsuv[4]; Standard_Real DuvBuf[4]; TColStd_Array1OfReal Duv (DuvBuf[0], 1, 4); Standard_Real UVd[4],UVf[4]; IntImp_ConstIsoparametric ChoixIso[4]; IntImp_ConstIsoparametric BestChoix = ChoixRef(0); const ThePSurface& Caro1 = myZerParFunc.AuxillarSurface1(); const ThePSurface& Caro2 = myZerParFunc.AuxillarSurface2(); ThePSurfaceTool::D1(Caro1, Param(1),Param(2),P1,DPUV[0],DPUV[1]); ThePSurfaceTool::D1(Caro2, Param(3),Param(4),P2,DPUV[2],DPUV[3]); Epsuv[0] = ThePSurfaceTool::UResolution(Caro1,Precision::Confusion()); Epsuv[1] = ThePSurfaceTool::VResolution(Caro1,Precision::Confusion()); Epsuv[2] = ThePSurfaceTool::UResolution(Caro2,Precision::Confusion()); Epsuv[3] = ThePSurfaceTool::VResolution(Caro2,Precision::Confusion()); for (Standard_Integer j=0;j<=3;j++) UVd[j] = Param(j+1); empty = Standard_True; Standard_Boolean Tangent = IntImp_ComputeTangence(DPUV,Epsuv,UVd,ChoixIso); if (Tangent) return BestChoix; Standard_Integer i=0; IntImp_ConstIsoparametric CurrentChoix = BestChoix; //-- Modif 17 Mai 93 while (empty && i<= 3) { CurrentChoix = Perform(Param,Rsnld,ChoixIso[i]); if(!empty) { BestChoix = CurrentChoix; } i++; } if (!empty) { // verifier que l on ne deborde pas les frontieres pint.Parameters(Duv(1),Duv(2),Duv(3),Duv(4)); UVd[0] = ua0; //-- ThePSurfaceTool::UIntervalFirst(Caro1); UVd[1] = va0; //-- ThePSurfaceTool::VIntervalFirst(Caro1); UVf[0] = ua1; //-- ThePSurfaceTool::UIntervalLast(Caro1); UVf[1] = va1; //-- ThePSurfaceTool::VIntervalLast(Caro1); UVd[2] = ub0; //-- ThePSurfaceTool::UIntervalFirst(Caro2); UVd[3] = vb0; //-- ubThePSurfaceTool::VIntervalFirst(Caro2); UVf[2] = ub1; //-- ThePSurfaceTool::UIntervalLast(Caro2); UVf[3] = vb1; //-- ThePSurfaceTool::VIntervalLast(Caro2); Standard_Integer Nc,Iiso; if (Duv(1) <= UVd[0]-Epsuv[0]) { Duv(1) = UVd[0]; Nc =0; Iiso=0; } else if (Duv(1) >= UVf[0]+ Epsuv[0]) { Duv(1) = UVf[0]; Nc =0; Iiso=0; } else if (Duv(2) <= UVd[1] -Epsuv[1]) { Duv(2) = UVd[1]; Nc =0; Iiso=1; } else if (Duv(2) >= UVf[1]+Epsuv[1]) { Duv(2) = UVf[1]; Nc =0; Iiso=1; } else if (Duv(3) <= UVd[2]-Epsuv[2]) { Duv(3) = UVd[2]; Nc =2; Iiso = 2; } else if (Duv(3) >= UVf[2]+Epsuv[2]) { Duv(3) = UVf[2]; Nc =2; Iiso = 2; } else if (Duv(4) <= UVd[3]-Epsuv[3]) { Duv(4) = UVd[3]; Nc =2; Iiso = 3; } else if (Duv(4) >= UVf[3]+Epsuv[3]) { Duv(4) = UVf[3]; Nc =2; Iiso =3; } else return BestChoix; // on a gagne empty = Standard_True; BestChoix = ChoixRef(Iiso); //en attendant BestChoix = Perform(Duv,Rsnld,BestChoix); if (!empty) { // verification si l on ne deborde pas sur le carreau // reciproque Nc =3-Nc; if (Duv(Nc) <= UVd[Nc-1]-Epsuv[Nc-1]) Duv(Nc)=UVd[Nc-1]; else if (Duv(Nc) >=UVf[Nc-1]+Epsuv[Nc-1]) Duv(Nc)=UVf[Nc-1]; else if (Duv(Nc+1) <= UVd[Nc]) { Nc = Nc + 1; Duv(Nc)=UVd[Nc-1]; } else if (Duv(Nc+1) >=UVf[Nc]) { Nc = Nc + 1; Duv(Nc)=UVf[Nc-1]; } else return BestChoix; empty = Standard_True; if(Nc == 4) Nc = 0; BestChoix = ChoixRef(Nc); //en attendant BestChoix = Perform(Duv,Rsnld,BestChoix); } } return BestChoix; }