Commit 2c34353a authored by Douglas Serson's avatar Douglas Serson
Browse files

Start fixing MultiLevelStaticCond with periodic bc

parent 333366c8
......@@ -692,6 +692,7 @@ namespace Nektar
if (i == pIt->second.size())
{
boost::add_vertex(boostGraphObj);
tempGraph[0][meshVertId] = tempGraphVertId++;
m_numNonDirVertexModes++;
}
......@@ -807,6 +808,7 @@ namespace Nektar
if (i == pIt->second.size())
{
boost::add_vertex(boostGraphObj);
tempGraph[1][meshEdgeId] = tempGraphVertId++;
m_numNonDirEdgeModes += EdgeSize[meshEdgeId];
m_numNonDirEdges++;
......@@ -857,6 +859,7 @@ namespace Nektar
meshFaceId = pIt->first;
ASSERTL0(graph[2].count(meshFaceId) == 0,
"This periodic boundary edge has been specified before");
boost::add_vertex(boostGraphObj);
tempGraph[2][meshFaceId] = tempGraphVertId++;
nFaceIntCoeffs = FaceSize[meshFaceId];
m_numNonDirFaceModes+=nFaceIntCoeffs;
......@@ -869,6 +872,7 @@ namespace Nektar
ASSERTL0(graph[2].count(pIt->second[0].id) == 0,
"This periodic boundary face has been specified before");
boost::add_vertex(boostGraphObj);
tempGraph[2][pIt->first] = tempGraphVertId;
tempGraph[2][pIt->second[0].id] = tempGraphVertId++;
nFaceIntCoeffs = FaceSize[pIt->first];
......@@ -1167,6 +1171,46 @@ namespace Nektar
}
}
}
for (pIt = periodicVerts.begin(); pIt != periodicVerts.end(); ++pIt)
{
if (graph[0].count(pIt->first) == 0)
{
for (i = 0; i < pIt->second.size(); ++i)
{
if (!pIt->second[i].isLocal)
{
partVerts.insert(tempGraph[0][pIt->first]);
}
}
}
}
for (pIt = periodicEdges.begin(); pIt != periodicEdges.end(); ++pIt)
{
if (graph[1].count(pIt->first) == 0)
{
for (i = 0; i < pIt->second.size(); ++i)
{
if (!pIt->second[i].isLocal)
{
partVerts.insert(tempGraph[1][pIt->first]);
}
}
}
}
for (pIt = periodicFaces.begin(); pIt != periodicFaces.end(); ++pIt)
{
if (graph[2].count(pIt->first) == 0)
{
for (i = 0; i < pIt->second.size(); ++i)
{
if (!pIt->second[0].isLocal)
{
partVerts.insert(tempGraph[2][pIt->first]);
}
}
}
}
}
int nGraphVerts = tempGraphVertId;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment