0029846: Data Exchange - STEP validation properties support is obsolete
[occt.git] / src / XDEDRAW / XDEDRAW_Props.cxx
1 // Created on: 2000-08-04
2 // Created by: Pavel TELKOV
3 // Copyright (c) 2000-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 // OCC532 sln 24.07.2002. Add epsilon parameter to SetProps function
17
18 #include <BRep_Tool.hxx>
19 #include <BRepGProp.hxx>
20 #include <BRepMesh_IncrementalMesh.hxx>
21 #include <DBRep.hxx>
22 #include <DDocStd.hxx>
23 #include <Draw.hxx>
24 #include <gp_Pln.hxx>
25 #include <GProp_GProps.hxx>
26 #include <Poly_Triangulation.hxx>
27 #include <Precision.hxx>
28 #include <Standard_ErrorHandler.hxx>
29 #include <Standard_Failure.hxx>
30 #include <Standard_Stream.hxx>
31 #include <TColgp_Array1OfPnt.hxx>
32 #include <TColgp_Array1OfXYZ.hxx>
33 #include <TCollection_AsciiString.hxx>
34 #include <TCollection_HAsciiString.hxx>
35 #include <TColStd_Array1OfReal.hxx>
36 #include <TDataStd_Name.hxx>
37 #include <TDataStd_TreeNode.hxx>
38 #include <TDF_Label.hxx>
39 #include <TDF_LabelSequence.hxx>
40 #include <TDF_Tool.hxx>
41 #include <TDocStd_Document.hxx>
42 #include <TNaming_NamedShape.hxx>
43 #include <TopExp_Explorer.hxx>
44 #include <TopoDS.hxx>
45 #include <TopoDS_Compound.hxx>
46 #include <TopoDS_Face.hxx>
47 #include <TopoDS_Shape.hxx>
48 #include <TopoDS_Vertex.hxx>
49 #include <XCAFDoc.hxx>
50 #include <XCAFDoc_Area.hxx>
51 #include <XCAFDoc_Centroid.hxx>
52 #include <XCAFDoc_DocumentTool.hxx>
53 #include <XCAFDoc_Location.hxx>
54 #include <XCAFDoc_MaterialTool.hxx>
55 #include <XCAFDoc_ShapeTool.hxx>
56 #include <XCAFDoc_Volume.hxx>
57 #include <XDEDRAW_Props.hxx>
58
59 #include <stdio.h>
60 // --------------------- VolumeFix Begin ---
61 //=======================================================================
62 //function : TetraVol
63 //purpose  : auxilary
64 //=======================================================================
65 static double TetraVol(gp_Pnt RefPoint, gp_Pnt Som1, gp_Pnt Som2, gp_Pnt Som3)
66 {
67   double curVolume = 0;
68   gp_Dir Line12;
69   gp_Pln Plane123;
70   gp_Vec N;
71   
72   {
73     try
74       {
75       OCC_CATCH_SIGNALS
76         Line12=gp_Dir( gp_Vec(Som1, Som2));
77         gp_Vec v1(Som1, Som2);
78         gp_Vec v2(Som2, Som3);
79         N=v1^v2;
80         Plane123=gp_Pln( Som1, gp_Dir( N ) );
81       }
82     catch(Standard_Failure) {return(0.);}
83   }
84   double L1, L2, L3;
85   L1 = Som1.Distance(Som2);
86   L2 = gp_Lin(Som1, Line12).Distance(Som3);
87   L3 = Plane123.Distance(RefPoint);
88   
89   curVolume = ((L1 * L2)/2) * (L3/3);
90   
91   gp_Vec Rad(RefPoint, Som1);
92   
93   if( (Rad*N)>0 )
94     return (curVolume);
95   else
96     return (-curVolume);
97 }
98
99
100 //=======================================================================
101 //function : TetraCen
102 //purpose  : auxilary
103 //=======================================================================
104 static gp_XYZ TetraCen(gp_Pnt RefPoint, gp_Pnt Som1, gp_Pnt Som2, gp_Pnt Som3)
105 {
106   gp_XYZ curCentr, plnPnt;
107   plnPnt = ( Som1.XYZ() + Som2.XYZ() + Som3.XYZ() )/3;
108   curCentr = plnPnt + (RefPoint.XYZ() - plnPnt)/4;
109   return curCentr;
110 }
111
112
113 //=======================================================================
114 //function : CalculVolume
115 //purpose  : auxilary
116 //=======================================================================
117 static Standard_Real CalculVolume(const TopoDS_Shape& So,
118                                   gp_Pnt&       aRefPoint,
119                                   Standard_Real       tol,
120                                   Standard_Boolean withForce,
121                                   Draw_Interpretor& di)
122 {
123   Standard_Real myVolume = 0, curVolume = 0;
124   gp_XYZ localCentroid(0,0,0), curCentroid(0,0,0);
125   Standard_Boolean haveVertex = Standard_False;
126   for (TopExp_Explorer ex(So,TopAbs_FACE) ; ex.More(); ex.Next())
127     {
128       TopoDS_Face F =TopoDS::Face(ex.Current());
129       TopLoc_Location L;
130       if ( ! haveVertex )
131         for (TopExp_Explorer Vex(F, TopAbs_VERTEX); Vex.More(); Vex.Next() )
132           {
133             TopoDS_Vertex v = TopoDS::Vertex(Vex.Current());
134             if ( ! v.IsNull() ) {
135               aRefPoint = BRep_Tool::Pnt(v);
136               haveVertex = Standard_True;
137               break;
138             }
139           }
140       
141       Handle (Poly_Triangulation) facing = BRep_Tool::Triangulation(F,L);
142       if(facing.IsNull() || withForce)
143         {
144           BRepMesh_IncrementalMesh MESH(F, tol);
145           
146           facing = BRep_Tool::Triangulation(F,L);
147         }
148
149       TColgp_Array1OfPnt tab(1,(facing->NbNodes()));
150       tab = facing->Nodes();
151       Poly_Array1OfTriangle tri(1,facing->NbTriangles());
152       tri = facing->Triangles();
153       for (Standard_Integer i=1;i<=(facing->NbTriangles());i++)
154         {
155           
156           Poly_Triangle trian = tri.Value(i);
157           Standard_Integer index1,index2,index3;//M,N;
158           if( F.Orientation() == TopAbs_REVERSED )
159             trian.Get(index1,index3,index2);
160           else
161             trian.Get(index1,index2,index3);
162           curVolume = TetraVol(aRefPoint, tab.Value(index1),
163                                tab.Value(index2), tab.Value(index3));
164           myVolume += curVolume;
165           curCentroid = TetraCen(aRefPoint, tab.Value(index1),
166                                  tab.Value(index2), tab.Value(index3));
167           
168           localCentroid = localCentroid + curCentroid*curVolume;
169         }
170     }
171   
172   localCentroid = localCentroid * (1./ myVolume);
173   
174   di << "Centroid:\n";
175   di << "X=\t" << localCentroid.X() << "\n";
176   di << "Y=\t" << localCentroid.Y() << "\n";
177   di << "Z=\t" << localCentroid.Z() << "\n";
178   return (myVolume);
179 }
180
181
182 // --------------------- VolumeFix End   ---
183
184
185 //=======================================================================
186 // Section: Work with val props
187 //=======================================================================
188
189 //=======================================================================
190 //function : SetProps
191 //purpose  : 
192 //=======================================================================
193
194 static Standard_Integer SetProps (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
195 {
196   if (argc < 3) {
197     di<<"Use: "<<argv[0]<<" DocName {Shape|Label} [epsilon = 0.001]\n";
198     return 1;
199   }
200   Handle(TDocStd_Document) Doc;
201   DDocStd::GetDocument(argv[1], Doc);
202   if ( Doc.IsNull() ) { di << argv[1] << " is not a document\n"; return 1; }
203   
204   Standard_Real Vres, Ares;
205
206   TDF_Label aLabel;
207   TDF_Tool::Label(Doc->GetData(), argv[2], aLabel);
208   TopoDS_Shape aShape;
209   if ( aLabel.IsNull() ) {
210     aShape= DBRep::Get(argv[2]);
211     if ( !aShape.IsNull() ) {
212       Handle(XCAFDoc_ShapeTool) STool = XCAFDoc_DocumentTool::ShapeTool(Doc->Main());
213       aLabel = STool->FindShape(aShape);
214     }
215   }
216   else {
217     aShape = XCAFDoc_ShapeTool::GetShape ( aLabel );
218   }
219   if ( !aLabel.IsNull() ) {
220     
221     // retrieve epsilon
222     Standard_Real anEps;
223     if(argc > 3 ) anEps = Draw::Atof(argv[3]);
224     else anEps = 0.001;
225     
226     GProp_GProps G;
227     BRepGProp::VolumeProperties(aShape,G,anEps,Standard_True);
228     Vres = G.Mass();
229     Handle(XCAFDoc_Volume) aVolume = new XCAFDoc_Volume;
230     if (!aLabel.FindAttribute (XCAFDoc_Volume::GetID(), aVolume)) aLabel.AddAttribute(aVolume);
231     aVolume->Set(Vres);
232     
233     gp_Pnt aPoint = G.CentreOfMass();
234     Handle(XCAFDoc_Centroid) aCentroid = new XCAFDoc_Centroid;
235     if (!aLabel.FindAttribute (XCAFDoc_Centroid::GetID(), aCentroid)) aLabel.AddAttribute(aCentroid);
236     aCentroid->Set(aPoint);
237
238     BRepGProp::SurfaceProperties(aShape,G,anEps);
239     Ares = G.Mass();
240     Handle(XCAFDoc_Area) aArea = new XCAFDoc_Area;
241     if (!aLabel.FindAttribute (XCAFDoc_Area::GetID(), aArea)) aLabel.AddAttribute(aArea);
242     aArea->Set(Ares);
243     
244     di << argv[2] << ": Volume = " << Vres << ", Area = " << Ares << 
245       ", Centroid is (" << aPoint.X() << ", " << aPoint.Y() << ", " << aPoint.Z() << ")";
246   }
247   return 0;
248 }
249
250
251 //=======================================================================
252 //function : SetVolume
253 //purpose  : 
254 //=======================================================================
255
256 static Standard_Integer SetVolume (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
257 {
258   if (argc!=4) {
259     di<<"Use: "<<argv[0]<<" DocName {Label|Shape} volume\n";
260     return 1;
261   }
262   Handle(TDocStd_Document) Doc;
263   DDocStd::GetDocument(argv[1], Doc);
264   if ( Doc.IsNull() ) { di << argv[1] << " is not a document\n"; return 1; }
265   
266   Standard_Real res=0.;
267
268   TDF_Label aLabel;
269   TDF_Tool::Label(Doc->GetData(), argv[2], aLabel);
270   if ( aLabel.IsNull() ) {
271     TopoDS_Shape aShape= DBRep::Get(argv[2]);
272     if ( !aShape.IsNull() ) {
273       Handle(XCAFDoc_ShapeTool) STool = XCAFDoc_DocumentTool::ShapeTool(Doc->Main());
274       aLabel = STool->FindShape(aShape);
275     }
276   }
277   if ( !aLabel.IsNull() ) {
278     res = Draw::Atof(argv[3]);
279     Handle(XCAFDoc_Volume) aVolume = new XCAFDoc_Volume;
280     if (!aLabel.FindAttribute (XCAFDoc_Volume::GetID(), aVolume)) aLabel.AddAttribute(aVolume);
281     aVolume->Set(res);
282   }
283   
284   di << res;
285   return 0;
286 }
287
288
289 //=======================================================================
290 //function : SetArea
291 //purpose  : 
292 //=======================================================================
293
294 static Standard_Integer SetArea (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
295 {
296   if (argc!=4) {
297     di<<"Use: "<<argv[0]<<" DocName {Label|Shape} area\n";
298     return 1;
299   }
300   Handle(TDocStd_Document) Doc;
301   DDocStd::GetDocument(argv[1], Doc);
302   if ( Doc.IsNull() ) { di << argv[1] << " is not a document\n"; return 1; }
303   
304   Standard_Real res=0.;
305
306   TDF_Label aLabel;
307   TDF_Tool::Label(Doc->GetData(), argv[2], aLabel);
308   if ( aLabel.IsNull() ) {
309     TopoDS_Shape aShape= DBRep::Get(argv[2]);
310     if ( !aShape.IsNull() ) {
311       Handle(XCAFDoc_ShapeTool) STool = XCAFDoc_DocumentTool::ShapeTool(Doc->Main());
312       aLabel = STool->FindShape(aShape);
313     }
314   }
315   if ( !aLabel.IsNull() ) {
316     res = Draw::Atof(argv[3]);
317     Handle(XCAFDoc_Area) aArea = new XCAFDoc_Area;
318     if (!aLabel.FindAttribute (XCAFDoc_Area::GetID(), aArea)) aLabel.AddAttribute(aArea);
319     aArea->Set(res);
320   }
321   di << res;
322   return 0;
323 }
324
325
326 //=======================================================================
327 //function : SetCentroid
328 //purpose  : 
329 //=======================================================================
330
331 static Standard_Integer SetCentroid (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
332 {
333   if (argc!=6) {
334     di<<"Use: "<<argv[0]<<" DocName {Label|Shape} x y z\n";
335     return 1;
336   }
337   Handle(TDocStd_Document) Doc;
338   DDocStd::GetDocument(argv[1], Doc);
339   if ( Doc.IsNull() ) { di << argv[1] << " is not a document\n"; return 1; }
340   
341   gp_Pnt aPoint;
342
343   TDF_Label aLabel;
344   TDF_Tool::Label(Doc->GetData(), argv[2], aLabel);
345   if ( aLabel.IsNull() ) {
346     TopoDS_Shape aShape= DBRep::Get(argv[2]);
347     if ( !aShape.IsNull() ) {
348       Handle(XCAFDoc_ShapeTool) STool = XCAFDoc_DocumentTool::ShapeTool(Doc->Main());
349       aLabel = STool->FindShape(aShape);
350     }
351   }
352   if ( !aLabel.IsNull() ) {
353     aPoint.SetX(Draw::Atof(argv[3]));
354     aPoint.SetY(Draw::Atof(argv[4]));
355     aPoint.SetZ(Draw::Atof(argv[5]));
356     Handle(XCAFDoc_Centroid) aCentroid = new XCAFDoc_Centroid;
357     if (!aLabel.FindAttribute (XCAFDoc_Centroid::GetID(), aCentroid)) aLabel.AddAttribute(aCentroid);
358     aCentroid->Set(aPoint);
359     di << Draw::Atof(argv[3])<<" "<<Draw::Atof(argv[4])<<" "<<Draw::Atof(argv[5]);
360   }
361   return 0;
362 }
363
364
365 //=======================================================================
366 //function : GetVolume
367 //purpose  : 
368 //=======================================================================
369
370 static Standard_Integer GetVolume (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
371 {
372   if (argc!=3) {
373     di<<"Use: "<<argv[0]<<" DocName {Shape|Label}\n";
374     return 1;
375   }
376   Handle(TDocStd_Document) Doc;
377   DDocStd::GetDocument(argv[1], Doc);
378   if ( Doc.IsNull() ) { di << argv[1] << " is not a document\n"; return 1; }
379   
380   TDF_Label aLabel;
381   TDF_Tool::Label(Doc->GetData(), argv[2], aLabel);
382   if ( aLabel.IsNull() ) {
383     TopoDS_Shape aShape= DBRep::Get(argv[2]);
384     if ( !aShape.IsNull() ) {
385       Handle(XCAFDoc_ShapeTool) STool = XCAFDoc_DocumentTool::ShapeTool(Doc->Main());
386       aLabel = STool->FindShape(aShape);
387     }
388   }
389   if ( !aLabel.IsNull() ) {
390 //    Handle(XCAFDoc_Volume) aVolume = new (XCAFDoc_Volume);
391 //    if (aLabel.FindAttribute (XCAFDoc_Volume::GetID(), aVolume)) di << aVolume->Get();
392     // another case
393     Standard_Real aVol;
394     if(XCAFDoc_Volume::Get(aLabel, aVol))
395       di << aVol;
396   }
397   return 0;
398 }
399
400
401 //=======================================================================
402 //function : GetArea
403 //purpose  : 
404 //=======================================================================
405
406 static Standard_Integer GetArea (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
407 {
408   if (argc!=3) {
409     di<<"Use: "<<argv[0]<<" DocName {Shape|Label}\n";
410     return 1;
411   }
412   Handle(TDocStd_Document) Doc;
413   DDocStd::GetDocument(argv[1], Doc);
414   if ( Doc.IsNull() ) { di << argv[1] << " is not a document\n"; return 1; }
415   
416   TDF_Label aLabel;
417   TDF_Tool::Label(Doc->GetData(), argv[2], aLabel);
418   if ( aLabel.IsNull() ) {
419     TopoDS_Shape aShape= DBRep::Get(argv[2]);
420     if ( !aShape.IsNull() ) {
421       Handle(XCAFDoc_ShapeTool) STool = XCAFDoc_DocumentTool::ShapeTool(Doc->Main());
422       aLabel = STool->FindShape(aShape);
423     }
424   }
425   if ( !aLabel.IsNull() ) {
426 //     Handle(XCAFDoc_Area) aArea = new (XCAFDoc_Area);
427 //     if (aLabel.FindAttribute (XCAFDoc_Area::GetID(), aArea)) di << aArea->Get();
428     // another case
429     Standard_Real anA;
430     if(XCAFDoc_Area::Get(aLabel, anA))
431       di << anA;
432   }
433   return 0;
434 }
435
436
437 //=======================================================================
438 //function : GetCentroid
439 //purpose  : 
440 //=======================================================================
441
442 static Standard_Integer GetCentroid (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
443 {
444   if (argc!=3) {
445     di<<"Use: "<<argv[0]<<" DocName {Shape|Label} \n";
446     return 1;
447   }
448   Handle(TDocStd_Document) Doc;
449   DDocStd::GetDocument(argv[1], Doc);
450   if ( Doc.IsNull() ) { di << argv[1] << " is not a document\n"; return 1; }
451   
452   gp_Pnt aPoint;
453   
454   TDF_Label aLabel;
455   TDF_Tool::Label(Doc->GetData(), argv[2], aLabel);
456   if ( aLabel.IsNull() ) {
457     TopoDS_Shape aShape= DBRep::Get(argv[2]);
458     if ( !aShape.IsNull() ) {
459       Handle(XCAFDoc_ShapeTool) STool = XCAFDoc_DocumentTool::ShapeTool(Doc->Main());
460       aLabel = STool->FindShape(aShape);
461     }
462   }
463   if ( !aLabel.IsNull() ) {
464     Handle(XCAFDoc_Centroid) aCentroid = new (XCAFDoc_Centroid);
465     if (aLabel.FindAttribute (XCAFDoc_Centroid::GetID(), aCentroid)) {
466 //       aPoint = aCentroid->Get();
467 //       di << aPoint.X()<<" "<<aPoint.Y()<<" "<<aPoint.Z();
468       // another case
469       if(XCAFDoc_Centroid::Get(aLabel, aPoint))
470         di << aPoint.X()<<" "<<aPoint.Y()<<" "<<aPoint.Z();
471     }
472   }
473   return 0;
474 }
475
476 //=======================================================================
477 //function : CheckProps
478 //purpose  : 
479 //=======================================================================
480
481 static Standard_Integer CheckProps (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
482 {
483   if (argc <2) {
484     di << "Use: "<<argv[0]<<" DocName [ 0|deflection [Shape|Label] ]\n";
485     di << "     If second argument is 0, standard CADCADE tool is used for\n";
486     di << "     computation of volume and CG.\n";
487     di << "     If second argument is not 0, it is treated as deflection\n";
488     di << "     and computation is done by triangulations\n";
489     di << "     If the second argument is negative, meshing is forced\n";
490     return 1;
491   }
492   Handle(TDocStd_Document) Doc;
493   DDocStd::GetDocument(argv[1], Doc);
494   if ( Doc.IsNull() ) { di << argv[1] << " is not a document\n"; return 1; }
495   Standard_Boolean withVolFix = Standard_False;
496   if ( argc >2 && Draw::Atof(argv[2]) != 0 ) withVolFix = Standard_True;
497   Standard_Boolean wholeDoc = ( argc <4 );
498   TDF_LabelSequence seq;
499   if ( ! wholeDoc ) {
500     TDF_Label aLabel;
501     TDF_Tool::Label(Doc->GetData(), argv[3], aLabel);
502     TopoDS_Shape aShape;
503     if ( aLabel.IsNull() ) {
504       aShape= DBRep::Get(argv[3]);
505       if ( !aShape.IsNull() ) {
506         Handle(XCAFDoc_ShapeTool) STool = XCAFDoc_DocumentTool::ShapeTool(Doc->Main());
507         aLabel = STool->FindShape(aShape);
508       }
509     }
510     if ( aLabel.IsNull() ) return 1;
511     seq.Append ( aLabel );
512   }
513   else {
514     Handle(XCAFDoc_ShapeTool) STool = XCAFDoc_DocumentTool::ShapeTool(Doc->Main());
515     STool->GetShapes(seq);
516   }
517   if ( wholeDoc ) {
518     di << "Label            Area defect   Volume defect    dX      dY      dZ    Name\n";
519   }
520   for ( Standard_Integer i=1; i <= seq.Length(); i++ ) {
521     TDF_Label aLabel = seq(i);
522     
523     // add instance labels to sequence to process them as well
524     if ( XCAFDoc_ShapeTool::IsAssembly ( aLabel ) ) {
525       TDF_LabelSequence comp;
526       XCAFDoc_ShapeTool::GetComponents ( aLabel, comp );
527       Standard_Integer m=i;
528       for ( Standard_Integer k=1; k <= comp.Length(); k++ ) {
529         TDF_Label lab = comp(k);
530         Handle(XCAFDoc_Volume) aVolume;
531         Handle(XCAFDoc_Area) aArea;
532         Handle(XCAFDoc_Centroid) aCentroid;
533         if ( ! lab.FindAttribute (XCAFDoc_Volume::GetID(), aVolume) &&
534              ! lab.FindAttribute (XCAFDoc_Area::GetID(), aArea) &&
535              ! lab.FindAttribute (XCAFDoc_Centroid::GetID(), aCentroid) ) continue;
536         seq.InsertAfter ( m++, lab );
537       }
538     }
539
540     TCollection_AsciiString str;
541     TDF_Tool::Entry ( aLabel, str );
542     //printf ( "%s%-12.12s", ( wholeDoc ? "" : "Label " ), str.ToCString() );
543     //fflush ( stdout );
544     char string1[260];
545     Sprintf (string1, "%s%-12.12s", ( wholeDoc ? "" : "Label " ), str.ToCString() );
546     di << string1;
547     Handle(TDataStd_Name) N;
548     if ( aLabel.FindAttribute ( TDataStd_Name::GetID(), N ) && ! wholeDoc ) {
549       TCollection_AsciiString name(N->Get(), '?');
550       di << " \"" << name.ToCString() << "\"";
551     }
552     if ( ! wholeDoc ) di << "\n";
553     
554     Handle(XCAFDoc_Volume) aVolume;
555     Handle(XCAFDoc_Area) aArea;
556     Handle(XCAFDoc_Centroid) aCentroid;
557     aLabel.FindAttribute (XCAFDoc_Volume::GetID(), aVolume);
558     aLabel.FindAttribute (XCAFDoc_Area::GetID(), aArea);
559     aLabel.FindAttribute (XCAFDoc_Centroid::GetID(), aCentroid);
560     GProp_GProps G;
561
562     TopoDS_Shape aShape = XCAFDoc_ShapeTool::GetShape ( aLabel );
563     if ( ! aArea.IsNull() ) {
564       try { 
565         OCC_CATCH_SIGNALS
566         BRepGProp::SurfaceProperties(aShape,G,0.001);
567         //printf ("%s%9.1f (%3d%%)%s", ( wholeDoc ? "" : "  Area defect:   " ),
568         //      aArea->Get() - G.Mass(), 
569         //      (Standard_Integer)( Abs ( G.Mass() ) > 1e-10 ? 100. * ( aArea->Get() - G.Mass() ) / G.Mass() : 999. ),
570         //      ( wholeDoc ? "" : "\n" ));
571         char string2[260];
572         Sprintf (string2, "%s%9.1f (%3d%%)%s", ( wholeDoc ? "" : "  Area defect:   " ),
573                  aArea->Get() - G.Mass(), 
574                  (Standard_Integer)( Abs ( G.Mass() ) > 1e-10 ? 100. * ( aArea->Get() - G.Mass() ) / G.Mass() : 999. ),
575                  ( wholeDoc ? "" : "\n" ));
576         di << string2;
577       }
578       catch (Standard_Failure) {
579         //printf ( "%-16.16s", "exception" );
580         char string3[260];
581         Sprintf (string3, "%-16.16s", "exception" );
582         di << string3;
583       }
584     }
585     else if ( wholeDoc ) {
586       //printf ( "%16.16s", "" );
587       char string4[260];
588       Sprintf (string4, "%16.16s", "" );
589       di << string4;
590     }
591
592     if ( ! aVolume.IsNull() || ! aCentroid.IsNull() ) {
593       try {
594         OCC_CATCH_SIGNALS
595         // Added for check Volume. PTV 08 Nov 2000.     
596         Standard_Real localVolume;
597         gp_Pnt pcg(0,0,0);
598         if ( withVolFix ) {
599           Standard_Real tol = Draw::Atof(argv[2]);
600           Standard_Boolean withForce = Standard_False;
601           if ( tol < 0 ) {
602             withForce = Standard_True;
603             tol = -tol;
604           }
605           localVolume = CalculVolume(aShape, pcg, tol, withForce, di);
606         }
607         else {
608           BRepGProp::VolumeProperties(aShape,G,0.001,Standard_True);
609           localVolume = G.Mass();
610           pcg = G.CentreOfMass();
611         }
612         
613         if ( ! aVolume.IsNull() ) {
614           //printf ("%s%9.1f (%3d%%)%s", ( wholeDoc ? "" : "  Volume defect: " ),
615                 //  aVolume->Get() - localVolume,
616                 //  (Standard_Integer)( Abs ( localVolume ) > 1e-10 ? 100. * ( aVolume->Get() - localVolume ) / localVolume : 999. ),
617                 //  ( wholeDoc ? "" : "\n" ));
618           char string5[260];
619           Sprintf (string5, "%s%9.1f (%3d%%)%s", ( wholeDoc ? "" : "  Volume defect: " ),
620                    aVolume->Get() - localVolume,
621                    (Standard_Integer)( Abs ( localVolume ) > 1e-10 ? 100. * ( aVolume->Get() - localVolume ) / localVolume : 999. ),
622                    ( wholeDoc ? "" : "\n" ));
623           di << string5;
624         }
625         else if ( wholeDoc ) {
626           //printf ( "%16.16s", "" );
627           char string6[260];
628           Sprintf (string6, "%16.16s", "" );
629           di << string6;
630         }
631
632         if ( ! aCentroid.IsNull() ) {
633           gp_Pnt p = aCentroid->Get();
634           char string7[260];
635           if ( wholeDoc ) {
636             //printf ( " %7.2f %7.2f %7.2f", 
637                 //    p.X() - pcg.X(), p.Y() - pcg.Y(), p.Z() - pcg.Z() );
638             Sprintf (string7, " %7.2f %7.2f %7.2f", 
639                     p.X() - pcg.X(), p.Y() - pcg.Y(), p.Z() - pcg.Z() );
640           } else {
641             //printf ( "  CG defect: dX=%.3f, dY=%.3f, dZ=%.3f\n", 
642                 //    p.X() - pcg.X(), p.Y() - pcg.Y(), p.Z() - pcg.Z() );
643             Sprintf (string7, "  CG defect: dX=%.3f, dY=%.3f, dZ=%.3f\n", 
644                     p.X() - pcg.X(), p.Y() - pcg.Y(), p.Z() - pcg.Z() );
645           }
646           di << string7;
647         }
648         else if ( wholeDoc ) {
649           //printf ( "%24.24s", "" );
650           char string8[260];
651           Sprintf (string8, "%24.24s", "" );
652           di << string8;
653         }
654       }
655       catch (Standard_Failure const& anException) {
656         //printf ( "%40.40s", "exception" );
657         char string9[260];
658         Sprintf (string9, "%40.40s", "exception" );
659         di << string9;
660 #ifdef OCCT_DEBUG
661         //fflush ( stdout );
662         di << ": ";
663         di << anException.GetMessageString();
664         di<<" ** Skip\n";
665 #endif
666         (void)anException;
667       }
668     }
669     else if ( wholeDoc ) {
670       //printf ( "%40.40s", "" );
671       char string10[260];
672       Sprintf (string10, "%40.40s", "" );
673       di << string10;
674     }
675     //fflush ( stdout );
676
677     if ( wholeDoc ) {
678       if ( ! N.IsNull() ) {
679         TCollection_AsciiString name(N->Get(),'?');
680         di << "  \"" << name.ToCString() << "\"";
681       }
682       di << "\n";
683     }
684   }
685   return 0;
686 }
687
688
689 //=======================================================================
690 //function : ShapeVolume
691 //purpose  : 
692 //=======================================================================
693
694 static Standard_Integer ShapeVolume (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
695 {
696   if (argc!=3) {
697     di<<"Use: "<<argv[0]<<" Shape deflection \n";
698     return 1;
699   }
700   TopoDS_Shape aShape = DBRep::Get(argv[1]);
701   if (aShape.IsNull()) return 1;
702   gp_Pnt aPoint(0,0,0);
703   Standard_Real localVolume;
704   Standard_Real tol = Draw::Atof(argv[2]);
705   Standard_Boolean withForce = Standard_False;
706   if ( tol < 0 ) {
707     withForce = Standard_True;
708     tol = -tol;
709   }
710   localVolume = CalculVolume(aShape, aPoint, tol, withForce, di);
711   //cout << "Volume : " << setw(15) << localVolume << "\n" << endl;
712   Standard_SStream aSStream;
713   aSStream << "Volume : " << setw(15) << localVolume << "\n";
714   di << aSStream;
715   return 0;
716 }
717
718
719 //=======================================================================
720 //function : GetMassProps
721 //purpose  : auxilary for ShapeMassProps
722 //=======================================================================
723
724 static Standard_Boolean GetMassProps(const TDF_Label& aLabel, gp_XYZ& theCenterGravity,
725                                      Standard_Real& theMassVal, const Standard_Real thetol)
726 {
727   Standard_Real aDensity = XCAFDoc_MaterialTool::GetDensityForShape(aLabel);
728
729   if(aDensity >0) {
730     TopoDS_Shape aShape = XCAFDoc_ShapeTool::GetShape(aLabel);
731     GProp_GProps G;
732     BRepGProp::VolumeProperties(aShape,G,0.001,Standard_True);
733     Standard_Real localVolume = G.Mass();
734     theMassVal = aDensity*localVolume;
735     theCenterGravity = G.CentreOfMass().XYZ();
736     return Standard_True;
737   }
738
739   if(aDensity==0) {
740     Handle(TNaming_NamedShape) NS;
741     if(aLabel.FindAttribute(TNaming_NamedShape::GetID(),NS)) {
742       //S = TNaming_Tool::GetShape(NS);
743       TopoDS_Shape aSh = NS->Get();
744       if(aSh.ShapeType()==TopAbs_SOLID) return Standard_False;
745     }
746
747     //TopoDS_Shape aSh = XCAFDoc_ShapeTool::GetShape(aLabel);
748     //if(aSh.ShapeType()==TopAbs_SOLID) return Standard_False;
749
750     Handle(TDataStd_TreeNode) Node;
751     if( aLabel.FindAttribute(XCAFDoc::ShapeRefGUID(),Node) && Node->HasFather() ) {
752       TDF_Label SubL = Node->Father()->Label();
753       if(GetMassProps(SubL,theCenterGravity,theMassVal,thetol)) {
754         Handle(XCAFDoc_Location) LocationAttribute;
755         if(aLabel.FindAttribute(XCAFDoc_Location::GetID(),LocationAttribute)) {
756           gp_XYZ tmp = LocationAttribute->Get().Transformation().TranslationPart();
757           theCenterGravity += tmp;
758         }
759         return Standard_True;
760       }
761       else
762         return Standard_False;
763     }
764     else {
765       // calculate for components
766       TDF_LabelSequence comp;
767       XCAFDoc_ShapeTool::GetComponents ( aLabel, comp );
768       if(!comp.Length())
769         return Standard_False;
770   
771       TColgp_Array1OfXYZ anArrCentres(1,comp.Length());
772       TColStd_Array1OfReal anArrMass(1,comp.Length());
773       anArrMass.Init(0.0);
774       Standard_Real aTotalMass =0.0;
775       Standard_Integer k=1;
776       for ( ; k <= comp.Length(); k++ ) {
777         TDF_Label lab = comp(k);
778         gp_XYZ aCenterGravity(0.0,0.0,0.0);
779         Standard_Real aMassVal =0.0;
780         if(GetMassProps(lab,aCenterGravity,aMassVal,thetol)) {
781           anArrCentres.SetValue(k,aCenterGravity);
782           anArrMass.SetValue(k,aMassVal);
783           aTotalMass += aMassVal;
784         }
785       }
786       if(aTotalMass>0) {
787         Standard_Real x= 0.0, y =0.0, z= 0.0;
788         for( k=1; k <= anArrMass.Length(); k++) {
789           Standard_Real koef = anArrMass.Value(k)/aTotalMass;
790           x+= (anArrCentres.Value(k).X()*koef);
791           y+= (anArrCentres.Value(k).Y()*koef);
792           z+= (anArrCentres.Value(k).Z()*koef);
793         }
794         theMassVal = aTotalMass;
795         theCenterGravity.SetCoord(x,y,z);
796       }
797     }
798   }
799   return Standard_True;
800 }
801
802
803 //=======================================================================
804 //function : ShapeMassProps
805 //purpose  : 
806 //=======================================================================
807
808 static Standard_Integer ShapeMassProps (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
809 {
810   
811   if (argc <2) {
812     di << "Use: "<<argv[0]<<" DocName [deflection [Shape|Label] ]\n";
813     di << "     If second argument is 0, standard CADCADE tool is used for\n";
814     di << "     computation of volume and CG.\n";
815     di << "     If second argument is not 0, it is treated as deflection\n";
816     di << "     and computation is done by triangulations\n";
817     di << "     If the second argument is negative, meshing is forced\n";
818     return 1;
819   }
820   Handle(TDocStd_Document) Doc;
821   DDocStd::GetDocument(argv[1], Doc);
822   Standard_Real atol = Precision::Confusion();
823   if(argc >2)
824     atol  = Draw::Atof(argv[2]);
825   if ( Doc.IsNull() ) { di << argv[1] << " is not a document\n"; return 1; }
826   Standard_Boolean wholeDoc = ( argc <4 );
827   TDF_LabelSequence seq;
828   if ( ! wholeDoc ) {
829     TDF_Label aLabel;
830     TDF_Tool::Label(Doc->GetData(), argv[3], aLabel);
831     TopoDS_Shape aShape;
832     if ( aLabel.IsNull() ) {
833       aShape= DBRep::Get(argv[3]);
834       if ( !aShape.IsNull() ) {
835         Handle(XCAFDoc_ShapeTool) STool = XCAFDoc_DocumentTool::ShapeTool(Doc->Main());
836         aLabel = STool->FindShape(aShape);
837       }
838     }
839     if ( aLabel.IsNull() ) return 1;
840     seq.Append ( aLabel );
841   }
842   else {
843     Handle(XCAFDoc_ShapeTool) STool = XCAFDoc_DocumentTool::ShapeTool(Doc->Main());
844     STool->GetShapes(seq);
845   }
846   //if ( wholeDoc ) {
847   //  di << "Label            Area defect   Volume defect    dX      dY      dZ    Name\n";
848  // }
849   gp_XYZ aCenterGravity(0.0,0.0,0.0);
850   Standard_Real aMassVal =0.0;
851   for ( Standard_Integer i=1; i <= seq.Length(); i++ ) {
852     TDF_Label aLabel = seq(i);
853     GetMassProps(aLabel,aCenterGravity,aMassVal,atol);
854 //    if(GetMassProps(aLabel,aCenterGravity,aMassVal,atol))
855 //    {
856       TCollection_AsciiString str;
857       TDF_Tool::Entry ( aLabel, str );
858     if(aMassVal>0) {
859       di<<"Shape from label : "<<str.ToCString()<<"\n";
860       di<<"Mass = "<<aMassVal<<"\n";
861       di<<"CenterOfGravity X = "<<aCenterGravity.X()<<",Y = "<<aCenterGravity.Y()<<",Z = "<<aCenterGravity.Z()<<"\n";
862       di<<"\n";
863     }
864     else {
865 //      di<<"For one component density is absent\n";
866       di<<"Shape from label : "<<str.ToCString()<<" not have a mass\n";
867     }
868   }
869   return 0;
870 }
871
872
873 //=======================================================================
874 //function : SetMaterial
875 //purpose  : 
876 //=======================================================================
877
878 static Standard_Integer SetMaterial (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
879 {
880   if (argc<5) {
881     di<<"Use: "<<argv[0]<<" Doc {Label|Shape} name density(g/cu sm) \n";
882     return 1;
883   }
884   Handle(TDocStd_Document) Doc;   
885   DDocStd::GetDocument(argv[1], Doc);
886   if ( Doc.IsNull() ) { di << argv[1] << " is not a document\n"; return 1; }
887   
888   TDF_Label aLabel;
889   TDF_Tool::Label(Doc->GetData(), argv[2], aLabel);
890
891   Handle(XCAFDoc_MaterialTool) MatTool = XCAFDoc_DocumentTool::MaterialTool(Doc->Main());
892
893   MatTool->SetMaterial(aLabel, new TCollection_HAsciiString(argv[3]),
894                        new TCollection_HAsciiString(""), Draw::Atof(argv[4]),
895                        new TCollection_HAsciiString("density measure"),
896                        new TCollection_HAsciiString("POSITIVE_RATIO_MEASURE"));
897
898   return 0;
899 }
900
901 static Standard_Integer GetValidationProps(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
902 {
903   if (argc< 2) {
904     di<<"Use: "<<argv[0]<<" Doc {Label|Shape}\n";
905     return 1;
906   }
907   Handle(TDocStd_Document) Doc;   
908   DDocStd::GetDocument(argv[1], Doc);
909   if ( Doc.IsNull() ) { di << argv[1] << " is not a document\n"; return 1; }
910   Handle(XCAFDoc_ShapeTool) STool = XCAFDoc_DocumentTool::ShapeTool(Doc->Main());
911   TDF_LabelSequence aLabels;
912   if( argc == 2)
913     STool->GetShapes (aLabels);
914   else
915   {
916     TDF_Label aLabel;
917     TopoDS_Shape aShape = DBRep::Get(argv[2]);
918     if( aShape.IsNull())
919     {
920       TDF_Label aLabel;
921       TDF_Tool::Label(Doc->GetData(), argv[2], aLabel);
922     }
923     else
924     {
925       aLabel = STool->FindShape(aShape);
926     }
927     if( !aLabel.IsNull())
928       aLabels.Append(aLabel);
929   }
930    
931   
932   enum
933   {
934     Vol =0,
935     Area = 1,
936     Centroid = 2
937   };
938   Standard_Integer nbProps[3];
939   Standard_Integer j =0;
940   for( ; j <= Centroid; j++)
941     nbProps[j] = 0;
942   Standard_Integer i = 1;
943   for( ; i <= aLabels.Length() ; i++)
944   {
945     TDF_Label aLabel = aLabels(i);
946     
947     Standard_Real aProp[2];
948     
949     aProp[0] = 0.;
950     aProp[1] = 0.;
951     
952
953     gp_Pnt aP(Precision::Infinite(), Precision::Infinite(), Precision::Infinite());
954     XCAFDoc_Volume::Get(aLabel, aProp[Vol]);
955     XCAFDoc_Area::Get(aLabel, aProp[Area]);
956
957     Handle(XCAFDoc_Centroid) aCentroid = new (XCAFDoc_Centroid);
958     if (aLabel.FindAttribute (XCAFDoc_Centroid::GetID(), aCentroid)) 
959       XCAFDoc_Centroid::Get(aLabel, aP);
960         
961     if(aProp[Vol] > 0 || aProp[Area] > 0 || !Precision::IsInfinite(aP.X()) )
962     {
963       TCollection_AsciiString str;
964       TDF_Tool::Entry ( aLabel, str );
965       di<<"Label : "<<str;
966       
967       for(j = 0 ; j <= Area; j++)
968       {
969         if( aProp[j] > 0)
970         {
971           di<<(j == Vol ? "; Volume - " : "; Area - ")<<aProp[j];
972           nbProps[j]++;
973         }
974       }
975
976       if( !Precision::IsInfinite(aP.X()) )
977       {      
978         di<< "; Centroid -  "<< aP.X()<<" "<<aP.Y()<<" "<<aP.Z();
979         nbProps[Centroid]++;
980       }
981        di<<"\n";
982     }
983   }
984   di<<"========================================================="<<"\n";
985   di<< "Number of the validation properties : ";
986   for( i = Vol; i <= Centroid; i++)
987     di<< ( i == Vol ? "Volume" : (i == Area ? "Area" : "Centroid"))<<" : "<<nbProps[i]<<" ; "; 
988     
989   di<<"\n";
990   return 0;
991 }
992
993 //=======================================================================
994 //function : InitCommands
995 //purpose  : 
996 //=======================================================================
997
998 void XDEDRAW_Props::InitCommands(Draw_Interpretor& di) 
999 {
1000   static Standard_Boolean initactor = Standard_False;
1001   if (initactor)
1002   {
1003     return;
1004   }
1005   initactor = Standard_True;
1006
1007   Standard_CString g = "XDE property's commands";
1008   
1009   di.Add ("XSetVolume","DocName {Label|Shape} volume \t: Seting volume to shape",
1010                    __FILE__, SetVolume, g);
1011
1012   di.Add ("XGetVolume","DocName {Shape|Label} \t: Getting volume of shape",
1013                    __FILE__, GetVolume, g);
1014
1015   di.Add ("XSetArea","DocName {Label|Shape} area \t: Seting area to shape",
1016                    __FILE__, SetArea, g);
1017
1018   di.Add ("XGetArea","DocName {Shape|Label} \t: Getting area of shape",
1019                    __FILE__, GetArea, g);
1020
1021   di.Add ("XSetCentroid","DocName  {Label|Shape} x y z \t: Seting centroid to shape",
1022                    __FILE__, SetCentroid, g);
1023
1024   di.Add ("XGetCentroid","DocName {Shape|Label} \t: Getting centroid of shape ",
1025                    __FILE__, GetCentroid, g);
1026
1027   di.Add ("XSetProps","DocName {Shape|Label} [epsilon = 0.001] \t: Compute properties for shape ",
1028                    __FILE__, SetProps, g);
1029
1030   di.Add ("XCheckProps","DocName [ 0|deflection [Shape|Label] ]\t: Check properties recorded for shape ",
1031                    __FILE__, CheckProps, g);
1032   di.Add ("XShapeVolume","Shape \t: Calculating volume of shape",
1033                    __FILE__, ShapeVolume, g);
1034   di.Add ("XShapeMassProps","XShapeMassProps DocName [deflection [Shape|Label] ]\t: Get mass value and center of gravity for shape ",
1035                    __FILE__,ShapeMassProps , g);
1036   di.Add ("XSetMaterial","Doc {Label|Shape} name density(g/cu sm) \t: Set material to shape given by Label",
1037                    __FILE__, SetMaterial, g);
1038   di.Add ("XGetValProps","Doc {Label|Shape}",__FILE__, GetValidationProps, g);
1039  
1040 }