Commit baf49f73 authored by Douglas Serson's avatar Douglas Serson
Browse files

Be more strict about using m_fielddef and m_exp at wrong place

parent 6045e6d7
......@@ -263,6 +263,10 @@ void ProcessCreateExp::Process(po::variables_map &vm)
m_f->m_exp[s * nfields + j]->UpdatePhys());
}
}
// Clear fielddef and data
// (they should not be used after running this module)
m_f->m_fielddef = vector<LibUtilities::FieldDefinitionsSharedPtr>();
m_f->m_data = vector<std::vector<NekDouble> >();
}
}
......
......@@ -409,6 +409,9 @@ void ProcessEquiSpacedOutput::Process(po::variables_map &vm)
{
SetHomogeneousConnectivity();
}
// Clear m_exp
m_f->m_exp = vector<MultiRegions::ExpListSharedPtr>();
}
void ProcessEquiSpacedOutput::SetHomogeneousConnectivity(void)
......
......@@ -173,12 +173,7 @@ GlobalMapping::MappingSharedPtr ProcessMapping::GetMapping(FieldSharedPtr f)
// Get field information
int npoints = f->m_exp[0]->GetNpoints();
int expdim = f->m_graph->GetMeshDimension();
int spacedim = expdim;
if ((f->m_fielddef[0]->m_numHomogeneousDir) == 1 ||
(f->m_fielddef[0]->m_numHomogeneousDir) == 2)
{
spacedim = 3;
}
int spacedim = expdim + f->m_numHomogeneousDir;
// Declare coordinates storage
Array<OneD, Array<OneD, NekDouble> > coords_new(3);
......
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