64 const TiXmlElement *pMapping)
71 "Mapping X = X(x,z) needs 3 velocity components.");
78 int physTot =
m_fields[0]->GetTotPoints();
96 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();
122 Vmath::Vsub(physTot, outarray[0], 1, wk, 1, outarray[0], 1);
135 int physTot =
m_fields[0]->GetTotPoints();
145 Vmath::Vadd(physTot, inarray[2], 1, outarray[2], 1, outarray[2], 1);
150 int physTot =
m_fields[0]->GetTotPoints();
158 int physTot =
m_fields[0]->GetTotPoints();
168 int physTot =
m_fields[0]->GetTotPoints();
172 for (
int i = 0; i < nvel * nvel; i++)
177 for (
int i = 1; i < nvel; i++)
179 Vmath::Sadd(physTot, 1.0, outarray[i + nvel * i], 1,
180 outarray[i + nvel * i], 1);
186 Vmath::Vcopy(physTot, wk, 1, outarray[0 * nvel + 2], 1);
187 Vmath::Vcopy(physTot, wk, 1, outarray[2 * nvel + 0], 1);
191 outarray[0 * nvel + 0], 1);
196 Vmath::Vadd(physTot, wk, 1, outarray[2 * nvel + 2], 1,
197 outarray[2 * nvel + 2], 1);
203 int physTot =
m_fields[0]->GetTotPoints();
207 for (
int i = 0; i < nvel * nvel; i++)
212 for (
int i = 0; i < nvel; i++)
214 Vmath::Sadd(physTot, 1.0, outarray[i + nvel * i], 1,
215 outarray[i + nvel * i], 1);
222 Vmath::Vcopy(physTot, wk, 1, outarray[0 * nvel + 2], 1);
223 Vmath::Vcopy(physTot, wk, 1, outarray[2 * nvel + 0], 1);
228 Vmath::Vadd(physTot, wk, 1, outarray[0 * nvel + 0], 1,
229 outarray[0 * nvel + 0], 1);
233 Vmath::Vdiv(physTot, outarray[0 * nvel + 0], 1, wk, 1,
234 outarray[0 * nvel + 0], 1);
241 int physTot =
m_fields[0]->GetTotPoints();
247 Vmath::Vmul(physTot, wk, 1, inarray[2], 1, outarray[0], 1);
248 Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[2], 1);
254 Vmath::Vadd(physTot, outarray[0], 1, wk, 1, outarray[0], 1);
265 Vmath::Vadd(physTot, wk, 1, outarray[2], 1, outarray[2], 1);
272 int physTot =
m_fields[0]->GetTotPoints();
278 Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[2], 1);
279 Vmath::Vsub(physTot, inarray[2], 1, outarray[2], 1, outarray[2], 1);
282 Vmath::Vmul(physTot, wk, 1, inarray[2], 1, outarray[0], 1);
288 Vmath::Vsub(physTot, wk, 1, outarray[0], 1, outarray[0], 1);
298 int physTot =
m_fields[0]->GetTotPoints();
302 for (
int i = 0; i < nvel; i++)
304 for (
int j = 0; j < nvel; j++)
314 Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[0 * nvel + 0], 1);
316 Vmath::Vvtvp(physTot, wk, 1, inarray[2], 1, outarray[0 * nvel + 0], 1,
317 outarray[0 * nvel + 0], 1);
320 Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[0 * nvel + 2], 1);
322 Vmath::Vvtvp(physTot, wk, 1, inarray[2], 1, outarray[0 * nvel + 2], 1,
323 outarray[0 * nvel + 2], 1);
330 int physTot =
m_fields[0]->GetTotPoints();
334 for (
int i = 0; i < nvel; i++)
336 for (
int j = 0; j < nvel; j++)
346 Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[0 * nvel + 0], 1);
350 Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[0 * nvel + 2], 1);
351 Vmath::Vcopy(physTot, outarray[0 * nvel + 2], 1, outarray[2 * nvel + 0], 1);
355 Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[2 * nvel + 2], 1);
360 int phystot =
m_fields[0]->GetTotPoints();
368 bool waveSpace =
m_fields[0]->GetWaveSpace();
384 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_ApplyChristoffelCovar(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) override
GLOBAL_MAPPING_EXPORT void v_UpdateGeomInfo() override
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_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_ContravarFromCartesian(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) override
static std::string className
Name of the class.
GLOBAL_MAPPING_EXPORT void v_DotGradJacobian(const Array< OneD, Array< OneD, NekDouble > > &inarray, 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_CovarToCartesian(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
MappingXofXZ(const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields)
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.
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_GetMetricTensor(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 Vvtvp(int n, const T *w, const int incw, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
vvtvp (vector times vector plus vector): z = w*x + y
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 Vdiv(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 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.