-// Copyright (C) 2007-2012 CEA/DEN, EDF R&D, OPEN CASCADE
+// Copyright (C) 2007-2016 CEA/DEN, EDF R&D, OPEN CASCADE
//
// Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
// CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
-// version 2.1 of the License.
+// version 2.1 of the License, or (at your option) any later version.
//
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// File : SMESH_Measurements_i.cxx
// Author : Pavel TELKOV, Open CASCADE S.A.S. (pavel.telkov@opencascade.com)
-#include "SMESH_Measurements_i.hxx"
+#ifdef WIN32
+#define NOMINMAX
+#endif
-#include "SMESH_Gen_i.hxx"
-#include "SMESH_PythonDump.hxx"
+#include "SMESH_Measurements_i.hxx"
+#include "SMDS_ElemIterator.hxx"
#include "SMDS_Mesh.hxx"
-#include "SMDS_MeshNode.hxx"
#include "SMDS_MeshElement.hxx"
-#include "SMDS_ElemIterator.hxx"
-
+#include "SMDS_MeshNode.hxx"
#include "SMESHDS_Mesh.hxx"
+#include "SMESH_Filter_i.hxx"
+#include "SMESH_Gen_i.hxx"
+#include "SMESH_MeshAlgos.hxx"
+#include "SMESH_PythonDump.hxx"
+#include <cmath>
using namespace SMESH;
return theTypes->length() > 0 && theTypes[0] == SMESH::NODE;
}
+static double getNumericalValue(SMESH::SMESH_IDSource_ptr theSource,
+ SMESH::Controls::NumericalFunctorPtr theFunctor)
+{
+ double value = 0;
+
+ if ( !CORBA::is_nil( theSource ) ) {
+ const SMESHDS_Mesh* aMesh = getMesh( theSource );
+ if ( aMesh ) {
+ theFunctor->SetMesh( aMesh );
+
+ SMESH::long_array_var anElementsId = theSource->GetIDs();
+ for ( CORBA::ULong i = 0; i < anElementsId->length(); i++) {
+ value += theFunctor->GetValue( anElementsId[i] );
+ }
+ }
+ }
+ return value;
+}
+
//=======================================================================
// name : MinDistance
// Purpose : minimal distance between two given entities
theMeasure.node1 = theNode->GetID();
}
else {
- theMeasure.minX = min( theMeasure.minX, theNode->X() );
- theMeasure.maxX = max( theMeasure.maxX, theNode->X() );
- theMeasure.minY = min( theMeasure.minY, theNode->Y() );
- theMeasure.maxY = max( theMeasure.maxY, theNode->Y() );
- theMeasure.minZ = min( theMeasure.minZ, theNode->Z() );
- theMeasure.maxZ = max( theMeasure.maxZ, theNode->Z() );
+ theMeasure.minX = std::min( theMeasure.minX, theNode->X() );
+ theMeasure.maxX = std::max( theMeasure.maxX, theNode->X() );
+ theMeasure.minY = std::min( theMeasure.minY, theNode->Y() );
+ theMeasure.maxY = std::max( theMeasure.maxY, theNode->Y() );
+ theMeasure.minZ = std::min( theMeasure.minZ, theNode->Z() );
+ theMeasure.maxZ = std::max( theMeasure.maxZ, theNode->Z() );
}
}
enlargeBoundingBox( aMesh->FindNode( aElementsId[i] ), theMeasure);
else
{
- const SMDS_MeshElement* elem = aMesh->FindElement( aElementsId[i] );
- if (!elem)
- continue;
- SMDS_ElemIteratorPtr aNodeIter = elem->nodesIterator();
- while( aNodeIter->more() )
- enlargeBoundingBox( dynamic_cast<const SMDS_MeshNode*>( aNodeIter->next() ), theMeasure);
+ if ( const SMDS_MeshElement* elem = aMesh->FindElement( aElementsId[i] ))
+ for (SMDS_NodeIteratorPtr aNodeIter = elem->nodeIterator(); aNodeIter->more(); )
+ enlargeBoundingBox( aNodeIter->next(), theMeasure);
}
}
}
-
+
//=======================================================================
// name : BoundingBox
// Purpose : compute common bounding box of entities
return aMeasure;
}
+
+//=======================================================================
+// name : Length
+// Purpose : sum of length of 1D elements of the source
+//=======================================================================
+double Measurements_i::Length(SMESH::SMESH_IDSource_ptr theSource)
+{
+ return getNumericalValue( theSource, SMESH::Controls::NumericalFunctorPtr(new SMESH::Controls::Length()) );
+}
+
+//=======================================================================
+// name : Area
+// Purpose : sum of area of 2D elements of the source
+//=======================================================================
+double Measurements_i::Area(SMESH::SMESH_IDSource_ptr theSource)
+{
+ return getNumericalValue( theSource, SMESH::Controls::NumericalFunctorPtr(new SMESH::Controls::Area()) );
+}
+
+//=======================================================================
+// name : Volume
+// Purpose : sum of volume of 3D elements of the source
+//=======================================================================
+double Measurements_i::Volume(SMESH::SMESH_IDSource_ptr theSource)
+{
+ return getNumericalValue( theSource, SMESH::Controls::NumericalFunctorPtr(new SMESH::Controls::Volume()) );
+}
+
+//=======================================================================
+//function : GravityCenter
+//purpose : return gravity center of the source: average coordinates of all nodes
+//=======================================================================
+
+SMESH::PointStruct Measurements_i::GravityCenter(SMESH::SMESH_IDSource_ptr theSource)
+{
+ SMESH::PointStruct grCenter = { 0.,0.,0. };
+ const SMESHDS_Mesh* mesh = getMesh( theSource );
+ if ( !mesh )
+ return grCenter;
+
+ // unmark all nodes; visited nodes will be marked
+ SMESH_MeshAlgos::MarkElems( mesh->nodesIterator(), /*isMarked=*/false );
+
+ gp_XYZ sumCoord( 0,0,0 );
+ int nodeCount = 0;
+
+ SMDS_ElemIteratorPtr eIt = SMESH_Mesh_i::GetElements( theSource, SMESH::ALL );
+ while ( eIt->more() )
+ {
+ const SMDS_MeshElement* elem = eIt->next();
+ for ( SMDS_NodeIteratorPtr nIt = elem->nodeIterator(); nIt->more(); )
+ {
+ const SMDS_MeshNode* n = nIt->next();
+ if ( !n->isMarked() )
+ {
+ sumCoord += SMESH_NodeXYZ( n );
+ ++nodeCount;
+ n->setIsMarked( true );
+ }
+ }
+ }
+ sumCoord /= nodeCount;
+
+ grCenter.x = sumCoord.X();
+ grCenter.y = sumCoord.Y();
+ grCenter.z = sumCoord.Z();
+
+ return grCenter;
+}
+
+//=======================================================================
+//function : Angle
+//purpose : Return angle in radians defined by 3 points <(p1,p2,p3)
+//=======================================================================
+
+CORBA::Double Measurements_i::Angle(const SMESH::PointStruct& p1,
+ const SMESH::PointStruct& p2,
+ const SMESH::PointStruct& p3 )
+{
+ gp_Vec v1( p1.x - p2.x, p1.y - p2.y, p1.z - p2.z );
+ gp_Vec v2( p3.x - p2.x, p3.y - p2.y, p3.z - p2.z );
+
+ double angle = -1;
+
+ try
+ {
+ angle = v1.Angle( v2 );
+ }
+ catch(...)
+ {
+ }
+ if ( std::isnan( angle ))
+ angle = -1;
+
+ return angle;
+}
