Commit c799ae0d authored by Kilian Lackhove's avatar Kilian Lackhove
Browse files

Forcing: added Smooth

parent d92e034c
......@@ -118,6 +118,19 @@ namespace Nektar
return vForceList;
}
void Forcing::Smooth(const MultiRegions::ExpListSharedPtr &field)
{
Array<OneD, NekDouble> tmpC(m_Forcing[0].num_elements());
for (int i = 0; i < m_NumVariable; ++i)
{
field->IProductWRTBase(m_Forcing[i], tmpC);
field->MultiplyByElmtInvMass(tmpC, tmpC);
field->LocalToGlobal(tmpC, tmpC);
field->GlobalToLocal(tmpC, tmpC);
field->BwdTrans(tmpC, m_Forcing[i]);
}
}
void Forcing::EvaluateTimeFunction(
LibUtilities::SessionReaderSharedPtr pSession,
std::string pFieldName,
......
......@@ -90,6 +90,10 @@ namespace SolverUtils
const Array<OneD, MultiRegions::ExpListSharedPtr>& pFields,
const unsigned int& pNumForcingFields = 0);
SOLVER_UTILS_EXPORT void Smooth(
const MultiRegions::ExpListSharedPtr &field);
protected:
/// Session reader
LibUtilities::SessionReaderSharedPtr m_session;
......
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