mirror of
https://git.dev.opencascade.org/repos/occt.git
synced 2025-04-10 18:51:21 +03:00
a4ed7309ff
Method BRepGProp::LinearProperties() is corrected to treat not geometric edges (their length is zero). To avoid wrong working of command nexplode global properties are calculated as point properties of vertexes. Small value Epsilon(1.) is used as density of points. This prevents exception on shapes that contain such edges, e.g. copy of a wire containing degenerated edge. Method BRep_Tool::IsGeometric() is optimized to avoid nested iteration for check of 3D curve for Null. Small bug is fixed in GProp_PGProps.cxx
472 lines
14 KiB
C++
472 lines
14 KiB
C++
// Copyright (c) 1995-1999 Matra Datavision
|
|
// Copyright (c) 1999-2014 OPEN CASCADE SAS
|
|
//
|
|
// This file is part of Open CASCADE Technology software library.
|
|
//
|
|
// This library is free software; you can redistribute it and/or modify it under
|
|
// the terms of the GNU Lesser General Public License version 2.1 as published
|
|
// by the Free Software Foundation, with special exception defined in the file
|
|
// OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
|
|
// distribution for complete text of the license and disclaimer of any warranty.
|
|
//
|
|
// Alternatively, this file may be used under the terms of Open CASCADE
|
|
// commercial license or contractual agreement.
|
|
|
|
#include <BRepGProp.hxx>
|
|
#include <BRepGProp_Cinert.hxx>
|
|
#include <BRepGProp_Sinert.hxx>
|
|
#include <BRepGProp_Vinert.hxx>
|
|
#include <BRepGProp_VinertGK.hxx>
|
|
#include <GProp_PGProps.hxx>
|
|
#include <BRepGProp_Face.hxx>
|
|
#include <BRepGProp_Domain.hxx>
|
|
#include <TopoDS.hxx>
|
|
#include <BRepAdaptor_Curve.hxx>
|
|
|
|
#include <TopTools.hxx>
|
|
#include <BRep_Tool.hxx>
|
|
#include <TopTools_ListOfShape.hxx>
|
|
#include <BRepCheck_Shell.hxx>
|
|
#include <TopTools_ListIteratorOfListOfShape.hxx>
|
|
#ifdef OCCT_DEBUG
|
|
static Standard_Integer AffichEps = 0;
|
|
#endif
|
|
|
|
static gp_Pnt roughBaryCenter(const TopoDS_Shape& S){
|
|
Standard_Integer i; TopExp_Explorer ex;
|
|
gp_XYZ xyz(0,0,0);
|
|
for (ex.Init(S,TopAbs_VERTEX), i = 0; ex.More(); ex.Next(), i++)
|
|
xyz += BRep_Tool::Pnt(TopoDS::Vertex(ex.Current())).XYZ();
|
|
if ( i > 0 ) xyz /= i;
|
|
return gp_Pnt(xyz);
|
|
}
|
|
|
|
void BRepGProp::LinearProperties(const TopoDS_Shape& S, GProp_GProps& SProps){
|
|
// find the origin
|
|
gp_Pnt P(0,0,0);
|
|
P.Transform(S.Location());
|
|
SProps = GProp_GProps(P);
|
|
|
|
BRepAdaptor_Curve BAC;
|
|
Standard_Real eps = Epsilon(1.);
|
|
TopExp_Explorer ex;
|
|
for (ex.Init(S,TopAbs_EDGE); ex.More(); ex.Next()) {
|
|
const TopoDS_Edge& aE = TopoDS::Edge(ex.Current());
|
|
if(!BRep_Tool::IsGeometric(aE))
|
|
{
|
|
GProp_PGProps aPProps;
|
|
TopoDS_Iterator anIter(aE);
|
|
for(; anIter.More(); anIter.Next())
|
|
{
|
|
const TopoDS_Vertex& aV = TopoDS::Vertex(anIter.Value());
|
|
aPProps.AddPoint(BRep_Tool::Pnt(aV), eps);
|
|
}
|
|
SProps.Add(aPProps);
|
|
}
|
|
else
|
|
{
|
|
BAC.Initialize(aE);
|
|
BRepGProp_Cinert CG(BAC,P);
|
|
SProps.Add(CG);
|
|
}
|
|
}
|
|
}
|
|
|
|
static Standard_Real surfaceProperties(const TopoDS_Shape& S, GProp_GProps& Props, const Standard_Real Eps){
|
|
Standard_Integer i;
|
|
#ifdef OCCT_DEBUG
|
|
Standard_Integer iErrorMax = 0;
|
|
#endif
|
|
Standard_Real ErrorMax = 0.0, Error;
|
|
TopExp_Explorer ex;
|
|
gp_Pnt P(roughBaryCenter(S));
|
|
BRepGProp_Sinert G; G.SetLocation(P);
|
|
|
|
BRepGProp_Face BF;
|
|
BRepGProp_Domain BD;
|
|
|
|
for (ex.Init(S,TopAbs_FACE), i = 1; ex.More(); ex.Next(), i++) {
|
|
const TopoDS_Face& F = TopoDS::Face(ex.Current());
|
|
BF.Load(F);
|
|
TopoDS_Iterator aWIter(F);
|
|
Standard_Boolean IsNatRestr = !aWIter.More();
|
|
if(!IsNatRestr) BD.Init(F);
|
|
if(Eps < 1.0) {
|
|
G.Perform(BF, BD, Eps);
|
|
Error = G.GetEpsilon();
|
|
if(ErrorMax < Error) {
|
|
ErrorMax = Error;
|
|
#ifdef OCCT_DEBUG
|
|
iErrorMax = i;
|
|
#endif
|
|
}
|
|
} else {
|
|
if(IsNatRestr) G.Perform(BF);
|
|
else G.Perform(BF, BD);
|
|
}
|
|
Props.Add(G);
|
|
#ifdef OCCT_DEBUG
|
|
if(AffichEps) cout<<"\n"<<i<<":\tEpsArea = "<< G.GetEpsilon();
|
|
#endif
|
|
}
|
|
#ifdef OCCT_DEBUG
|
|
if(AffichEps) cout<<"\n-----------------\n"<<iErrorMax<<":\tMaxError = "<<ErrorMax<<"\n";
|
|
#endif
|
|
return ErrorMax;
|
|
}
|
|
void BRepGProp::SurfaceProperties(const TopoDS_Shape& S, GProp_GProps& Props){
|
|
// find the origin
|
|
gp_Pnt P(0,0,0);
|
|
P.Transform(S.Location());
|
|
Props = GProp_GProps(P);
|
|
surfaceProperties(S,Props,1.0);
|
|
}
|
|
Standard_Real BRepGProp::SurfaceProperties(const TopoDS_Shape& S, GProp_GProps& Props, const Standard_Real Eps){
|
|
// find the origin
|
|
gp_Pnt P(0,0,0); P.Transform(S.Location());
|
|
Props = GProp_GProps(P);
|
|
Standard_Real ErrorMax = surfaceProperties(S,Props,Eps);
|
|
return ErrorMax;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : volumeProperties
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
static Standard_Real volumeProperties(const TopoDS_Shape& S, GProp_GProps& Props, const Standard_Real Eps){
|
|
Standard_Integer i;
|
|
#ifdef OCCT_DEBUG
|
|
Standard_Integer iErrorMax = 0;
|
|
#endif
|
|
Standard_Real ErrorMax = 0.0, Error = 0.0;
|
|
TopExp_Explorer ex;
|
|
gp_Pnt P(roughBaryCenter(S));
|
|
BRepGProp_Vinert G; G.SetLocation(P);
|
|
|
|
BRepGProp_Face BF;
|
|
BRepGProp_Domain BD;
|
|
|
|
for (ex.Init(S,TopAbs_FACE), i = 1; ex.More(); ex.Next(), i++) {
|
|
const TopoDS_Face& F = TopoDS::Face(ex.Current());
|
|
if ((F.Orientation() == TopAbs_FORWARD) || (F.Orientation() == TopAbs_REVERSED)){
|
|
BF.Load(F);
|
|
TopoDS_Iterator aWIter(F);
|
|
Standard_Boolean IsNatRestr = !aWIter.More();
|
|
if(!IsNatRestr) BD.Init(F);
|
|
if(Eps < 1.0) {
|
|
G.Perform(BF, BD, Eps);
|
|
Error = G.GetEpsilon();
|
|
if(ErrorMax < Error) {
|
|
ErrorMax = Error;
|
|
#ifdef OCCT_DEBUG
|
|
iErrorMax = i;
|
|
#endif
|
|
}
|
|
}
|
|
else {
|
|
if(IsNatRestr) G.Perform(BF);
|
|
else G.Perform(BF, BD);
|
|
}
|
|
Props.Add(G);
|
|
#ifdef OCCT_DEBUG
|
|
if(AffichEps) cout<<"\n"<<i<<":\tEpsVolume = "<< G.GetEpsilon();
|
|
#endif
|
|
}
|
|
}
|
|
#ifdef OCCT_DEBUG
|
|
if(AffichEps) cout<<"\n-----------------\n"<<iErrorMax<<":\tMaxError = "<<ErrorMax<<"\n";
|
|
#endif
|
|
return ErrorMax;
|
|
}
|
|
void BRepGProp::VolumeProperties(const TopoDS_Shape& S, GProp_GProps& Props, const Standard_Boolean OnlyClosed){
|
|
// find the origin
|
|
gp_Pnt P(0,0,0); P.Transform(S.Location());
|
|
Props = GProp_GProps(P);
|
|
if(OnlyClosed){
|
|
TopExp_Explorer ex(S,TopAbs_SHELL);
|
|
for (; ex.More(); ex.Next()) {
|
|
const TopoDS_Shape& Sh = ex.Current();
|
|
if(BRep_Tool::IsClosed(Sh)) volumeProperties(Sh,Props,1.0);
|
|
}
|
|
} else volumeProperties(S,Props,1.0);
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : VolumeProperties
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
Standard_Real BRepGProp::VolumeProperties(const TopoDS_Shape& S, GProp_GProps& Props,
|
|
const Standard_Real Eps, const Standard_Boolean OnlyClosed)
|
|
{
|
|
// find the origin
|
|
gp_Pnt P(0,0,0); P.Transform(S.Location());
|
|
Props = GProp_GProps(P);
|
|
Standard_Integer i;
|
|
#ifdef OCCT_DEBUG
|
|
Standard_Integer iErrorMax = 0;
|
|
#endif
|
|
Standard_Real ErrorMax = 0.0, Error = 0.0;
|
|
if(OnlyClosed){
|
|
TopExp_Explorer ex(S,TopAbs_SHELL);
|
|
for (i = 1; ex.More(); ex.Next(), i++) {
|
|
const TopoDS_Shape& Sh = ex.Current();
|
|
if(BRep_Tool::IsClosed(Sh)) {
|
|
Error = volumeProperties(Sh,Props,Eps);
|
|
if(ErrorMax < Error) {
|
|
ErrorMax = Error;
|
|
#ifdef OCCT_DEBUG
|
|
iErrorMax = i;
|
|
#endif
|
|
}
|
|
}
|
|
}
|
|
} else ErrorMax = volumeProperties(S,Props,Eps);
|
|
#ifdef OCCT_DEBUG
|
|
if(AffichEps) cout<<"\n\n==================="<<iErrorMax<<":\tMaxEpsVolume = "<<ErrorMax<<"\n";
|
|
#endif
|
|
return ErrorMax;
|
|
}
|
|
|
|
//===========================================================================================//
|
|
// Volume properties by Gauss-Kronrod integration
|
|
//===========================================================================================//
|
|
//=======================================================================
|
|
//function : VolumePropertiesGK
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
static Standard_Real volumePropertiesGK(const TopoDS_Shape &theShape,
|
|
GProp_GProps &theProps,
|
|
const Standard_Real theTol,
|
|
const Standard_Boolean IsUseSpan,
|
|
const Standard_Boolean CGFlag,
|
|
const Standard_Boolean IFlag)
|
|
{
|
|
TopExp_Explorer anExp;
|
|
anExp.Init(theShape, TopAbs_FACE);
|
|
|
|
Standard_Real aTol = theTol;
|
|
|
|
// Compute properties.
|
|
gp_Pnt aLoc(roughBaryCenter(theShape));
|
|
BRepGProp_VinertGK aVProps;
|
|
BRepGProp_Face aPropFace(IsUseSpan);
|
|
BRepGProp_Domain aPropDomain;
|
|
Standard_Real aLocalError;
|
|
Standard_Real anError = 0.;
|
|
|
|
aVProps.SetLocation(aLoc);
|
|
|
|
for (; anExp.More(); anExp.Next()) {
|
|
TopoDS_Face aFace = TopoDS::Face(anExp.Current());
|
|
|
|
if (aFace.Orientation() == TopAbs_FORWARD ||
|
|
aFace.Orientation() == TopAbs_REVERSED) {
|
|
aPropFace.Load(aFace);
|
|
|
|
TopoDS_Iterator aWIter(aFace);
|
|
Standard_Boolean IsNatRestr = !aWIter.More();
|
|
if(IsNatRestr)
|
|
aLocalError = aVProps.Perform(aPropFace, aTol, CGFlag, IFlag);
|
|
else {
|
|
aPropDomain.Init(aFace);
|
|
aLocalError = aVProps.Perform(aPropFace, aPropDomain, aTol, CGFlag, IFlag);
|
|
}
|
|
|
|
if (aLocalError < 0.)
|
|
return aLocalError;
|
|
|
|
anError += aLocalError;
|
|
theProps.Add(aVProps);
|
|
}
|
|
}
|
|
|
|
return anError;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : VolumePropertiesGK
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
Standard_Real BRepGProp::VolumePropertiesGK(const TopoDS_Shape &S,
|
|
GProp_GProps &Props,
|
|
const Standard_Real Eps,
|
|
const Standard_Boolean OnlyClosed,
|
|
const Standard_Boolean IsUseSpan,
|
|
const Standard_Boolean CGFlag,
|
|
const Standard_Boolean IFlag)
|
|
{
|
|
gp_Pnt P(0,0,0);
|
|
Standard_Real anError = 0.;
|
|
|
|
P.Transform(S.Location());
|
|
Props = GProp_GProps(P);
|
|
|
|
if(OnlyClosed) {
|
|
// To select closed shells.
|
|
TopExp_Explorer anExp;
|
|
TopTools_ListOfShape aClosedShells;
|
|
|
|
anExp.Init(S, TopAbs_SHELL);
|
|
|
|
for (; anExp.More(); anExp.Next()) {
|
|
const TopoDS_Shape &aShell = anExp.Current();
|
|
|
|
BRepCheck_Shell aChecker(TopoDS::Shell(aShell));
|
|
BRepCheck_Status aStatus = aChecker.Closed(Standard_False);
|
|
|
|
if(aStatus == BRepCheck_NoError)
|
|
aClosedShells.Append(aShell);
|
|
|
|
}
|
|
|
|
if (aClosedShells.IsEmpty())
|
|
return -1.;
|
|
|
|
// Compute the properties for each closed shell.
|
|
Standard_Real aTol = Eps;
|
|
Standard_Real aLocalError;
|
|
TopTools_ListIteratorOfListOfShape anIter(aClosedShells);
|
|
|
|
for (; anIter.More(); anIter.Next()) {
|
|
const TopoDS_Shape &aShell = anIter.Value();
|
|
|
|
aLocalError = volumePropertiesGK(aShell, Props, aTol, IsUseSpan, CGFlag, IFlag);
|
|
|
|
if (aLocalError < 0)
|
|
return aLocalError;
|
|
|
|
anError += aLocalError;
|
|
}
|
|
|
|
} else
|
|
anError = volumePropertiesGK(S, Props, Eps, IsUseSpan, CGFlag, IFlag);
|
|
|
|
Standard_Real vol = Props.Mass();
|
|
if(vol > Epsilon(1.)) anError /= vol;
|
|
return anError;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : VolumeProperties
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
static Standard_Real volumePropertiesGK(const TopoDS_Shape &theShape,
|
|
GProp_GProps &theProps,
|
|
const gp_Pln &thePln,
|
|
const Standard_Real theTol,
|
|
const Standard_Boolean IsUseSpan,
|
|
const Standard_Boolean CGFlag,
|
|
const Standard_Boolean IFlag)
|
|
{
|
|
TopExp_Explorer anExp;
|
|
anExp.Init(theShape, TopAbs_FACE);
|
|
|
|
Standard_Real aTol = theTol;
|
|
|
|
// Compute properties.
|
|
gp_Pnt aLoc(roughBaryCenter(theShape));
|
|
BRepGProp_VinertGK aVProps;
|
|
BRepGProp_Face aPropFace(IsUseSpan);
|
|
BRepGProp_Domain aPropDomain;
|
|
Standard_Real aLocalError;
|
|
Standard_Real anError = 0.;
|
|
|
|
aVProps.SetLocation(aLoc);
|
|
|
|
for (; anExp.More(); anExp.Next()) {
|
|
TopoDS_Face aFace = TopoDS::Face(anExp.Current());
|
|
|
|
if (aFace.Orientation() == TopAbs_FORWARD ||
|
|
aFace.Orientation() == TopAbs_REVERSED) {
|
|
aPropFace.Load(aFace);
|
|
|
|
TopoDS_Iterator aWIter(aFace);
|
|
Standard_Boolean IsNatRestr = !aWIter.More();
|
|
if(IsNatRestr)
|
|
aLocalError = aVProps.Perform(aPropFace, thePln, aTol, CGFlag, IFlag);
|
|
else {
|
|
aPropDomain.Init(aFace);
|
|
aLocalError = aVProps.Perform(aPropFace, aPropDomain, thePln, aTol, CGFlag, IFlag);
|
|
}
|
|
|
|
if (aLocalError < 0.)
|
|
return aLocalError;
|
|
|
|
anError += aLocalError;
|
|
theProps.Add(aVProps);
|
|
}
|
|
}
|
|
|
|
return anError;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : VolumeProperties
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
Standard_Real BRepGProp::VolumePropertiesGK(const TopoDS_Shape &S,
|
|
GProp_GProps &Props,
|
|
const gp_Pln &thePln,
|
|
const Standard_Real Eps,
|
|
const Standard_Boolean OnlyClosed,
|
|
const Standard_Boolean IsUseSpan,
|
|
const Standard_Boolean CGFlag,
|
|
const Standard_Boolean IFlag)
|
|
{
|
|
gp_Pnt P(0,0,0);
|
|
Standard_Real anError = 0.;
|
|
|
|
P.Transform(S.Location());
|
|
Props = GProp_GProps(P);
|
|
|
|
if(OnlyClosed) {
|
|
// To select closed shells.
|
|
TopExp_Explorer anExp;
|
|
TopTools_ListOfShape aClosedShells;
|
|
|
|
anExp.Init(S, TopAbs_SHELL);
|
|
|
|
for (; anExp.More(); anExp.Next()) {
|
|
const TopoDS_Shape &aShell = anExp.Current();
|
|
|
|
BRepCheck_Shell aChecker(TopoDS::Shell(aShell));
|
|
BRepCheck_Status aStatus = aChecker.Closed(Standard_False);
|
|
|
|
if(aStatus == BRepCheck_NoError)
|
|
aClosedShells.Append(aShell);
|
|
|
|
}
|
|
|
|
if (aClosedShells.IsEmpty())
|
|
return -1.;
|
|
|
|
// Compute the properties for each closed shell.
|
|
Standard_Real aTol = Eps;
|
|
Standard_Real aLocalError;
|
|
TopTools_ListIteratorOfListOfShape anIter(aClosedShells);
|
|
|
|
for (; anIter.More(); anIter.Next()) {
|
|
const TopoDS_Shape &aShell = anIter.Value();
|
|
|
|
aLocalError = volumePropertiesGK(aShell, Props, thePln, aTol, IsUseSpan, CGFlag, IFlag);
|
|
|
|
if (aLocalError < 0)
|
|
return aLocalError;
|
|
|
|
anError += aLocalError;
|
|
}
|
|
} else
|
|
anError = volumePropertiesGK(S, Props, thePln, Eps, IsUseSpan, CGFlag, IFlag);
|
|
|
|
Standard_Real vol = Props.Mass();
|
|
if(vol > Epsilon(1.)) anError /= vol;
|
|
|
|
return anError;
|
|
}
|