+++ /dev/null
-// 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.
-
-//========================================================================
-// circulaire tangent a un element de type : - Cercle. +
-// - Ligne. +
-// - Point. +
-// centre sur un deuxieme element de type : - Cercle. +
-// - Ligne. +
-// de rayon donne : Radius. +
-//========================================================================
-
-#include <ElCLib.hxx>
-#include <math_DirectPolynomialRoots.hxx>
-#include <TColStd_Array1OfReal.hxx>
-#include <Standard_NegativeValue.hxx>
-#include <gp_Dir2d.hxx>
-#include <Standard_OutOfRange.hxx>
-#include <StdFail_NotDone.hxx>
-#include <GccEnt_BadQualifier.hxx>
-#include <IntRes2d_Domain.hxx>
-#include <IntRes2d_IntersectionPoint.hxx>
-
-//=========================================================================
-// Cercle tangent : a un cercle Qualified1 (C1). +
-// centre : sur une droite OnLine. +
-// de rayon : Radius. +
-// +
-// On initialise le tableau de solutions cirsol ainsi que tous les +
-// champs. +
-// On elimine en fonction du qualifieur les cas ne presentant pas de +
-// solutions. +
-// On resoud l equation du second degre indiquant que le point de centre +
-// recherche (xc,yc) est a une distance Radius du cercle C1 et +
-// sur la droite OnLine. +
-// Les solutions sont representees par les cercles : +
-// - de centre Pntcen(xc,yc) +
-// - de rayon Radius. +
-//=========================================================================
-
-GccGeo_Circ2dTanOnRad::
- GccGeo_Circ2dTanOnRad (const TheQCurve& Qualified1,
- const gp_Lin2d& OnLine ,
- const Standard_Real Radius ,
- const Standard_Real Tolerance ):
-
-//=========================================================================
-// Initialisation des champs. +
-//=========================================================================
-
- cirsol(1,8) ,
- qualifier1(1,8) ,
- TheSame1(1,8) ,
- pnttg1sol(1,8) ,
- pntcen3(1,8) ,
- par1sol(1,8) ,
- pararg1(1,8) ,
- parcen3(1,8)
-{
-
-//=========================================================================
-// Traitement. +
-//=========================================================================
-
- gp_Dir2d dirx(1.0,0.0);
- Standard_Real Tol = Abs(Tolerance);
- Standard_Real thefirst = -100000.;
- Standard_Real thelast = 100000.;
- Standard_Real firstparam;
- Standard_Real lastparam;
- WellDone = Standard_False;
- NbrSol = 0;
- if (!(Qualified1.IsEnclosed() || Qualified1.IsEnclosing() ||
- Qualified1.IsOutside() || Qualified1.IsUnqualified())) {
- GccEnt_BadQualifier::Raise();
- return;
- }
- Standard_Integer nbrcote1 = 0;
- TColStd_Array1OfReal Coef(1,2);
- TheCurve Cu1 = Qualified1.Qualified();
-
- if (Radius < 0.0) { Standard_NegativeValue::Raise(); }
- else {
- if (Qualified1.IsEnclosed()) {
-// ===========================
- nbrcote1 = 1;
- Coef(1) = Radius;
- }
- else if(Qualified1.IsOutside()) {
-// ===============================
- nbrcote1 = 1;
- Coef(1) = -Radius;
- }
- else if(Qualified1.IsUnqualified()) {
-// ===================================
- nbrcote1 = 2;
- Coef(1) = Radius;
- Coef(2) = -Radius;
- }
- IntRes2d_Domain D1;
- TheIntConicCurve Intp;
- for (Standard_Integer jcote1 = 1 ; jcote1 <= nbrcote1 ; jcote1++) {
- Handle(TheHParGenCurve) HCu1 = new TheHParGenCurve(Cu1);
- TheParGenCurve C2(HCu1,Coef(jcote1));
- firstparam = Max(TheCurvePGTool::FirstParameter(C2),thefirst);
- lastparam = Min(TheCurvePGTool::LastParameter(C2),thelast);
- IntRes2d_Domain D2(TheCurvePGTool::Value(C2,firstparam),firstparam,Tol,
- TheCurvePGTool::Value(C2,lastparam),lastparam,Tol);
- Intp.Perform(OnLine,D1,C2,D2,Tol,Tol);
- if (Intp.IsDone()) {
- if (!Intp.IsEmpty()) {
- for (Standard_Integer i = 1 ; i <= Intp.NbPoints() ; i++) {
- NbrSol++;
- gp_Pnt2d Center(Intp.Point(i).Value());
- cirsol(NbrSol) = gp_Circ2d(gp_Ax2d(Center,dirx),Radius);
-// =======================================================
- qualifier1(NbrSol) = Qualified1.Qualifier();
- TheSame1(NbrSol) = 0;
- pararg1(NbrSol) = Intp.Point(i).ParamOnSecond();
- parcen3(NbrSol) = Intp.Point(i).ParamOnFirst();
- par1sol(NbrSol)=ElCLib::Parameter(cirsol(NbrSol),
- pnttg1sol(NbrSol));
- pnttg1sol(NbrSol) = gp_Pnt2d(TheTool::Value(Cu1,pararg1(NbrSol)));
- pntcen3(NbrSol) = Center;
- }
- }
- WellDone = Standard_True;
- }
- }
- }
- }
-
-//=========================================================================
-// Cercle tangent : a un cercle Qualified1 (C1). +
-// centre : sur une droite OnLine. +
-// de rayon : Radius. +
-// +
-// On initialise le tableau de solutions cirsol ainsi que tous les +
-// champs. +
-// On elimine en fonction du qualifieur les cas ne presentant pas de +
-// solutions. +
-// On resoud l equation du second degre indiquant que le point de centre +
-// recherche (xc,yc) est a une distance Radius du cercle C1 et +
-// sur la droite OnLine. +
-// Les solutions sont representees par les cercles : +
-// - de centre Pntcen(xc,yc) +
-// - de rayon Radius. +
-//=========================================================================
-
-GccGeo_Circ2dTanOnRad::
- GccGeo_Circ2dTanOnRad (const TheQCurve& Qualified1,
- const gp_Circ2d& OnCirc ,
- const Standard_Real Radius ,
- const Standard_Real Tolerance ):
-
-//=========================================================================
-// Initialisation des champs. +
-//=========================================================================
-
- cirsol(1,8) ,
- qualifier1(1,8) ,
- TheSame1(1,8) ,
- pnttg1sol(1,8) ,
- pntcen3(1,8) ,
- par1sol(1,8) ,
- pararg1(1,8) ,
- parcen3(1,8)
-{
-
-//=========================================================================
-// Traitement. +
-//=========================================================================
-
- gp_Dir2d dirx(1.0,0.0);
- Standard_Real thefirst = -100000.;
- Standard_Real thelast = 100000.;
- Standard_Real firstparam;
- Standard_Real lastparam;
- Standard_Real Tol = Abs(Tolerance);
- Standard_Integer nbrcote1=0;
- WellDone = Standard_False;
- NbrSol = 0;
- if (!(Qualified1.IsEnclosed() || Qualified1.IsEnclosing() ||
- Qualified1.IsOutside() || Qualified1.IsUnqualified())) {
- GccEnt_BadQualifier::Raise();
- return;
- }
- TColStd_Array1OfReal cote1(1,2);
- TheCurve Cu1 = Qualified1.Qualified();
-
- if (Radius < 0.0) {
- Standard_NegativeValue::Raise();
- }
- else {
- if (Qualified1.IsEnclosed()) {
-// ===========================
- nbrcote1 = 1;
- cote1(1) = Radius;
- }
- else if(Qualified1.IsOutside()) {
-// ===============================
- nbrcote1 = 1;
- cote1(1) = -Radius;
- }
- else if(Qualified1.IsUnqualified()) {
-// ===================================
- nbrcote1 = 2;
- cote1(1) = Radius;
- cote1(2) = -Radius;
- }
- IntRes2d_Domain D1(ElCLib::Value(0.,OnCirc), 0.,Tol,
- ElCLib::Value(2.*M_PI,OnCirc),2.*M_PI,Tol);
- D1.SetEquivalentParameters(0.,2.*M_PI);
- TheIntConicCurve Intp;
- for (Standard_Integer jcote1 = 1 ; jcote1 <= nbrcote1 ; jcote1++) {
- Handle(TheHParGenCurve) HCu1 = new TheHParGenCurve(Cu1);
- TheParGenCurve C2(HCu1,cote1(jcote1));
- firstparam = Max(TheCurvePGTool::FirstParameter(C2),thefirst);
- lastparam = Min(TheCurvePGTool::LastParameter(C2),thelast);
- IntRes2d_Domain D2(TheCurvePGTool::Value(C2,firstparam),firstparam,Tol,
- TheCurvePGTool::Value(C2,lastparam),lastparam,Tol);
- Intp.Perform(OnCirc,D1,C2,D2,Tol,Tol);
- if (Intp.IsDone()) {
- if (!Intp.IsEmpty()) {
- for (Standard_Integer i = 1 ; i <= Intp.NbPoints() ; i++) {
- NbrSol++;
- gp_Pnt2d Center(Intp.Point(i).Value());
- cirsol(NbrSol) = gp_Circ2d(gp_Ax2d(Center,dirx),Radius);
-// =======================================================
- qualifier1(NbrSol) = Qualified1.Qualifier();
- TheSame1(NbrSol) = 0;
- pararg1(NbrSol) = Intp.Point(i).ParamOnSecond();
- parcen3(NbrSol) = Intp.Point(i).ParamOnFirst();
- par1sol(NbrSol)=ElCLib::Parameter(cirsol(NbrSol),
- pnttg1sol(NbrSol));
- pnttg1sol(NbrSol) = gp_Pnt2d(TheTool::Value(Cu1,pararg1(NbrSol)));
- pntcen3(NbrSol) = Center;
- }
- }
- WellDone = Standard_True;
- }
- }
- }
- }
-
-//=========================================================================
-// Cercle tangent : a un cercle Qualified1 (C1). +
-// centre : sur une droite OnLine. +
-// de rayon : Radius. +
-// +
-// On initialise le tableau de solutions cirsol ainsi que tous les +
-// champs. +
-// On elimine en fonction du qualifieur les cas ne presentant pas de +
-// solutions. +
-// On resoud l equation du second degre indiquant que le point de centre +
-// recherche (xc,yc) est a une distance Radius du cercle C1 et +
-// sur la droite OnLine. +
-// Les solutions sont representees par les cercles : +
-// - de centre Pntcen(xc,yc) +
-// - de rayon Radius. +
-//=========================================================================
-
-GccGeo_Circ2dTanOnRad::
- GccGeo_Circ2dTanOnRad (const GccEnt_QualifiedCirc& Qualified1,
- const TheCurve& OnCurv ,
- const Standard_Real Radius ,
- const Standard_Real Tolerance ):
-
-//=========================================================================
-// Initialisation des champs. +
-//=========================================================================
-
- cirsol(1,8) ,
- qualifier1(1,8) ,
- TheSame1(1,8) ,
- pnttg1sol(1,8) ,
- pntcen3(1,8) ,
- par1sol(1,8) ,
- pararg1(1,8) ,
- parcen3(1,8)
-{
-
-//=========================================================================
-// Traitement. +
-//=========================================================================
-
- gp_Dir2d dirx(1.0,0.0);
- Standard_Real thefirst = -100000.;
- Standard_Real thelast = 100000.;
- Standard_Real firstparam;
- Standard_Real lastparam;
- Standard_Real Tol = Abs(Tolerance);
- Standard_Integer nbrcote1=0;
- WellDone = Standard_False;
- NbrSol = 0;
- if (!(Qualified1.IsEnclosed() || Qualified1.IsEnclosing() ||
- Qualified1.IsOutside() || Qualified1.IsUnqualified())) {
- GccEnt_BadQualifier::Raise();
- return;
- }
- TColStd_Array1OfReal cote1(1,2);
- gp_Circ2d C1 = Qualified1.Qualified();
- gp_Pnt2d center1(C1.Location());
- Standard_Real R1 = C1.Radius();
-
- if (Radius < 0.0) {
- Standard_NegativeValue::Raise();
- }
- else {
- if (Qualified1.IsEnclosed()) {
-// ===========================
- nbrcote1 = 1;
- cote1(1) = Radius;
- }
- else if(Qualified1.IsOutside()) {
-// ===============================
- nbrcote1 = 1;
- cote1(1) = -Radius;
- }
- else if(Qualified1.IsUnqualified()) {
-// ===================================
- nbrcote1 = 2;
- cote1(1) = Radius;
- cote1(2) = -Radius;
- }
- TheIntConicCurve Intp;
- for (Standard_Integer jcote1 = 1 ; jcote1 <= nbrcote1 ; jcote1++) {
- gp_Circ2d Circ(C1.XAxis(),R1 + cote1(jcote1));
- IntRes2d_Domain D1(ElCLib::Value(0.,Circ), 0.,Tol,
- ElCLib::Value(2.*M_PI,Circ),2.*M_PI,Tol);
- D1.SetEquivalentParameters(0.,2.*M_PI);
- firstparam = Max(TheTool::FirstParameter(OnCurv),thefirst);
- lastparam = Min(TheTool::LastParameter(OnCurv),thelast);
- IntRes2d_Domain D2(TheTool::Value(OnCurv,firstparam),firstparam,Tol,
- TheTool::Value(OnCurv,lastparam),lastparam,Tol);
- Intp.Perform(Circ,D1,OnCurv,D2,Tol,Tol);
- if (Intp.IsDone()) {
- if (!Intp.IsEmpty()) {
- for (Standard_Integer i = 1 ; i <= Intp.NbPoints() ; i++) {
- NbrSol++;
- gp_Pnt2d Center(Intp.Point(i).Value());
- cirsol(NbrSol) = gp_Circ2d(gp_Ax2d(Center,dirx),Radius);
-// =======================================================
- Standard_Real distcc1 = Center.Distance(center1);
- if (!Qualified1.IsUnqualified()) {
- qualifier1(NbrSol) = Qualified1.Qualifier();
- }
- else if (Abs(distcc1+Radius-R1) < Tol) {
- qualifier1(NbrSol) = GccEnt_enclosed;
- }
- else if (Abs(distcc1-R1-Radius) < Tol) {
- qualifier1(NbrSol) = GccEnt_outside;
- }
- else { qualifier1(NbrSol) = GccEnt_enclosing; }
- TheSame1(NbrSol) = 0;
- pararg1(NbrSol) = Intp.Point(i).ParamOnFirst();
- parcen3(NbrSol) = Intp.Point(i).ParamOnSecond();
- par1sol(NbrSol)=ElCLib::Parameter(cirsol(NbrSol),
- pnttg1sol(NbrSol));
- pnttg1sol(NbrSol) = ElCLib::Value(pararg1(NbrSol),C1);
- pntcen3(NbrSol) = Center;
- }
- }
- WellDone = Standard_True;
- }
- }
- }
- }
-
-//=========================================================================
-// Cercle tangent : a un cercle Qualified1 (C1). +
-// centre : sur une droite OnLine. +
-// de rayon : Radius. +
-// +
-// On initialise le tableau de solutions cirsol ainsi que tous les +
-// champs. +
-// On elimine en fonction du qualifieur les cas ne presentant pas de +
-// solutions. +
-// On resoud l equation du second degre indiquant que le point de centre +
-// recherche (xc,yc) est a une distance Radius du cercle C1 et +
-// sur la droite OnLine. +
-// Les solutions sont representees par les cercles : +
-// - de centre Pntcen(xc,yc) +
-// - de rayon Radius. +
-//=========================================================================
-
-GccGeo_Circ2dTanOnRad::
- GccGeo_Circ2dTanOnRad (const GccEnt_QualifiedLin& Qualified1,
- const TheCurve& OnCurv ,
- const Standard_Real Radius ,
- const Standard_Real Tolerance ):
-
-//=========================================================================
-// Initialisation des champs. +
-//=========================================================================
-
- cirsol(1,8) ,
- qualifier1(1,8) ,
- TheSame1(1,8) ,
- pnttg1sol(1,8) ,
- pntcen3(1,8) ,
- par1sol(1,8) ,
- pararg1(1,8) ,
- parcen3(1,8)
-{
-
-//=========================================================================
-// Traitement. +
-//=========================================================================
-
- gp_Dir2d dirx(1.0,0.0);
- Standard_Real thefirst = -100000.;
- Standard_Real thelast = 100000.;
- Standard_Real firstparam;
- Standard_Real lastparam;
- Standard_Real Tol = Abs(Tolerance);
- WellDone = Standard_False;
- NbrSol = 0;
- if (!(Qualified1.IsEnclosed() ||
- Qualified1.IsOutside() || Qualified1.IsUnqualified())) {
- GccEnt_BadQualifier::Raise();
- return;
- }
- Standard_Integer nbrcote1=0;
- TColStd_Array1OfReal cote1(1,2);
- gp_Lin2d L1 = Qualified1.Qualified();
- gp_Pnt2d origin1(L1.Location());
- gp_Dir2d dir1(L1.Direction());
- gp_Dir2d norm1(-dir1.Y(),dir1.X());
-
- if (Radius < 0.0) {
- Standard_NegativeValue::Raise();
- }
- else {
- if (Qualified1.IsEnclosed()) {
-// ===========================
- nbrcote1 = 1;
- cote1(1) = Radius;
- }
- else if(Qualified1.IsOutside()) {
-// ===============================
- nbrcote1 = 1;
- cote1(1) = -Radius;
- }
- else if(Qualified1.IsUnqualified()) {
-// ===================================
- nbrcote1 = 2;
- cote1(1) = Radius;
- cote1(2) = -Radius;
- }
- TheIntConicCurve Intp;
- for (Standard_Integer jcote1 = 1 ; jcote1 <= nbrcote1 ; jcote1++) {
- gp_Pnt2d Point(dir1.XY()+cote1(jcote1)*norm1.XY());
- gp_Lin2d Line(Point,dir1); // ligne avec deport.
- IntRes2d_Domain D1;
- firstparam = Max(TheTool::FirstParameter(OnCurv),thefirst);
- lastparam = Min(TheTool::LastParameter(OnCurv),thelast);
- IntRes2d_Domain D2(TheTool::Value(OnCurv,firstparam),firstparam,Tol,
- TheTool::Value(OnCurv,lastparam),lastparam,Tol);
- Intp.Perform(Line,D1,OnCurv,D2,Tol,Tol);
- if (Intp.IsDone()) {
- if (!Intp.IsEmpty()) {
- for (Standard_Integer i = 1 ; i <= Intp.NbPoints() ; i++) {
- NbrSol++;
- gp_Pnt2d Center(Intp.Point(i).Value());
- cirsol(NbrSol) = gp_Circ2d(gp_Ax2d(Center,dirx),Radius);
-// =======================================================
- gp_Dir2d dc1(origin1.XY()-Center.XY());
- if (!Qualified1.IsUnqualified()) {
- qualifier1(NbrSol) = Qualified1.Qualifier();
- }
- else if (dc1.Dot(norm1) > 0.0) {
- qualifier1(NbrSol) = GccEnt_outside;
- }
- else { qualifier1(NbrSol) = GccEnt_enclosed; }
- TheSame1(NbrSol) = 0;
- pararg1(NbrSol) = Intp.Point(i).ParamOnFirst();
- parcen3(NbrSol) = Intp.Point(i).ParamOnSecond();
- par1sol(NbrSol)=ElCLib::Parameter(cirsol(NbrSol),
- pnttg1sol(NbrSol));
- pnttg1sol(NbrSol) = ElCLib::Value(pararg1(NbrSol),L1);
- pntcen3(NbrSol) = Center;
- }
- }
- WellDone = Standard_True;
- }
- }
- }
- }
-
-//=========================================================================
-// Cercle tangent : a un cercle Qualified1 (C1). +
-// centre : sur une droite OnLine. +
-// de rayon : Radius. +
-// +
-// On initialise le tableau de solutions cirsol ainsi que tous les +
-// champs. +
-// On elimine en fonction du qualifieur les cas ne presentant pas de +
-// solutions. +
-// On resoud l equation du second degre indiquant que le point de centre +
-// recherche (xc,yc) est a une distance Radius du cercle C1 et +
-// sur la droite OnLine. +
-// Les solutions sont representees par les cercles : +
-// - de centre Pntcen(xc,yc) +
-// - de rayon Radius. +
-//=========================================================================
-
-GccGeo_Circ2dTanOnRad::
- GccGeo_Circ2dTanOnRad (const TheQCurve& Qualified1,
- const TheCurve& OnCurv ,
- const Standard_Real Radius ,
- const Standard_Real Tolerance ):
-
-//=========================================================================
-// Initialisation des champs. +
-//=========================================================================
-
- cirsol(1,8) ,
- qualifier1(1,8) ,
- TheSame1(1,8) ,
- pnttg1sol(1,8) ,
- pntcen3(1,8) ,
- par1sol(1,8) ,
- pararg1(1,8) ,
- parcen3(1,8)
-{
-
-//=========================================================================
-// Traitement. +
-//=========================================================================
-
- gp_Dir2d dirx(1.0,0.0);
- Standard_Real thefirst = -100000.;
- Standard_Real thelast = 100000.;
- Standard_Real firstparam;
- Standard_Real lastparam;
- Standard_Real Tol = Abs(Tolerance);
- Standard_Integer nbrcote1=0;
- WellDone = Standard_False;
- NbrSol = 0;
- if (!(Qualified1.IsEnclosed() || Qualified1.IsEnclosing() ||
- Qualified1.IsOutside() || Qualified1.IsUnqualified())) {
- GccEnt_BadQualifier::Raise();
- return;
- }
- TColStd_Array1OfReal cote1(1,2);
- TheCurve Cu1 = Qualified1.Qualified();
-
- if (Radius < 0.0) {
- Standard_NegativeValue::Raise();
- }
- else {
- if (Qualified1.IsEnclosed()) {
-// ===========================
- nbrcote1 = 1;
- cote1(1) = Radius;
- }
- else if(Qualified1.IsOutside()) {
-// ===============================
- nbrcote1 = 1;
- cote1(1) = -Radius;
- }
- else if(Qualified1.IsUnqualified()) {
-// ===================================
- nbrcote1 = 2;
- cote1(1) = Radius;
- cote1(2) = -Radius;
- }
- TheIntCurveCurve Intp;
- for (Standard_Integer jcote1 = 1 ; jcote1 <= nbrcote1 ; jcote1++) {
- Handle(TheHParGenCurve) HCu1 = new TheHParGenCurve(Cu1);
- TheParGenCurve C1(HCu1,cote1(jcote1));
- firstparam = Max(TheCurvePGTool::FirstParameter(C1),thefirst);
- lastparam = Min(TheCurvePGTool::LastParameter(C1),thelast);
- IntRes2d_Domain D1(TheCurvePGTool::Value(C1,firstparam),firstparam,Tol,
- TheCurvePGTool::Value(C1,lastparam),lastparam,Tol);
- Handle(TheHParGenCurve) HOnCurv = new TheHParGenCurve(OnCurv);
- TheParGenCurve C2(HOnCurv);
- firstparam = Max(TheCurvePGTool::FirstParameter(C2),thefirst);
- lastparam = Min(TheCurvePGTool::LastParameter(C2),thelast);
- IntRes2d_Domain D2(TheCurvePGTool::Value(C2,firstparam),firstparam,Tol,
- TheCurvePGTool::Value(C2,lastparam),lastparam,Tol);
- Intp.Perform(C1,D1,C2,D2,Tol,Tol);
- if (Intp.IsDone()) {
- if (!Intp.IsEmpty()) {
- for (Standard_Integer i = 1 ; i <= Intp.NbPoints() ; i++) {
- NbrSol++;
- gp_Pnt2d Center(Intp.Point(i).Value());
- cirsol(NbrSol) = gp_Circ2d(gp_Ax2d(Center,dirx),Radius);
-// =======================================================
- qualifier1(NbrSol) = Qualified1.Qualifier();
- TheSame1(NbrSol) = 0;
- pararg1(NbrSol) = Intp.Point(i).ParamOnFirst();
- parcen3(NbrSol) = Intp.Point(i).ParamOnSecond();
- par1sol(NbrSol)=ElCLib::Parameter(cirsol(NbrSol),
- pnttg1sol(NbrSol));
- pnttg1sol(NbrSol) = gp_Pnt2d(TheTool::Value(Cu1,pararg1(NbrSol)));
- pntcen3(NbrSol) = Center;
- }
- }
- WellDone = Standard_True;
- }
- }
- }
- }
-
-//=========================================================================
-// Cercle tangent : a un cercle Qualified1 (C1). +
-// centre : sur une droite OnLine. +
-// de rayon : Radius. +
-// +
-// On initialise le tableau de solutions cirsol ainsi que tous les +
-// champs. +
-// On elimine en fonction du qualifieur les cas ne presentant pas de +
-// solutions. +
-// On resoud l equation du second degre indiquant que le point de centre +
-// recherche (xc,yc) est a une distance Radius du cercle C1 et +
-// sur la droite OnLine. +
-// Les solutions sont representees par les cercles : +
-// - de centre Pntcen(xc,yc) +
-// - de rayon Radius. +
-//=========================================================================
-
-GccGeo_Circ2dTanOnRad::
- GccGeo_Circ2dTanOnRad (const gp_Pnt2d& Point1 ,
- const TheCurve& OnCurv ,
- const Standard_Real Radius ,
- const Standard_Real Tolerance ):
-
-//=========================================================================
-// Initialisation des champs. +
-//=========================================================================
-
- cirsol(1,8) ,
- qualifier1(1,8) ,
- TheSame1(1,8) ,
- pnttg1sol(1,8) ,
- pntcen3(1,8) ,
- par1sol(1,8) ,
- pararg1(1,8) ,
- parcen3(1,8)
-{
-
-//=========================================================================
-// Traitement. +
-//=========================================================================
-
- gp_Dir2d dirx(1.0,0.0);
- Standard_Real thefirst = -100000.;
- Standard_Real thelast = 100000.;
- Standard_Real firstparam;
- Standard_Real lastparam;
- Standard_Real Tol = Abs(Tolerance);
- WellDone = Standard_False;
- NbrSol = 0;
-
- if (Radius < 0.0) {
- Standard_NegativeValue::Raise();
- }
- else {
-// gp_Dir2d Dir(-y1dir,x1dir);
- gp_Circ2d Circ(gp_Ax2d(Point1,gp_Dir2d(1.,0.)),Radius);
- IntRes2d_Domain D1(ElCLib::Value(0.,Circ), 0.,Tol,
- ElCLib::Value(2.*M_PI,Circ),2*M_PI,Tol);
- D1.SetEquivalentParameters(0.,2.*M_PI);
- firstparam = Max(TheTool::FirstParameter(OnCurv),thefirst);
- lastparam = Min(TheTool::LastParameter(OnCurv),thelast);
- IntRes2d_Domain D2(TheTool::Value(OnCurv,firstparam),firstparam,Tol,
- TheTool::Value(OnCurv,lastparam),lastparam,Tol);
- TheIntConicCurve Intp(Circ,D1,OnCurv,D2,Tol,Tol);
- if (Intp.IsDone()) {
- if (!Intp.IsEmpty()) {
- for (Standard_Integer i = 1 ; i <= Intp.NbPoints() ; i++) {
- NbrSol++;
- gp_Pnt2d Center(Intp.Point(i).Value());
- cirsol(NbrSol) = gp_Circ2d(gp_Ax2d(Center,dirx),Radius);
-// =======================================================
- qualifier1(NbrSol) = GccEnt_noqualifier;
- TheSame1(NbrSol) = 0;
- pararg1(NbrSol) = Intp.Point(i).ParamOnFirst();
- parcen3(NbrSol) = Intp.Point(i).ParamOnSecond();
- par1sol(NbrSol)=ElCLib::Parameter(cirsol(NbrSol),
- pnttg1sol(NbrSol));
- pnttg1sol(NbrSol) = Point1;
- pntcen3(NbrSol) = Center;
- }
- WellDone = Standard_True;
- }
- }
- }
- }
-
-//=========================================================================
-
-Standard_Boolean GccGeo_Circ2dTanOnRad::
- IsDone () const { return WellDone; }
-
-Standard_Integer GccGeo_Circ2dTanOnRad::
- NbSolutions () const { return NbrSol; }
-
-gp_Circ2d GccGeo_Circ2dTanOnRad::
- ThisSolution (const Standard_Integer Index) const
-{
-
- if (Index > NbrSol || Index <= 0)
- Standard_OutOfRange::Raise();
-
- return cirsol(Index);
-}
-
-void GccGeo_Circ2dTanOnRad::
- WhichQualifier(const Standard_Integer Index ,
- GccEnt_Position& Qualif1 ) const
-{
- if (!WellDone) { StdFail_NotDone::Raise(); }
- else if (Index <= 0 ||Index > NbrSol) { Standard_OutOfRange::Raise(); }
- else {
- Qualif1 = qualifier1(Index);
- }
-}
-
-void GccGeo_Circ2dTanOnRad::
- Tangency1 (const Standard_Integer Index,
- Standard_Real& ParSol,
- Standard_Real& ParArg,
- gp_Pnt2d& PntSol) const{
- if (!WellDone) {
- StdFail_NotDone::Raise();
- }
- else if (Index <= 0 ||Index > NbrSol) {
- Standard_OutOfRange::Raise();
- }
- else {
- ParSol = par1sol(Index);
- ParArg = pararg1(Index);
- PntSol = gp_Pnt2d(pnttg1sol(Index));
- }
- }
-
-void GccGeo_Circ2dTanOnRad::
- CenterOn3 (const Standard_Integer Index,
- Standard_Real& ParArg,
- gp_Pnt2d& PntSol) const {
- if (!WellDone) {
- StdFail_NotDone::Raise();
- }
- else if (Index <= 0 ||Index > NbrSol) {
- Standard_OutOfRange::Raise();
- }
- else {
- ParArg = parcen3(Index);
- PntSol = pnttg1sol(Index);
- }
- }
-
-Standard_Boolean GccGeo_Circ2dTanOnRad::
- IsTheSame1 (const Standard_Integer Index) const
-{
- if (!WellDone) StdFail_NotDone::Raise();
- if (Index <= 0 ||Index > NbrSol) Standard_OutOfRange::Raise();
-
- if (TheSame1(Index) == 0)
- return Standard_False;
-
- return Standard_True;
-}
projects / occt.git / blobdiff