mirror of
https://git.salome-platform.org/gitpub/modules/smesh.git
synced 2025-03-15 18:31:27 +05:00
2270 lines
94 KiB
C++
2270 lines
94 KiB
C++
// Copyright (C) 2007-2020 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, 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
|
|
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
// Lesser General Public License for more details.
|
|
//
|
|
// You should have received a copy of the GNU Lesser General Public
|
|
// License along with this library; if not, write to the Free Software
|
|
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
|
|
//
|
|
// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
|
|
//
|
|
|
|
// SMESH SMESHDS : management of mesh data and SMESH document
|
|
// File : SMESH_Mesh.cxx
|
|
// Author : Yves FRICAUD, OCC
|
|
// Module : SMESH
|
|
//
|
|
#include "SMESHDS_Mesh.hxx"
|
|
|
|
#include "SMDS_Downward.hxx"
|
|
#include "SMDS_EdgePosition.hxx"
|
|
#include "SMDS_FacePosition.hxx"
|
|
#include "SMDS_SpacePosition.hxx"
|
|
#include "SMDS_VertexPosition.hxx"
|
|
#include "SMESHDS_Group.hxx"
|
|
#include "SMESHDS_GroupOnGeom.hxx"
|
|
#include "SMESHDS_Script.hxx"
|
|
#include "SMESHDS_TSubMeshHolder.hxx"
|
|
|
|
#include <Standard_ErrorHandler.hxx>
|
|
#include <Standard_OutOfRange.hxx>
|
|
#include <TopExp.hxx>
|
|
#include <TopExp_Explorer.hxx>
|
|
#include <TopoDS_Edge.hxx>
|
|
#include <TopoDS_Face.hxx>
|
|
#include <TopoDS_Iterator.hxx>
|
|
#include <TopoDS_Shell.hxx>
|
|
#include <TopoDS_Solid.hxx>
|
|
#include <TopoDS_Vertex.hxx>
|
|
|
|
#include "utilities.h"
|
|
|
|
class SMESHDS_Mesh::SubMeshHolder : public SMESHDS_TSubMeshHolder< const SMESHDS_SubMesh >
|
|
{
|
|
};
|
|
|
|
//=======================================================================
|
|
//function : Create
|
|
//purpose :
|
|
//=======================================================================
|
|
SMESHDS_Mesh::SMESHDS_Mesh(int theMeshID, bool theIsEmbeddedMode):
|
|
myMeshID(theMeshID),
|
|
mySubMeshHolder( new SubMeshHolder ),
|
|
myIsEmbeddedMode(theIsEmbeddedMode)
|
|
{
|
|
myScript = new SMESHDS_Script(theIsEmbeddedMode);
|
|
SetPersistentId(theMeshID);
|
|
}
|
|
|
|
//=======================================================================
|
|
bool SMESHDS_Mesh::IsEmbeddedMode()
|
|
{
|
|
return myIsEmbeddedMode;
|
|
}
|
|
|
|
//================================================================================
|
|
/*!
|
|
* \brief Store ID persistent during lifecycle
|
|
*
|
|
* Initially it was used to have a persistent reference to the mesh from the hypothesis
|
|
*/
|
|
//================================================================================
|
|
|
|
void SMESHDS_Mesh::SetPersistentId(int id)
|
|
{
|
|
if (NbNodes() == 0)
|
|
myPersistentID = id;
|
|
}
|
|
//================================================================================
|
|
/*!
|
|
* \brief Return ID persistent during lifecycle
|
|
*/
|
|
//================================================================================
|
|
|
|
int SMESHDS_Mesh::GetPersistentId() const
|
|
{
|
|
return myPersistentID;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : ShapeToMesh
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::ShapeToMesh(const TopoDS_Shape & S)
|
|
{
|
|
if ( !myShape.IsNull() && S.IsNull() ) // case: "save study" after geometry removal
|
|
{
|
|
// removal of a shape to mesh, delete ...
|
|
// - hypotheses
|
|
myShapeToHypothesis.Clear();
|
|
// - shape indices in SMDS_Position of nodes
|
|
SMESHDS_SubMeshIteratorPtr smIt = SubMeshes();
|
|
while ( SMESHDS_SubMesh* sm = const_cast< SMESHDS_SubMesh* >( smIt->next() )) {
|
|
if ( !sm->IsComplexSubmesh() ) {
|
|
SMDS_NodeIteratorPtr nIt = sm->GetNodes();
|
|
while ( nIt->more() )
|
|
sm->RemoveNode(nIt->next());
|
|
}
|
|
}
|
|
// - sub-meshes
|
|
mySubMeshHolder->DeleteAll();
|
|
|
|
myIndexToShape.Clear();
|
|
// - groups on geometry
|
|
std::set<SMESHDS_GroupBase*>::iterator gr = myGroups.begin();
|
|
while ( gr != myGroups.end() ) {
|
|
if ( dynamic_cast<SMESHDS_GroupOnGeom*>( *gr ))
|
|
myGroups.erase( gr++ );
|
|
else
|
|
gr++;
|
|
}
|
|
}
|
|
else {
|
|
myShape = S;
|
|
if ( !S.IsNull() )
|
|
TopExp::MapShapes(myShape, myIndexToShape);
|
|
}
|
|
|
|
SMDS_Mesh::setNbShapes( MaxShapeIndex() );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddHypothesis
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
bool SMESHDS_Mesh::AddHypothesis(const TopoDS_Shape & SS,
|
|
const SMESHDS_Hypothesis * H)
|
|
{
|
|
if (!myShapeToHypothesis.IsBound(SS/*.Oriented(TopAbs_FORWARD)*/)) {
|
|
std::list<const SMESHDS_Hypothesis *> aList;
|
|
myShapeToHypothesis.Bind(SS/*.Oriented(TopAbs_FORWARD)*/, aList);
|
|
}
|
|
std::list<const SMESHDS_Hypothesis *>& alist =
|
|
myShapeToHypothesis(SS/*.Oriented(TopAbs_FORWARD)*/); // ignore orientation of SS
|
|
|
|
//Check if the Hypothesis is still present
|
|
std::list<const SMESHDS_Hypothesis*>::iterator ith = find(alist.begin(),alist.end(), H );
|
|
|
|
if (alist.end() != ith) return false;
|
|
|
|
alist.push_back(H);
|
|
return true;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : RemoveHypothesis
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
bool SMESHDS_Mesh::RemoveHypothesis(const TopoDS_Shape & S,
|
|
const SMESHDS_Hypothesis * H)
|
|
{
|
|
if( myShapeToHypothesis.IsBound( S/*.Oriented(TopAbs_FORWARD)*/ ) )
|
|
{
|
|
std::list<const SMESHDS_Hypothesis *>& alist=myShapeToHypothesis.ChangeFind( S/*.Oriented(TopAbs_FORWARD)*/ );
|
|
std::list<const SMESHDS_Hypothesis*>::iterator ith=find(alist.begin(),alist.end(), H );
|
|
if (ith != alist.end())
|
|
{
|
|
alist.erase(ith);
|
|
return true;
|
|
}
|
|
}
|
|
return false;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddNode
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshNode* SMESHDS_Mesh::AddNode(double x, double y, double z){
|
|
SMDS_MeshNode* node = SMDS_Mesh::AddNode(x, y, z);
|
|
if(node!=NULL) myScript->AddNode(node->GetID(), x, y, z);
|
|
return node;
|
|
}
|
|
|
|
SMDS_MeshNode* SMESHDS_Mesh::AddNodeWithID(double x, double y, double z, smIdType ID){
|
|
SMDS_MeshNode* node = SMDS_Mesh::AddNodeWithID(x,y,z,ID);
|
|
if(node!=NULL) myScript->AddNode(node->GetID(), x, y, z);
|
|
return node;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : MoveNode
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
void SMESHDS_Mesh::MoveNode(const SMDS_MeshNode *n, double x, double y, double z)
|
|
{
|
|
SMDS_Mesh::MoveNode( n, x, y, z );
|
|
myScript->MoveNode(n->GetID(), x, y, z);
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : ChangeElementNodes
|
|
//purpose : Changed nodes of an element provided that nb of nodes does not change
|
|
//=======================================================================
|
|
|
|
bool SMESHDS_Mesh::ChangeElementNodes(const SMDS_MeshElement * elem,
|
|
const SMDS_MeshNode * nodes[],
|
|
const smIdType nbnodes)
|
|
{
|
|
if ( ! SMDS_Mesh::ChangeElementNodes( elem, nodes, nbnodes ))
|
|
return false;
|
|
|
|
std::vector<int> IDs( nbnodes );
|
|
for ( int i = 0; i < nbnodes; i++ )
|
|
IDs [ i ] = nodes[ i ]->GetID();
|
|
myScript->ChangeElementNodes( elem->GetID(), &IDs[0], nbnodes);
|
|
|
|
return true;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : ChangePolygonNodes
|
|
//purpose :
|
|
//=======================================================================
|
|
bool SMESHDS_Mesh::ChangePolygonNodes
|
|
(const SMDS_MeshElement * elem,
|
|
std::vector<const SMDS_MeshNode*> nodes)
|
|
{
|
|
ASSERT(nodes.size() > 3);
|
|
|
|
return ChangeElementNodes(elem, &nodes[0], nodes.size());
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : ChangePolyhedronNodes
|
|
//purpose :
|
|
//=======================================================================
|
|
bool SMESHDS_Mesh
|
|
::ChangePolyhedronNodes (const SMDS_MeshElement * elem,
|
|
const std::vector<const SMDS_MeshNode*>& nodes,
|
|
const std::vector<int> & quantities)
|
|
{
|
|
ASSERT(nodes.size() > 3);
|
|
|
|
if ( !SMDS_Mesh::ChangePolyhedronNodes( elem, nodes, quantities ))
|
|
return false;
|
|
|
|
int i, len = nodes.size();
|
|
std::vector<int> nodes_ids (len);
|
|
for (i = 0; i < len; i++) {
|
|
nodes_ids[i] = nodes[i]->GetID();
|
|
}
|
|
myScript->ChangePolyhedronNodes(elem->GetID(), nodes_ids, quantities);
|
|
|
|
return true;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : Renumber
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
void SMESHDS_Mesh::Renumber (const bool /*isNodes*/, const int /*startID*/, const int /*deltaID*/)
|
|
{
|
|
// TODO not possible yet to have node numbers not starting to O and continuous.
|
|
if ( !this->IsCompacted() )
|
|
this->CompactMesh();
|
|
// SMDS_Mesh::Renumber( isNodes, startID, deltaID );
|
|
// myScript->Renumber( isNodes, startID, deltaID );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : Add0DElement
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_Mesh0DElement* SMESHDS_Mesh::Add0DElementWithID(smIdType nodeID, smIdType ID)
|
|
{
|
|
SMDS_Mesh0DElement* anElem = SMDS_Mesh::Add0DElementWithID(nodeID, ID);
|
|
if (anElem) myScript->Add0DElement(ID, nodeID);
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_Mesh0DElement* SMESHDS_Mesh::Add0DElementWithID
|
|
(const SMDS_MeshNode * node, smIdType ID)
|
|
{
|
|
return Add0DElementWithID(node->GetID(), ID);
|
|
}
|
|
|
|
SMDS_Mesh0DElement* SMESHDS_Mesh::Add0DElement(const SMDS_MeshNode * node)
|
|
{
|
|
SMDS_Mesh0DElement* anElem = SMDS_Mesh::Add0DElement(node);
|
|
if (anElem) myScript->Add0DElement(anElem->GetID(), node->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function :AddBallWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
SMDS_BallElement* SMESHDS_Mesh::AddBallWithID(smIdType node, double diameter, smIdType ID)
|
|
{
|
|
SMDS_BallElement* anElem = SMDS_Mesh::AddBallWithID(node,diameter,ID);
|
|
if (anElem) myScript->AddBall(anElem->GetID(), node, diameter);
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_BallElement* SMESHDS_Mesh::AddBallWithID(const SMDS_MeshNode * node,
|
|
double diameter,
|
|
smIdType ID)
|
|
{
|
|
SMDS_BallElement* anElem = SMDS_Mesh::AddBallWithID(node,diameter,ID);
|
|
if (anElem) myScript->AddBall(anElem->GetID(), node->GetID(), diameter);
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_BallElement* SMESHDS_Mesh::AddBall (const SMDS_MeshNode * node,
|
|
double diameter)
|
|
{
|
|
SMDS_BallElement* anElem = SMDS_Mesh::AddBall(node,diameter);
|
|
if (anElem) myScript->AddBall(anElem->GetID(), node->GetID(), diameter);
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function :AddEdgeWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
SMDS_MeshEdge* SMESHDS_Mesh::AddEdgeWithID(smIdType n1, smIdType n2, smIdType ID)
|
|
{
|
|
SMDS_MeshEdge* anElem = SMDS_Mesh::AddEdgeWithID(n1,n2,ID);
|
|
if(anElem) myScript->AddEdge(ID,n1,n2);
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshEdge* SMESHDS_Mesh::AddEdgeWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
smIdType ID)
|
|
{
|
|
return AddEdgeWithID(n1->GetID(),
|
|
n2->GetID(),
|
|
ID);
|
|
}
|
|
|
|
SMDS_MeshEdge* SMESHDS_Mesh::AddEdge(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2)
|
|
{
|
|
SMDS_MeshEdge* anElem = SMDS_Mesh::AddEdge(n1,n2);
|
|
if(anElem) myScript->AddEdge(anElem->GetID(),
|
|
n1->GetID(),
|
|
n2->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function :AddFace
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFaceWithID(smIdType n1, smIdType n2, smIdType n3, smIdType ID)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddFaceWithID(n1, n2, n3, ID);
|
|
if(anElem) myScript->AddFace(ID,n1,n2,n3);
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFaceWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
smIdType ID)
|
|
{
|
|
return AddFaceWithID(n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID(),
|
|
ID);
|
|
}
|
|
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFace( const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddFace(n1, n2, n3);
|
|
if(anElem) myScript->AddFace(anElem->GetID(),
|
|
n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function :AddFace
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFaceWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4, smIdType ID)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddFaceWithID(n1, n2, n3, n4, ID);
|
|
if(anElem) myScript->AddFace(ID, n1, n2, n3, n4);
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFaceWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
smIdType ID)
|
|
{
|
|
return AddFaceWithID(n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID(),
|
|
n4->GetID(),
|
|
ID);
|
|
}
|
|
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFace(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddFace(n1, n2, n3, n4);
|
|
if(anElem) myScript->AddFace(anElem->GetID(),
|
|
n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID(),
|
|
n4->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function :AddVolume
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4, smIdType ID)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolumeWithID(n1, n2, n3, n4, ID);
|
|
if(anElem) myScript->AddVolume(ID, n1, n2, n3, n4);
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
smIdType ID)
|
|
{
|
|
return AddVolumeWithID(n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID(),
|
|
n4->GetID(),
|
|
ID);
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolume(n1, n2, n3, n4);
|
|
if(anElem) myScript->AddVolume(anElem->GetID(),
|
|
n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID(),
|
|
n4->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function :AddVolume
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4, smIdType n5, smIdType ID)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolumeWithID(n1, n2, n3, n4, n5, ID);
|
|
if(anElem) myScript->AddVolume(ID, n1, n2, n3, n4, n5);
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
smIdType ID)
|
|
{
|
|
return AddVolumeWithID(n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID(),
|
|
n4->GetID(),
|
|
n5->GetID(),
|
|
ID);
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolume(n1, n2, n3, n4, n5);
|
|
if(anElem) myScript->AddVolume(anElem->GetID(),
|
|
n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID(),
|
|
n4->GetID(),
|
|
n5->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function :AddVolume
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4, smIdType n5, smIdType n6, smIdType ID)
|
|
{
|
|
SMDS_MeshVolume *anElem= SMDS_Mesh::AddVolumeWithID(n1, n2, n3, n4, n5, n6, ID);
|
|
if(anElem) myScript->AddVolume(ID, n1, n2, n3, n4, n5, n6);
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6,
|
|
smIdType ID)
|
|
{
|
|
return AddVolumeWithID(n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID(),
|
|
n4->GetID(),
|
|
n5->GetID(),
|
|
n6->GetID(),
|
|
ID);
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolume(n1, n2, n3, n4, n5, n6);
|
|
if(anElem) myScript->AddVolume(anElem->GetID(),
|
|
n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID(),
|
|
n4->GetID(),
|
|
n5->GetID(),
|
|
n6->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function :AddVolume
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4, smIdType n5, smIdType n6, smIdType n7, smIdType n8, smIdType ID)
|
|
{
|
|
SMDS_MeshVolume *anElem= SMDS_Mesh::AddVolumeWithID(n1, n2, n3, n4, n5, n6, n7, n8, ID);
|
|
if(anElem) myScript->AddVolume(ID, n1, n2, n3, n4, n5, n6, n7, n8);
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6,
|
|
const SMDS_MeshNode * n7,
|
|
const SMDS_MeshNode * n8,
|
|
smIdType ID)
|
|
{
|
|
return AddVolumeWithID(n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID(),
|
|
n4->GetID(),
|
|
n5->GetID(),
|
|
n6->GetID(),
|
|
n7->GetID(),
|
|
n8->GetID(),
|
|
ID);
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6,
|
|
const SMDS_MeshNode * n7,
|
|
const SMDS_MeshNode * n8)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolume(n1, n2, n3, n4, n5, n6, n7, n8);
|
|
if(anElem) myScript->AddVolume(anElem->GetID(),
|
|
n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID(),
|
|
n4->GetID(),
|
|
n5->GetID(),
|
|
n6->GetID(),
|
|
n7->GetID(),
|
|
n8->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
|
|
//=======================================================================
|
|
//function :AddVolume
|
|
//purpose : add hexagonal prism
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
|
|
smIdType n5, smIdType n6, smIdType n7, smIdType n8,
|
|
smIdType n9, smIdType n10, smIdType n11, smIdType n12,
|
|
smIdType ID)
|
|
{
|
|
SMDS_MeshVolume *anElem= SMDS_Mesh::AddVolumeWithID(n1, n2, n3, n4, n5, n6, n7, n8, n9, n10, n11, n12, ID);
|
|
if(anElem) myScript->AddVolume(ID, n1, n2, n3, n4, n5, n6, n7, n8, n9, n10, n11, n12);
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6,
|
|
const SMDS_MeshNode * n7,
|
|
const SMDS_MeshNode * n8,
|
|
const SMDS_MeshNode * n9,
|
|
const SMDS_MeshNode * n10,
|
|
const SMDS_MeshNode * n11,
|
|
const SMDS_MeshNode * n12,
|
|
smIdType ID)
|
|
{
|
|
return AddVolumeWithID(n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID(),
|
|
n4->GetID(),
|
|
n5->GetID(),
|
|
n6->GetID(),
|
|
n7->GetID(),
|
|
n8->GetID(),
|
|
n9->GetID(),
|
|
n10->GetID(),
|
|
n11->GetID(),
|
|
n12->GetID(),
|
|
ID);
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6,
|
|
const SMDS_MeshNode * n7,
|
|
const SMDS_MeshNode * n8,
|
|
const SMDS_MeshNode * n9,
|
|
const SMDS_MeshNode * n10,
|
|
const SMDS_MeshNode * n11,
|
|
const SMDS_MeshNode * n12)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolume(n1, n2, n3, n4, n5, n6, n7, n8, n9, n10, n11, n12);
|
|
if(anElem) myScript->AddVolume(anElem->GetID(),
|
|
n1->GetID(),
|
|
n2->GetID(),
|
|
n3->GetID(),
|
|
n4->GetID(),
|
|
n5->GetID(),
|
|
n6->GetID(),
|
|
n7->GetID(),
|
|
n8->GetID(),
|
|
n9->GetID(),
|
|
n10->GetID(),
|
|
n11->GetID(),
|
|
n12->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
|
|
//=======================================================================
|
|
//function : AddPolygonalFace
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddPolygonalFaceWithID (const std::vector<smIdType>& nodes_ids,
|
|
const smIdType ID)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddPolygonalFaceWithID(nodes_ids, ID);
|
|
if (anElem) {
|
|
myScript->AddPolygonalFace(ID, nodes_ids);
|
|
}
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshFace*
|
|
SMESHDS_Mesh::AddPolygonalFaceWithID (const std::vector<const SMDS_MeshNode*>& nodes,
|
|
const smIdType ID)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddPolygonalFaceWithID(nodes, ID);
|
|
if (anElem) {
|
|
smIdType i, len = nodes.size();
|
|
std::vector<smIdType> nodes_ids (len);
|
|
for (i = 0; i < len; i++) {
|
|
nodes_ids[i] = nodes[i]->GetID();
|
|
}
|
|
myScript->AddPolygonalFace(ID, nodes_ids);
|
|
}
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshFace*
|
|
SMESHDS_Mesh::AddPolygonalFace (const std::vector<const SMDS_MeshNode*>& nodes)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddPolygonalFace(nodes);
|
|
if (anElem) {
|
|
smIdType i, len = nodes.size();
|
|
std::vector<smIdType> nodes_ids (len);
|
|
for (i = 0; i < len; i++) {
|
|
nodes_ids[i] = nodes[i]->GetID();
|
|
}
|
|
myScript->AddPolygonalFace(anElem->GetID(), nodes_ids);
|
|
}
|
|
return anElem;
|
|
}
|
|
|
|
|
|
//=======================================================================
|
|
//function : AddQuadPolygonalFace
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddQuadPolygonalFaceWithID (const std::vector<smIdType>& nodes_ids,
|
|
const smIdType ID)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddQuadPolygonalFaceWithID(nodes_ids, ID);
|
|
if (anElem) {
|
|
myScript->AddQuadPolygonalFace(ID, nodes_ids);
|
|
}
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshFace*
|
|
SMESHDS_Mesh::AddQuadPolygonalFaceWithID (const std::vector<const SMDS_MeshNode*>& nodes,
|
|
const smIdType ID)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddQuadPolygonalFaceWithID(nodes, ID);
|
|
if (anElem) {
|
|
smIdType i, len = nodes.size();
|
|
std::vector<smIdType> nodes_ids (len);
|
|
for (i = 0; i < len; i++) {
|
|
nodes_ids[i] = nodes[i]->GetID();
|
|
}
|
|
myScript->AddQuadPolygonalFace(ID, nodes_ids);
|
|
}
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshFace*
|
|
SMESHDS_Mesh::AddQuadPolygonalFace (const std::vector<const SMDS_MeshNode*>& nodes)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddQuadPolygonalFace(nodes);
|
|
if (anElem) {
|
|
smIdType i, len = nodes.size();
|
|
std::vector<smIdType> nodes_ids (len);
|
|
for (i = 0; i < len; i++) {
|
|
nodes_ids[i] = nodes[i]->GetID();
|
|
}
|
|
myScript->AddQuadPolygonalFace(anElem->GetID(), nodes_ids);
|
|
}
|
|
return anElem;
|
|
}
|
|
|
|
|
|
//=======================================================================
|
|
//function : AddPolyhedralVolume
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddPolyhedralVolumeWithID (const std::vector<smIdType>& nodes_ids,
|
|
const std::vector<smIdType>& quantities,
|
|
const smIdType ID)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddPolyhedralVolumeWithID(nodes_ids, quantities, ID);
|
|
if (anElem) {
|
|
myScript->AddPolyhedralVolume(ID, nodes_ids, quantities);
|
|
}
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddPolyhedralVolumeWithID
|
|
(const std::vector<const SMDS_MeshNode*>& nodes,
|
|
const std::vector<smIdType>& quantities,
|
|
const smIdType ID)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddPolyhedralVolumeWithID(nodes, quantities, ID);
|
|
if (anElem) {
|
|
smIdType i, len = nodes.size();
|
|
std::vector<smIdType> nodes_ids (len);
|
|
for (i = 0; i < len; i++) {
|
|
nodes_ids[i] = nodes[i]->GetID();
|
|
}
|
|
myScript->AddPolyhedralVolume(ID, nodes_ids, quantities);
|
|
}
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddPolyhedralVolume
|
|
(const std::vector<const SMDS_MeshNode*>& nodes,
|
|
const std::vector<smIdType>& quantities)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddPolyhedralVolume(nodes, quantities);
|
|
if (anElem) {
|
|
smIdType i, len = nodes.size();
|
|
std::vector<smIdType> nodes_ids (len);
|
|
for (i = 0; i < len; i++) {
|
|
nodes_ids[i] = nodes[i]->GetID();
|
|
}
|
|
myScript->AddPolyhedralVolume(anElem->GetID(), nodes_ids, quantities);
|
|
}
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : removeFromContainers
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
static void removeFromContainers (SMESHDS_Mesh* /*theMesh*/,
|
|
std::set<SMESHDS_GroupBase*>& theGroups,
|
|
std::vector<const SMDS_MeshElement*>& theElems)
|
|
{
|
|
if ( theElems.empty() )
|
|
return;
|
|
|
|
// Rm from group
|
|
// Element can belong to several groups
|
|
if ( !theGroups.empty() )
|
|
{
|
|
std::set<SMESHDS_GroupBase*>::iterator GrIt = theGroups.begin();
|
|
for ( ; GrIt != theGroups.end(); GrIt++ )
|
|
{
|
|
SMESHDS_Group* group = dynamic_cast<SMESHDS_Group*>( *GrIt );
|
|
if ( !group || group->IsEmpty() ) continue;
|
|
|
|
std::vector<const SMDS_MeshElement *>::iterator elIt = theElems.begin();
|
|
for ( ; elIt != theElems.end(); elIt++ )
|
|
{
|
|
group->SMDSGroup().Remove( *elIt );
|
|
if ( group->IsEmpty() ) break;
|
|
}
|
|
}
|
|
}
|
|
|
|
//const bool deleted=true;
|
|
|
|
// Rm from sub-meshes
|
|
// Element should belong to only one sub-mesh
|
|
// if ( theMesh->SubMeshes()->more() )
|
|
// {
|
|
// std::list<const SMDS_MeshElement *>::iterator elIt = theElems.begin();
|
|
// if ( isNode ) {
|
|
// for ( ; elIt != theElems.end(); ++elIt )
|
|
// if ( SMESHDS_SubMesh* sm = theMesh->MeshElements( (*elIt)->getshapeId() ))
|
|
// sm->RemoveNode( static_cast<const SMDS_MeshNode*> (*elIt), deleted );
|
|
// }
|
|
// else {
|
|
// for ( ; elIt != theElems.end(); ++elIt )
|
|
// if ( SMESHDS_SubMesh* sm = theMesh->MeshElements( (*elIt)->getshapeId() ))
|
|
// sm->RemoveElement( *elIt, deleted );
|
|
// }
|
|
// }
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : RemoveNode
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::RemoveNode(const SMDS_MeshNode * n)
|
|
{
|
|
if ( RemoveFreeNode( n, 0, true ))
|
|
return;
|
|
|
|
myScript->RemoveNode(n->GetID());
|
|
|
|
// remove inverse elements from the sub-meshes
|
|
for ( SMDS_ElemIteratorPtr eIt = n->GetInverseElementIterator(); eIt->more() ; )
|
|
{
|
|
const SMDS_MeshElement* e = eIt->next();
|
|
if ( SMESHDS_SubMesh * sm = MeshElements( e->getshapeId() ))
|
|
sm->RemoveElement( e );
|
|
}
|
|
if ( SMESHDS_SubMesh * sm = MeshElements( n->getshapeId() ))
|
|
sm->RemoveNode( n );
|
|
|
|
|
|
std::vector<const SMDS_MeshElement *> removedElems;
|
|
std::vector<const SMDS_MeshElement *> removedNodes;
|
|
|
|
SMDS_Mesh::RemoveElement( n, removedElems, removedNodes, true );
|
|
|
|
removeFromContainers( this, myGroups, removedElems );
|
|
removeFromContainers( this, myGroups, removedNodes );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : RemoveFreeNode
|
|
//purpose :
|
|
//=======================================================================
|
|
bool SMESHDS_Mesh::RemoveFreeNode(const SMDS_MeshNode * n,
|
|
SMESHDS_SubMesh * subMesh,
|
|
bool fromGroups)
|
|
{
|
|
if ( n->NbInverseElements() > 0 )
|
|
return false;
|
|
|
|
myScript->RemoveNode(n->GetID());
|
|
|
|
// Rm from group
|
|
// Node can belong to several groups
|
|
if ( fromGroups && !myGroups.empty() ) {
|
|
std::set<SMESHDS_GroupBase*>::iterator GrIt = myGroups.begin();
|
|
for (; GrIt != myGroups.end(); GrIt++) {
|
|
SMESHDS_Group* group = dynamic_cast<SMESHDS_Group*>(*GrIt);
|
|
if (group && group->GetType() == SMDSAbs_Node )
|
|
group->SMDSGroup().Remove(n);
|
|
}
|
|
}
|
|
|
|
// Rm from sub-mesh
|
|
// Node should belong to only one sub-mesh
|
|
if ( !subMesh || !subMesh->RemoveNode( n ))
|
|
if (( subMesh = MeshElements( n->getshapeId() )))
|
|
subMesh->RemoveNode(n);
|
|
|
|
SMDS_Mesh::RemoveFreeElement(n);
|
|
|
|
return true;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : RemoveElement
|
|
//purpose :
|
|
//========================================================================
|
|
void SMESHDS_Mesh::RemoveElement(const SMDS_MeshElement * elt)
|
|
{
|
|
if (elt->GetType() == SMDSAbs_Node)
|
|
{
|
|
RemoveNode( static_cast<const SMDS_MeshNode*>( elt ));
|
|
return;
|
|
}
|
|
//if (!hasConstructionEdges() && !hasConstructionFaces())
|
|
{
|
|
SMESHDS_SubMesh* subMesh=0;
|
|
if ( elt->getshapeId() > 0 )
|
|
subMesh = MeshElements( elt->getshapeId() );
|
|
|
|
RemoveFreeElement( elt, subMesh, true );
|
|
return;
|
|
}
|
|
|
|
myScript->RemoveElement(elt->GetID());
|
|
|
|
std::vector<const SMDS_MeshElement *> removedElems;
|
|
std::vector<const SMDS_MeshElement *> removedNodes;
|
|
|
|
SMDS_Mesh::RemoveElement(elt, removedElems, removedNodes );
|
|
|
|
removeFromContainers( this, myGroups, removedElems );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : RemoveFreeElement
|
|
//purpose :
|
|
//========================================================================
|
|
void SMESHDS_Mesh::RemoveFreeElement(const SMDS_MeshElement * elt,
|
|
SMESHDS_SubMesh * subMesh,
|
|
bool fromGroups)
|
|
{
|
|
if (elt->GetType() == SMDSAbs_Node) {
|
|
RemoveFreeNode( static_cast<const SMDS_MeshNode*>(elt), subMesh, fromGroups);
|
|
return;
|
|
}
|
|
|
|
// if (hasConstructionEdges() || hasConstructionFaces())
|
|
// // this methods is only for meshes without descendants
|
|
// return;
|
|
|
|
myScript->RemoveElement(elt->GetID());
|
|
|
|
// Rm from group
|
|
// Element can belong to several groups
|
|
if ( fromGroups && !myGroups.empty() ) {
|
|
std::set<SMESHDS_GroupBase*>::iterator GrIt = myGroups.begin();
|
|
for (; GrIt != myGroups.end(); GrIt++) {
|
|
SMESHDS_Group* group = dynamic_cast<SMESHDS_Group*>(*GrIt);
|
|
if (group && !group->IsEmpty())
|
|
group->SMDSGroup().Remove(elt);
|
|
}
|
|
}
|
|
|
|
// Rm from sub-mesh
|
|
// Element should belong to only one sub-mesh
|
|
if ( !subMesh && elt->getshapeId() > 0 )
|
|
subMesh = MeshElements( elt->getshapeId() );
|
|
if ( subMesh )
|
|
subMesh->RemoveElement( elt );
|
|
|
|
SMDS_Mesh::RemoveFreeElement( elt );
|
|
}
|
|
|
|
//================================================================================
|
|
/*!
|
|
* \brief Remove all data from the mesh
|
|
*/
|
|
//================================================================================
|
|
|
|
void SMESHDS_Mesh::ClearMesh()
|
|
{
|
|
myScript->ClearMesh();
|
|
SMDS_Mesh::Clear();
|
|
|
|
// clear submeshes
|
|
SMESHDS_SubMeshIteratorPtr smIt = SubMeshes();
|
|
while ( SMESHDS_SubMesh* sm = const_cast< SMESHDS_SubMesh* >( smIt->next() ))
|
|
sm->Clear();
|
|
|
|
// clear groups
|
|
TGroups::iterator group, groupEnd = myGroups.end();
|
|
for ( group = myGroups.begin(); group != groupEnd; ++group ) {
|
|
if ( SMESHDS_Group* g = dynamic_cast<SMESHDS_Group*>(*group)) {
|
|
SMDSAbs_ElementType groupType = g->GetType();
|
|
g->Clear();
|
|
g->SetType( groupType );
|
|
}
|
|
else
|
|
{
|
|
(*group)->Extent(); // to free cached elements in GroupOnFilter's
|
|
}
|
|
}
|
|
}
|
|
|
|
//================================================================================
|
|
/*!
|
|
* \brief return submesh by shape
|
|
* \param shape - the sub-shape
|
|
* \retval SMESHDS_SubMesh* - the found submesh
|
|
*/
|
|
//================================================================================
|
|
|
|
SMESHDS_SubMesh* SMESHDS_Mesh::getSubmesh( const TopoDS_Shape & shape )
|
|
{
|
|
if ( shape.IsNull() )
|
|
return 0;
|
|
|
|
return NewSubMesh( ShapeToIndex( shape ));
|
|
}
|
|
|
|
//================================================================================
|
|
/*!
|
|
* \brief Add element or node to submesh
|
|
* \param elem - element to add
|
|
* \param subMesh - submesh to be filled in
|
|
*/
|
|
//================================================================================
|
|
|
|
int SMESHDS_Mesh::add(const SMDS_MeshElement* elem, SMESHDS_SubMesh* subMesh )
|
|
{
|
|
if ( elem && subMesh ) {
|
|
subMesh->AddElement( elem );
|
|
return subMesh->GetID();
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : SetNodeOnVolume
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::SetNodeInVolume(const SMDS_MeshNode* aNode,
|
|
const TopoDS_Shell & S)
|
|
{
|
|
if ( int shapeID = add( aNode, getSubmesh( S )))
|
|
const_cast< SMDS_MeshNode* >
|
|
( aNode )->SetPosition( SMDS_SpacePosition::originSpacePosition(), shapeID );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : SetNodeOnVolume
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::SetNodeInVolume(const SMDS_MeshNode * aNode,
|
|
const TopoDS_Solid & S)
|
|
{
|
|
if ( int shapeID = add( aNode, getSubmesh( S )))
|
|
const_cast< SMDS_MeshNode* >
|
|
( aNode )->SetPosition( SMDS_SpacePosition::originSpacePosition(), shapeID );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : SetNodeOnFace
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::SetNodeOnFace(const SMDS_MeshNode * aNode,
|
|
const TopoDS_Face & S,
|
|
double u,
|
|
double v)
|
|
{
|
|
if ( int shapeID = add( aNode, getSubmesh( S )))
|
|
const_cast< SMDS_MeshNode* >
|
|
( aNode )->SetPosition(SMDS_PositionPtr( new SMDS_FacePosition( u, v )), shapeID );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : SetNodeOnEdge
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::SetNodeOnEdge(const SMDS_MeshNode * aNode,
|
|
const TopoDS_Edge & S,
|
|
double u)
|
|
{
|
|
if ( int shapeID = add( aNode, getSubmesh( S )))
|
|
const_cast< SMDS_MeshNode* >
|
|
( aNode )->SetPosition(SMDS_PositionPtr( new SMDS_EdgePosition( u )), shapeID );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : SetNodeOnVertex
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::SetNodeOnVertex(const SMDS_MeshNode * aNode,
|
|
const TopoDS_Vertex & S)
|
|
{
|
|
if ( int shapeID = add( aNode, getSubmesh( S )))
|
|
const_cast< SMDS_MeshNode* >
|
|
( aNode )->SetPosition(SMDS_PositionPtr( new SMDS_VertexPosition()), shapeID );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : UnSetNodeOnShape
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::UnSetNodeOnShape(const SMDS_MeshNode* aNode)
|
|
{
|
|
int shapeId = aNode->getshapeId();
|
|
if (shapeId > 0)
|
|
if ( SMESHDS_SubMesh* sm = MeshElements( shapeId ))
|
|
sm->RemoveNode(aNode);
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : SetMeshElementOnShape
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::SetMeshElementOnShape(const SMDS_MeshElement * anElement,
|
|
const TopoDS_Shape & S)
|
|
{
|
|
add( anElement, getSubmesh(S) );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : UnSetMeshElementOnShape
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::UnSetMeshElementOnShape(const SMDS_MeshElement * elem,
|
|
const TopoDS_Shape & S)
|
|
{
|
|
if ( SMESHDS_SubMesh* sm = MeshElements( S ))
|
|
sm->RemoveElement(elem);
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : ShapeToMesh
|
|
//purpose :
|
|
//=======================================================================
|
|
TopoDS_Shape SMESHDS_Mesh::ShapeToMesh() const
|
|
{
|
|
return myShape;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : IsGroupOfSubShapes
|
|
//purpose : return true if at least one sub-shape of theShape is a sub-shape
|
|
// of myShape or theShape == myShape
|
|
//=======================================================================
|
|
|
|
bool SMESHDS_Mesh::IsGroupOfSubShapes (const TopoDS_Shape& theShape) const
|
|
{
|
|
if ( myIndexToShape.Contains(theShape) )
|
|
return true;
|
|
|
|
for ( TopoDS_Iterator it( theShape ); it.More(); it.Next() )
|
|
if ( IsGroupOfSubShapes( it.Value() ))
|
|
return true;
|
|
|
|
return false;
|
|
}
|
|
|
|
///////////////////////////////////////////////////////////////////////////////
|
|
/// Return the sub mesh linked to the a given TopoDS_Shape or NULL if the given
|
|
/// TopoDS_Shape is unknown
|
|
///////////////////////////////////////////////////////////////////////////////
|
|
SMESHDS_SubMesh * SMESHDS_Mesh::MeshElements(const TopoDS_Shape & S) const
|
|
{
|
|
int Index = ShapeToIndex(S);
|
|
return (SMESHDS_SubMesh *) ( Index ? mySubMeshHolder->Get( Index ) : 0 );
|
|
}
|
|
|
|
///////////////////////////////////////////////////////////////////////////////
|
|
/// Return the sub mesh by Id of shape it is linked to
|
|
///////////////////////////////////////////////////////////////////////////////
|
|
SMESHDS_SubMesh * SMESHDS_Mesh::MeshElements(const int Index) const
|
|
{
|
|
return const_cast< SMESHDS_SubMesh* >( mySubMeshHolder->Get( Index ));
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : SubMeshIndices
|
|
//purpose :
|
|
//=======================================================================
|
|
std::list<int> SMESHDS_Mesh::SubMeshIndices() const
|
|
{
|
|
std::list<int> anIndices;
|
|
SMESHDS_SubMeshIteratorPtr smIt = SubMeshes();
|
|
while ( const SMESHDS_SubMesh* sm = smIt->next() )
|
|
anIndices.push_back( sm->GetID() );
|
|
|
|
return anIndices;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : SubMeshes
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
SMESHDS_SubMeshIteratorPtr SMESHDS_Mesh::SubMeshes() const
|
|
{
|
|
return SMESHDS_SubMeshIteratorPtr( mySubMeshHolder->GetIterator() );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : GetHypothesis
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
const std::list<const SMESHDS_Hypothesis*>&
|
|
SMESHDS_Mesh::GetHypothesis(const TopoDS_Shape & S) const
|
|
{
|
|
if ( myShapeToHypothesis.IsBound( S/*.Oriented(TopAbs_FORWARD)*/ ) ) // ignore orientation of S
|
|
return myShapeToHypothesis.Find( S/*.Oriented(TopAbs_FORWARD)*/ );
|
|
|
|
static std::list<const SMESHDS_Hypothesis*> empty;
|
|
return empty;
|
|
}
|
|
|
|
//================================================================================
|
|
/*!
|
|
* \brief returns true if the hypothesis is assigned to any sub-shape
|
|
*/
|
|
//================================================================================
|
|
|
|
bool SMESHDS_Mesh::IsUsedHypothesis(const SMESHDS_Hypothesis * H) const
|
|
{
|
|
ShapeToHypothesis::Iterator s2h( myShapeToHypothesis );
|
|
for ( ; s2h.More(); s2h.Next() )
|
|
if ( std::find( s2h.Value().begin(), s2h.Value().end(), H ) != s2h.Value().end() )
|
|
return true;
|
|
return false;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : GetScript
|
|
//purpose :
|
|
//=======================================================================
|
|
SMESHDS_Script* SMESHDS_Mesh::GetScript()
|
|
{
|
|
return myScript;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : ClearScript
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::ClearScript()
|
|
{
|
|
myScript->Clear();
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : HasMeshElements
|
|
//purpose :
|
|
//=======================================================================
|
|
bool SMESHDS_Mesh::HasMeshElements(const TopoDS_Shape & S) const
|
|
{
|
|
int Index = myIndexToShape.FindIndex(S);
|
|
return mySubMeshHolder->Get( Index );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : HasHypothesis
|
|
//purpose :
|
|
//=======================================================================
|
|
bool SMESHDS_Mesh::HasHypothesis(const TopoDS_Shape & S)
|
|
{
|
|
return myShapeToHypothesis.IsBound(S/*.Oriented(TopAbs_FORWARD)*/);
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : NewSubMesh
|
|
//purpose :
|
|
//=======================================================================
|
|
SMESHDS_SubMesh * SMESHDS_Mesh::NewSubMesh(int Index)
|
|
{
|
|
SMESHDS_SubMesh* SM = MeshElements( Index );
|
|
if ( !SM )
|
|
{
|
|
SM = new SMESHDS_SubMesh(this, Index);
|
|
mySubMeshHolder->Add( Index, SM );
|
|
}
|
|
return SM;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddCompoundSubmesh
|
|
//purpose :
|
|
//=======================================================================
|
|
|
|
int SMESHDS_Mesh::AddCompoundSubmesh(const TopoDS_Shape& S,
|
|
TopAbs_ShapeEnum type)
|
|
{
|
|
int aMainIndex = 0;
|
|
if ( IsGroupOfSubShapes( S ))
|
|
{
|
|
aMainIndex = myIndexToShape.Add( S );
|
|
bool all = ( type == TopAbs_SHAPE );
|
|
if ( all ) // corresponding simple submesh may exist
|
|
aMainIndex = -aMainIndex;
|
|
SMESHDS_SubMesh * aNewSub = NewSubMesh( aMainIndex );
|
|
if ( !aNewSub->IsComplexSubmesh() ) // is empty
|
|
{
|
|
int shapeType = Max( TopAbs_SOLID, all ? myShape.ShapeType() : type );
|
|
int typeLimit = all ? TopAbs_VERTEX : type;
|
|
for ( ; shapeType <= typeLimit; shapeType++ )
|
|
{
|
|
TopExp_Explorer exp( S, TopAbs_ShapeEnum( shapeType ));
|
|
for ( ; exp.More(); exp.Next() )
|
|
{
|
|
int index = myIndexToShape.FindIndex( exp.Current() );
|
|
if ( index )
|
|
aNewSub->AddSubMesh( NewSubMesh( index ));
|
|
}
|
|
}
|
|
}
|
|
}
|
|
return aMainIndex;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : IndexToShape
|
|
//purpose :
|
|
//=======================================================================
|
|
const TopoDS_Shape& SMESHDS_Mesh::IndexToShape(int ShapeIndex) const
|
|
{
|
|
try
|
|
{
|
|
if ( ShapeIndex > 0 )
|
|
return myIndexToShape.FindKey(ShapeIndex);
|
|
}
|
|
catch ( ... )
|
|
{
|
|
}
|
|
static TopoDS_Shape nullShape;
|
|
return nullShape;
|
|
}
|
|
|
|
//================================================================================
|
|
/*!
|
|
* \brief Return max index of sub-mesh
|
|
*/
|
|
//================================================================================
|
|
|
|
int SMESHDS_Mesh::MaxSubMeshIndex() const
|
|
{
|
|
return mySubMeshHolder->GetMaxID();
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : ShapeToIndex
|
|
//purpose :
|
|
//=======================================================================
|
|
int SMESHDS_Mesh::ShapeToIndex(const TopoDS_Shape & S) const
|
|
{
|
|
int index = myIndexToShape.FindIndex(S);
|
|
return index;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : SetNodeOnVolume
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::SetNodeInVolume(const SMDS_MeshNode* aNode, int Index)
|
|
{
|
|
if ( int shapeID = add( aNode, NewSubMesh( Index )))
|
|
((SMDS_MeshNode*) aNode)->SetPosition( SMDS_SpacePosition::originSpacePosition(), shapeID );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : SetNodeOnFace
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::SetNodeOnFace(const SMDS_MeshNode* aNode, int Index, double u, double v)
|
|
{
|
|
//Set Position on Node
|
|
if ( int shapeID = add( aNode, NewSubMesh( Index )))
|
|
const_cast< SMDS_MeshNode* >
|
|
( aNode )->SetPosition(SMDS_PositionPtr(new SMDS_FacePosition( u, v )), shapeID );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : SetNodeOnEdge
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::SetNodeOnEdge(const SMDS_MeshNode* aNode,
|
|
int Index,
|
|
double u)
|
|
{
|
|
//Set Position on Node
|
|
if ( int shapeID = add( aNode, NewSubMesh( Index )))
|
|
const_cast< SMDS_MeshNode* >
|
|
( aNode )->SetPosition(SMDS_PositionPtr(new SMDS_EdgePosition( u )), shapeID );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : SetNodeOnVertex
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::SetNodeOnVertex(const SMDS_MeshNode* aNode, int Index)
|
|
{
|
|
//Set Position on Node
|
|
if ( int shapeID = add( aNode, NewSubMesh( Index )))
|
|
const_cast< SMDS_MeshNode* >
|
|
( aNode )->SetPosition(SMDS_PositionPtr(new SMDS_VertexPosition()), shapeID );
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : SetMeshElementOnShape
|
|
//purpose :
|
|
//=======================================================================
|
|
void SMESHDS_Mesh::SetMeshElementOnShape(const SMDS_MeshElement* anElement,
|
|
int Index)
|
|
{
|
|
add( anElement, NewSubMesh( Index ));
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : ~SMESHDS_Mesh
|
|
//purpose :
|
|
//=======================================================================
|
|
SMESHDS_Mesh::~SMESHDS_Mesh()
|
|
{
|
|
// myScript
|
|
delete myScript;
|
|
// submeshes
|
|
delete mySubMeshHolder;
|
|
}
|
|
|
|
|
|
//********************************************************************
|
|
//********************************************************************
|
|
//******** *********
|
|
//***** Methods for addition of quadratic elements ******
|
|
//******** *********
|
|
//********************************************************************
|
|
//********************************************************************
|
|
|
|
//=======================================================================
|
|
//function : AddEdgeWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshEdge* SMESHDS_Mesh::AddEdgeWithID(smIdType n1, smIdType n2, smIdType n12, smIdType ID)
|
|
{
|
|
SMDS_MeshEdge* anElem = SMDS_Mesh::AddEdgeWithID(n1,n2,n12,ID);
|
|
if(anElem) myScript->AddEdge(ID,n1,n2,n12);
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddEdge
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshEdge* SMESHDS_Mesh::AddEdge(const SMDS_MeshNode* n1,
|
|
const SMDS_MeshNode* n2,
|
|
const SMDS_MeshNode* n12)
|
|
{
|
|
SMDS_MeshEdge* anElem = SMDS_Mesh::AddEdge(n1,n2,n12);
|
|
if(anElem) myScript->AddEdge(anElem->GetID(),
|
|
n1->GetID(),
|
|
n2->GetID(),
|
|
n12->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddEdgeWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshEdge* SMESHDS_Mesh::AddEdgeWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n12,
|
|
smIdType ID)
|
|
{
|
|
return AddEdgeWithID(n1->GetID(),
|
|
n2->GetID(),
|
|
n12->GetID(),
|
|
ID);
|
|
}
|
|
|
|
|
|
//=======================================================================
|
|
//function : AddFace
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFace(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n31)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddFace(n1,n2,n3,n12,n23,n31);
|
|
if(anElem) myScript->AddFace(anElem->GetID(),
|
|
n1->GetID(), n2->GetID(), n3->GetID(),
|
|
n12->GetID(), n23->GetID(), n31->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddFaceWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFaceWithID(smIdType n1, smIdType n2, smIdType n3,
|
|
smIdType n12,smIdType n23,smIdType n31, smIdType ID)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddFaceWithID(n1,n2,n3,n12,n23,n31,ID);
|
|
if(anElem) myScript->AddFace(ID,n1,n2,n3,n12,n23,n31);
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddFaceWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFaceWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n31,
|
|
smIdType ID)
|
|
{
|
|
return AddFaceWithID(n1->GetID(), n2->GetID(), n3->GetID(),
|
|
n12->GetID(), n23->GetID(), n31->GetID(),
|
|
ID);
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddFace
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFace(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n31,
|
|
const SMDS_MeshNode * nCenter)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddFace(n1,n2,n3,n12,n23,n31,nCenter);
|
|
if(anElem) myScript->AddFace(anElem->GetID(),
|
|
n1->GetID(), n2->GetID(), n3->GetID(),
|
|
n12->GetID(), n23->GetID(), n31->GetID(),
|
|
nCenter->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddFaceWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFaceWithID(smIdType n1, smIdType n2, smIdType n3,
|
|
smIdType n12,smIdType n23,smIdType n31, smIdType nCenter, smIdType ID)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddFaceWithID(n1,n2,n3,n12,n23,n31,nCenter,ID);
|
|
if(anElem) myScript->AddFace(ID,n1,n2,n3,n12,n23,n31,nCenter);
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddFaceWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFaceWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n31,
|
|
const SMDS_MeshNode * nCenter,
|
|
smIdType ID)
|
|
{
|
|
return AddFaceWithID(n1->GetID(), n2->GetID(), n3->GetID(),
|
|
n12->GetID(), n23->GetID(), n31->GetID(),
|
|
nCenter->GetID(), ID);
|
|
}
|
|
|
|
|
|
//=======================================================================
|
|
//function : AddFace
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFace(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n34,
|
|
const SMDS_MeshNode * n41)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddFace(n1,n2,n3,n4,n12,n23,n34,n41);
|
|
if(anElem) myScript->AddFace(anElem->GetID(),
|
|
n1->GetID(), n2->GetID(), n3->GetID(), n4->GetID(),
|
|
n12->GetID(), n23->GetID(), n34->GetID(), n41->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddFaceWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFaceWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
|
|
smIdType n12,smIdType n23,smIdType n34,smIdType n41, smIdType ID)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddFaceWithID(n1,n2,n3,n4,n12,n23,n34,n41,ID);
|
|
if(anElem) myScript->AddFace(ID,n1,n2,n3,n4,n12,n23,n34,n41);
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddFaceWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFaceWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n34,
|
|
const SMDS_MeshNode * n41,
|
|
smIdType ID)
|
|
{
|
|
return AddFaceWithID(n1->GetID(), n2->GetID(), n3->GetID(), n4->GetID(),
|
|
n12->GetID(), n23->GetID(), n34->GetID(), n41->GetID(),
|
|
ID);
|
|
}
|
|
|
|
|
|
//=======================================================================
|
|
//function : AddFace
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFace(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n34,
|
|
const SMDS_MeshNode * n41,
|
|
const SMDS_MeshNode * nCenter)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddFace(n1,n2,n3,n4,n12,n23,n34,n41,nCenter);
|
|
if(anElem) myScript->AddFace(anElem->GetID(),
|
|
n1->GetID(), n2->GetID(), n3->GetID(), n4->GetID(),
|
|
n12->GetID(), n23->GetID(), n34->GetID(), n41->GetID(),
|
|
nCenter->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddFaceWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFaceWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
|
|
smIdType n12,smIdType n23,smIdType n34,smIdType n41,
|
|
smIdType nCenter, smIdType ID)
|
|
{
|
|
SMDS_MeshFace *anElem = SMDS_Mesh::AddFaceWithID(n1,n2,n3,n4,n12,n23,n34,n41,nCenter,ID);
|
|
if(anElem) myScript->AddFace(ID,n1,n2,n3,n4,n12,n23,n34,n41,nCenter);
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddFaceWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshFace* SMESHDS_Mesh::AddFaceWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n34,
|
|
const SMDS_MeshNode * n41,
|
|
const SMDS_MeshNode * nCenter,
|
|
smIdType ID)
|
|
{
|
|
return AddFaceWithID(n1->GetID(), n2->GetID(), n3->GetID(), n4->GetID(),
|
|
n12->GetID(), n23->GetID(), n34->GetID(), n41->GetID(),
|
|
nCenter->GetID(), ID);
|
|
}
|
|
|
|
|
|
//=======================================================================
|
|
//function : AddVolume
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n31,
|
|
const SMDS_MeshNode * n14,
|
|
const SMDS_MeshNode * n24,
|
|
const SMDS_MeshNode * n34)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolume(n1,n2,n3,n4,n12,n23,n31,n14,n24,n34);
|
|
if(anElem) myScript->AddVolume(anElem->GetID(),
|
|
n1->GetID(), n2->GetID(), n3->GetID(), n4->GetID(),
|
|
n12->GetID(), n23->GetID(), n31->GetID(),
|
|
n14->GetID(), n24->GetID(), n34->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddVolumeWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
|
|
smIdType n12,smIdType n23,smIdType n31,
|
|
smIdType n14,smIdType n24,smIdType n34, smIdType ID)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolumeWithID(n1,n2,n3,n4,n12,n23,
|
|
n31,n14,n24,n34,ID);
|
|
if(anElem) myScript->AddVolume(ID,n1,n2,n3,n4,n12,n23,n31,n14,n24,n34);
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddVolumeWithID
|
|
//purpose : 2d order tetrahedron of 10 nodes
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n31,
|
|
const SMDS_MeshNode * n14,
|
|
const SMDS_MeshNode * n24,
|
|
const SMDS_MeshNode * n34,
|
|
smIdType ID)
|
|
{
|
|
return AddVolumeWithID(n1->GetID(), n2->GetID(), n3->GetID(), n4->GetID(),
|
|
n12->GetID(), n23->GetID(), n31->GetID(),
|
|
n14->GetID(), n24->GetID(), n34->GetID(), ID);
|
|
}
|
|
|
|
|
|
//=======================================================================
|
|
//function : AddVolume
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n34,
|
|
const SMDS_MeshNode * n41,
|
|
const SMDS_MeshNode * n15,
|
|
const SMDS_MeshNode * n25,
|
|
const SMDS_MeshNode * n35,
|
|
const SMDS_MeshNode * n45)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolume(n1,n2,n3,n4,n5,n12,n23,n34,n41,
|
|
n15,n25,n35,n45);
|
|
if(anElem)
|
|
myScript->AddVolume(anElem->GetID(), n1->GetID(), n2->GetID(),
|
|
n3->GetID(), n4->GetID(), n5->GetID(),
|
|
n12->GetID(), n23->GetID(), n34->GetID(), n41->GetID(),
|
|
n15->GetID(), n25->GetID(), n35->GetID(), n45->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddVolumeWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4, smIdType n5,
|
|
smIdType n12,smIdType n23,smIdType n34,smIdType n41,
|
|
smIdType n15,smIdType n25,smIdType n35,smIdType n45, smIdType ID)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolumeWithID(n1,n2,n3,n4,n5,
|
|
n12,n23,n34,n41,
|
|
n15,n25,n35,n45,ID);
|
|
if(anElem) myScript->AddVolume(ID,n1,n2,n3,n4,n5,n12,n23,n34,n41,
|
|
n15,n25,n35,n45);
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddVolumeWithID
|
|
//purpose : 2d order pyramid of 13 nodes
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n34,
|
|
const SMDS_MeshNode * n41,
|
|
const SMDS_MeshNode * n15,
|
|
const SMDS_MeshNode * n25,
|
|
const SMDS_MeshNode * n35,
|
|
const SMDS_MeshNode * n45,
|
|
smIdType ID)
|
|
{
|
|
return AddVolumeWithID(n1->GetID(), n2->GetID(), n3->GetID(),
|
|
n4->GetID(), n5->GetID(),
|
|
n12->GetID(), n23->GetID(), n34->GetID(), n41->GetID(),
|
|
n15->GetID(), n25->GetID(), n35->GetID(), n45->GetID(),
|
|
ID);
|
|
}
|
|
|
|
|
|
//=======================================================================
|
|
//function : AddVolume
|
|
//purpose : 2nd order pentahedron (prism) with 15 nodes
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n31,
|
|
const SMDS_MeshNode * n45,
|
|
const SMDS_MeshNode * n56,
|
|
const SMDS_MeshNode * n64,
|
|
const SMDS_MeshNode * n14,
|
|
const SMDS_MeshNode * n25,
|
|
const SMDS_MeshNode * n36)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolume(n1,n2,n3,n4,n5,n6,n12,n23,n31,
|
|
n45,n56,n64,n14,n25,n36);
|
|
if(anElem)
|
|
myScript->AddVolume(anElem->GetID(), n1->GetID(), n2->GetID(),
|
|
n3->GetID(), n4->GetID(), n5->GetID(), n6->GetID(),
|
|
n12->GetID(), n23->GetID(), n31->GetID(),
|
|
n45->GetID(), n56->GetID(), n64->GetID(),
|
|
n14->GetID(), n25->GetID(), n36->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddVolumeWithID
|
|
//purpose : 2nd order pentahedron (prism) with 15 nodes
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3,
|
|
smIdType n4, smIdType n5, smIdType n6,
|
|
smIdType n12,smIdType n23,smIdType n31,
|
|
smIdType n45,smIdType n56,smIdType n64,
|
|
smIdType n14,smIdType n25,smIdType n36, smIdType ID)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolumeWithID(n1,n2,n3,n4,n5,n6,
|
|
n12,n23,n31,
|
|
n45,n56,n64,
|
|
n14,n25,n36,ID);
|
|
if(anElem) myScript->AddVolume(ID,n1,n2,n3,n4,n5,n6,n12,n23,n31,
|
|
n45,n56,n64,n14,n25,n36);
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddVolumeWithID
|
|
//purpose : 2d order Pentahedron (prism) with 15 nodes
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n31,
|
|
const SMDS_MeshNode * n45,
|
|
const SMDS_MeshNode * n56,
|
|
const SMDS_MeshNode * n64,
|
|
const SMDS_MeshNode * n14,
|
|
const SMDS_MeshNode * n25,
|
|
const SMDS_MeshNode * n36,
|
|
smIdType ID)
|
|
{
|
|
return AddVolumeWithID(n1->GetID(), n2->GetID(), n3->GetID(),
|
|
n4->GetID(), n5->GetID(), n6->GetID(),
|
|
n12->GetID(), n23->GetID(), n31->GetID(),
|
|
n45->GetID(), n56->GetID(), n64->GetID(),
|
|
n14->GetID(), n25->GetID(), n36->GetID(),
|
|
ID);
|
|
}
|
|
//=======================================================================
|
|
//function : AddVolume
|
|
//purpose : 2nd order pentahedron (prism) with 18 nodes
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n31,
|
|
const SMDS_MeshNode * n45,
|
|
const SMDS_MeshNode * n56,
|
|
const SMDS_MeshNode * n64,
|
|
const SMDS_MeshNode * n14,
|
|
const SMDS_MeshNode * n25,
|
|
const SMDS_MeshNode * n36,
|
|
const SMDS_MeshNode * n1245,
|
|
const SMDS_MeshNode * n2356,
|
|
const SMDS_MeshNode * n1346)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolume(n1,n2,n3,n4,n5,n6,n12,n23,n31,
|
|
n45,n56,n64,n14,n25,n36,
|
|
n1245, n2356, n1346);
|
|
if(anElem)
|
|
myScript->AddVolume(anElem->GetID(), n1->GetID(), n2->GetID(),
|
|
n3->GetID(), n4->GetID(), n5->GetID(), n6->GetID(),
|
|
n12->GetID(), n23->GetID(), n31->GetID(),
|
|
n45->GetID(), n56->GetID(), n64->GetID(),
|
|
n14->GetID(), n25->GetID(), n36->GetID(),
|
|
n1245->GetID(), n2356->GetID(), n1346->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddVolumeWithID
|
|
//purpose : 2nd order pentahedron (prism) with 18 nodes
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3,
|
|
smIdType n4, smIdType n5, smIdType n6,
|
|
smIdType n12,smIdType n23,smIdType n31,
|
|
smIdType n45,smIdType n56,smIdType n64,
|
|
smIdType n14,smIdType n25,smIdType n36,
|
|
smIdType n1245, smIdType n2356, smIdType n1346,
|
|
smIdType ID)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolumeWithID(n1,n2,n3,n4,n5,n6,
|
|
n12,n23,n31,
|
|
n45,n56,n64,
|
|
n14,n25,n36,
|
|
n1245, n2356, n1346, ID);
|
|
if(anElem) myScript->AddVolume(ID,n1,n2,n3,n4,n5,n6,n12,n23,n31,
|
|
n45,n56,n64,n14,n25,n36, n1245, n2356, n1346);
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddVolumeWithID
|
|
//purpose : 2d order Pentahedron with 18 nodes
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n31,
|
|
const SMDS_MeshNode * n45,
|
|
const SMDS_MeshNode * n56,
|
|
const SMDS_MeshNode * n64,
|
|
const SMDS_MeshNode * n14,
|
|
const SMDS_MeshNode * n25,
|
|
const SMDS_MeshNode * n36,
|
|
const SMDS_MeshNode * n1245,
|
|
const SMDS_MeshNode * n2356,
|
|
const SMDS_MeshNode * n1346,
|
|
smIdType ID)
|
|
{
|
|
return AddVolumeWithID(n1->GetID(), n2->GetID(), n3->GetID(),
|
|
n4->GetID(), n5->GetID(), n6->GetID(),
|
|
n12->GetID(), n23->GetID(), n31->GetID(),
|
|
n45->GetID(), n56->GetID(), n64->GetID(),
|
|
n14->GetID(), n25->GetID(), n36->GetID(),
|
|
n1245->GetID(), n2356->GetID(), n1346->GetID(), ID);
|
|
}
|
|
|
|
|
|
//=======================================================================
|
|
//function : AddVolume
|
|
//purpose : add quadratic hexahedron
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6,
|
|
const SMDS_MeshNode * n7,
|
|
const SMDS_MeshNode * n8,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n34,
|
|
const SMDS_MeshNode * n41,
|
|
const SMDS_MeshNode * n56,
|
|
const SMDS_MeshNode * n67,
|
|
const SMDS_MeshNode * n78,
|
|
const SMDS_MeshNode * n85,
|
|
const SMDS_MeshNode * n15,
|
|
const SMDS_MeshNode * n26,
|
|
const SMDS_MeshNode * n37,
|
|
const SMDS_MeshNode * n48)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolume(n1,n2,n3,n4,n5,n6,n7,n8,
|
|
n12,n23,n34,n41,
|
|
n56,n67,n78,n85,
|
|
n15,n26,n37,n48);
|
|
if(anElem)
|
|
myScript->AddVolume(anElem->GetID(), n1->GetID(), n2->GetID(),
|
|
n3->GetID(), n4->GetID(), n5->GetID(),
|
|
n6->GetID(), n7->GetID(), n8->GetID(),
|
|
n12->GetID(), n23->GetID(), n34->GetID(), n41->GetID(),
|
|
n56->GetID(), n67->GetID(), n78->GetID(), n85->GetID(),
|
|
n15->GetID(), n26->GetID(), n37->GetID(), n48->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddVolumeWithID
|
|
//purpose :
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
|
|
smIdType n5, smIdType n6, smIdType n7, smIdType n8,
|
|
smIdType n12,smIdType n23,smIdType n34,smIdType n41,
|
|
smIdType n56,smIdType n67,smIdType n78,smIdType n85,
|
|
smIdType n15,smIdType n26,smIdType n37,smIdType n48, smIdType ID)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolumeWithID(n1,n2,n3,n4,n5,n6,n7,n8,
|
|
n12,n23,n34,n41,
|
|
n56,n67,n78,n85,
|
|
n15,n26,n37,n48,ID);
|
|
if(anElem) myScript->AddVolume(ID,n1,n2,n3,n4,n5,n6,n7,n8,n12,n23,n34,n41,
|
|
n56,n67,n78,n85,n15,n26,n37,n48);
|
|
return anElem;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddVolumeWithID
|
|
//purpose : 2d order Hexahedrons with 20 nodes
|
|
//=======================================================================
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6,
|
|
const SMDS_MeshNode * n7,
|
|
const SMDS_MeshNode * n8,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n34,
|
|
const SMDS_MeshNode * n41,
|
|
const SMDS_MeshNode * n56,
|
|
const SMDS_MeshNode * n67,
|
|
const SMDS_MeshNode * n78,
|
|
const SMDS_MeshNode * n85,
|
|
const SMDS_MeshNode * n15,
|
|
const SMDS_MeshNode * n26,
|
|
const SMDS_MeshNode * n37,
|
|
const SMDS_MeshNode * n48,
|
|
smIdType ID)
|
|
{
|
|
return AddVolumeWithID(n1->GetID(), n2->GetID(), n3->GetID(), n4->GetID(),
|
|
n5->GetID(), n6->GetID(), n7->GetID(), n8->GetID(),
|
|
n12->GetID(), n23->GetID(), n34->GetID(), n41->GetID(),
|
|
n56->GetID(), n67->GetID(), n78->GetID(), n85->GetID(),
|
|
n15->GetID(), n26->GetID(), n37->GetID(), n48->GetID(),
|
|
ID);
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : AddVolume
|
|
//purpose : add tri-quadratic hexahedron of 27 nodes
|
|
//=======================================================================
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6,
|
|
const SMDS_MeshNode * n7,
|
|
const SMDS_MeshNode * n8,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n34,
|
|
const SMDS_MeshNode * n41,
|
|
const SMDS_MeshNode * n56,
|
|
const SMDS_MeshNode * n67,
|
|
const SMDS_MeshNode * n78,
|
|
const SMDS_MeshNode * n85,
|
|
const SMDS_MeshNode * n15,
|
|
const SMDS_MeshNode * n26,
|
|
const SMDS_MeshNode * n37,
|
|
const SMDS_MeshNode * n48,
|
|
const SMDS_MeshNode * n1234,
|
|
const SMDS_MeshNode * n1256,
|
|
const SMDS_MeshNode * n2367,
|
|
const SMDS_MeshNode * n3478,
|
|
const SMDS_MeshNode * n1458,
|
|
const SMDS_MeshNode * n5678,
|
|
const SMDS_MeshNode * nCenter)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolume(n1,n2,n3,n4,n5,n6,n7,n8,
|
|
n12,n23,n34,n41,
|
|
n56,n67,n78,n85,
|
|
n15,n26,n37,n48,
|
|
n1234,n1256,n2367,n3478,n1458,n5678,nCenter);
|
|
if(anElem)
|
|
myScript->AddVolume(anElem->GetID(), n1->GetID(), n2->GetID(),
|
|
n3->GetID(), n4->GetID(), n5->GetID(),
|
|
n6->GetID(), n7->GetID(), n8->GetID(),
|
|
n12->GetID(), n23->GetID(), n34->GetID(), n41->GetID(),
|
|
n56->GetID(), n67->GetID(), n78->GetID(), n85->GetID(),
|
|
n15->GetID(), n26->GetID(), n37->GetID(), n48->GetID(),
|
|
n1234->GetID(),n1256->GetID(),n2367->GetID(),n3478->GetID(),
|
|
n1458->GetID(),n5678->GetID(),nCenter->GetID());
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
|
|
smIdType n5, smIdType n6, smIdType n7, smIdType n8,
|
|
smIdType n12,smIdType n23,smIdType n34,smIdType n41,
|
|
smIdType n56,smIdType n67,smIdType n78,smIdType n85,
|
|
smIdType n15,smIdType n26,smIdType n37,smIdType n48,
|
|
smIdType n1234,smIdType n1256,smIdType n2367,smIdType n3478,
|
|
smIdType n1458,smIdType n5678,smIdType nCenter,
|
|
smIdType ID)
|
|
{
|
|
SMDS_MeshVolume *anElem = SMDS_Mesh::AddVolumeWithID(n1,n2,n3,n4,n5,n6,n7,n8,
|
|
n12,n23,n34,n41,
|
|
n56,n67,n78,n85,
|
|
n15,n26,n37,n48,
|
|
n1234, n1256, n2367, n3478,
|
|
n1458, n5678, nCenter,
|
|
ID);
|
|
if(anElem) myScript->AddVolume(ID,n1,n2,n3,n4,n5,n6,n7,n8,n12,n23,n34,n41,
|
|
n56,n67,n78,n85,n15,n26,n37,n48,
|
|
n1234, n1256, n2367, n3478,
|
|
n1458, n5678, nCenter);
|
|
return anElem;
|
|
}
|
|
|
|
SMDS_MeshVolume* SMESHDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
|
|
const SMDS_MeshNode * n2,
|
|
const SMDS_MeshNode * n3,
|
|
const SMDS_MeshNode * n4,
|
|
const SMDS_MeshNode * n5,
|
|
const SMDS_MeshNode * n6,
|
|
const SMDS_MeshNode * n7,
|
|
const SMDS_MeshNode * n8,
|
|
const SMDS_MeshNode * n12,
|
|
const SMDS_MeshNode * n23,
|
|
const SMDS_MeshNode * n34,
|
|
const SMDS_MeshNode * n41,
|
|
const SMDS_MeshNode * n56,
|
|
const SMDS_MeshNode * n67,
|
|
const SMDS_MeshNode * n78,
|
|
const SMDS_MeshNode * n85,
|
|
const SMDS_MeshNode * n15,
|
|
const SMDS_MeshNode * n26,
|
|
const SMDS_MeshNode * n37,
|
|
const SMDS_MeshNode * n48,
|
|
const SMDS_MeshNode * n1234,
|
|
const SMDS_MeshNode * n1256,
|
|
const SMDS_MeshNode * n2367,
|
|
const SMDS_MeshNode * n3478,
|
|
const SMDS_MeshNode * n1458,
|
|
const SMDS_MeshNode * n5678,
|
|
const SMDS_MeshNode * nCenter,
|
|
smIdType ID)
|
|
{
|
|
return AddVolumeWithID(n1->GetID(), n2->GetID(), n3->GetID(), n4->GetID(),
|
|
n5->GetID(), n6->GetID(), n7->GetID(), n8->GetID(),
|
|
n12->GetID(), n23->GetID(), n34->GetID(), n41->GetID(),
|
|
n56->GetID(), n67->GetID(), n78->GetID(), n85->GetID(),
|
|
n15->GetID(), n26->GetID(), n37->GetID(), n48->GetID(),
|
|
n1234->GetID(),n1256->GetID(),n2367->GetID(),n3478->GetID(),
|
|
n1458->GetID(),n5678->GetID(),nCenter->GetID(), ID);
|
|
}
|
|
|
|
void SMESHDS_Mesh::CompactMesh()
|
|
{
|
|
if ( IsCompacted() )
|
|
return;
|
|
|
|
SMDS_Mesh::CompactMesh();
|
|
|
|
this->myScript->SetModified(true); // notify GUI client for buildPrs when update
|
|
}
|
|
|
|
void SMESHDS_Mesh::CleanDownWardConnectivity()
|
|
{
|
|
myGrid->CleanDownwardConnectivity();
|
|
}
|
|
|
|
void SMESHDS_Mesh::BuildDownWardConnectivity(bool withEdges)
|
|
{
|
|
myGrid->BuildDownwardConnectivity(withEdges);
|
|
}
|
|
|
|
/*! change some nodes in cell without modifying type or internal connectivity.
|
|
* Nodes inverse connectivity is maintained up to date.
|
|
* @param vtkVolId vtk id of the cell.
|
|
* @param localClonedNodeIds map old node id to new node id.
|
|
* @return ok if success.
|
|
*/
|
|
bool SMESHDS_Mesh::ModifyCellNodes(int vtkVolId, std::map<int,int> localClonedNodeIds)
|
|
{
|
|
myGrid->ModifyCellNodes(vtkVolId, localClonedNodeIds);
|
|
return true;
|
|
}
|