0031501: Foundation Classes, Message_Printer - remove theToPutEndl argument -- prepar...
[occt.git] / src / IGESGeom / IGESGeom_ToolBSplineCurve.cxx
1 // Created by: CKY / Contract Toubro-Larsen
2 // Copyright (c) 1993-1999 Matra Datavision
3 // Copyright (c) 1999-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 //--------------------------------------------------------------------
17 //--------------------------------------------------------------------
18 //#54 rln 24.12.98 CCI60005
19
20 #include <gp_XYZ.hxx>
21 #include <IGESData_DirChecker.hxx>
22 #include <IGESData_Dump.hxx>
23 #include <IGESData_IGESDumper.hxx>
24 #include <IGESData_IGESReaderData.hxx>
25 #include <IGESData_IGESWriter.hxx>
26 #include <IGESData_ParamCursor.hxx>
27 #include <IGESData_ParamReader.hxx>
28 #include <IGESGeom_BSplineCurve.hxx>
29 #include <IGESGeom_ToolBSplineCurve.hxx>
30 #include <Interface_Check.hxx>
31 #include <Interface_CopyTool.hxx>
32 #include <Interface_EntityIterator.hxx>
33 #include <Interface_Macros.hxx>
34 #include <Interface_ShareTool.hxx>
35 #include <Message_Messenger.hxx>
36 #include <Message_Msg.hxx>
37 #include <Standard_DomainError.hxx>
38 #include <TColgp_HArray1OfXYZ.hxx>
39 #include <TColStd_HArray1OfReal.hxx>
40
41 // MGE 29/07/98
42 //=======================================================================
43 //function : IGESGeom_ToolBSplineCurve
44 //purpose  : 
45 //=======================================================================
46 IGESGeom_ToolBSplineCurve::IGESGeom_ToolBSplineCurve ()
47 {
48 }
49
50
51 //=======================================================================
52 //function : ReadOwnParams
53 //purpose  : 
54 //=======================================================================
55
56 void IGESGeom_ToolBSplineCurve::ReadOwnParams(const Handle(IGESGeom_BSplineCurve)& ent,
57                                               const Handle(IGESData_IGESReaderData)& /* IR */,
58                                               IGESData_ParamReader& PR) const
59 {
60   // MGE 29/07/98
61   // Building of messages
62   //========================================
63   Message_Msg Msg99("XSTEP_99");
64   Message_Msg Msg100("XSTEP_100");
65   Message_Msg Msg101("XSTEP_101");
66   Message_Msg Msg102("XSTEP_102");
67   Message_Msg Msg103("XSTEP_103");
68   //========================================
69
70   Standard_Integer anIndex, aDegree;
71   Standard_Boolean aPlanar, aClosed, aPolynomial, aPeriodic;
72   Standard_Real aUmin, aUmax, normX,normY,normZ;
73   gp_XYZ aNorm (0.,0.,0.);
74   Handle(TColStd_HArray1OfReal) allKnots;
75   Handle(TColStd_HArray1OfReal) allWeights;
76   Handle(TColgp_HArray1OfXYZ)   allPoles;
77
78   //Standard_Boolean st; //szv#4:S4163:12Mar99 moved down
79
80   //st = PR.ReadInteger(PR.Current(), Msg97, anIndex); //szv#4:S4163:12Mar99 moved in if
81   //st = PR.ReadInteger(PR.Current(), "Upper Index Of Sum", anIndex);
82
83   //szv#4:S4163:12Mar99 optimized
84   /*if (st && anIndex >= 0) {
85     allPoles   = new TColgp_HArray1OfXYZ(0, anIndex);
86     // allWeights = new TColStd_HArray1OfReal(1, anIndex+1);  done by ReadReals
87   }
88   
89   if (st && anIndex < 0)
90   { 
91     PR.SendFail(Msg97);
92     anIndex = 0;
93   }*/
94   if (PR.ReadInteger(PR.Current(), anIndex)) {
95     if (anIndex < 0) {
96       Message_Msg Msg97("XSTEP_97");
97       PR.SendFail(Msg97);
98       anIndex = 0;
99     }
100     else {
101       allPoles   = new TColgp_HArray1OfXYZ(0, anIndex);
102       // allWeights = new TColStd_HArray1OfReal(1, anIndex+1);  done by ReadReals
103     }
104   }
105   else{
106     Message_Msg Msg97("XSTEP_97");
107     PR.SendFail(Msg97);
108   }
109
110   //st = PR.ReadInteger(PR.Current(), Msg98, aDegree); //szv#4:S4163:12Mar99 moved in if
111 //  if (st && ! allWeights.IsNull() )   done by ReadReals
112 //    allKnots = new TColStd_HArray1OfReal(-aDegree, anIndex+1);
113   if (!PR.ReadInteger(PR.Current(), aDegree)){
114     aDegree = 0; //szv#4:S4163:12Mar99 `st=` not needed
115     Message_Msg Msg98("XSTEP_98");
116     PR.SendFail(Msg98);
117   }
118   //szv#4:S4163:12Mar99 `st=` not needed
119   PR.ReadBoolean(PR.Current(), Msg99, aPlanar);
120   PR.ReadBoolean(PR.Current(), Msg100, aClosed);
121   PR.ReadBoolean(PR.Current(), Msg101, aPolynomial);
122   PR.ReadBoolean(PR.Current(), Msg102, aPeriodic);
123
124 //st = PR.ReadBoolean(PR.Current(), "Planar/Non Planar Flag", aPlanar);
125 //st = PR.ReadBoolean(PR.Current(), "Open/Closed Flag", aClosed);
126 //st = PR.ReadBoolean(PR.Current(), "Rational/Polynomial Flag", aPolynomial);
127 //st = PR.ReadBoolean(PR.Current(), "NonPeriodic/Periodic Flag", aPeriodic);
128
129   Standard_Integer nbKnots = anIndex + aDegree + 2;
130   // Reading all the knot sequences 
131
132   PR.ReadReals(PR.CurrentList(nbKnots), Msg103 , allKnots, -aDegree); //szv#4:S4163:12Mar99 `st=` not needed
133
134 //st = PR.ReadReals
135 //  (PR.CurrentList(nbKnots), "Knot sequence values", allKnots, -aDegree);
136
137   if (! allPoles.IsNull() )
138     {
139       Message_Msg Msg104("XSTEP_104");
140       Message_Msg Msg105("XSTEP_105");
141       PR.ReadReals(PR.CurrentList(anIndex+1), Msg104, allWeights,0); //szv#4:S4163:12Mar99 `st=` not needed
142       //st = PR.ReadReals(PR.CurrentList(anIndex+1), "Weights", allWeights,0);
143
144       for (Standard_Integer I = 0; I <= anIndex; I ++) 
145         {
146           gp_XYZ tempPole;
147           //st = PR.ReadXYZ(PR.CurrentList(1, 3), Msg105, tempPole); //szv#4:S4163:12Mar99 moved down
148           //st = PR.ReadXYZ(PR.CurrentList(1, 3), "Control Points", tempPole);
149           if (PR.ReadXYZ(PR.CurrentList(1, 3), Msg105, tempPole)) allPoles->SetValue(I, tempPole);
150         }
151     }
152
153   if (!PR.ReadReal(PR.Current(), aUmin)){
154     Message_Msg Msg106("XSTEP_106");
155     PR.SendFail(Msg106);
156   } //szv#4:S4163:12Mar99 `st=` not needed
157   if (!PR.ReadReal(PR.Current(), aUmax)){
158     Message_Msg Msg107("XSTEP_107");
159     PR.SendFail(Msg107);
160   } //szv#4:S4163:12Mar99 `st=` not needed
161 /*
162   st = PR.ReadReal(PR.Current(), "Starting Parameter Value", aUmin);
163   st = PR.ReadReal(PR.Current(), "Ending Parameter Value", aUmax);
164 */
165   Standard_Boolean st = Standard_False;
166   if (PR.DefinedElseSkip()){
167     st = PR.ReadReal(PR.Current(),  normX);
168     if(!st){
169       Message_Msg Msg108("XSTEP_108");
170       PR.SendFail(Msg108);
171     }
172   }
173     //st = PR.ReadReal(PR.Current(), "Unit Normal X", normX);
174   else normX = 0.;
175   if (PR.DefinedElseSkip()){
176     st = PR.ReadReal(PR.Current(),  normY);
177     if(!st){
178       Message_Msg Msg108("XSTEP_108");
179       PR.SendFail(Msg108);
180     }
181   }
182     //st = PR.ReadReal(PR.Current(), "Unit Normal Y", normY);
183   else normY = 0.;
184   if (PR.DefinedElseSkip()){
185     st = PR.ReadReal(PR.Current(), normZ);
186     if(!st){
187       Message_Msg Msg108("XSTEP_108");
188       PR.SendFail(Msg108);
189     }
190   }
191     //st = PR.ReadReal(PR.Current(), "Unit Normal Z", normZ);
192   else normZ = 0.;
193   if (st) aNorm.SetCoord(normX,normY,normZ);
194
195   DirChecker(ent).CheckTypeAndForm(PR.CCheck(),ent);
196   ent->Init
197     (anIndex,aDegree, aPlanar, aClosed, aPolynomial, aPeriodic,
198      allKnots, allWeights, allPoles, aUmin, aUmax, aNorm);
199 }
200
201
202 //=======================================================================
203 //function : WriteOwnParams
204 //purpose  : 
205 //=======================================================================
206
207 void IGESGeom_ToolBSplineCurve::WriteOwnParams
208   (const Handle(IGESGeom_BSplineCurve)& ent, IGESData_IGESWriter& IW)  const
209 {
210   Standard_Integer low, up ;
211   Standard_Integer I ;
212   Standard_Integer index  = ent->UpperIndex();
213   Standard_Integer degree = ent->Degree();
214   IW.Send(index);
215   IW.Send(degree);
216   IW.SendBoolean(ent->IsPlanar());
217   IW.SendBoolean(ent->IsClosed());
218   IW.SendBoolean(ent->IsPolynomial());
219   IW.SendBoolean(ent->IsPeriodic());
220
221   low = -degree;
222   up  = index + 1;
223   for (I = low; I <= up; I ++)
224     IW.Send(ent->Knot(I));
225
226   low = 0;
227   up  = index;
228   for (I = low; I <= up; I ++)
229     IW.Send(ent->Weight(I));
230
231   for (I = low; I <= up; I ++) {
232     IW.Send((ent->Pole(I)).X());
233     IW.Send((ent->Pole(I)).Y());
234     IW.Send((ent->Pole(I)).Z());
235   }
236   IW.Send(ent->UMin());
237   IW.Send(ent->UMax());
238   IW.Send(ent->Normal().X());
239   IW.Send(ent->Normal().Y());
240   IW.Send(ent->Normal().Z());
241 }
242
243
244 //=======================================================================
245 //function : OwnShared
246 //purpose  : 
247 //=======================================================================
248
249 void IGESGeom_ToolBSplineCurve::OwnShared(const Handle(IGESGeom_BSplineCurve)& /* ent */,
250                                           Interface_EntityIterator& /* iter */) const
251 {
252 }
253
254
255 //=======================================================================
256 //function : OwnCopy
257 //purpose  : 
258 //=======================================================================
259
260 void IGESGeom_ToolBSplineCurve::OwnCopy 
261   (const Handle(IGESGeom_BSplineCurve)& another,
262    const Handle(IGESGeom_BSplineCurve)& ent, Interface_CopyTool& /* TC */) const
263 {
264   Standard_Integer I;
265   Standard_Integer low, up;
266   Standard_Integer anIndex, aDegree;
267   Standard_Boolean aPlanar, aClosed, aPolynomial, aPeriodic;
268   Handle(TColStd_HArray1OfReal) allKnots, allWeights;
269   Handle(TColgp_HArray1OfXYZ) allPoles;
270   Standard_Real aUmin, aUmax;
271   gp_XYZ aNorm;
272
273   anIndex = another->UpperIndex();
274   aDegree = another->Degree();
275   aPlanar = another->IsPlanar();
276   aClosed = another->IsClosed();
277   aPolynomial = another->IsPolynomial();
278   aPeriodic = another->IsPeriodic();
279
280   allKnots = new TColStd_HArray1OfReal(-aDegree, anIndex+1);
281
282   low = -aDegree;
283   up  = anIndex + 1;
284   for (I = low; I <= up; I++)
285     allKnots->SetValue(I, another->Knot(I));
286
287   allWeights = new TColStd_HArray1OfReal(0, anIndex);
288
289   low = 0;
290   up  = anIndex;
291   for (I = low; I <= up; I++)
292     allWeights->SetValue(I, another->Weight(I));
293
294   allPoles = new TColgp_HArray1OfXYZ(0, anIndex);
295
296   for (I = low; I <= up; I++) 
297     allPoles->SetValue(I, (another->Pole(I)).XYZ());
298
299   aUmin = another->UMin();
300   aUmax = another->UMax();
301   aNorm = another->Normal();
302
303   ent->Init (anIndex,aDegree, aPlanar, aClosed, aPolynomial, aPeriodic,
304              allKnots, allWeights, allPoles, aUmin, aUmax, aNorm);
305 }
306
307
308 //=======================================================================
309 //function : DirChecker
310 //purpose  : 
311 //=======================================================================
312
313 IGESData_DirChecker IGESGeom_ToolBSplineCurve::DirChecker
314   (const Handle(IGESGeom_BSplineCurve)& /* ent */ )   const
315 {
316   IGESData_DirChecker DC(126, 0, 5);
317   DC.Structure(IGESData_DefVoid);
318   DC.LineFont(IGESData_DefAny);
319 //  DC.LineWeight(IGESData_DefValue);
320   DC.Color(IGESData_DefAny);
321   DC.HierarchyStatusIgnored();
322   return DC;
323 }
324
325
326 //=======================================================================
327 //function : OwnCheck
328 //purpose  : 
329 //=======================================================================
330
331 void IGESGeom_ToolBSplineCurve::OwnCheck(const Handle(IGESGeom_BSplineCurve)& ent,
332                                          const Interface_ShareTool&,
333                                          Handle(Interface_Check)& ach)  const
334 {
335   // MGE 29/07/98
336   // Building of messages
337   //========================================
338   //Message_Msg Msg104("XSTEP_104");
339   //Message_Msg Msg109("XSTEP_109");
340   //========================================
341
342   Standard_Real eps = 1.E-04;    // Tolerance des tests ??
343 //  Standard_Real norm = ent->Normal().SquareModulus();
344
345   //modified by rln 17/12/97 check of flag PROP2 according to IGES Standard
346   //It is possible to compare V(0) and V(1) only if StartingParameter = FirstKnot
347   //and EndingParameter = LastKnot (else we must build real geometrical curve)
348   //The fail is replaced with warning because it is not a serious problem
349   //if (ent->UMin() == ent->Knot(-ent->Degree()        ) &&
350   //  ent->UMax() == ent->Knot( ent->UpperIndex() + 1)   ) {
351   //  Standard_Real udif = ent->UMax() - ent->UMin();
352   //  if (udif < 0) udif = -udif;
353   //    Standard_Real udif = ent->Pole(0).SquareDistance (ent->Pole(ent->UpperIndex()));
354   //   if (udif <  eps * eps && !ent->IsClosed())
355   //     ach.AddWarning("V(0) == V(1) for an Open Curve (PROP2 = 0)");
356   //   if (udif >= eps * eps &&  ent->IsClosed())
357   //     ach.AddWarning("V(0) != V(1) for a Closed Curve (PROP2 = 1)");
358   //}
359
360   Standard_Integer lower = 0;
361   Standard_Integer upper = ent->UpperIndex();
362   Standard_Boolean Flag  = Standard_True;
363   
364   Standard_Integer I; // svv Jan 11 2000 : porting on DEC
365   for (I = 0; ((I < upper) && (Flag)); I++)
366     Flag &= (ent->Weight(I) > 0);
367   
368   if (!Flag) {
369     Message_Msg Msg104("XSTEP_104");
370     ach->SendFail(Msg104);
371   }
372   
373   Flag = Standard_True;
374   Standard_Real tempVal = ent->Weight(lower);
375   
376   for (I = lower; ((I < upper) && (Flag)); I++)
377     Flag &= (ent->Weight(I) == tempVal);
378 /*
379   if (Flag && !ent->IsPolynomial(Standard_True))
380     ach.AddWarning("All weights equal & PROP3 != 1 (Curve Not Polynomial)");
381   if (!Flag && ent->IsPolynomial(Standard_True))
382     ach.AddWarning("All weights not equal & PROP3 != 0 (Curve Not Rational)");
383 */
384
385   if (ent->IsPlanar()) {
386     gp_XYZ aNorm = ent->Normal();
387     Standard_Real epsn = eps * 10.;    // Tolerance ?? ici large
388     Standard_Real normod = aNorm.SquareModulus();
389     if (normod < epsn) {
390       Message_Msg Msg109("XSTEP_109");
391       ach->AddWarning(Msg109);
392     }
393   }
394 }
395
396
397 //=======================================================================
398 //function : OwnDump
399 //purpose  : 
400 //=======================================================================
401
402 void IGESGeom_ToolBSplineCurve::OwnDump(const Handle(IGESGeom_BSplineCurve)& ent,
403                                         const IGESData_IGESDumper& /* dumper */,
404                                         Standard_OStream& S,
405                                         const Standard_Integer level)  const
406 {
407   Standard_Integer upind = ent->UpperIndex();
408   S << "BSplineCurve from IGESGeom\n"
409     << "Sum UpperIndex : " << upind
410     << "   Degree : " << ent->Degree() << "  "
411     << (ent->IsPlanar() ? "Planar" : "NonPlanar") << "\n"
412     << (ent->IsClosed() ? "Closed" : "Open") << "  "
413     << (ent->IsPeriodic() ? "Periodic" : "NonPeriodic") << "  " //#54 rln 24.12.98 CCI60005
414     << (ent->IsPolynomial(Standard_True) ? "Polynomial" : "Rational")
415     << "\nKnots : ";
416   IGESData_DumpVals(S,level,-ent->Degree(), upind+1,ent->Knot);
417   S << "\nWeights : ";
418   IGESData_DumpVals(S,level,0, upind,ent->Weight);
419   S << "\nControl Points (Poles) : ";
420   IGESData_DumpListXYZL(S,level,0, upind, ent->Pole, ent->Location());
421   S << "\nStarting Parameter Value : " << ent->UMin()
422     << "  Ending Parameter Value : " << ent->UMax() << "\n"
423     << "Unit Normal : ";
424   IGESData_DumpXYZL(S,level, ent->Normal(), ent->Location());
425   S << std::endl;
426 }