Commit 6c69665e authored by Spencer Sherwin's avatar Spencer Sherwin

Modified how the m_expansion is set up to be consistent with the previous...

Modified how the m_expansion is set up to be consistent with the previous implementation before MeshPartition was moved. The issue related to looping over all composites rather than just the ones deinfed in the EXPANSION section
parent 49db729e
......@@ -369,15 +369,16 @@ void MeshPartition::ReadExpansions()
// construct mapping (elmt id, field name) -> nummodes
map<int, CompositeSharedPtr> &compMap =
m_meshgraph->GetComposites();
for (auto &i : compMap)
{
for (int i = 0; i < composite.size(); ++i)
{
for (int j = 0; j < fieldName.size(); j++)
{
for (int k = 0; k < i.second->m_geomVec.size(); k++)
for (int k = 0; k < compMap[i]->m_geomVec.size(); k++)
{
int elid = i.second->m_geomVec[k]->GetGlobalID();
int elid = compMap[i]->m_geomVec[k]->GetGlobalID();
m_expansions[elid][fieldName[j]] = nummodes;
m_shape[elid] = i.second->m_geomVec[k]->GetShapeType();
m_shape[elid] = compMap[i]->m_geomVec[k]->GetShapeType();
}
}
}
......
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