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

parallel refinement

Browse Source
This commit is contained in:
Joachim Schöberl 2016-02-28 08:11:34 +01:00
parent 3e82f61424
commit fbb7071939
1 changed files with 16 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
libsrc/meshing/paralleltop.cpp
View File

@ -206,7 +206,7 @@ namespace netgen
void ParallelMeshTopology :: UpdateCoarseGrid () void ParallelMeshTopology :: UpdateCoarseGrid ()
{ {
cout << "updatecoarsegrid called, is_updated = " << is_updated << endl; cout << "updatecoarsegrid called, is_updated = " << is_updated << endl;
if (is_updated) return; // if (is_updated) return;
Reset(); Reset();
static int timer = NgProfiler::CreateTimer ("UpdateCoarseGrid"); static int timer = NgProfiler::CreateTimer ("UpdateCoarseGrid");
@ -239,15 +239,16 @@ namespace netgen
cout << "update refined vertices" << endl; cout << "update refined vertices" << endl;
// update new vertices after mesh-refinement // update new vertices after mesh-refinement
if (loc2distvert.Size() < mesh.mlbetweennodes.Size()) if (mesh.mlbetweennodes.Size() > 0)
{ {
int oldnv = loc2distvert.Size(); // int oldnv = loc2distvert.Size();
int newnv = mesh.mlbetweennodes.Size(); int newnv = mesh.mlbetweennodes.Size();
loc2distvert.ChangeSize(mesh.mlbetweennodes.Size()); loc2distvert.ChangeSize(mesh.mlbetweennodes.Size());
for (PointIndex pi = oldnv+PointIndex::BASE; pi < newnv+PointIndex::BASE; pi++) for (PointIndex pi = PointIndex::BASE; pi < newnv+PointIndex::BASE; pi++)
{ {
PointIndex v1 = mesh.mlbetweennodes[pi][0]; PointIndex v1 = mesh.mlbetweennodes[pi][0];
PointIndex v2 = mesh.mlbetweennodes[pi][1]; PointIndex v2 = mesh.mlbetweennodes[pi][1];
if (mesh.mlbetweennodes[pi][0] != PointIndex::BASE-1)
for (int dest = 1; dest < ntasks; dest++) for (int dest = 1; dest < ntasks; dest++)
if (IsExchangeVert (dest, v1) && IsExchangeVert (dest, v2)) if (IsExchangeVert (dest, v1) && IsExchangeVert (dest, v2))
SetDistantPNum(dest, pi); SetDistantPNum(dest, pi);
@ -291,19 +292,21 @@ namespace netgen
for (int edge = 1; edge <= ned; edge++) for (int edge = 1; edge <= ned; edge++)
{ {
topology.GetEdgeVertices (edge, v1, v2); topology.GetEdgeVertices (edge, v1, v2);
INDEX_2 es(GetGlobalPNum(v1), GetGlobalPNum(v2)); // INDEX_2 es(GetGlobalPNum(v1), GetGlobalPNum(v2));
es.Sort(); // es.Sort();
gv2e.Set (es, edge); // gv2e.Set (es, edge);
for (int dest = 1; dest < ntasks; dest++) for (int dest = 1; dest < ntasks; dest++)
if (IsExchangeVert (dest, v1) && IsExchangeVert (dest, v2)) if (IsExchangeVert (dest, v1) && IsExchangeVert (dest, v2))
{ {
send_edges.Add (dest-1, es[0]); SetDistantEdgeNum(dest, edge);
send_edges.Add (dest-1, es[1]); // send_edges.Add (dest-1, es[0]);
// send_edges.Add (dest-1, es[1]);
} }
} }
/*
TABLE<int> recv_edges(ntasks-1); TABLE<int> recv_edges(ntasks-1);
MyMPI_ExchangeTable (send_edges, recv_edges, MPI_TAG_MESH+9, MPI_LocalComm); MyMPI_ExchangeTable (send_edges, recv_edges, MPI_TAG_MESH+9, MPI_LocalComm);
@ -318,7 +321,7 @@ namespace netgen
SetDistantEdgeNum (sender, gv2e.Get(gv12)); SetDistantEdgeNum (sender, gv2e.Get(gv12));
} }
} }
*/
NgProfiler::StopTimer (timere); NgProfiler::StopTimer (timere);