0024053: Section between plane and sphere is not correct
[occt.git] / src / QABugs / QABugs_19.cxx
CommitLineData
b311480e 1// Created on: 2002-05-21
2// Created by: QA Admin
3// Copyright (c) 2002-2012 OPEN CASCADE SAS
4//
5// The content of this file is subject to the Open CASCADE Technology Public
6// License Version 6.5 (the "License"). You may not use the content of this file
7// except in compliance with the License. Please obtain a copy of the License
8// at http://www.opencascade.org and read it completely before using this file.
9//
10// The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
11// main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
12//
13// The Original Code and all software distributed under the License is
14// distributed on an "AS IS" basis, without warranty of any kind, and the
15// Initial Developer hereby disclaims all such warranties, including without
16// limitation, any warranties of merchantability, fitness for a particular
17// purpose or non-infringement. Please see the License for the specific terms
18// and conditions governing the rights and limitations under the License.
19
7fd59977 20
1cd84fee 21#include <QABugs.hxx>
7fd59977 22
23#include <Draw_Interpretor.hxx>
24#include <DBRep.hxx>
25#include <DrawTrSurf.hxx>
26#include <AIS_InteractiveContext.hxx>
27#include <ViewerTest.hxx>
28#include <AIS_Shape.hxx>
29#include <TopoDS_Shape.hxx>
30
31#include <gp_Pnt2d.hxx>
f6f03db9 32#include <gp_Ax1.hxx>
7fd59977 33#include <GCE2d_MakeSegment.hxx>
34#include <Geom2d_TrimmedCurve.hxx>
35#include <DrawTrSurf.hxx>
36
f6f03db9 37#include <Precision.hxx>
38
7fd59977 39#include <PCollection_HAsciiString.hxx>
40
c2ae831c 41#include <cstdio>
42#include <cmath>
43#include <iostream>
44#include <OSD_PerfMeter.hxx>
45#include <OSD_Timer.hxx>
46#include <BRepPrimAPI_MakeBox.hxx>
47#include <BRepPrimAPI_MakeSphere.hxx>
48#include <BRepAlgo_Cut.hxx>
49
d2786814 50
7fd59977 51static Standard_Integer OCC230 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
52{
53 if ( argc != 4) {
54 di << "ERROR OCC230: Usage : " << argv[0] << " TrimmedCurve Pnt2d Pnt2d" << "\n";
55 return 1;
56 }
57
58 gp_Pnt2d P1, P2;
59 if ( !DrawTrSurf::GetPoint2d(argv[2],P1)) {
60 di << "ERROR OCC230: " << argv[2] << " is not Pnt2d" << "\n";
61 return 1;
62 }
63 if ( !DrawTrSurf::GetPoint2d(argv[3],P2)) {
64 di << "ERROR OCC230: " << argv[3] << " is not Pnt2d" << "\n";
65 return 1;
66 }
67
68 GCE2d_MakeSegment MakeSegment(P1,P2);
69 Handle(Geom2d_TrimmedCurve) TrimmedCurve = MakeSegment.Value();
70 DrawTrSurf::Set(argv[1], TrimmedCurve);
71 return 0;
72}
73
74static Standard_Integer OCC142 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
75{
76 for(Standard_Integer i= 0;i <= 20;i++){
77 Handle(PCollection_HAsciiString) pstr = new PCollection_HAsciiString("TEST");
78 pstr->Clear();
79 }
80 di << "OCC142: OK" << "\n";
81 return 0;
82}
83
f6f03db9 84static Standard_Integer OCC23361 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
85{
86 gp_Pnt p(0, 0, 2);
87
88 gp_Trsf t1, t2;
89 t1.SetRotation(gp_Ax1(p, gp_Dir(0, 1, 0)), -0.49328285294022267);
90 t2.SetRotation(gp_Ax1(p, gp_Dir(0, 0, 1)), 0.87538474718473880);
91
92 gp_Trsf tComp = t2 * t1;
93
94 gp_Pnt p1(10, 3, 4);
95 gp_Pnt p2 = p1.Transformed(tComp);
96 gp_Pnt p3 = p1.Transformed(t1);
97 p3.Transform(t2);
98
99 // points must be equal
100 if ( ! p2.IsEqual(p3, Precision::Confusion()) )
101 di << "ERROR OCC23361: equivalent transformations does not produce equal points" << "\n";
102 else
103 di << "OCC23361: OK" << "\n";
104
105 return 0;
106}
107
c2ae831c 108static Standard_Integer OCC23237 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char** /*argv*/)
109{
110 OSD_PerfMeter aPM("TestMeter",0);
111 OSD_Timer aTM;
112
113 // run some operation in cycle for about 2 seconds to have good values of times to compare
114 int count = 0;
115 printf("OSD_PerfMeter test.\nRunning Boolean operation on solids in loop.\n");
116 for (; aTM.ElapsedTime() < 2.; count++)
117 {
118 aPM.Start();
119 aTM.Start();
120
121 // do some operation that will take considerable time compared with time or starting / stopping timers
122 BRepPrimAPI_MakeBox aBox (10., 10., 10.);
123 BRepPrimAPI_MakeSphere aSphere (10.);
124 BRepAlgo_Cut (aBox.Shape(), aSphere.Shape());
125
126 aTM.Stop();
127 aPM.Stop();
128 }
129
130 int aNbEnters = 0;
131 Standard_Real aPerfMeter_CPUtime = 0., aTimer_ElapsedTime = aTM.ElapsedTime();
132
133 perf_get_meter("TestMeter", &aNbEnters, &aPerfMeter_CPUtime);
134
135 Standard_Real aTimeDiff = (fabs(aTimer_ElapsedTime - aPerfMeter_CPUtime) / aTimer_ElapsedTime);
136
137 printf("\nMeasurement results (%d cycles):\n", count);
138 printf("\nOSD_PerfMeter CPU time: %lf\nOSD_Timer elapsed time: %lf\n", aPerfMeter_CPUtime, aTimer_ElapsedTime);
139 printf("Time delta is: %.3lf %%\n", aTimeDiff * 100);
140
141 if (aTimeDiff > 0.2)
142 di << "OCC23237: Error: too much difference between CPU and elapsed times";
143 else if (aNbEnters != count)
144 di << "OCC23237: Error: counter reported by PerfMeter (" << aNbEnters << ") does not correspond to actual number of cycles";
145 else
146 di << "OCC23237: OK";
147
148 return 0;
149}
150
1365140b
RL
151#ifdef HAVE_TBB
152
153#include <Standard_Atomic.hxx>
154#include <tbb/blocked_range.h>
155#include <tbb/parallel_for.h>
156
157class IncrementerDecrementer
158{
159public:
160 IncrementerDecrementer (Standard_Integer* theVal, Standard_Boolean thePositive) : myVal (theVal), myPositive (thePositive)
161 {}
162 void operator() (const tbb::blocked_range<size_t>& r) const
163 {
164 if (myPositive)
165 for (size_t i = r.begin(); i != r.end(); ++i)
166 Standard_Atomic_Increment (myVal);
167 else
168 for (size_t i = r.begin(); i != r.end(); ++i)
169 Standard_Atomic_Decrement (myVal);
170 }
171private:
172 Standard_Integer* myVal;
173 Standard_Boolean myPositive;
174};
175#endif
176
177#define QCOMPARE(val1, val2) \
178 di << "Checking " #val1 " == " #val2 << \
179 ((val1) == (val2) ? ": OK\n" : ": Error\n")
180
181#ifdef HAVE_TBB
182static Standard_Integer OCC22980 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
183{
184 int aSum = 0;
185
186 //check returned value
187 QCOMPARE (Standard_Atomic_Decrement (&aSum), -1);
188 QCOMPARE (Standard_Atomic_Increment (&aSum), 0);
189 QCOMPARE (Standard_Atomic_Increment (&aSum), 1);
190 QCOMPARE (Standard_Atomic_Increment (&aSum), 2);
191// QCOMPARE (Standard_Atomic_DecrementTest (&aSum), 0);
192// QCOMPARE (Standard_Atomic_DecrementTest (&aSum), 1);
193
194 //check atomicity
195 aSum = 0;
196 const int N = 1 << 24; //big enough to ensure concurrency
197
198 //increment
199 tbb::parallel_for (tbb::blocked_range<size_t> (0, N), IncrementerDecrementer (&aSum, true));
200 QCOMPARE (aSum, N);
201
202 //decrement
203 tbb::parallel_for (tbb::blocked_range<size_t> (0, N), IncrementerDecrementer (&aSum, false));
204 QCOMPARE (aSum, 0);
205
206 return 0;
207}
208
209#else /* HAVE_TBB */
210
211static Standard_Integer OCC22980 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char **argv)
212{
213 di << "Test skipped: command " << argv[0] << " requires TBB library\n";
214 return 0;
215}
216
217#endif /* HAVE_TBB */
218
e6aad0ee
RL
219#include <TDocStd_Application.hxx>
220#include <XCAFApp_Application.hxx>
221#include <TDocStd_Document.hxx>
222#include <XCAFDoc_ShapeTool.hxx>
223#include <XCAFDoc_DocumentTool.hxx>
224#include <TDF_Label.hxx>
225#include <TDataStd_Name.hxx>
226
227static Standard_Integer OCC23595 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char **argv)
228{
229 const Handle(TDocStd_Application)& anApp = XCAFApp_Application::GetApplication();
230 Handle(TDocStd_Document) aDoc;
231 anApp->NewDocument ("XmlXCAF", aDoc);
232 QCOMPARE (!aDoc.IsNull(), Standard_True);
233
234 Handle(XCAFDoc_ShapeTool) aShTool = XCAFDoc_DocumentTool::ShapeTool (aDoc->Main());
235
236 //check default value
237 Standard_Boolean aValue = XCAFDoc_ShapeTool::AutoNaming();
238 QCOMPARE (aValue, Standard_True);
239
240 //true
241 XCAFDoc_ShapeTool::SetAutoNaming (Standard_True);
242 TopoDS_Shape aShape = BRepPrimAPI_MakeBox (100., 200., 300.).Shape();
243 TDF_Label aLabel = aShTool->AddShape (aShape);
244 Handle(TDataStd_Name) anAttr;
245 QCOMPARE (aLabel.FindAttribute (TDataStd_Name::GetID(), anAttr), Standard_True);
246
247 //false
248 XCAFDoc_ShapeTool::SetAutoNaming (Standard_False);
249 aShape = BRepPrimAPI_MakeBox (300., 200., 100.).Shape();
250 aLabel = aShTool->AddShape (aShape);
251 QCOMPARE (!aLabel.FindAttribute (TDataStd_Name::GetID(), anAttr), Standard_True);
252
253 //restore
254 XCAFDoc_ShapeTool::SetAutoNaming (aValue);
255
256 return 0;
257}
258
5d1833ef 259#include <ExprIntrp_GenExp.hxx>
260Standard_Integer OCC22611 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
261{
262
263 if (argc != 3) {
264 di << "Usage : " << argv[0] << " string nb\n";
265 return 1;
266 }
267
268 TCollection_AsciiString aToken = argv[1];
269 Standard_Integer aNb = atoi(argv[2]);
270
271 Handle(ExprIntrp_GenExp) aGen = ExprIntrp_GenExp::Create();
272 for (Standard_Integer i=0; i < aNb; i++)
273 {
274 aGen->Process(aToken);
275 Handle(Expr_GeneralExpression) aExpr = aGen->Expression();
276 }
277
278 return 0;
279}
280
e187cc0a 281Standard_Integer OCC22595 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
282{
283 gp_Mat M0;
284 di << "M0 = "
285 << " {" << M0(1,1) << "} {" << M0(1,2) << "} {" << M0(1,3) <<"}"
286 << " {" << M0(2,1) << "} {" << M0(2,2) << "} {" << M0(2,3) <<"}"
287 << " {" << M0(1,1) << "} {" << M0(1,2) << "} {" << M0(1,3) <<"}";
288 return 0;
289}
290
bead40f2 291#include <TopoDS_Face.hxx>
292#include <TopoDS_Face.hxx>
293#include <TopoDS.hxx>
294#include <BRepBuilderAPI_Transform.hxx>
295#include <BRepExtrema_DistShapeShape.hxx>
296#include <BRepTools.hxx>
297
298Standard_Boolean static OCC23774Test(const TopoDS_Face& grossPlateFace, const TopoDS_Shape& originalWire, Draw_Interpretor& di)
299{
300 BRepExtrema_DistShapeShape distShapeShape(grossPlateFace,originalWire,Extrema_ExtFlag_MIN);
301 if(!distShapeShape.IsDone()) {
302 di <<"Distance ShapeShape is Not Done\n";
303 return Standard_False;
304 }
305
306 if(distShapeShape.Value() > 0.01) {
307 di << "Wrong Dist = " <<distShapeShape.Value() << "\n";
308 return Standard_False;
309 } else
310 di << "Dist0 = " <<distShapeShape.Value() <<"\n";
311
312 //////////////////////////////////////////////////////////////////////////
313 /// First Flip Y
314 const gp_Pnt2d axis1P1(1474.8199035519228,1249.9995745636970);
315 const gp_Pnt2d axis1P2(1474.8199035519228,1250.9995745636970);
316
317 gp_Vec2d mirrorVector1(axis1P1,axis1P2);
318
319 gp_Trsf2d mirror1;
320 mirror1.SetMirror(gp_Ax2d(axis1P1,mirrorVector1));
321
322 BRepBuilderAPI_Transform transformer1(mirror1);
323 transformer1.Perform(originalWire);
324 if(!transformer1.IsDone()) {
325 di << "Not Done1 " << "\n";
326 return Standard_False;
327 }
328 const TopoDS_Shape& step1ModifiedShape = transformer1.ModifiedShape(originalWire);
329
330 BRepExtrema_DistShapeShape distShapeShape1(grossPlateFace,step1ModifiedShape,Extrema_ExtFlag_MIN);
331 if(!distShapeShape1.IsDone())
332 return Standard_False;
333 if(distShapeShape1.Value() > 0.01) {
334 di << "Dist = " <<distShapeShape1.Value() <<"\n";
335 return Standard_False;
336 } else
337 di << "Dist1 = " <<distShapeShape1.Value() <<"\n";
338
339 //////////////////////////////////////////////////////////////////////////
340 /// Second flip Y
341 transformer1.Perform(step1ModifiedShape);
342 if(!transformer1.IsDone()) {
343 di << "Not Done1 \n";
344 return Standard_False;
345 }
346 const TopoDS_Shape& step2ModifiedShape = transformer1.ModifiedShape(step1ModifiedShape);
bead40f2 347
348 //This is identity matrix for values but for type is gp_Rotation ?!
349 gp_Trsf2d mirror11 = mirror1;
350 mirror11.PreMultiply(mirror1);
351
352 BRepExtrema_DistShapeShape distShapeShape2(grossPlateFace,step2ModifiedShape);//,Extrema_ExtFlag_MIN);
353 if(!distShapeShape2.IsDone())
354 return Standard_False;
355
356 //This last test case give error (the value is 1008.8822038689706)
357 if(distShapeShape2.Value() > 0.01) {
358 di << "Wrong Dist2 = " <<distShapeShape2.Value() <<"\n";
359 Standard_Integer N = distShapeShape2.NbSolution();
360 di << "Nb = " <<N <<"\n";
361 for (Standard_Integer i=1;i <= N;i++)
362 di <<"Sol(" <<i<<") = " <<distShapeShape2.PointOnShape1(i).Distance(distShapeShape2.PointOnShape2(i)) <<"\n";
363 return Standard_False;
364 }
365 di << "Distance2 = " <<distShapeShape2.Value() <<"\n";
366
367 return Standard_True;
368}
369static Standard_Integer OCC23774(Draw_Interpretor& di, Standard_Integer n, const char** a)
370{
371
372 if (n != 3) {
373 di <<"OCC23774: invalid number of input parameters\n";
374 return 1;
375 }
376
377 const char *ns1 = (a[1]), *ns2 = (a[2]);
378 TopoDS_Shape S1(DBRep::Get(ns1)), S2(DBRep::Get(ns2));
379 if (S1.IsNull() || S2.IsNull()) {
380 di <<"OCC23774: Null input shapes\n";
381 return 1;
382 }
383 const TopoDS_Face& aFace = TopoDS::Face(S1);
384 if(!OCC23774Test(aFace, S2, di))
385 di << "Something is wrong\n";
386
387 return 0;
388}
389
1ef32e96
RL
390#include <GeomConvert_ApproxSurface.hxx>
391#include <Geom_BSplineSurface.hxx>
392#include <Draw.hxx>
393#include <OSD_Thread.hxx>
394static void GeomConvertTest (Draw_Interpretor& di, Standard_Integer theTargetNbUPoles, Standard_CString theFileName)
395{
396 Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(theFileName);
397 GeomConvert_ApproxSurface aGAS (aSurf, 1e-4, GeomAbs_C1, GeomAbs_C1, 9, 9, 100, 1);
398 if (!aGAS.IsDone()) {
399 di << "ApproxSurface is not done!" << "\n";
400 return;
401 }
402 const Handle(Geom_BSplineSurface)& aBSurf = aGAS.Surface();
403 if (aBSurf.IsNull()) {
404 di << "BSplineSurface is not created!" << "\n";
405 return;
406 }
407 di << "Number of UPoles:" << aBSurf->NbUPoles() << "\n";
408 QCOMPARE (aBSurf->NbUPoles(), theTargetNbUPoles);
409}
410
411struct aData {
412 Draw_Interpretor* di;
413 Standard_Integer nbupoles;
414 Standard_CString filename;
415};
416
417Standard_EXPORT Standard_Address convert(Standard_Address data)
418{
419 aData* info = (aData*) data;
420 GeomConvertTest(*(info->di),info->nbupoles,info->filename);
421 return NULL;
422}
423
424static Standard_Integer OCC23952sweep (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
425{
426 if (argc != 3) {
427 di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
428 return 1;
429 }
430 struct aData aStorage;
431 aStorage.di = &di;
432 aStorage.nbupoles = Draw::Atoi(argv[1]);
433 aStorage.filename = argv[2];
434
435 OSD_Thread aThread1(convert);
436 aThread1.Run(&aStorage);
437 GeomConvertTest(di,aStorage.nbupoles,aStorage.filename);
85c44a05 438 cout << "result of thread: " << aThread1.Wait() << endl;
1ef32e96
RL
439
440 return 0;
441}
d2786814 442
1ef32e96
RL
443#include <GeomInt_IntSS.hxx>
444static void GeomIntSSTest (Draw_Interpretor& di, Standard_Integer theNbSol, Standard_CString theFileName1, Standard_CString theFileName2)
445{
446 Handle(Geom_Surface) aSurf1 = DrawTrSurf::GetSurface(theFileName1);
447 Handle(Geom_Surface) aSurf2 = DrawTrSurf::GetSurface(theFileName2);
448 GeomInt_IntSS anInter;
449 anInter.Perform(aSurf1, aSurf2, Precision::Confusion(), Standard_True);
450 if (!anInter.IsDone()) {
451 di << "An intersection is not done!" << "\n";
452 return;
453 }
454
455 di << "Number of Lines:" << anInter.NbLines() << "\n";
456 QCOMPARE (anInter.NbLines(), theNbSol);
457}
458
459struct aNewData {
460 Draw_Interpretor* di;
461 Standard_Integer nbsol;
462 Standard_CString filename1;
463 Standard_CString filename2;
464};
465Standard_EXPORT Standard_Address convert_inter(Standard_Address data)
466{
467 aNewData* info = (aNewData*) data;
468 GeomIntSSTest(*(info->di),info->nbsol,info->filename1,info->filename2);
469 return NULL;
470}
471
472static Standard_Integer OCC23952intersect (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
473{
474 if (argc != 4) {
475 di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
476 return 1;
477 }
478 struct aNewData aStorage;
479 aStorage.di = &di;
480 aStorage.nbsol = Draw::Atoi(argv[1]);
481 aStorage.filename1 = argv[2];
482 aStorage.filename2 = argv[3];
483
484 OSD_Thread aThread1(convert_inter);
485 aThread1.Run(&aStorage);
486 GeomIntSSTest(di,aStorage.nbsol,aStorage.filename1,aStorage.filename2);
85c44a05 487 cout << "result of thread: " << aThread1.Wait() << endl;
1ef32e96
RL
488
489 return 0;
490}
491
492#include <Geom_SurfaceOfRevolution.hxx>
06be28a4
RL
493static Standard_Integer OCC23683 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
494{
495 if (argc < 2) {
496 di<<"Usage: " << argv[0] << " invalid number of arguments"<<"\n";
497 return 1;
498 }
499
500 Standard_Integer ucontinuity = 1;
501 Standard_Integer vcontinuity = 1;
502 Standard_Boolean iscnu = false;
503 Standard_Boolean iscnv = false;
504
505 Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(argv[1]);
506
507 QCOMPARE (aSurf->IsCNu (ucontinuity), iscnu);
508 QCOMPARE (aSurf->IsCNv (vcontinuity), iscnv);
509
510 return 0;
511}
512
1939140c 513#include <gp_Ax1.hxx>
514#include <gp_Ax22d.hxx>
515#include <Geom_Plane.hxx>
516#include <Geom2d_Circle.hxx>
517#include <Geom2d_TrimmedCurve.hxx>
518#include <BRepBuilderAPI_MakeEdge.hxx>
519#include <BRepPrimAPI_MakeRevol.hxx>
520#include <Geom2d_OffsetCurve.hxx>
521
522static int test_offset(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
523{
524 // Check the command arguments
525 if ( argc != 1 )
526 {
527 di << "Error: " << argv[0] << " - invalid number of arguments" << "\n";
528 di << "Usage: type help " << argv[0] << "\n";
529 return 1; // TCL_ERROR
530 }
531
532 gp_Ax1 RotoAx( gp::Origin(), gp::DZ() );
533 gp_Ax22d Ax2( gp::Origin2d(), gp::DY2d(), gp::DX2d() );
534 Handle(Geom_Surface) Plane = new Geom_Plane( gp::YOZ() );
535
536 di << "<<<< Preparing sample surface of revolution based on trimmed curve >>>>" << "\n";
537 di << "-----------------------------------------------------------------------" << "\n";
538
539 Handle(Geom2d_Circle) C2d1 = new Geom2d_Circle(Ax2, 1.0);
540 Handle(Geom2d_TrimmedCurve) C2d1Trimmed = new Geom2d_TrimmedCurve(C2d1, 0.0, M_PI/2.0);
541 TopoDS_Edge E1 = BRepBuilderAPI_MakeEdge(C2d1Trimmed, Plane);
542
543 DBRep::Set("e1", E1);
544
545 BRepPrimAPI_MakeRevol aRevolBuilder1(E1, RotoAx);
546 TopoDS_Face F1 = TopoDS::Face( aRevolBuilder1.Shape() );
547
548 DBRep::Set("f1", F1);
549
550 di << "Result: f1" << "\n";
551
552 di << "<<<< Preparing sample surface of revolution based on offset curve >>>>" << "\n";
553 di << "-----------------------------------------------------------------------" << "\n";
554
555 Handle(Geom2d_OffsetCurve) C2d2Offset = new Geom2d_OffsetCurve(C2d1Trimmed, -0.5);
556 TopoDS_Edge E2 = BRepBuilderAPI_MakeEdge(C2d2Offset, Plane);
557
558 DBRep::Set("e2", E2);
559
560 BRepPrimAPI_MakeRevol aRevolBuilder2(E2, RotoAx);
561 TopoDS_Face F2 = TopoDS::Face( aRevolBuilder2.Shape() );
562
563 DBRep::Set("f2", F2);
564
565 di << "Result: f2" << "\n";
566
567 return 0;
568}
569
470ebb43 570#include <GeomAdaptor_Surface.hxx>
571#include <Draw.hxx>
572//=======================================================================
573//function : OCC23945
574//purpose :
575//=======================================================================
576
577static Standard_Integer OCC23945 (Draw_Interpretor& di,Standard_Integer n, const char** a)
578{
579 if (n < 5) return 1;
580
581 Handle(Geom_Surface) aS = DrawTrSurf::GetSurface(a[1]);
582 if (aS.IsNull()) return 1;
583
584 GeomAdaptor_Surface GS(aS);
585
586 Standard_Real U = Draw::Atof(a[2]);
587 Standard_Real V = Draw::Atof(a[3]);
588
589 Standard_Boolean DrawPoint = ( n%3 == 2);
590 if ( DrawPoint) n--;
591
592 gp_Pnt P;
593 if (n >= 13) {
594 gp_Vec DU,DV;
595 if (n >= 22) {
596 gp_Vec D2U,D2V,D2UV;
597 GS.D2(U,V,P,DU,DV,D2U,D2V,D2UV);
598 Draw::Set(a[13],D2U.X());
599 Draw::Set(a[14],D2U.Y());
600 Draw::Set(a[15],D2U.Z());
601 Draw::Set(a[16],D2V.X());
602 Draw::Set(a[17],D2V.Y());
603 Draw::Set(a[18],D2V.Z());
604 Draw::Set(a[19],D2UV.X());
605 Draw::Set(a[20],D2UV.Y());
606 Draw::Set(a[21],D2UV.Z());
607 }
608 else
609 GS.D1(U,V,P,DU,DV);
610
611 Draw::Set(a[7],DU.X());
612 Draw::Set(a[8],DU.Y());
613 Draw::Set(a[9],DU.Z());
614 Draw::Set(a[10],DV.X());
615 Draw::Set(a[11],DV.Y());
616 Draw::Set(a[12],DV.Z());
617 }
618 else
619 GS.D0(U,V,P);
620
621 if ( n > 6) {
622 Draw::Set(a[4],P.X());
623 Draw::Set(a[5],P.Y());
624 Draw::Set(a[6],P.Z());
625 }
626 if ( DrawPoint) {
627 DrawTrSurf::Set(a[n],P);
628 }
629
630 return 0;
631}
632
1cd84fee 633void QABugs::Commands_19(Draw_Interpretor& theCommands) {
1365140b 634 const char *group = "QABugs";
7fd59977 635
636 theCommands.Add ("OCC230", "OCC230 TrimmedCurve Pnt2d Pnt2d", __FILE__, OCC230, group);
637 theCommands.Add ("OCC142", "OCC142", __FILE__, OCC142, group);
f6f03db9 638 theCommands.Add ("OCC23361", "OCC23361", __FILE__, OCC23361, group);
1365140b
RL
639 theCommands.Add ("OCC23237", "OCC23237", __FILE__, OCC23237, group);
640 theCommands.Add ("OCC22980", "OCC22980", __FILE__, OCC22980, group);
e6aad0ee 641 theCommands.Add ("OCC23595", "OCC23595", __FILE__, OCC23595, group);
5d1833ef 642 theCommands.Add ("OCC22611", "OCC22611 string nb", __FILE__, OCC22611, group);
e187cc0a 643 theCommands.Add ("OCC22595", "OCC22595", __FILE__, OCC22595, group);
bead40f2 644 theCommands.Add ("OCC23774", "OCC23774 shape1 shape2", __FILE__, OCC23774, group);
06be28a4 645 theCommands.Add ("OCC23683", "OCC23683 shape", __FILE__, OCC23683, group);
1ef32e96
RL
646 theCommands.Add ("OCC23952sweep", "OCC23952sweep nbupoles shape", __FILE__, OCC23952sweep, group);
647 theCommands.Add ("OCC23952intersect", "OCC23952intersect nbsol shape1 shape2", __FILE__, OCC23952intersect, group);
1939140c 648 theCommands.Add ("test_offset", "test_offset", __FILE__, test_offset, group);
470ebb43 649 theCommands.Add("OCC23945", "OCC23945 surfname U V X Y Z [DUX DUY DUZ DVX DVY DVZ [D2UX D2UY D2UZ D2VX D2VY D2VZ D2UVX D2UVY D2UVZ]]", __FILE__, OCC23945,group);
7fd59977 650 return;
651}
a9926dd8 652
653