Salome HOME
Merge from V6_main 01/04/2013
[modules/geom.git] / src / GEOMImpl / GEOMImpl_LineDriver.cxx
1 // Copyright (C) 2007-2013  CEA/DEN, EDF R&D, OPEN CASCADE
2 //
3 // Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
5 //
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.
10 //
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.
15 //
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
19 //
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 //
22
23 #include <Standard_Stream.hxx>
24
25 #include <GEOMImpl_LineDriver.hxx>
26 #include <GEOMImpl_ILine.hxx>
27 #include <GEOMImpl_Types.hxx>
28 #include <GEOM_Function.hxx>
29
30 #include <BRep_Tool.hxx>
31 #include <BRepBuilderAPI_MakeEdge.hxx>
32 #include <BRepBuilderAPI_MakeVertex.hxx>
33 #include <BRepAlgoAPI_Section.hxx>
34 #include <TopAbs.hxx>
35 #include <TopExp.hxx>
36 #include <TopoDS.hxx>
37 #include <TopoDS_Edge.hxx>
38 #include <TopoDS_Shape.hxx>
39 #include <TopoDS_Vertex.hxx>
40 #include <TopExp_Explorer.hxx>
41 #include <TopTools_MapOfShape.hxx>
42
43 #include <gp_Pnt.hxx>
44 #include <Precision.hxx>
45 #include <Standard_ConstructionError.hxx>
46 #include <Standard_NullObject.hxx>
47
48 //=======================================================================
49 //function : GetID
50 //purpose  :
51 //=======================================================================
52 const Standard_GUID& GEOMImpl_LineDriver::GetID()
53 {
54   static Standard_GUID aLineDriver("FF1BBB06-5D14-4df2-980B-3A668264EA16");
55   return aLineDriver;
56 }
57
58
59 //=======================================================================
60 //function : GEOMImpl_LineDriver
61 //purpose  :
62 //=======================================================================
63 GEOMImpl_LineDriver::GEOMImpl_LineDriver()
64 {
65 }
66
67 //=======================================================================
68 //function : Execute
69 //purpose  :
70 //=======================================================================
71 Standard_Integer GEOMImpl_LineDriver::Execute(TFunction_Logbook& log) const
72 {
73   if (Label().IsNull())  return 0;
74   Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
75
76   GEOMImpl_ILine aPI (aFunction);
77   Standard_Integer aType = aFunction->GetType();
78
79   TopoDS_Shape aShape;
80
81   if (aType == LINE_TWO_PNT) {
82     Handle(GEOM_Function) aRefPnt1 = aPI.GetPoint1();
83     Handle(GEOM_Function) aRefPnt2 = aPI.GetPoint2();
84     TopoDS_Shape aShape1 = aRefPnt1->GetValue();
85     TopoDS_Shape aShape2 = aRefPnt2->GetValue();
86     if (aShape1.ShapeType() != TopAbs_VERTEX ||
87         aShape2.ShapeType() != TopAbs_VERTEX) {
88       Standard_ConstructionError::Raise("Wrong arguments: two points must be given");
89     }
90     if (aShape1.IsSame(aShape2)) {
91       Standard_ConstructionError::Raise("The end points must be different");
92     }
93     TopoDS_Vertex V1 = TopoDS::Vertex(aShape1);
94     TopoDS_Vertex V2 = TopoDS::Vertex(aShape2);
95     gp_Pnt P1 = BRep_Tool::Pnt(V1);
96     gp_Pnt P2 = BRep_Tool::Pnt(V2);
97     if (P1.Distance(P2) < Precision::Confusion()) {
98       Standard_ConstructionError::Raise("The end points are too close");
99     }
100     aShape = BRepBuilderAPI_MakeEdge(V1, V2).Shape();
101
102   } else if (aType == LINE_TWO_FACES) {
103     Handle(GEOM_Function) aRefFace1 = aPI.GetFace1();
104     Handle(GEOM_Function) aRefFace2 = aPI.GetFace2();
105     TopoDS_Shape aShape1 = aRefFace1->GetValue();
106     TopoDS_Shape aShape2 = aRefFace2->GetValue();
107     if (aShape1.ShapeType() != TopAbs_FACE ||
108         aShape2.ShapeType() != TopAbs_FACE) {
109       Standard_ConstructionError::Raise("Wrong arguments: two faces must be given");
110     }
111     if (aShape1.IsSame(aShape2)) {
112       Standard_ConstructionError::Raise("The end faces must be different");
113     }
114     BRepAlgoAPI_Section E (aShape1, aShape2, Standard_False);
115     E.Approximation(Standard_True);
116     E.Build();
117     if (!E.IsDone()) {
118       Standard_ConstructionError::Raise("Line can not be performed on the given faces");
119     }
120     else
121     {
122         TopExp_Explorer Exp (E, TopAbs_EDGE);
123         if ( Exp.More() ){
124             aShape = Exp.Current();
125             Exp.Next();
126         }
127         else
128           {
129             Standard_ConstructionError::Raise("Faces not have intersection line");
130             aShape = E.Shape();
131           }
132         if ( Exp.More() )
133           aShape = E.Shape();
134     }
135
136   } else if (aType == LINE_PNT_DIR) {
137     Handle(GEOM_Function) aRefPnt = aPI.GetPoint1();
138     Handle(GEOM_Function) aRefDir = aPI.GetPoint2();
139     TopoDS_Shape aShape1 = aRefPnt->GetValue();
140     TopoDS_Shape aShape2 = aRefDir->GetValue();
141     if (aShape1.ShapeType() != TopAbs_VERTEX) {
142       Standard_ConstructionError::Raise("Wrong first argument: must be point");
143     }
144     if (aShape2.ShapeType() != TopAbs_EDGE) {
145       Standard_ConstructionError::Raise("Wrong second argument: must be vector");
146     }
147     if (aShape1.IsSame(aShape2)) {
148       Standard_ConstructionError::Raise("The end points must be different");
149     }
150     TopoDS_Vertex Vstart = TopoDS::Vertex(aShape1);
151     gp_Pnt P1 = BRep_Tool::Pnt(Vstart);
152
153     TopoDS_Edge anE = TopoDS::Edge(aShape2);
154     TopoDS_Vertex V1, V2;
155     TopExp::Vertices(anE, V1, V2, Standard_True);
156     if (V1.IsNull() || V2.IsNull()) {
157       Standard_NullObject::Raise("Line creation aborted: vector is not defined");
158     }
159     gp_Pnt PV1 = BRep_Tool::Pnt(V1);
160     gp_Pnt PV2 = BRep_Tool::Pnt(V2);
161     if (PV1.Distance(PV2) < Precision::Confusion()) {
162       Standard_ConstructionError::Raise("Vector with null magnitude");
163     }
164     TopoDS_Vertex Vend;
165     if ( V1.IsSame( Vstart ) )
166       Vend = V2;
167     else if ( V2.IsSame( Vstart ) )
168       Vend = V1;
169     else
170       Vend = BRepBuilderAPI_MakeVertex( gp_Pnt(P1.XYZ() + PV2.XYZ() - PV1.XYZ()) );
171     aShape = BRepBuilderAPI_MakeEdge(Vstart, Vend).Shape();
172   } else {
173   }
174
175   if (aShape.IsNull()) return 0;
176   //aShape.Infinite(true); // VSR: 05/04/2010: Fix 20668 (Fit All for points & lines)
177
178   aFunction->SetValue(aShape);
179
180   log.SetTouched(Label());
181
182   return 1;
183 }
184
185
186 //=======================================================================
187 //function :  GEOMImpl_LineDriver_Type_
188 //purpose  :
189 //=======================================================================
190 Standard_EXPORT Handle_Standard_Type& GEOMImpl_LineDriver_Type_()
191 {
192
193   static Handle_Standard_Type aType1 = STANDARD_TYPE(TFunction_Driver);
194   if ( aType1.IsNull()) aType1 = STANDARD_TYPE(TFunction_Driver);
195   static Handle_Standard_Type aType2 = STANDARD_TYPE(MMgt_TShared);
196   if ( aType2.IsNull()) aType2 = STANDARD_TYPE(MMgt_TShared);
197   static Handle_Standard_Type aType3 = STANDARD_TYPE(Standard_Transient);
198   if ( aType3.IsNull()) aType3 = STANDARD_TYPE(Standard_Transient);
199
200
201   static Handle_Standard_Transient _Ancestors[]= {aType1,aType2,aType3,NULL};
202   static Handle_Standard_Type _aType = new Standard_Type("GEOMImpl_LineDriver",
203                                                          sizeof(GEOMImpl_LineDriver),
204                                                          1,
205                                                          (Standard_Address)_Ancestors,
206                                                          (Standard_Address)NULL);
207
208   return _aType;
209 }
210
211 //=======================================================================
212 //function : DownCast
213 //purpose  :
214 //=======================================================================
215 const Handle(GEOMImpl_LineDriver) Handle(GEOMImpl_LineDriver)::DownCast
216        (const Handle(Standard_Transient)& AnObject)
217 {
218   Handle(GEOMImpl_LineDriver) _anOtherObject;
219
220   if (!AnObject.IsNull()) {
221      if (AnObject->IsKind(STANDARD_TYPE(GEOMImpl_LineDriver))) {
222        _anOtherObject = Handle(GEOMImpl_LineDriver)((Handle(GEOMImpl_LineDriver)&)AnObject);
223      }
224   }
225
226   return _anOtherObject ;
227 }