65 const TiXmlElement *pMapping)
72 "Mapping X = x + f(z) needs 3 velocity components.");
79 int physTot =
m_fields[0]->GetTotPoints();
84 Vmath::Vadd(physTot, wk, 1, inarray[0], 1, outarray[0], 1);
97 int physTot =
m_fields[0]->GetTotPoints();
108 Vmath::Vsub(physTot, inarray[2], 1, wk, 1, outarray[2], 1);
115 int physTot =
m_fields[0]->GetTotPoints();
120 Vmath::Vsub(physTot, inarray[0], 1, wk, 1, outarray[0], 1);
133 int physTot =
m_fields[0]->GetTotPoints();
143 Vmath::Vadd(physTot, inarray[2], 1, outarray[2], 1, outarray[2], 1);
148 int physTot =
m_fields[0]->GetTotPoints();
156 int physTot =
m_fields[0]->GetTotPoints();
164 int physTot =
m_fields[0]->GetTotPoints();
167 for (
int i = 0; i < nvel * nvel; i++)
172 for (
int i = 0; i < nvel; i++)
174 Vmath::Sadd(physTot, 1.0, outarray[i + nvel * i], 1,
175 outarray[i + nvel * i], 1);
184 outarray[2 * nvel + 2], 1);
190 int physTot =
m_fields[0]->GetTotPoints();
194 for (
int i = 0; i < nvel * nvel; i++)
199 for (
int i = 0; i < nvel; i++)
201 Vmath::Sadd(physTot, 1.0, outarray[i + nvel * i], 1,
202 outarray[i + nvel * i], 1);
208 Vmath::Vcopy(physTot, wk, 1, outarray[0 * nvel + 2], 1);
209 Vmath::Vcopy(physTot, wk, 1, outarray[2 * nvel + 0], 1);
213 outarray[0 * nvel + 0], 1);
220 int physTot =
m_fields[0]->GetTotPoints();
229 Vmath::Vadd(physTot, outarray[0], 1, inarray[0], 1, outarray[0], 1);
237 Vmath::Vadd(physTot, wk, 1, outarray[2], 1, outarray[2], 1);
244 int physTot =
m_fields[0]->GetTotPoints();
252 Vmath::Vsub(physTot, inarray[2], 1, outarray[2], 1, outarray[2], 1);
260 Vmath::Vsub(physTot, wk, 1, outarray[0], 1, outarray[0], 1);
267 int physTot =
m_fields[0]->GetTotPoints();
270 for (
int i = 0; i < nvel; i++)
272 for (
int j = 0; j < nvel; j++)
282 outarray[0 * nvel + 2], 1);
289 int physTot =
m_fields[0]->GetTotPoints();
292 for (
int i = 0; i < nvel; i++)
294 for (
int j = 0; j < nvel; j++)
304 outarray[2 * nvel + 2], 1);
309 int phystot =
m_fields[0]->GetTotPoints();
317 bool waveSpace =
m_fields[0]->GetWaveSpace();
328 m_fields[0]->SetWaveSpace(waveSpace);
#define ASSERTL0(condition, msg)
Base class for mapping to be applied to the coordinate system.
int m_nConvectiveFields
Number of velocity components.
Array< OneD, Array< OneD, NekDouble > > m_GeometricInfo
Array with metric terms of the mapping.
Array< OneD, Array< OneD, NekDouble > > m_coords
Array with the Cartesian coordinates.
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
virtual GLOBAL_MAPPING_EXPORT void v_InitObject(const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const TiXmlElement *pMapping)
bool m_constantJacobian
Flag defining if the Jacobian is constant.
GLOBAL_MAPPING_EXPORT void v_GetMetricTensor(Array< OneD, Array< OneD, NekDouble > > &outarray) override
MappingXofZ(const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields)
GLOBAL_MAPPING_EXPORT void v_RaiseIndex(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) override
GLOBAL_MAPPING_EXPORT void v_CovarToCartesian(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) override
static GLOBAL_MAPPING_EXPORT MappingSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const TiXmlElement *pMapping)
Creates an instance of this class.
static std::string className
Name of the class.
GLOBAL_MAPPING_EXPORT void v_UpdateGeomInfo() override
GLOBAL_MAPPING_EXPORT void v_ApplyChristoffelCovar(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) override
GLOBAL_MAPPING_EXPORT void v_LowerIndex(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) override
GLOBAL_MAPPING_EXPORT void v_ApplyChristoffelContravar(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) override
GLOBAL_MAPPING_EXPORT void v_ContravarFromCartesian(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) override
GLOBAL_MAPPING_EXPORT void v_GetInvMetricTensor(Array< OneD, Array< OneD, NekDouble > > &outarray) override
GLOBAL_MAPPING_EXPORT void v_InitObject(const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const TiXmlElement *pMapping) override
GLOBAL_MAPPING_EXPORT void v_DotGradJacobian(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, NekDouble > &outarray) override
GLOBAL_MAPPING_EXPORT void v_GetJacobian(Array< OneD, NekDouble > &outarray) override
GLOBAL_MAPPING_EXPORT void v_ContravarToCartesian(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) override
GLOBAL_MAPPING_EXPORT void v_CovarFromCartesian(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) override
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
MappingFactory & GetMappingFactory()
Declaration of the mapping factory singleton.
std::shared_ptr< SessionReader > SessionReaderSharedPtr
MultiRegions::Direction const DirCartesianMap[]
void Vmul(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Multiply vector z = x*y.
void Neg(int n, T *x, const int incx)
Negate x = -x.
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 Zero(int n, T *x, const int incx)
Zero vector.
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
void Sadd(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Add vector y = alpha + x.
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
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.