ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-01-13 06:30:34 +05:00
Code Issues Packages Projects Releases Wiki Activity

update parallel topology

Browse Source
This commit is contained in:
Joachim Schöberl 2016-02-28 05:50:15 +01:00
parent 9ef2028729
commit bcbc4294bc
1 changed files with 20 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
libsrc/meshing/paralleltop.cpp
View File

@ -231,6 +231,26 @@ namespace netgen
if (id == 0) return; if (id == 0) return;
const MeshTopology & topology = mesh.GetTopology();
cout << "update refined vertices" << endl;
// update new vertices after mesh-refinement
if (loc2distvert.Size() < mesh.mlbetweennodes)
{
int oldnv = loc2distvert.Size();
int newnv = mesh.mlbetweennodes.Size();
loc2distvert.ChangeSize(mesh.mlbetweennodes.Size());
for (PointIndex pi = oldnv+PointIndex::BASE; nr < newnv+PointIndex::Base; nr++)
{
PointIndex v1 = mesh.mlbetweennodex[pi][0];
PointIndex v2 = mesh.mlbetweennodex[pi][1];
for (int dest = 1; dest < ntasks; dest++)
if (IsExchangeVert (dest, v1) && IsExchangeVert (dest, v2))
SetDistantPNum(dest, pi);
}
}
cout << "update refined vertices done" << endl;
Array<int> sendarray, recvarray; Array<int> sendarray, recvarray;
@ -246,7 +266,6 @@ namespace netgen
NgProfiler::StartTimer (timere); NgProfiler::StartTimer (timere);
const MeshTopology & topology = mesh.GetTopology();
int nfa = topology . GetNFaces(); int nfa = topology . GetNFaces();
int ned = topology . GetNEdges(); int ned = topology . GetNEdges();