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

0-dim elements in topology

Browse Source
This commit is contained in:
Joachim Schöberl 2016-08-17 14:51:02 +02:00
parent 982fee589e
commit f5a9991d24
2 changed files with 25 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
libsrc/meshing/topology.cpp
View File

@ -55,6 +55,7 @@ namespace netgen
delete vert2element;
delete vert2surfelement;
delete vert2segment;
delete vert2pointelement;
}
void MeshTopology :: Update()
@ -90,6 +91,7 @@ namespace netgen
delete vert2element;
delete vert2surfelement;
delete vert2segment;
delete vert2pointelement;
Array<int,PointIndex::BASE> cnt(nv);
Array<int> vnums;
@ -150,6 +152,22 @@ namespace netgen
vert2segment->AddSave (seg[1], si);
}
cnt = 0;
for (int pi = 0; pi < mesh.pointelements.Size(); pi++)
{
const Element0d & pointel = mesh.pointelements[pi];
cnt[pointel.pnum]++;
}
vert2pointelement = new TABLE<int,PointIndex::BASE> (cnt);
for (int pi = 0; pi < mesh.pointelements.Size(); pi++)
{
const Element0d & pointel = mesh.pointelements[pi];
vert2pointelement->AddSave (pointel.pnum, pi);
}
if (buildedges)
{
static int timer1 = NgProfiler::CreateTimer ("topology::buildedges");
@ -734,11 +752,12 @@ namespace netgen
nfa = oldnfa;
INDEX_3_CLOSED_HASHTABLE<int> vert2face(2*max_face_on_vertex+10);
for (int v = PointIndex::BASE; v < nv+PointIndex::BASE; v++)
{
int first_fa = nfa;
INDEX_3_CLOSED_HASHTABLE<int> vert2face(2*max_face_on_vertex+10);
vert2face.DeleteData();
for (int j = 0; j < vert2oldface[v].Size(); j++)
{
@ -770,12 +789,8 @@ namespace netgen
}
for (int j = 0; j < (*vert2element)[v].Size(); j++)
for (ElementIndex elnr : (*vert2element)[v])
{
// NgProfiler::RegionTimer reg3 (timer2d);
ElementIndex elnr = (*vert2element)[v][j];
const Element & el = mesh[elnr];
int nelfaces = GetNFaces (el.GetType());

4
libsrc/meshing/topology.hpp
View File

@ -46,6 +46,7 @@ class MeshTopology
TABLE<ElementIndex,PointIndex::BASE> *vert2element;
TABLE<SurfaceElementIndex,PointIndex::BASE> *vert2surfelement;
TABLE<SegmentIndex,PointIndex::BASE> *vert2segment;
TABLE<int,PointIndex::BASE> *vert2pointelement = nullptr;
int timestamp;
public:
int GetNSurfedges() const {return surfedges.Size();}
@ -143,6 +144,9 @@ public:
FlatArray<SegmentIndex> GetVertexSegments (int vnr) const
{ return (*vert2segment)[vnr]; }
FlatArray<int> GetVertexPointElements (int vnr) const
{ return (*vert2pointelement)[vnr]; }
int GetVerticesEdge ( int v1, int v2) const;
void GetSegmentVolumeElements ( int segnr, Array<ElementIndex> & els ) const;
void GetSegmentSurfaceElements ( int segnr, Array<SurfaceElementIndex> & els ) const;