0023319: Typo in AIS_ConcentricRelation.cxx line 156
[occt.git] / src / AIS / AIS_ConcentricRelation.cxx
1 // Created on: 1996-12-05
2 // Created by: Flore Lantheaume/Odile Olivier
3 // Copyright (c) 1996-1999 Matra Datavision
4 // Copyright (c) 1999-2012 OPEN CASCADE SAS
5 //
6 // The content of this file is subject to the Open CASCADE Technology Public
7 // License Version 6.5 (the "License"). You may not use the content of this file
8 // except in compliance with the License. Please obtain a copy of the License
9 // at http://www.opencascade.org and read it completely before using this file.
10 //
11 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
12 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
13 //
14 // The Original Code and all software distributed under the License is
15 // distributed on an "AS IS" basis, without warranty of any kind, and the
16 // Initial Developer hereby disclaims all such warranties, including without
17 // limitation, any warranties of merchantability, fitness for a particular
18 // purpose or non-infringement. Please see the License for the specific terms
19 // and conditions governing the rights and limitations under the License.
20
21
22 #include <Standard_NotImplemented.hxx>
23
24 #include <AIS_ConcentricRelation.ixx>
25
26 #include <SelectMgr_EntityOwner.hxx>
27 #include <Select3D_SensitiveCircle.hxx>
28 #include <Select3D_SensitiveSegment.hxx>
29
30 #include <DsgPrs_ConcentricPresentation.hxx>
31
32 #include <TopoDS.hxx>
33
34 #include <BRepAdaptor_Curve.hxx>
35
36 #include <GeomAbs_CurveType.hxx>
37 #include <Geom_Circle.hxx>
38 #include <gp_Dir.hxx>
39 #include <gp_Vec.hxx>
40 #include <gp_Ax2.hxx>
41 #include <gp_Ax1.hxx>
42 #include <gp_Pnt.hxx>
43 #include <gp_Pln.hxx>
44
45 #include <TopoDS_Vertex.hxx>
46 #include <AIS.hxx>
47
48 //=======================================================================
49 //function : Constructor
50 //purpose  : 
51 //=======================================================================
52
53 AIS_ConcentricRelation::AIS_ConcentricRelation(
54         const TopoDS_Shape& aFShape, 
55         const TopoDS_Shape& aSShape, 
56         const Handle(Geom_Plane)& aPlane)
57 {
58   myFShape = aFShape;
59   mySShape = aSShape;
60   myPlane = aPlane;
61   myDir = aPlane->Pln().Axis().Direction();
62 }
63
64 //=======================================================================
65 //function : Compute
66 //purpose  : 
67 //=======================================================================
68 void AIS_ConcentricRelation::Compute(const Handle(PrsMgr_PresentationManager3d)&, 
69                                      const Handle(Prs3d_Presentation)& aPresentation, 
70                                      const Standard_Integer)
71 {
72   aPresentation->Clear();
73
74   TopAbs_ShapeEnum type2(mySShape.ShapeType());
75   aPresentation->SetInfiniteState(Standard_True);
76   switch (myFShape.ShapeType()) {
77   case TopAbs_EDGE: 
78     {
79       if (type2 == TopAbs_EDGE) ComputeTwoEdgesConcentric(aPresentation);
80       else if (type2 == TopAbs_VERTEX) ComputeEdgeVertexConcentric(aPresentation);
81     }
82   break;
83   
84   case TopAbs_VERTEX: 
85     {
86       if (type2 == TopAbs_VERTEX) ComputeTwoVerticesConcentric(aPresentation);
87       else if (type2 == TopAbs_EDGE) ComputeEdgeVertexConcentric(aPresentation);      
88     }
89   break;
90   default: {return;}
91   }  
92 }
93
94 void AIS_ConcentricRelation::Compute(const Handle_Prs3d_Projector& aProjector, const Handle_Geom_Transformation& aTransformation, const Handle_Prs3d_Presentation& aPresentation)
95 {
96 // Standard_NotImplemented::Raise("AIS_ConcentricRelation::Compute(const Handle_Prs3d_Projector&, const Handle_Geom_Transformation&, const Handle_Prs3d_Presentation&)");
97   PrsMgr_PresentableObject::Compute( aProjector , aTransformation , aPresentation ) ;
98 }
99
100 //=======================================================================
101 //function : ComputeTwoEdgesConcentric
102 //purpose  : 
103 //=======================================================================
104 void AIS_ConcentricRelation::ComputeEdgeVertexConcentric(const Handle(Prs3d_Presentation)& aPresentation)
105 {
106   TopoDS_Edge E;
107   TopoDS_Vertex V;
108   if (myFShape.ShapeType() == TopAbs_EDGE) {
109     E = TopoDS::Edge(myFShape);
110     V = TopoDS::Vertex(mySShape);
111   }
112   else {
113     E = TopoDS::Edge(mySShape);
114     V = TopoDS::Vertex(myFShape);
115   }
116   gp_Pnt p1,p2;
117   Handle(Geom_Curve) C;
118   Handle(Geom_Curve) extCurv;
119   Standard_Boolean isInfinite;
120   Standard_Boolean isOnPlanEdge, isOnPlanVertex;
121   if (!AIS::ComputeGeometry(E,C,p1,p2,extCurv,isInfinite,isOnPlanEdge,myPlane)) return;
122   gp_Pnt P;
123   AIS::ComputeGeometry(V,P, myPlane, isOnPlanVertex);
124
125   const Handle(Geom_Circle)& CIRCLE = (Handle(Geom_Circle)&) C;
126   myCenter = CIRCLE->Location();
127   myRad = Min(CIRCLE->Radius()/5.,15.);
128   gp_Dir vec(p1.XYZ() - myCenter.XYZ() );
129   gp_Vec vectrans(vec);
130   myPnt = myCenter.Translated(vectrans.Multiplied(myRad));
131   DsgPrs_ConcentricPresentation::Add(aPresentation,myDrawer,myCenter,myRad,myDir,myPnt);
132   if (!isOnPlanEdge) AIS::ComputeProjEdgePresentation(aPresentation,myDrawer,E,CIRCLE,p1,p2);
133   if (!isOnPlanVertex) AIS::ComputeProjVertexPresentation(aPresentation,myDrawer,V,P);
134 }
135
136 //=======================================================================
137 //function : ComputeTwoEdgesConcentric
138 //purpose  : 
139 //=======================================================================
140 void AIS_ConcentricRelation::ComputeTwoVerticesConcentric(const Handle(Prs3d_Presentation)& aPresentation)
141 {
142   TopoDS_Vertex V1,V2;
143   V1 = TopoDS::Vertex(myFShape);
144   V2 = TopoDS::Vertex(myFShape);  
145   Standard_Boolean isOnPlanVertex1(Standard_True),isOnPlanVertex2(Standard_True);
146   gp_Pnt P1,P2;
147   AIS::ComputeGeometry(V1,P1, myPlane,isOnPlanVertex1);
148   AIS::ComputeGeometry(V2,P2, myPlane,isOnPlanVertex2);
149   myCenter = P1;
150   myRad    = 15.;
151   gp_Dir vec(myPlane->Pln().Position().XDirection());
152   gp_Vec vectrans(vec);
153   myPnt = myCenter.Translated(vectrans.Multiplied(myRad));
154   DsgPrs_ConcentricPresentation::Add(aPresentation,myDrawer,myCenter,myRad,myDir,myPnt);
155   if (!isOnPlanVertex1) AIS::ComputeProjVertexPresentation(aPresentation,myDrawer,V1,P1);
156   if (!isOnPlanVertex2) AIS::ComputeProjVertexPresentation(aPresentation,myDrawer,V2,P2);
157 }
158
159 //=======================================================================
160 //function : ComputeTwoEdgesConcentric
161 //purpose  : 
162 //=======================================================================
163 void AIS_ConcentricRelation::ComputeTwoEdgesConcentric(const Handle(Prs3d_Presentation)& aPresentation)
164 {
165   BRepAdaptor_Curve curv1(TopoDS::Edge(myFShape));
166   BRepAdaptor_Curve curv2(TopoDS::Edge(mySShape));
167   
168   gp_Pnt ptat11,ptat12,ptat21,ptat22;
169   Handle(Geom_Curve) geom1,geom2;
170   Standard_Boolean isInfinite1,isInfinite2;
171   Handle(Geom_Curve) extCurv;
172   if (!AIS::ComputeGeometry(TopoDS::Edge(myFShape),
173                             TopoDS::Edge(mySShape),
174                             myExtShape,
175                             geom1,
176                             geom2,
177                             ptat11,
178                             ptat12,
179                             ptat21,
180                             ptat22,
181                             extCurv,
182                             isInfinite1,isInfinite2,
183                             myPlane)) {
184     return;
185   }
186   
187   const Handle(Geom_Circle)& gcirc1 = (Handle(Geom_Circle)&) geom1;
188   const Handle(Geom_Circle)& gcirc2 = (Handle(Geom_Circle)&) geom2;
189   
190   myCenter = gcirc1->Location();
191   
192   // choose the radius equal to 1/5 of the smallest radius of 
193   // 2 circles. Limit is imposed ( 0.02 by chance)
194   Standard_Real rad1 = gcirc1->Radius();
195   Standard_Real rad2 = gcirc2->Radius();
196   myRad = (rad1 > rad2 ) ? rad2 : rad1;
197   myRad /= 5;
198   if (myRad > 15.) myRad =15.;
199   
200   
201   //Calculate a point of circle of radius myRad
202   gp_Dir vec(ptat11.XYZ() - myCenter.XYZ() );
203   gp_Vec vectrans(vec);
204   myPnt = myCenter.Translated(vectrans.Multiplied(myRad));
205   
206   DsgPrs_ConcentricPresentation::Add(aPresentation,
207                                      myDrawer,
208                                      myCenter,
209                                      myRad,
210                                      myDir,
211                                      myPnt);
212   if ( (myExtShape != 0) &&  !extCurv.IsNull()) {
213     gp_Pnt pf, pl;
214     if ( myExtShape == 1 ) {
215       if (!isInfinite1) {
216         pf = ptat11; 
217         pl = ptat12;
218       }
219       ComputeProjEdgePresentation(aPresentation,TopoDS::Edge(myFShape),gcirc1,pf,pl);
220     }
221     else {
222       if (!isInfinite2) {
223         pf = ptat21; 
224         pl = ptat22;
225       }
226       ComputeProjEdgePresentation(aPresentation,TopoDS::Edge(mySShape),gcirc2,pf,pl);
227     }
228   }
229 }
230
231 //=======================================================================
232 //function : Compute
233 //purpose  : to avoid warning
234 //=======================================================================
235
236 void AIS_ConcentricRelation::Compute(const Handle(Prs3d_Projector)&, 
237                                      const Handle(Prs3d_Presentation)&)
238 {
239 }
240
241 //=======================================================================
242 //function : Compute
243 //purpose  : to avoid warning
244 //=======================================================================
245
246 void AIS_ConcentricRelation::Compute(const Handle(PrsMgr_PresentationManager2d)&, 
247                                      const Handle(Graphic2d_GraphicObject)&,
248                                      const Standard_Integer)
249 {
250 }
251
252 //=======================================================================
253 //function : ComputeSelection
254 //purpose  : 
255 //=======================================================================
256
257 void AIS_ConcentricRelation::ComputeSelection(const Handle(SelectMgr_Selection)& aSelection, 
258                                               const Standard_Integer)
259 {
260   Handle(SelectMgr_EntityOwner) own = new SelectMgr_EntityOwner(this,7);
261   
262   //Creation of 2 sensitive circles
263      // the greater
264   gp_Ax2 ax(myCenter, myDir);
265   Handle(Geom_Circle) Circ = new Geom_Circle(ax, myRad) ;
266   Handle(Select3D_SensitiveCircle) 
267     sensit = new Select3D_SensitiveCircle (own,
268                                            Circ);
269   aSelection->Add(sensit);
270      // the smaller
271   Circ->SetRadius(myRad/2);
272   sensit = new Select3D_SensitiveCircle (own,
273                                          Circ);
274   aSelection->Add(sensit);
275
276   //Creation of 2 segments sensitive for the cross
277   Handle(Select3D_SensitiveSegment) seg;
278   gp_Pnt otherPnt = myPnt.Mirrored(myCenter);
279   seg = new Select3D_SensitiveSegment(own,
280                                       otherPnt,
281                                       myPnt);
282   aSelection->Add(seg);
283
284   gp_Ax1 RotateAxis(myCenter, myDir);
285   gp_Pnt FPnt = myCenter.Rotated(RotateAxis, M_PI/2);
286   gp_Pnt SPnt = myCenter.Rotated(RotateAxis, -M_PI/2);
287   seg = new Select3D_SensitiveSegment(own,
288                                       FPnt,
289                                       SPnt);
290   aSelection->Add(seg);
291
292 }
293