mirror of
https://github.com/Open-Cascade-SAS/OCCT.git
synced 2026-05-10 09:30:48 +08:00
127 lines
3.6 KiB
C++
Executable File
127 lines
3.6 KiB
C++
Executable File
// File: BRepClass_Intersector.cxx
|
|
// Created: Thu Nov 19 15:04:21 1992
|
|
// Author: Remi LEQUETTE
|
|
// <rle@phylox>
|
|
|
|
|
|
#include <BRepClass_Intersector.ixx>
|
|
#include <BRep_Tool.hxx>
|
|
#include <BRepAdaptor_Curve2d.hxx>
|
|
#include <BRepAdaptor_Surface.hxx>
|
|
#include <TopoDS_Vertex.hxx>
|
|
#include <TopExp.hxx>
|
|
#include <IntRes2d_Domain.hxx>
|
|
#include <Geom2dLProp_CLProps2d.hxx>
|
|
#include <Geom2d_Curve.hxx>
|
|
#include <ElCLib.hxx>
|
|
#include <Precision.hxx>
|
|
|
|
#include <Geom2d_Line.hxx>
|
|
|
|
#include <Geom2dInt_GInter.hxx>
|
|
|
|
//=======================================================================
|
|
//function : BRepClass_Intersector
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
BRepClass_Intersector::BRepClass_Intersector()
|
|
{
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : Perform
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
void BRepClass_Intersector::Perform(const gp_Lin2d& L,
|
|
const Standard_Real P,
|
|
const Standard_Real Tol,
|
|
const BRepClass_Edge& E)
|
|
{
|
|
|
|
Standard_Real pfbid,plbid;
|
|
if (BRep_Tool::CurveOnSurface(E.Edge(),E.Face(),pfbid,plbid).IsNull()) {
|
|
done = Standard_False; // !IsDone()
|
|
}
|
|
else {
|
|
IntRes2d_Domain DL;
|
|
if(P!=RealLast())
|
|
DL.SetValues(L.Location(),0.,Tol,ElCLib::Value(P,L),P,Tol);
|
|
else
|
|
DL.SetValues(L.Location(),0.,Tol,Standard_True);
|
|
|
|
const TopoDS_Edge& EE = E.Edge();
|
|
const TopoDS_Face& F = E.Face();
|
|
TopoDS_Vertex Vdeb, Vfin;
|
|
TopExp::Vertices(EE, Vdeb, Vfin);
|
|
BRepAdaptor_Curve2d C(EE,F);
|
|
Standard_Real deb = C.FirstParameter(), fin = C.LastParameter();
|
|
gp_Pnt2d pdeb,pfin;
|
|
C.D0(deb,pdeb);
|
|
C.D0(fin,pfin);
|
|
Standard_Real toldeb = 1.e-5, tolfin = 1.e-5;
|
|
#if 0
|
|
// essai de calcul juste des tolerances du domaine
|
|
// qui ne couche pas avec les modeles pourris de
|
|
// styler !!
|
|
BRepAdaptor_Surface S(F);
|
|
gp_Vec2d vdeb,vfin;
|
|
C.D1(deb,pdeb,vdeb);
|
|
C.D1(fin,pfin,vfin);
|
|
gp_Pnt P; gp_Vec DU, DV;
|
|
S.D1(pdeb.X(),pdeb.Y(),P,DU,DV);
|
|
Standard_Real scaldeb = (vdeb.X()*DU + vdeb.Y()*DV).Magnitude();
|
|
scaldeb = Max(scaldeb, 1.e-5);
|
|
toldeb = BRep_Tool::Tolerance(Vdeb)/scaldeb;
|
|
S.D1(pfin.X(),pfin.Y(),P,DU,DV);
|
|
Standard_Real scalfin = (vfin.X()*DU + vfin.Y()*DV).Magnitude();
|
|
scalfin = Max(scalfin, 1.e-5);
|
|
tolfin = BRep_Tool::Tolerance(Vfin)/scalfin;
|
|
#endif
|
|
|
|
IntRes2d_Domain DE(pdeb,deb,toldeb,pfin,fin,tolfin);
|
|
// temporary periodic domain
|
|
if (C.Curve()->IsPeriodic()) {
|
|
DE.SetEquivalentParameters(C.FirstParameter(),
|
|
C.FirstParameter() +
|
|
C.Curve()->LastParameter() -
|
|
C.Curve()->FirstParameter());
|
|
}
|
|
|
|
Handle(Geom2d_Line) GL= new Geom2d_Line(L);
|
|
Geom2dAdaptor_Curve CGA(GL);
|
|
Geom2dInt_GInter Inter(CGA,DL,C,DE,
|
|
Precision::PConfusion(),
|
|
Precision::PIntersection());
|
|
this->SetValues(Inter);
|
|
}
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : LocalGeometry
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
void BRepClass_Intersector::LocalGeometry(const BRepClass_Edge& E,
|
|
const Standard_Real U,
|
|
gp_Dir2d& Tang,
|
|
gp_Dir2d& Norm,
|
|
Standard_Real& C) const
|
|
{
|
|
Standard_Real f,l;
|
|
Geom2dLProp_CLProps2d Prop(BRep_Tool::CurveOnSurface(E.Edge(),E.Face(),f,l),
|
|
U,2,Precision::PConfusion());
|
|
Prop.Tangent(Tang);
|
|
C = Prop.Curvature();
|
|
if (C > Precision::PConfusion())
|
|
Prop.Normal(Norm);
|
|
else
|
|
Norm.SetCoord(Tang.Y(),-Tang.X());
|
|
}
|
|
|
|
|
|
|
|
|
|
|