0028643: Coding rules - eliminate GCC compiler warnings -Wmisleading-indentation
[occt.git] / src / SWDRAW / SWDRAW_ShapeCustom.cxx
1 // Created on: 1999-03-09
2 // Created by: Roman LYGIN
3 // Copyright (c) 1999-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 <BRep_Tool.hxx>
19 #include <DBRep.hxx>
20 #include <Draw.hxx>
21 #include <DrawTrSurf.hxx>
22 #include <Geom2d_BezierCurve.hxx>
23 #include <Geom2d_BSplineCurve.hxx>
24 #include <Geom2d_Curve.hxx>
25 #include <Geom2d_OffsetCurve.hxx>
26 #include <Geom2d_TrimmedCurve.hxx>
27 #include <Geom_BezierCurve.hxx>
28 #include <Geom_BezierSurface.hxx>
29 #include <Geom_BSplineCurve.hxx>
30 #include <Geom_BSplineSurface.hxx>
31 #include <Geom_Curve.hxx>
32 #include <Geom_OffsetCurve.hxx>
33 #include <Geom_OffsetSurface.hxx>
34 #include <Geom_Plane.hxx>
35 #include <Geom_RectangularTrimmedSurface.hxx>
36 #include <Geom_Surface.hxx>
37 #include <Geom_SweptSurface.hxx>
38 #include <Geom_TrimmedCurve.hxx>
39 #include <GeomAbs_Shape.hxx>
40 #include <ShapeCustom.hxx>
41 #include <ShapeCustom_BSplineRestriction.hxx>
42 #include <ShapeCustom_RestrictionParameters.hxx>
43 #include <ShapeFix.hxx>
44 #include <SWDRAW.hxx>
45 #include <SWDRAW_ShapeCustom.hxx>
46 #include <TColStd_Array1OfReal.hxx>
47 #include <TColStd_Array2OfInteger.hxx>
48 #include <TopAbs.hxx>
49 #include <TopExp_Explorer.hxx>
50 #include <TopLoc_Location.hxx>
51 #include <TopoDS.hxx>
52 #include <TopoDS_Edge.hxx>
53 #include <TopoDS_Face.hxx>
54 #include <TopoDS_Shape.hxx>
55
56 static Standard_Integer ContToInteger( const GeomAbs_Shape Cont)
57 {
58   Standard_Integer result =0;
59   switch(Cont) {
60     case GeomAbs_C0:
61     case GeomAbs_G1: result = 0; break;
62     case GeomAbs_C1:
63     case GeomAbs_G2: result = 1; break;
64     case GeomAbs_C2: result = 2; break;  
65     case GeomAbs_C3: result = 3; break;
66     default : result = 4; break;
67   }
68   return result;
69 }
70
71 static Standard_Integer directfaces(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
72 {
73   if (argc < 3) { di<<"Donner un nom de SHAPE + un nom de RESULTAT\n"; return 1 /* Error */; }
74   Standard_CString arg1 = argv[1];
75   Standard_CString arg2 = argv[2];
76   TopoDS_Shape Shape = DBRep::Get(arg2);
77   if (Shape.IsNull()) { di<<"Shape unknown : "<<arg2<<"\n"; return 1 /* Error */; }
78
79   TopoDS_Shape result = ShapeCustom::DirectFaces (Shape);
80   if (result.IsNull()) { di<<"NO RESULT\n"; return 1; }
81   else if (result == Shape) { di<<"No modif\n"; return 0; }
82   di<<"DirectFaces -> Result : "<<arg1<<"\n";
83   DBRep::Set (arg1,result);
84   return 0; // Done
85 }
86
87 static Standard_Integer ckeckKnots(const TColStd_Array1OfReal& theKnots,
88                                    Standard_Real theFirst,
89                                    Standard_Real theLast)
90 {
91   Standard_Integer i =1, nb = theKnots.Length();
92   Standard_Integer aNum =0;
93   for( ; i< nb; i++)
94   {
95     if(!aNum)
96     {
97       Standard_Real aDF = theKnots(i) - theFirst;
98       if(fabs(aDF) <= Precision::PConfusion() ||
99       aDF > Precision::PConfusion())
100       {
101         aNum++;
102         continue;
103       }
104     }
105     else 
106     {
107       Standard_Real aDL = theKnots(i) - theLast;
108       if(fabs(aDL) <= Precision::PConfusion() ||
109         aDL > Precision::PConfusion())
110         break;
111       aNum++;
112
113     }
114
115   }
116   return aNum;
117 }
118
119 static void expcurv2d(const Handle(Geom2d_Curve)& aCurve, TColStd_Array2OfInteger& Nb, const Standard_Integer Degree,
120                       const Standard_Integer MaxSeg,
121                       const Standard_Integer theCont,
122                       Standard_Real theFirst,
123                       Standard_Real theLast)   
124 {
125   if(aCurve.IsNull()) return;
126   
127   if (aCurve->IsKind(STANDARD_TYPE(Geom2d_TrimmedCurve))) {
128     Handle(Geom2d_TrimmedCurve) tmp = Handle(Geom2d_TrimmedCurve)::DownCast (aCurve);
129     Handle(Geom2d_Curve) BasCurve = tmp->BasisCurve();
130     expcurv2d(BasCurve,Nb,Degree,MaxSeg,theCont,theFirst,theLast);
131     return;
132   }
133   if (aCurve->IsKind(STANDARD_TYPE(Geom2d_OffsetCurve))) {
134     Handle(Geom2d_OffsetCurve) tmp = Handle(Geom2d_OffsetCurve)::DownCast (aCurve);
135     Handle(Geom2d_Curve) BasCurve = tmp->BasisCurve();
136     expcurv2d(BasCurve,Nb,Degree,MaxSeg,theCont,theFirst,theLast);
137     return;
138   }
139   Standard_Integer aCont = ContToInteger(aCurve->Continuity());
140   if(aCurve->IsKind(STANDARD_TYPE(Geom2d_BSplineCurve))) {
141     Handle(Geom2d_BSplineCurve) Bs = Handle(Geom2d_BSplineCurve)::DownCast(aCurve);
142     if(Bs->Degree() > Degree) Nb.ChangeValue(1,1)++; 
143     if(Bs->NbKnots() -1 > MaxSeg)  Nb.ChangeValue(1,2)++;
144     if(Bs->IsRational())  Nb.ChangeValue(1,3)++;
145     if(aCont < theCont && Bs->NbKnots() > 2)
146     {
147       TColStd_Array1OfReal aKnots(1,Bs->NbKnots());
148       Bs->Knots(aKnots);
149       Standard_Integer nbInt = ckeckKnots(aKnots,theFirst,theLast);
150       if(nbInt >1)
151         Nb.ChangeValue(1,4)++;
152     }
153     return;  
154   }
155   if (aCurve->IsKind(STANDARD_TYPE(Geom2d_BezierCurve))) {
156     Handle(Geom2d_BezierCurve) Bs = Handle(Geom2d_BezierCurve)::DownCast(aCurve); 
157     if(Bs->Degree() > Degree) Nb.ChangeValue(2,1)++;
158     if(Bs->IsRational())  Nb.ChangeValue(2,3)++;
159     if(aCont < theCont)
160     {
161
162       Nb.ChangeValue(2,4)++;
163     }
164     return;
165   }
166   else
167     Nb.ChangeValue(1,5)++;
168   return;
169 }
170
171 static void expcurv(const Handle(Geom_Curve)& aCurve, TColStd_Array2OfInteger& Nb, const Standard_Integer Degree,
172                     const Standard_Integer MaxSeg,
173                     const Standard_Integer theCont,
174                     Standard_Real theFirst,
175                     Standard_Real theLast)   
176 {
177   if(aCurve.IsNull()) return;
178   if (aCurve->IsKind(STANDARD_TYPE(Geom_TrimmedCurve))) {
179     Handle(Geom_TrimmedCurve) tmp = Handle(Geom_TrimmedCurve)::DownCast (aCurve);
180     Handle(Geom_Curve) BasCurve = tmp->BasisCurve();
181     expcurv(BasCurve,Nb,Degree,MaxSeg,theCont,theFirst,theLast);
182     return;
183   }
184   if (aCurve->IsKind(STANDARD_TYPE(Geom_OffsetCurve))) {
185     Handle(Geom_OffsetCurve) tmp = Handle(Geom_OffsetCurve)::DownCast (aCurve);
186     Handle(Geom_Curve) BasCurve = tmp->BasisCurve();
187     expcurv(BasCurve,Nb,Degree,MaxSeg,theCont,theFirst,theLast);
188     return;
189   }
190
191   Standard_Integer aCont = ContToInteger(aCurve->Continuity()); 
192   if(aCurve->IsKind(STANDARD_TYPE(Geom_BSplineCurve))) {
193     Handle(Geom_BSplineCurve) Bs = Handle(Geom_BSplineCurve)::DownCast(aCurve);
194     if(Bs->Degree() > Degree) Nb.ChangeValue(1,1)++; 
195     if(Bs->NbKnots() -1 > MaxSeg)  Nb.ChangeValue(1,2)++;
196     if(Bs->IsRational())  Nb.ChangeValue(1,3)++;
197     if(aCont < theCont && Bs->NbKnots() > 2) 
198     {
199       TColStd_Array1OfReal aKnots(1,Bs->NbKnots());
200       Bs->Knots(aKnots);
201       Standard_Integer nbInt = ckeckKnots(aKnots,theFirst,theLast);
202       if(nbInt >1)
203         Nb.ChangeValue(1,4)++;
204     }
205     return;  
206   }
207   if (aCurve->IsKind(STANDARD_TYPE(Geom_BezierCurve))) {
208     Handle(Geom_BezierCurve) Bs = Handle(Geom_BezierCurve)::DownCast(aCurve); 
209     if(Bs->Degree() > Degree) Nb.ChangeValue(2,1)++;
210     if(Bs->IsRational())  Nb.ChangeValue(2,3)++;
211     if(aCont < theCont) 
212       Nb.ChangeValue(2,4)++;
213     return;
214   }
215   else
216     Nb.ChangeValue(1,5)++;
217   return;
218 }
219 static void expsurf(const Handle(Geom_Surface)& aSurface, TColStd_Array2OfInteger& NbSurf, const Standard_Integer Degree,
220                     const Standard_Integer MaxSeg,
221                     const Standard_Integer theCont)   
222 {
223   if (aSurface.IsNull()) return;
224   if (aSurface->IsKind(STANDARD_TYPE(Geom_SweptSurface))) {      
225     Handle(Geom_SweptSurface) aSurf = Handle(Geom_SweptSurface)::DownCast(aSurface);
226     Handle(Geom_Curve) BasCurve = aSurf->BasisCurve();
227     expcurv(BasCurve,NbSurf,Degree,MaxSeg,theCont,BasCurve->FirstParameter(),
228       BasCurve->LastParameter()); 
229     return;  
230   }
231   if (aSurface->IsKind(STANDARD_TYPE(Geom_RectangularTrimmedSurface))) {
232     Handle(Geom_RectangularTrimmedSurface) aSurf = Handle(Geom_RectangularTrimmedSurface)::DownCast(aSurface);
233     Handle(Geom_Surface) theSurf = aSurf->BasisSurface();
234     expsurf(theSurf,NbSurf,Degree,MaxSeg,theCont);
235     return;
236   }
237   if(aSurface->IsKind(STANDARD_TYPE(Geom_OffsetSurface))) { 
238     Handle(Geom_OffsetSurface) aSurf = Handle(Geom_OffsetSurface)::DownCast(aSurface);
239     Handle(Geom_Surface) theSurf = aSurf->BasisSurface();
240     expsurf(theSurf,NbSurf,Degree,MaxSeg,theCont);
241     return;
242
243   }
244   Standard_Integer aCont = ContToInteger(aSurface->Continuity());
245
246   if (aSurface->IsKind(STANDARD_TYPE(Geom_BSplineSurface))) {
247     Handle(Geom_BSplineSurface) BS = Handle(Geom_BSplineSurface)::DownCast(aSurface); 
248     if(BS->UDegree() > Degree || BS->VDegree() > Degree) NbSurf.ChangeValue(1,1)++;
249     if(((BS->NbUKnots()-1) * (BS->NbVKnots()-1)) >MaxSeg) NbSurf.ChangeValue(1,2)++;
250     if(BS->IsURational() || BS->IsVRational()) NbSurf.ChangeValue(1,3)++;
251     if(aCont < theCont)
252       NbSurf.ChangeValue(1,4)++;
253     return;
254   }
255   if (aSurface->IsKind(STANDARD_TYPE(Geom_BezierSurface))) { 
256     Handle(Geom_BezierSurface) BS = Handle(Geom_BezierSurface)::DownCast(aSurface);
257     if(BS->UDegree() > Degree || BS->VDegree() > Degree) NbSurf.ChangeValue(2,1)++;
258     if(BS->IsURational() || BS->IsVRational()) NbSurf.ChangeValue(2,3)++;
259     if(aCont < theCont)
260       NbSurf.ChangeValue(2,4)++;
261     return;
262   }
263   if(aSurface->IsKind(STANDARD_TYPE(Geom_Plane)))
264   {
265     NbSurf.ChangeValue(1,5)++;
266   }
267   else
268     NbSurf.ChangeValue(2,5)++;
269 }
270
271 static Standard_Integer expshape(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
272 {
273   if (argc <4) {
274     di << "Incorrect number of arguments. Must be 3\n";
275     return 1 /* Error */;
276   }
277   Standard_CString arg2 = argv[1];
278   TopoDS_Shape Shape = DBRep::Get(arg2);
279   Standard_Integer Degree = Draw::Atoi(argv[2]);
280   Standard_Integer MaxSeg = Draw::Atoi(argv[3]);
281   GeomAbs_Shape aCont3 = GeomAbs_C0;
282   Standard_Integer k = 4;
283   if( argc > k)
284   {
285     if(strcmp(argv[k],"C0") == 0)  aCont3 = GeomAbs_C0;
286     else if (strcmp(argv[k],"C1") == 0)  aCont3 = GeomAbs_C1;
287     else if (strcmp(argv[k],"C2") == 0)  aCont3 = GeomAbs_C2;
288     else if (strcmp(argv[k],"C3") == 0)  aCont3 = GeomAbs_C3;
289     else if (strcmp(argv[k],"CN") == 0)  aCont3 = GeomAbs_CN;
290     else if (strcmp(argv[k],"G1") == 0)  aCont3 = GeomAbs_C0;
291     else if (strcmp(argv[k],"G2") == 0)  aCont3 = GeomAbs_C1;
292     else {
293       di<< "Invalid argument Cont3e\n"; return 1; }
294
295
296   }
297
298   Standard_Integer aCont = ContToInteger(aCont3);
299   TColStd_Array2OfInteger NbSurf(1,2,1,5);
300   TColStd_Array2OfInteger NbCurv(1,2,1,5);
301   TColStd_Array2OfInteger NbCurv2d(1,2,1,5);
302   Standard_Integer nbSeam =0;
303   NbSurf.Init(0);
304   NbCurv.Init(0);
305   NbCurv2d.Init(0);
306   if (Shape.IsNull()) { di << "Shape unknown: " << arg2 << "\n"; return 1 /* Error */; }
307   TopExp_Explorer Ex;
308   
309   Standard_Integer nbF =1;
310   for (Ex.Init(Shape,TopAbs_FACE); Ex.More(); Ex.Next(),nbF++) {
311     TopoDS_Face F = TopoDS::Face(Ex.Current());
312     TopLoc_Location L;
313     Handle(Geom_Surface) aSurface = BRep_Tool::Surface(F,L);
314     expsurf(aSurface,NbSurf,Degree,MaxSeg,aCont);
315     TopExp_Explorer exp;
316
317     Standard_Integer nbE = 1;
318     for(exp.Init(F,TopAbs_EDGE); exp.More(); exp.Next(),nbE++) {
319       TopoDS_Edge E = TopoDS::Edge(exp.Current());
320       if(BRep_Tool::IsClosed(E,F))
321         nbSeam++;
322       Standard_Real First,Last;
323       Handle(Geom_Curve) aCurve = BRep_Tool::Curve(E,L,First, Last);
324       expcurv(aCurve,NbCurv,Degree,MaxSeg,aCont,First, Last);
325       Standard_Real First2d,Last2d;
326       Handle(Geom2d_Curve) aCurve2d = BRep_Tool::CurveOnSurface(E,F,First2d,Last2d);
327       expcurv2d(aCurve2d,NbCurv2d,Degree,MaxSeg,aCont,First2d,Last2d);
328     }
329   }
330   di<<"Number of seam edges - "<<nbSeam<<"\n";
331   //if(NbSurf.Value(1,1) !=0)
332   di<< "Number of BSpline surfaces with degree more then "<< Degree << " - " <<NbSurf.Value(1,1)<<"\n";
333   //if(NbSurf.Value(1,2) !=0)
334   di<< "Number of BSpline surfaces with number of spans more then "<< MaxSeg << " - " <<NbSurf.Value(1,2)<<"\n";
335   //if(NbSurf.Value(1,3) !=0)
336   di<< "Number of Rational BSpline surfaces " << NbSurf.Value(1,3)<<"\n";
337   //if(NbSurf.Value(1,4) !=0)
338   di<< "Number of BSpline surfaces with continuity less than specified - " << NbSurf.Value(1,4)<<"\n";
339   //if(NbSurf.Value(2,1) !=0)
340   di<< "Number of Bezier surfaces with degree more then "<< Degree << " - " <<NbSurf.Value(2,1)<<"\n";
341   //if(NbSurf.Value(2,3) !=0)
342   di<< "Number of Rational Bezier surfaces  - " <<NbSurf.Value(2,3)<<"\n";
343   //if(NbSurf.Value(2,4) !=0)
344   di<< "Number of Bezier surfaces with continuity less than specified - " << NbSurf.Value(2,4)<<"\n";
345
346   //if(NbSurf.Value(1,5) !=0)
347   di<< "Number of Planes - " << NbSurf.Value(1,5)<<"\n";
348   //if(NbSurf.Value(2,5) !=0)
349   di<< "Number of other surfaces - " << NbSurf.Value(2,5)<<"\n";
350
351   //if(NbCurv.Value(1,1) !=0)
352   di<< "Number of BSpline curves with degree more then "<< Degree << " - " <<NbCurv.Value(1,1)<<"\n";
353   //if(NbCurv.Value(1,2) !=0)
354   di<< "Number of BSpline curves with number of spans more then - "<< MaxSeg << " - " <<NbCurv.Value(1,2)<<"\n";
355   //if(NbCurv.Value(1,3) !=0)
356   di<< "Number of Rational BSpline curves " << NbCurv.Value(1,3)<<"\n";
357   //if(NbCurv.Value(1,4) !=0)
358   di<< "Number of  BSpline curves with less continuity - " << NbCurv.Value(1,4)<<"\n";
359   //if(NbCurv.Value(2,1) !=0)
360   di<< "Number of Bezier curves with degree more then - "<< Degree << " - " <<NbCurv.Value(2,1)<<"\n";
361   //if(NbCurv.Value(2,3) !=0)
362   di<< "Number of Rational Bezier curves  - " <<NbCurv.Value(2,3)<<"\n"; 
363   //if(NbCurv.Value(2,4) !=0)
364   di<< "Number of  Bezier curves with less continuity - " << NbCurv.Value(2,4)<<"\n";
365
366   //if(NbCurv.Value(1,5) !=0)
367   di<< "Number of  other curves - " << NbCurv.Value(1,5)<<"\n";
368
369   //if(NbCurv2d.Value(1,1) !=0)
370   di<< "Number of BSpline pcurves with degree more then - "<< Degree << " - " <<NbCurv2d.Value(1,1)<<"\n";
371   //if(NbCurv2d.Value(1,2) !=0)
372   di<< "Number of BSpline pcurves with number of spans more then "<< MaxSeg << " - " <<NbCurv2d.Value(1,2)<<"\n";
373   //if(NbCurv2d.Value(1,3) !=0)
374   di<< "Number of Rational BSpline pcurves - " << NbCurv2d.Value(1,3)<<"\n";
375   // if(NbCurv2d.Value(1,4) !=0)
376   di<< "Number of  BSpline pcurves with less continuity - " << NbCurv2d.Value(1,4)<<"\n";
377   //if(NbCurv2d.Value(2,1) !=0)
378   di<< "Number of Bezier pcurves with degree more then "<< Degree <<" -  - " <<NbCurv2d.Value(2,1)<<"\n";
379   //if(NbCurv2d.Value(2,3) !=0)
380   di<< "Number of Rational Bezier pcurves  - " <<NbCurv2d.Value(2,3)<<"\n";
381   //if(NbCurv2d.Value(2,4) !=0)
382   di<< "Number of  Bezier pcurves with less continuity - " << NbCurv2d.Value(2,4)<<"\n";
383
384   //if(NbCurv2d.Value(1,5) !=0)
385   di<< "Number of  other pcurves - " << NbCurv2d.Value(1,5)<<"\n";
386   return 0;
387 }
388
389 static Standard_Integer scaleshape(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
390 {
391   if (argc !=4) {
392     di << "Incorrect number of arguments. Must be 4\n";
393     return 1 /* Error */;
394   }
395   Standard_CString arg2 = argv[2];
396   TopoDS_Shape Shape = DBRep::Get(arg2);
397   if (Shape.IsNull()) { di << "Shape unknown: " << arg2 << "\n"; return 1 /* Error */; }
398
399   TopoDS_Shape result = ShapeCustom::ScaleShape (Shape, Draw::Atof(argv[3]));
400   if (result.IsNull()) { di << "NO RESULT\n"; return 1; }
401   else if (result == Shape) { di << "NO MODIFICATIONS\n"; return 0; }
402   DBRep::Set (argv[1],result);
403   return 0;
404 }
405 static Standard_Integer BSplRes(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
406 {
407   if (argc < 11) {
408     di << "Incorrect number of arguments. Must be 10\n";
409     return 1 /* Error */;
410   }
411   Standard_CString arg2 = argv[2];
412   TopoDS_Shape Shape = DBRep::Get(arg2);
413   if (Shape.IsNull()) { di << "Shape unknown: " << arg2 << "\n"; return 1 /* Error */; }
414   GeomAbs_Shape aCont3;
415   if(strcmp(argv[7],"C0") == 0)  aCont3 = GeomAbs_C0;
416   else if (strcmp(argv[7],"C1") == 0)  aCont3 = GeomAbs_C1;
417   else if (strcmp(argv[7],"C2") == 0)  aCont3 = GeomAbs_C2;
418   else if (strcmp(argv[7],"C3") == 0)  aCont3 = GeomAbs_C3;
419   else if (strcmp(argv[7],"CN") == 0)  aCont3 = GeomAbs_CN;
420   else if (strcmp(argv[7],"G1") == 0)  aCont3 = GeomAbs_C0;
421   else if (strcmp(argv[7],"G2") == 0)  aCont3 = GeomAbs_C1;
422   else {
423     di<< "Invalid argument Cont3e\n"; return 1; }
424   GeomAbs_Shape aCont2;
425   if(strcmp(argv[8],"C0") == 0)  aCont2 = GeomAbs_C0;
426   else if (strcmp(argv[8],"C1") == 0)  aCont2 = GeomAbs_C1;
427   else if (strcmp(argv[8],"C2") == 0)  aCont2 = GeomAbs_C2;
428   else if (strcmp(argv[8],"C3") == 0)  aCont2 = GeomAbs_C3;
429   else if (strcmp(argv[8],"CN") == 0)  aCont2 = GeomAbs_CN;
430   else if (strcmp(argv[8],"G1") == 0)  aCont2 = GeomAbs_C0;
431   else if (strcmp(argv[8],"G2") == 0)  aCont2 = GeomAbs_C1;
432   else {
433     di<< "Invalid argument Cont3e\n"; return 1; }
434
435   Handle(ShapeCustom_RestrictionParameters) aParameters = new ShapeCustom_RestrictionParameters;
436   TopoDS_Shape result = ShapeCustom::BSplineRestriction (Shape, Draw::Atof(argv[3]), Draw::Atof(argv[4]), Draw::Atoi(argv[5]),Draw::Atoi(argv[6]),aCont3,aCont2, Draw::Atoi(argv[9]) != 0, Draw::Atoi(argv[10]) != 0,aParameters);
437   if (result.IsNull()) { di << "NO RESULT\n"; return 1; }
438   else if (result == Shape) { 
439     di << "NO MODIFICATIONS\n";
440     DBRep::Set (argv[1],result); return 0; 
441   }
442   ShapeFix::SameParameter(result,Standard_False);
443
444   DBRep::Set (argv[1],result);
445   return 0;
446 }
447
448 static Standard_Integer convtorevol(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
449 {
450   if (argc < 3) { di<<"Convert surfaces to revolution: convtorevol result shape\n"; return 1; }
451
452   Standard_CString arg1 = argv[1];
453   Standard_CString arg2 = argv[2];
454   TopoDS_Shape Shape = DBRep::Get(arg2);
455   if (Shape.IsNull()) { di<<"Shape unknown : "<<arg2<<"\n"; return 1; }
456
457   TopoDS_Shape result = ShapeCustom::ConvertToRevolution (Shape);
458   if (result.IsNull()) { di<<"NO RESULT\n"; return 1; }
459   else if (result == Shape) { di<<"No modif\n"; return 0; }
460   di<<"ConvertToRevolution -> Result : "<<arg1<<"\n";
461   DBRep::Set (arg1,result);
462   return 0; // Done
463 }  
464
465 //======================================================================
466 //function : InitCommands
467 //purpose  : 
468 //=======================================================================
469
470 void SWDRAW_ShapeCustom::InitCommands(Draw_Interpretor& theCommands) 
471 {
472   static Standard_Integer initactor = 0;
473   if (initactor)
474   {
475     return;
476   }
477   initactor = 1;
478
479   Standard_CString g = SWDRAW::GroupName();
480
481   theCommands.Add ("directfaces", "directfaces result shape", __FILE__, directfaces, g);
482   theCommands.Add ("expshape", "expshape shape maxdegree maxseg [min_continuity]", __FILE__, expshape, g);
483   theCommands.Add ("scaleshape", "scaleshape result shape scale", __FILE__, scaleshape, g);
484   theCommands.Add ("bsplres", "BSplineRestriction result shape tol3d tol2d reqdegree reqnbsegments continuity3d continuity2d PriorDeg RationalConvert", __FILE__,BSplRes, g);
485   theCommands.Add ("convtorevol", "convtorevol result shape", __FILE__, convtorevol, g);
486 }