ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-04-08 22:37:28 +05:00
Code Issues Packages Projects Releases Wiki Activity

identify edges

Browse Source
This commit is contained in:
Joachim Schöberl 2016-02-28 14:00:41 +01:00
parent 2e6f7850d7
commit fcac89f793
1 changed files with 23 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
libsrc/meshing/paralleltop.cpp
View File

@ -287,8 +287,8 @@ namespace netgen
} }
TABLE<int> send_edges(cnt_send); TABLE<int> send_edges(cnt_send);
/*
INDEX_2_HASHTABLE<int> gv2e(2*ned); INDEX_2_HASHTABLE<int> gv2e(2*ned);
for (int edge = 1; edge <= ned; edge++) for (int edge = 1; edge <= ned; edge++)
{ {
topology.GetEdgeVertices (edge, v1, v2); topology.GetEdgeVertices (edge, v1, v2);
@ -305,6 +305,28 @@ namespace netgen
send_edges.Add (dest-1, es[1]); send_edges.Add (dest-1, es[1]);
} }
} }
*/
Array<int, PointIndex::BASE> loc2exchange(mesh.GetNV());
for (int dest = 1; dest < ntasks; dest++)
{
loc2exchange = -1;
cnt = 0;
for (PointIndex pi : mesh.Points().Range())
if (IsExchangeVert(dest, pi))
loc2exchange[pi] = cnt++;
for (int edge = 1; edge <= ned; edge++)
{
topology.GetEdgeVertices (edge, v1, v2);
if (IsExchangeVert (dest, v1) && IsExchangeVert (dest, v2))
{
INDEX_2 es(loc2exchange[v1], loc2exchange[v2]);
es.Sort();
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);