mirror of
https://git.dev.opencascade.org/repos/occt.git
synced 2025-04-21 10:13:43 +03:00
9f091145df
blend/complex/D4. Calculation of gravity center was removed for degenerated edge. It is revealed that function nexplode does not work on shapes with degenerated edges. Comments. Calculation of gravity center was made by vertices data in case of degenerated edge without representations especially for test case blend/complex/D4. The mass of the edge was set to 1. The mass near 0 does not satisfy the test case. As a result test case bugs/moddata/bug268 became failed. Because only the mass near 0 satisfies last test case.
440 lines
13 KiB
C++
Executable File
440 lines
13 KiB
C++
Executable File
// Copyright (c) 1995-1999 Matra Datavision
|
|
// Copyright (c) 1999-2012 OPEN CASCADE SAS
|
|
//
|
|
// The content of this file is subject to the Open CASCADE Technology Public
|
|
// License Version 6.5 (the "License"). You may not use the content of this file
|
|
// except in compliance with the License. Please obtain a copy of the License
|
|
// at http://www.opencascade.org and read it completely before using this file.
|
|
//
|
|
// The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
|
|
// main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
|
|
//
|
|
// The Original Code and all software distributed under the License is
|
|
// distributed on an "AS IS" basis, without warranty of any kind, and the
|
|
// Initial Developer hereby disclaims all such warranties, including without
|
|
// limitation, any warranties of merchantability, fitness for a particular
|
|
// purpose or non-infringement. Please see the License for the specific terms
|
|
// and conditions governing the rights and limitations under the License.
|
|
|
|
|
|
#include <BRepGProp.hxx>
|
|
#include <BRepGProp_Cinert.hxx>
|
|
#include <BRepGProp_Sinert.hxx>
|
|
#include <BRepGProp_Vinert.hxx>
|
|
#include <BRepGProp_VinertGK.hxx>
|
|
#include <BRepGProp_Face.hxx>
|
|
#include <BRepGProp_Domain.hxx>
|
|
#include <GProp_PGProps.hxx>
|
|
#include <TopoDS.hxx>
|
|
#include <TopoDS_Vertex.hxx>
|
|
#include <BRepAdaptor_Curve.hxx>
|
|
|
|
#include <TopTools.hxx>
|
|
#include <BRep_Tool.hxx>
|
|
#include <TopTools_ListOfShape.hxx>
|
|
#include <BRepCheck_Shell.hxx>
|
|
#include <TopTools_ListIteratorOfListOfShape.hxx>
|
|
#include <BRep_TEdge.hxx>
|
|
#include <BRep_TVertex.hxx>
|
|
#ifdef DEB
|
|
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_Integer n,i;
|
|
TopExp_Explorer ex;
|
|
for (ex.Init(S,TopAbs_EDGE); ex.More(); ex.Next()) {
|
|
TopoDS_Edge anES = TopoDS::Edge(ex.Current());
|
|
Handle_BRep_TEdge & anEG = (Handle_BRep_TEdge &)anES.TShape();
|
|
BRep_TEdge dsd;
|
|
if (!anEG->Degenerated())
|
|
{
|
|
BAC.Initialize(anES);
|
|
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, iErrorMax = 0;
|
|
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);
|
|
if(!BF.NaturalRestriction()) BD.Init(F);
|
|
if(Eps < 1.0) {
|
|
G.Perform(BF, BD, Eps);
|
|
Error = G.GetEpsilon();
|
|
if(ErrorMax < Error){ ErrorMax = Error; iErrorMax = i;}
|
|
} else {
|
|
if(BF.NaturalRestriction()) G.Perform(BF);
|
|
else G.Perform(BF, BD);
|
|
}
|
|
Props.Add(G);
|
|
#ifdef DEB
|
|
if(AffichEps) cout<<"\n"<<i<<":\tEpsArea = "<< G.GetEpsilon();
|
|
#endif
|
|
}
|
|
#ifdef DEB
|
|
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, iErrorMax = 0;
|
|
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);
|
|
if(!BF.NaturalRestriction()) BD.Init(F);
|
|
if(Eps < 1.0) {
|
|
G.Perform(BF, BD, Eps);
|
|
Error = G.GetEpsilon();
|
|
if(ErrorMax < Error){ ErrorMax = Error; iErrorMax = i;}
|
|
}
|
|
else {
|
|
if(BF.NaturalRestriction()) G.Perform(BF);
|
|
else G.Perform(BF, BD);
|
|
}
|
|
Props.Add(G);
|
|
#ifdef DEB
|
|
if(AffichEps) cout<<"\n"<<i<<":\tEpsVolume = "<< G.GetEpsilon();
|
|
#endif
|
|
}
|
|
}
|
|
#ifdef DEB
|
|
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);
|
|
Standard_Real Error = 0.0;
|
|
if(OnlyClosed){
|
|
TopExp_Explorer ex(S,TopAbs_SHELL);
|
|
for (; ex.More(); ex.Next()) {
|
|
const TopoDS_Shape& Sh = ex.Current();
|
|
if(BRep_Tool::IsClosed(Sh)) Error = volumeProperties(Sh,Props,1.0);
|
|
}
|
|
} else Error = 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, iErrorMax = 0;
|
|
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; iErrorMax = i;}
|
|
}
|
|
}
|
|
} else ErrorMax = volumeProperties(S,Props,Eps);
|
|
#ifdef DEB
|
|
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);
|
|
|
|
if(aPropFace.NaturalRestriction())
|
|
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);
|
|
|
|
if(aPropFace.NaturalRestriction())
|
|
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;
|
|
}
|