0027960: Configuration - fix compilation of OSD_Directory with MinGW-w64
[occt.git] / src / ChFiKPart / ChFiKPart_ComputeData_Rotule.cxx
1 // Created on: 1994-06-27
2 // Created by: Laurent BOURESCHE
3 // Copyright (c) 1994-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
5 //
6 // This file is part of Open CASCADE Technology software library.
7 //
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.
13 //
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
16
17
18 #include <Adaptor3d_HSurface.hxx>
19 #include <ChFiDS_Spine.hxx>
20 #include <ChFiDS_SurfData.hxx>
21 #include <ChFiKPart_ComputeData.hxx>
22 #include <ChFiKPart_ComputeData_Fcts.hxx>
23 #include <ElCLib.hxx>
24 #include <ElSLib.hxx>
25 #include <Geom2d_Circle.hxx>
26 #include <Geom2d_Line.hxx>
27 #include <Geom_Circle.hxx>
28 #include <Geom_Line.hxx>
29 #include <Geom_ToroidalSurface.hxx>
30 #include <gp.hxx>
31 #include <gp_Ax2.hxx>
32 #include <gp_Ax3.hxx>
33 #include <gp_Ax22d.hxx>
34 #include <gp_Dir.hxx>
35 #include <gp_Dir2d.hxx>
36 #include <gp_Lin.hxx>
37 #include <gp_Lin2d.hxx>
38 #include <gp_Pln.hxx>
39 #include <gp_Pnt.hxx>
40 #include <gp_Pnt2d.hxx>
41 #include <gp_Vec.hxx>
42 #include <IntAna_QuadQuadGeo.hxx>
43 #include <Precision.hxx>
44 #include <TopOpeBRepDS_DataStructure.hxx>
45
46 //=======================================================================
47 //function : MakeRotule
48 //Purpose  : cas plan/plan.
49 //=======================================================================
50 Standard_Boolean ChFiKPart_MakeRotule(TopOpeBRepDS_DataStructure& DStr,
51                                       const Handle(ChFiDS_SurfData)& Data, 
52                                       const gp_Pln& pl, 
53                                       const gp_Pln& pl1, 
54                                       const gp_Pln& pl2, 
55                                       const TopAbs_Orientation opl,
56                                       const TopAbs_Orientation opl1,
57                                       const TopAbs_Orientation opl2,
58                                       const Standard_Real r, 
59                                       const TopAbs_Orientation ofpl)
60 {
61
62   //calcul du tore.
63   //---------------
64   gp_Ax3 pos = pl.Position();
65   gp_Dir dpl = pos.XDirection().Crossed(pos.YDirection());
66   gp_Dir dfpl = dpl;
67   gp_Dir dplnat = dpl;
68   if (opl == TopAbs_REVERSED) { dpl.Reverse(); }
69   if (ofpl == TopAbs_REVERSED) { dfpl.Reverse(); }
70   pos = pl1.Position();
71   gp_Dir dpl1 = pos.XDirection().Crossed(pos.YDirection());
72   if (opl1 == TopAbs_REVERSED) { dpl1.Reverse(); }
73   pos = pl2.Position();
74   gp_Dir dpl2 = pos.XDirection().Crossed(pos.YDirection());
75   if (opl2 == TopAbs_REVERSED) { dpl2.Reverse(); }
76
77   Standard_Real alpha = dpl1.Angle(dpl2);
78
79   IntAna_QuadQuadGeo LInt (pl1,pl2,Precision::Angular(),
80                            Precision::Confusion());
81   gp_Pnt ptor,pcirc;
82   if (LInt.IsDone()) {
83
84     pcirc = ElCLib::Value(ElCLib::Parameter(LInt.Line(1),pl.Location()),
85                          LInt.Line(1));
86     ptor.SetCoord(pcirc.X()+r*dpl.X(),
87                   pcirc.Y()+r*dpl.Y(),
88                   pcirc.Z()+r*dpl.Z());
89   }
90   else { return Standard_False; }
91
92   gp_Ax3 ppos(ptor,dpl.Reversed(),dpl1);
93   if(ppos.YDirection().Dot(dpl2) < 0.) ppos.YReverse();
94   Handle(Geom_ToroidalSurface) 
95     gtor = new Geom_ToroidalSurface(ppos,r,r);
96   Data->ChangeSurf(ChFiKPart_IndexSurfaceInDS(gtor,DStr));
97
98   //on compare l orientation du tore a celle de la face en bout.
99   //------------------------------------------------------------
100   gp_Pnt pp;
101   gp_Vec du,dv;
102   ElSLib::TorusD1(0.,M_PI/2,ppos,r,r,pp,du,dv);
103   gp_Dir drot(du.Crossed(dv));
104   Standard_Boolean reversecur = ( drot.Dot(dplnat) <= 0. );
105   Standard_Boolean reversefil = ( drot.Dot(dfpl) <= 0. );
106   if (reversefil) { Data->ChangeOrientation() = TopAbs_REVERSED; }
107   else { Data->ChangeOrientation() = TopAbs_FORWARD; }
108
109   //on charge les FaceInterferences avec les pcurves et courbes 3d.
110   //-----------------------------------------------------------------
111   
112   //du cote du plan
113   //---------------
114   gp_Ax2 circAx2 = ppos.Ax2();
115   circAx2.SetLocation(pcirc);
116   Handle(Geom_Circle) GC = new Geom_Circle(circAx2,r);
117   Standard_Real u,v;
118   ElSLib::Parameters(pl,pcirc,u,v);
119   gp_Pnt2d p2dcirc(u,v);
120   gp_Dir2d dx2d(dpl1.Dot(pl.Position().XDirection()),
121                 dpl1.Dot(pl.Position().YDirection()));
122   gp_Dir2d dy2d(ppos.YDirection().Dot(pl.Position().XDirection()),
123                 ppos.YDirection().Dot(pl.Position().YDirection()));
124   gp_Ax22d circ2dax(p2dcirc,dx2d,dy2d);
125   Handle(Geom2d_Circle) GC2d = new Geom2d_Circle(circ2dax,r);
126   gp_Pnt2d p2dlin(0.,M_PI/2);
127   Handle(Geom2d_Line) GL2d = new Geom2d_Line(p2dlin,gp::DX2d());
128   TopAbs_Orientation trans = TopAbs_REVERSED; 
129   if (reversecur) trans = TopAbs_FORWARD; 
130   Data->ChangeInterferenceOnS1().
131     SetInterference(ChFiKPart_IndexCurveInDS(GC,DStr),trans,GC2d,GL2d);
132   
133   //du cote pointu
134   //--------------
135   Handle(Geom_Curve) bid;
136   Handle(Geom2d_Curve) bid2d;
137   p2dlin.SetCoord(0.,M_PI);
138   Handle(Geom2d_Line) GL2dcoin = new Geom2d_Line(p2dlin,gp::DX2d());
139   Data->ChangeInterferenceOnS2().
140     SetInterference(ChFiKPart_IndexCurveInDS(bid,DStr),trans,bid2d,GL2dcoin);
141
142   //et les points
143   //-------------
144   Data->ChangeVertexFirstOnS1().SetPoint(pp);
145   ElSLib::TorusD0(alpha,M_PI/2,ppos,r,r,pp);
146   Data->ChangeVertexLastOnS1().SetPoint(pp);
147   Data->ChangeInterferenceOnS1().SetFirstParameter(0.);
148   Data->ChangeInterferenceOnS1().SetLastParameter(alpha);
149   Data->ChangeInterferenceOnS2().SetFirstParameter(0.);
150   Data->ChangeInterferenceOnS2().SetLastParameter(alpha);
151
152   return Standard_True;
153 }