ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-12-26 13:50:33 +05:00
Code Issues Packages Projects Releases Wiki Activity

getvertexsurfaceelement for 2D/3D

Browse Source
This commit is contained in:
Joachim Schoeberl 2014-06-25 02:15:45 +00:00
parent b06e66cece
commit e0293b1bda
1 changed files with 41 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

50
libsrc/interface/nginterface.cpp
View File

@ -1932,14 +1932,29 @@ int Ng_GetVertex_Elements( int vnr, int* elems )
///// Added by Roman Stainko .... ///// Added by Roman Stainko ....
int Ng_GetVertex_SurfaceElements( int vnr, int* elems ) int Ng_GetVertex_SurfaceElements( int vnr, int* elems )
{ {
const MeshTopology& topology = mesh->GetTopology(); switch (mesh->GetDimension())
ArrayMem<int,4> indexArray; {
topology.GetVertexSurfaceElements( vnr, indexArray ); case 3:
{
const MeshTopology& topology = mesh->GetTopology();
ArrayMem<int,4> indexArray;
topology.GetVertexSurfaceElements( vnr, indexArray );
for( int i=0; i<indexArray.Size(); i++ ) for( int i=0; i<indexArray.Size(); i++ )
elems[i] = indexArray[i]; elems[i] = indexArray[i];
return indexArray.Size(); return indexArray.Size();
}
case 2:
{
int cnt = 0;
for (SegmentIndex i = 0; i < mesh->GetNSeg(); i++)
if ( ((*mesh)[i][0] == vnr) || ((*mesh)[i][1] == vnr) )
elems[cnt++] = i+1;
return cnt;
}
}
return 0;
} }
///// Added by Roman Stainko .... ///// Added by Roman Stainko ....
@ -1955,11 +1970,24 @@ int Ng_GetVertex_NElements( int vnr )
///// Added by Roman Stainko .... ///// Added by Roman Stainko ....
int Ng_GetVertex_NSurfaceElements( int vnr ) int Ng_GetVertex_NSurfaceElements( int vnr )
{ {
const MeshTopology& topology = mesh->GetTopology(); switch (mesh->GetDimension())
ArrayMem<int,4> indexArray; {
topology.GetVertexSurfaceElements( vnr, indexArray ); case 3:
{
return indexArray.Size(); const MeshTopology& topology = mesh->GetTopology();
ArrayMem<int,4> indexArray;
topology.GetVertexSurfaceElements( vnr, indexArray );
return indexArray.Size();
}
case 2:
{
int cnt = 0;
for (SegmentIndex i = 0; i < mesh->GetNSeg(); i++)
if ( ((*mesh)[i][0] == vnr) || ((*mesh)[i][1] == vnr) ) cnt++;
return cnt;
}
}
return 0;
} }