0024002: Overall code and build procedure refactoring -- automatic
[occt.git] / src / BRepLProp / BRepLProp.cxx
CommitLineData
b311480e 1// Created on: 1994-02-24
2// Created by: Laurent BOURESCHE
3// Copyright (c) 1994-1999 Matra Datavision
973c2be1 4// Copyright (c) 1999-2014 OPEN CASCADE SAS
b311480e 5//
973c2be1 6// This file is part of Open CASCADE Technology software library.
b311480e 7//
d5f74e42 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
973c2be1 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.
b311480e 13//
973c2be1 14// Alternatively, this file may be used under the terms of Open CASCADE
15// commercial license or contractual agreement.
7fd59977 16
42cf5bc1 17
18#include <BRepAdaptor_Curve.hxx>
19#include <BRepLProp.hxx>
7fd59977 20#include <BRepLProp_CLProps.hxx>
21#include <gp_Dir.hxx>
22#include <gp_Vec.hxx>
42cf5bc1 23#include <Precision.hxx>
7fd59977 24#include <TopAbs_Orientation.hxx>
42cf5bc1 25#include <TopoDS_Edge.hxx>
7fd59977 26
27//=======================================================================
28//function : Continuity
29//purpose :
30//=======================================================================
7fd59977 31GeomAbs_Shape BRepLProp::Continuity(const BRepAdaptor_Curve& C1,
32 const BRepAdaptor_Curve& C2,
33 const Standard_Real u1,
34 const Standard_Real u2,
35 const Standard_Real tl,
36 const Standard_Real ta)
37{
38 GeomAbs_Shape cont = GeomAbs_C0;
39 Standard_Boolean fini = Standard_False;
40 gp_Vec d1,d2;
41 gp_Dir dir1,dir2;
42 GeomAbs_Shape cont1 = C1.Continuity(), cont2 = C2.Continuity();
43 Standard_Integer n1 = 0, n2 = 0;
44 if (cont1 >= 5) n1 = 3;
45 else if(cont1 == 4) n1 = 2;
46 else if(cont1 == 2) n1 = 1;
47 if (cont2 >= 5) n2 = 3;
48 else if(cont2 == 4) n2 = 2;
49 else if(cont2 == 2) n2 = 1;
50 BRepLProp_CLProps clp1(C1,u1,n1,tl);
51 BRepLProp_CLProps clp2(C2,u2,n2,tl);
52 if(!(clp1.Value().IsEqual(clp2.Value(),tl))) {
53 Standard_Failure::Raise("Courbes non jointives");
54 }
55 Standard_Integer min = Min(n1,n2);
56 if ( min >= 1 ) {
57 d1 = clp1.D1();
58 d2 = clp2.D1();
59 if(C1.Edge().Orientation() == TopAbs_REVERSED) d1.Reverse();
60 if(C2.Edge().Orientation() == TopAbs_REVERSED) d2.Reverse();
61 if(d1.IsEqual(d2,tl,ta)) {
62 cont = GeomAbs_C1;
63 }
64 else if(clp1.IsTangentDefined() && clp2.IsTangentDefined()){
65 clp1.Tangent(dir1);
66 clp2.Tangent(dir2);
67 if(C1.Edge().Orientation() == TopAbs_REVERSED) dir1.Reverse();
68 if(C2.Edge().Orientation() == TopAbs_REVERSED) dir2.Reverse();
69 if(dir1.IsEqual(dir2,ta)){
70 cont = GeomAbs_G1;
71 }
72 fini = Standard_True;
73 }
74 else {fini = Standard_True; }
75 }
76 if ( min >= 2 && !fini ) {
77 d1 = clp1.D2();
78 d2 = clp2.D2();
79 if(d1.IsEqual(d2,tl,ta)){
80 cont = GeomAbs_C2;
81 }
82 }
83 const TopoDS_Edge& E1 = C1.Edge();
84 const TopoDS_Edge& E2 = C2.Edge();
85 if (E1.IsSame(E2) && C1.IsPeriodic() && cont >= GeomAbs_G1)
86 cont = GeomAbs_CN;
87 return cont;
88}
89
90
91//=======================================================================
92//function : Continuity
93//purpose :
94//=======================================================================
95
96GeomAbs_Shape BRepLProp::Continuity(const BRepAdaptor_Curve& C1,
97 const BRepAdaptor_Curve& C2,
98 const Standard_Real u1,
99 const Standard_Real u2)
100{
101 return Continuity(C1,C2,u1,u2,Precision::Confusion(),Precision::Angular());
102}
103
104