1 // Copyright (C) 2007-2008 CEA/DEN, EDF R&D, OPEN CASCADE
3 // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License.
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 // Lesser General Public License for more details.
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
22 // GEOM SKETCHER : basic sketcher
23 // File : Sketcher_Profile.cxx
24 // Author : Damien COQUERET
28 #include <Standard_Stream.hxx>
30 #include <Sketcher_Profile.hxx>
32 #include <TopoDS_Vertex.hxx>
33 #include <TopoDS_Face.hxx>
34 #include <BRepLib.hxx>
35 #include <BRepBuilderAPI_MakeVertex.hxx>
36 #include <BRepBuilderAPI_MakeEdge.hxx>
37 #include <BRepBuilderAPI_MakeWire.hxx>
38 #include <BRepBuilderAPI_MakeFace.hxx>
40 #include <GeomAPI.hxx>
41 #include <Geom2d_Line.hxx>
42 #include <Geom2d_Circle.hxx>
43 #include <Geom_Surface.hxx>
45 #include <Precision.hxx>
49 #include <TCollection_AsciiString.hxx>
50 #include <TColStd_Array1OfAsciiString.hxx>
52 #include "utilities.h"
54 //=======================================================================
56 // command to build a profile
57 //=======================================================================
58 Sketcher_Profile::Sketcher_Profile()
63 //=======================================================================
65 // command to build a profile
66 //=======================================================================
67 Sketcher_Profile::Sketcher_Profile(const char* aCmd)
69 enum {line, circle, point, none} move;
71 Standard_Integer i = 1;
72 Standard_Real x0, y0, x, y, dx, dy;
73 x0 = y0 = x = y = dy = 0;
76 Standard_Boolean first, stayfirst, face, close;
77 first = Standard_True;
78 stayfirst = face = close = Standard_False;
82 BRepBuilderAPI_MakeWire MW;
83 gp_Ax3 DummyHP(gp::XOY());
85 TopLoc_Location TheLocation;
86 Handle(Geom_Surface) Surface;
88 myOK = Standard_False;
90 //TCollection_AsciiString aCommand(CORBA::string_dup(aCmd));
91 TCollection_AsciiString aCommand ((char*)aCmd);
92 TCollection_AsciiString aToken = aCommand.Token(":", 1);
95 TColStd_Array1OfAsciiString aTab (0, aCommand.Length() - 1);
96 if ( aCommand.Length() )
98 while(aToken.Length() != 0) {
99 if(aCommand.Token(":", n + 1).Length() > 0)
100 aTab(n) = aCommand.Token(":", n + 1);
101 aToken = aCommand.Token(":", ++n);
105 if ( aTab.Length() && aTab(0).Length() )
107 Standard_Real length = 0, radius = 0, angle = 0;
111 TColStd_Array1OfAsciiString a (0, aTab(0).Length());
112 aToken = aTab(i).Token(" ", 1);
113 while (aToken.Length() != 0) {
114 if (aTab(i).Token(" ", n1 + 1).Length() > 0)
115 a(n1) = aTab(i).Token(" ", n1 + 1);
116 aToken = aTab(i).Token(" ", ++n1);
120 switch(a(0).Value(1))
124 if (n1 != 3) goto badargs;
126 MESSAGE("profile : The F instruction must precede all moves");
129 x0 = x = a(1).RealValue();
130 y0 = y = a(2).RealValue();
131 stayfirst = Standard_True;
136 if (n1 != 4) goto badargs;
137 P.SetLocation(gp_Pnt(a(1).RealValue(), a(2).RealValue(), a(3).RealValue()));
138 stayfirst = Standard_True;
143 if (n1 != 7) goto badargs;
144 gp_Vec vn(a(1).RealValue(), a(2).RealValue(), a(3).RealValue());
145 gp_Vec vx(a(4).RealValue(), a(5).RealValue(), a(6).RealValue());
146 if (vn.Magnitude() <= Precision::Confusion() || vx.Magnitude() <= Precision::Confusion()) {
147 MESSAGE("profile : null direction");
150 gp_Ax2 ax(P.Location(), vn, vx);
152 stayfirst = Standard_True;
157 if (n1 != 2) goto badargs;
158 length = a(1).RealValue();
167 if (n1 != 2) goto badargs;
168 length = a(1).RealValue();
177 if (n1 != 2) goto badargs;
178 length = a(1).RealValue();
179 if (Abs(length) > Precision::Confusion())
187 if (n1 != 3) goto badargs;
188 Standard_Real vx = a(1).RealValue();
189 Standard_Real vy = a(2).RealValue();
194 length = Sqrt(vx * vx + vy * vy);
195 if (length > Precision::Confusion()) {
206 if (n1 != 2) goto badargs;
207 angle = a(1).RealValue() * PI180;
213 Standard_Real c = Cos(angle);
214 Standard_Real s = Sin(angle);
215 Standard_Real t = c * dx - s * dy;
216 dy = s * dx + c * dy;
223 if (n1 != 3) goto badargs;
224 Standard_Real vx = a(1).RealValue();
225 Standard_Real vy = a(2).RealValue();
226 length = Sqrt(vx * vx + vy * vy);
227 if (length > Precision::Confusion()) {
237 if (n1 != 3) goto badargs;
238 radius = a(1).RealValue();
239 if (Abs(radius) > Precision::Confusion()) {
240 angle = a(2).RealValue() * PI180;
249 if (n1 != 2) goto badargs;
250 length = a(1).RealValue();
252 if (Abs(dx) < Precision::Confusion()) {
253 MESSAGE("profile : cannot intersect, arg "<<i-1);
256 length = (length - x) / dx;
258 else if (a(0) == "IY") {
259 if (Abs(dy) < Precision::Confusion()) {
260 MESSAGE("profile : cannot intersect, arg "<<i-1);
263 length = (length - y) / dy;
265 if (Abs(length) > Precision::Confusion())
274 close = Standard_True;
275 else if(a(0) == "WF") {
276 close = Standard_True;
277 face = Standard_True;
284 MESSAGE("profile : unknown code " << a(i));
299 Handle(Geom2d_Line) l = new Geom2d_Line(gp_Pnt2d(x,y),gp_Dir2d(dx,dy));
300 BRepBuilderAPI_MakeEdge ME (GeomAPI::To3d(l,P),0,length);
310 Standard_Boolean sense = Standard_True;
317 gp_Ax2d ax(gp_Pnt2d(x-radius*dy,y+radius*dx),gp_Dir2d(dy,-dx));
322 Handle(Geom2d_Circle) c = new Geom2d_Circle(ax,radius,sense);
323 BRepBuilderAPI_MakeEdge ME (GeomAPI::To3d(c,P),0,angle);
338 MP = BRepBuilderAPI_MakeVertex(gp_Pnt(x, y, 0.0));
350 stayfirst = Standard_False;
352 if(!(dx == 0 && dy == 0))
353 myLastDir.SetCoord(dx, dy, 0.0);
361 if ((i == n) && close) {
362 // the closing segment
365 length = Sqrt(dx * dx + dy * dy);
367 if (length > Precision::Confusion()) {
375 // get the result, face or wire
378 } else if (move == point) {
384 BRepBuilderAPI_MakeFace MF (P, MW.Wire());
396 if(!TheLocation.IsIdentity())
404 MESSAGE("profile : bad number of arguments");