48 const std::weak_ptr<EquationSystem> &pEquation)
49 : m_session(pSession), m_equ(pEquation)
59 if (&inarray != &outarray)
61 int nvar = std::min(inarray.size(), outarray.size());
62 for (
int i = 0; i < nvar; ++i)
64 if (&inarray[i] != &outarray[i])
66 int np = std::min(inarray[i].size(), outarray[i].size());
78 const std::weak_ptr<EquationSystem> &pEquation,
80 const unsigned int &pNumForcingFields)
82 std::vector<ForcingSharedPtr> vForceList;
84 if (!pSession->DefinesElement(
"Nektar/Forcing"))
89 TiXmlElement *vForcing = pSession->GetElement(
"Nektar/Forcing");
92 unsigned int vNumForcingFields = pNumForcingFields;
93 if (!pNumForcingFields)
95 vNumForcingFields = pFields.size();
98 TiXmlElement *vForce = vForcing->FirstChildElement(
"FORCE");
101 std::string vType = vForce->Attribute(
"TYPE");
103 TiXmlElement *vForceParam = vForce;
105 vForceParam, pSession->GetTimeLevel(),
false);
107 vType, pSession, pEquation, pFields, vNumForcingFields,
110 vForce = vForce->NextSiblingElement(
"FORCE");
131 ASSERTL0(pSession->DefinesFunction(pFunctionName),
132 "Function '" + pFunctionName +
"' does not exist.");
135 pSession->GetFunction(pFunctionName, pFieldName);
147 pEqn->Evaluate(x0, x0, x0, pTime, pArray);
163 pSession, pFields[0], pName, pCache);
181 ASSERTL0(
false,
"v_ApplyCoeff not defined");
#define ASSERTL0(condition, msg)
Provides a generic Factory class.
static void GetXMLElementTimeLevel(TiXmlElement *&element, const size_t timeLevel, const bool enableCheck=true)
Get XML elment time level (Parallel-in-Time)
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.
virtual SOLVER_UTILS_EXPORT void v_PreApply(const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble &time)
Array< OneD, Array< OneD, NekDouble > > m_Forcing
Evaluated forcing function.
static SOLVER_UTILS_EXPORT std::vector< ForcingSharedPtr > Load(const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const unsigned int &pNumForcingFields=0)
SOLVER_UTILS_EXPORT Array< OneD, Array< OneD, NekDouble > > & UpdateForces()
SOLVER_UTILS_EXPORT const Array< OneD, const Array< OneD, NekDouble > > & GetForces()
std::map< std::string, SolverUtils::SessionFunctionSharedPtr > m_sessionFunctions
Map of known SessionFunctions.
SOLVER_UTILS_EXPORT void EvaluateTimeFunction(LibUtilities::SessionReaderSharedPtr pSession, std::string pFieldName, Array< OneD, NekDouble > &pArray, const std::string &pFunctionName, NekDouble pTime=NekDouble(0))
SOLVER_UTILS_EXPORT SessionFunctionSharedPtr GetFunction(const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const LibUtilities::SessionReaderSharedPtr &pSession, std::string pName, bool pCache=false)
Get a SessionFunction by name.
SOLVER_UTILS_EXPORT Forcing(const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation)
Constructor.
virtual SOLVER_UTILS_EXPORT void v_ApplyCoeff(const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble &time)
std::shared_ptr< SessionReader > SessionReaderSharedPtr
std::shared_ptr< Equation > EquationSharedPtr
ForcingFactory & GetForcingFactory()
Declaration of the forcing factory singleton.
std::shared_ptr< SessionFunction > SessionFunctionSharedPtr
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)