54 const std::shared_ptr<EquationSystem> &pEquation,
const ParamMap &pParams)
55 :
Filter(pSession, pEquation)
58 std::string ext =
".mdl";
62 auto it = pParams.find(
"OutputFrequency");
63 if (it == pParams.end())
75 m_session->MatchSolverInfo(
"CalculatePerturbationEnergy",
"True",
85 it = pParams.find(
"OutputPlane");
86 if (it == pParams.end())
116 if (vComm->GetRank() == 0)
120 m_session->MatchSolverInfo(
"Driver",
"Adaptive", adaptive,
false);
166 int colrank = vComm->GetColumnComm()->GetRank();
167 int nproc = vComm->GetColumnComm()->GetSize();
183 ASSERTL0(
false,
"Stability analysis module not "
184 "implemented for the Compressible Flow "
185 "Solver. Please remove the function BaseFlow "
186 "from your .xml file");
194 string file =
m_session->GetFunctionFilename(
"BaseFlow", 0);
197 for (
int i = 0; i < pFields.size() - 1; ++i)
200 pFields[i]->GetCoeffs(), 1,
201 m_base[i]->GetCoeffs(), 1,
202 pFields[i]->UpdateCoeffs(), 1);
204 energy_tmp = pFields[i]->HomogeneousEnergy();
205 Vmath::Vadd(locsize, energy_tmp, 1, energy, 1, energy, 1);
208 pFields[i]->GetCoeffs(), 1,
209 m_base[i]->GetCoeffs(), 1,
210 pFields[i]->UpdateCoeffs(), 1);
222 for (
int i = 1; i < pFields.size() - 1; ++i)
224 energy_tmp = pFields[i]->HomogeneousEnergy();
225 Vmath::Vadd(locsize, energy_tmp, 1, energy, 1, energy, 1);
232 for (
int i = 0; i < pFields.size() - 1; ++i)
234 energy_tmp = pFields[i]->HomogeneousEnergy();
235 Vmath::Vadd(locsize, energy_tmp, 1, energy, 1, energy, 1);
245 for (j = 0; j < energy.size(); ++j, ++m)
248 << energy[j] << endl;
251 for (
int i = 1; i < nproc; ++i)
253 vComm->GetColumnComm()->Recv(i, energy);
255 for (j = 0; j < energy.size(); ++j, ++m)
258 << setw(18) << energy[j] << endl;
264 vComm->GetColumnComm()->Send(0, energy);
270 ASSERTL0(
false,
"3D Homogeneous 2D energy "
271 "dumping not implemented yet");
282 for (
int i = 1; i < pFields.size() - 1; ++i)
284 pFields[i]->SetPhysState(
true);
286 energy += norm * norm;
299 for (
int i = 0; i < pFields.size() - 1; ++i)
301 pFields[i]->SetPhysState(
true);
303 energy += norm * norm;
320 if (pFields[0]->GetComm()->GetRank() == 0)
339 if (pFields[
field]->GetPhysState() ==
false)
341 pFields[
field]->BwdTrans(pFields[
field]->GetCoeffs(),
342 pFields[
field]->UpdatePhys());
346 L2error = pFields[
field]->L2(pFields[
field]->GetPhys());
357 int m_expdim = graphShrPtr->GetMeshDimension();
360 if (
m_session->DefinesSolverInfo(
"PROJECTION"))
362 std::string ProjectStr =
m_session->GetSolverInfo(
"PROJECTION");
364 if ((ProjectStr ==
"Continuous") || (ProjectStr ==
"Galerkin") ||
365 (ProjectStr ==
"CONTINUOUS") || (ProjectStr ==
"GALERKIN"))
369 else if ((ProjectStr ==
"MixedCGDG") ||
370 (ProjectStr ==
"Mixed_CG_Discontinuous"))
374 else if (ProjectStr ==
"DisContinuous")
380 ASSERTL0(
false,
"PROJECTION value not recognised");
385 cerr <<
"Projection type not specified in SOLVERINFO,"
386 "defaulting to continuous Galerkin"
391 if (
m_session->DefinesSolverInfo(
"ModeType"))
405 if (
m_session->DefinesSolverInfo(
"ModeType"))
421 ASSERTL0(
false,
"SolverInfo ModeType not valid");
436 for (i = 0; i <
m_base.size(); i++)
455 for (i = 0; i <
m_base.size(); i++)
464 m_base[i]->SetWaveSpace(
true);
476 for (i = 0; i <
m_base.size(); i++)
485 m_base[i]->SetWaveSpace(
true);
495 for (i = 0; i <
m_base.size(); i++)
504 m_base[i]->SetWaveSpace(
false);
518 for (i = 1; i <
m_base.size(); i++)
533 for (i = 1; i <
m_base.size(); i++)
543 "Expansion dimension not recognised");
555 for (i = 0; i <
m_base.size(); i++)
565 for (i = 0; i <
m_base.size(); i++)
578 "Expansion dimension not recognised");
589 std::vector<LibUtilities::FieldDefinitionsSharedPtr> FieldDef;
590 std::vector<std::vector<NekDouble>> FieldData;
593 m_fld->Import(pInfile, FieldDef, FieldData);
595 int nvar =
m_session->GetVariables().size();
596 if (
m_session->DefinesSolverInfo(
"HOMOGENEOUS"))
598 std::string HomoStr =
m_session->GetSolverInfo(
"HOMOGENEOUS");
601 for (
int j = 0; j < nvar; ++j)
603 for (
int i = 0; i < FieldDef.size(); ++i)
605 bool flag = FieldDef[i]->m_fields[j] ==
m_session->GetVariable(j);
607 ASSERTL0(flag, (std::string(
"Order of ") + pInfile +
608 std::string(
" data and that defined in "
609 "m_boundaryconditions differs"))
612 m_base[j]->ExtractDataToCoeffs(FieldDef[i], FieldData[i],
613 FieldDef[i]->m_fields[j],
614 m_base[j]->UpdateCoeffs());
#define ASSERTL0(condition, msg)
#define NEKERROR(type, msg)
Assert Level 0 – Fundamental assert which is used whether in FULLDEBUG, DEBUG or OPT compilation mode...
Describes the specification for a Basis.
NekDouble Evaluate() const
static std::shared_ptr< FieldIO > CreateDefault(const LibUtilities::SessionReaderSharedPtr session)
Returns an object for the default FieldIO method.
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
Defines a specification for a set of points.
General purpose memory allocation routines with the ability to allocate from thread specific memory p...
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.
SOLVER_UTILS_EXPORT std::string SetupOutput(const std::string ext, const ParamMap &pParams)
LibUtilities::SessionReaderSharedPtr m_session
std::map< std::string, std::string > ParamMap
bool m_updateOnInitialise
static FilterSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession, const std::shared_ptr< EquationSystem > &pEquation, const ParamMap &pParams)
void v_Initialise(const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time) override
void v_Finalise(const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time) override
void v_Update(const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time) override
static std::string className
enum MultiRegions::ProjectionType m_projectionType
unsigned int m_outputFrequency
bool v_IsTimeDependent() override
SOLVER_UTILS_EXPORT FilterModalEnergy(const LibUtilities::SessionReaderSharedPtr &pSession, const std::shared_ptr< EquationSystem > &pEquation, const ParamMap &pParams)
bool m_homogen_dealiasing
SOLVER_UTILS_EXPORT ~FilterModalEnergy() override
LibUtilities::FieldIOSharedPtr m_fld
unsigned int m_outputPlane
std::ofstream m_outputStream
void ImportFldBase(std::string pInfile)
NekDouble L2Error(const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, unsigned int field, const NekDouble &time)
Array< OneD, MultiRegions::ExpListSharedPtr > m_base
void SetUpBaseFields(SpatialDomains::MeshGraphSharedPtr &mesh)
static MeshGraphSharedPtr Read(const LibUtilities::SessionReaderSharedPtr pSession, LibUtilities::DomainRangeShPtr rng=LibUtilities::NullDomainRangeShPtr, bool fillGraph=true, SpatialDomains::MeshGraphSharedPtr partitionedGraph=nullptr)
std::shared_ptr< SessionReader > SessionReaderSharedPtr
@ eFourierEvenlySpaced
1D Evenly-spaced points using Fourier Fit
@ eFourierSingleModeSpaced
1D Non Evenly-spaced points for Single Mode analysis
std::shared_ptr< Comm > CommSharedPtr
Pointer to a Communicator object.
@ eFourierHalfModeRe
Fourier Modified expansions with just the real part of the first mode .
@ eFourier
Fourier Expansion .
@ eMixed_CG_Discontinuous
std::shared_ptr< ContField > ContFieldSharedPtr
FilterFactory & GetFilterFactory()
std::shared_ptr< MeshGraph > MeshGraphSharedPtr
void Vadd(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Add vector z = x+y.
void Vsub(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Subtract vector z = x-y.