src/QABugs/QABugs_9.cxx

   1 // Created on: 2002-03-18
   2 // Created by: QA Admin
   3 // Copyright (c) 2002-2014 OPEN CASCADE SAS
   4 //
   5 // This file is part of Open CASCADE Technology software library.
   6 //
   7 // This library is free software; you can redistribute it and/or modify it under
   8 // the terms of the GNU Lesser General Public License version 2.1 as published
   9 // by the Free Software Foundation, with special exception defined in the file
  10 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
  11 // distribution for complete text of the license and disclaimer of any warranty.
  12 //
  13 // Alternatively, this file may be used under the terms of Open CASCADE
  14 // commercial license or contractual agreement.
  15
  16 #include <QABugs.hxx>
  17
  18 #include <Draw.hxx>
  19 #include <Draw_Interpretor.hxx>
  20 #include <DBRep.hxx>
  21 #include <DrawTrSurf.hxx>
  22 #include <AIS_InteractiveContext.hxx>
  23 #include <ViewerTest.hxx>
  24 #include <AIS_Shape.hxx>
  25 #include <TopoDS_Shape.hxx>
  26
  27 #include <gp_Ax2.hxx>
  28 #include <gp_Pnt.hxx>
  29 #include <Geom_RectangularTrimmedSurface.hxx>
  30 #include <BRepBuilderAPI_MakeFace.hxx>
  31 #include <GC_MakeTrimmedCone.hxx>
  32
  33 static Standard_Integer BUC60857 (Draw_Interpretor& di, Standard_Integer /*argc*/,const char ** argv)
  34 {
  35   gp_Ax2  Cone_Ax;
  36   double R1=8, R2=16;
  37   gp_Pnt P0(0,0,0), P1(0,0,20), P2(0,0,45);
  38
  39   Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
  40   if(aContext.IsNull()) {
  41     di << "Use vinit command before " << argv[0] << "\n";
  42     return 1;
  43   }
  44
  45   Handle(Geom_RectangularTrimmedSurface) S = GC_MakeTrimmedCone (P1, P2, R1, R2).Value();
  46   TopoDS_Shape myshape = BRepBuilderAPI_MakeFace(S, Precision::Confusion()).Shape();
  47   DBRep::Set("BUC60857_BLUE",myshape);
  48   Handle(AIS_Shape) ais1 = new AIS_Shape(myshape);
  49   aContext->Display (ais1, Standard_False);
  50   aContext->SetColor (ais1, Quantity_NOC_BLUE1, Standard_False);
  51
  52   Handle(Geom_RectangularTrimmedSurface) S2 = GC_MakeTrimmedCone (P1, P2, R1, 0).Value();
  53   TopoDS_Shape myshape2 = BRepBuilderAPI_MakeFace(S2, Precision::Confusion()).Shape();
  54   DBRep::Set("BUC60857_RED",myshape2);
  55   Handle(AIS_Shape) ais2 = new AIS_Shape(myshape2);
  56   aContext->Display (ais2, Standard_False);
  57   aContext->SetColor (ais2, Quantity_NOC_RED, Standard_False);
  58
  59   Handle(Geom_RectangularTrimmedSurface) S3 = GC_MakeTrimmedCone (P1, P2, R2, R1).Value();
  60   TopoDS_Shape myshape3 = BRepBuilderAPI_MakeFace(S3, Precision::Confusion()).Shape();
  61   DBRep::Set("BUC60857_GREEN",myshape3);
  62   Handle(AIS_Shape) ais3 = new AIS_Shape(myshape3);
  63   aContext->Display (ais3, Standard_False);
  64   aContext->SetColor (ais3, Quantity_NOC_GREEN, Standard_True);
  65
  66   return 0;
  67 }
  68
  69 #include <GccEnt_Position.hxx>
  70 #include <Geom2dGcc_QualifiedCurve.hxx>
  71 #include <Geom2dGcc_Circ2d2TanRad.hxx>
  72 #include <gp_Elips2d.hxx>
  73 #include <Geom2d_Ellipse.hxx>
  74 #include <Geom2d_Circle.hxx>
  75 #include <AIS_InteractiveObject.hxx>
  76
  77 static Standard_Integer OCC24303(Draw_Interpretor& di, Standard_Integer n, const char** a)
  78 {
  79   if(n < 2)
  80     return 1;
  81
  82   const Standard_Integer SolID = Draw::Atoi(a[1]);
  83
  84   //Ellipses
  85   Standard_Real majorRadius = 2.0;
  86   Standard_Real minorRadius = 1.0;
  87   gp_Pnt2d p0(gp::Origin2d());
  88   gp_Pnt2d p1(4.0,0.0);
  89
  90   gp_Elips2d ellipse1 = gp_Elips2d( gp_Ax2d(p0,gp::DX2d()),majorRadius, minorRadius,true);
  91   gp_Elips2d ellipse2 = gp_Elips2d( gp_Ax2d(p1,gp::DX2d()),majorRadius, minorRadius,true);
  92
  93   Handle(Geom2d_Curve) curve1 = new Geom2d_Ellipse(ellipse1);
  94   Handle(Geom2d_Curve) curve2 = new Geom2d_Ellipse(ellipse2);
  95   DrawTrSurf::Set("c1", curve1);
  96   DrawTrSurf::Set("c2", curve2);
  97   //Expected tangent
  98   gp_Pnt2d centre(5.0,0.0);
  99   Standard_Real radius = 3.0;
 100   gp_Circ2d theorical_tangent = gp_Circ2d(gp_Ax2d(centre,gp::DX2d()),radius);
 101
 102   //Calculate the tangent with Geom2dGcc_Circ2dTanRan
 103
 104   const Geom2dAdaptor_Curve  AdaptedCurve1 ( curve1 );
 105   const Geom2dAdaptor_Curve  AdaptedCurve2 ( curve2 );
 106
 107   GccEnt_Position  curveQualif1 = GccEnt_unqualified;
 108   GccEnt_Position  curveQualif2 = GccEnt_unqualified;
 109
 110   const Geom2dGcc_QualifiedCurve qualifiedCurve1 ( AdaptedCurve1, curveQualif1 );
 111   const Geom2dGcc_QualifiedCurve qualifiedCurve2 ( AdaptedCurve2, curveQualif2 );
 112
 113   const Geom2dGcc_Circ2d2TanRad circCalc(qualifiedCurve1,qualifiedCurve2,radius,/*Precision::Approximation()*/ 1.0e-9);
 114
 115   const Standard_Integer aNbSol = circCalc.NbSolutions();
 116   di << "Solutions " << aNbSol << "\n";
 117
 118   if((SolID < 1) || (SolID > aNbSol))
 119   {
 120     di << "Wrong SolID value\n";
 121     return 1;
 122   }
 123
 124   gp_Circ2d calculated_tangent = circCalc.ThisSolution(SolID);
 125
 126   char Buf[10];
 127   for (Standard_Integer i = 1; i <= aNbSol; i++)
 128   {
 129     gp_Circ2d ct = circCalc.ThisSolution(i);
 130     Handle (Geom2d_Circle) GSol = new Geom2d_Circle(ct);
 131     Sprintf(Buf, "Sol%d",i);
 132     DrawTrSurf::Set(Buf, GSol);
 133   }
 134
 135   //This distance is different in OC 6.5.4 and OC 6.6.0
 136   Standard_Real dist = theorical_tangent.Location().Distance(calculated_tangent.Location());
 137   di << "Distance = " << dist << "\n";
 138
 139   return 0;
 140 }
 141
 142 void QABugs::Commands_9(Draw_Interpretor& theCommands) {
 143   const char *group = "QABugs";
 144
 145   theCommands.Add ("BUC60857", "BUC60857", __FILE__, BUC60857, group);
 146   theCommands.Add("OCC24303", "OCC24303 SolID ",        __FILE__,       OCC24303,group);
 147
 148   return;
 149 }