Nektar++
Static Public Member Functions | Static Public Attributes | Protected Member Functions | Friends | List of all members
Nektar::GlobalMapping::MappingXYofZ Class Reference

#include <MappingXYofZ.h>

Inheritance diagram for Nektar::GlobalMapping::MappingXYofZ:
[legend]

Static Public Member Functions

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. More...
 
- Static Public Member Functions inherited from Nektar::GlobalMapping::Mapping
static GLOBAL_MAPPING_EXPORT MappingSharedPtr Load (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields)
 Return a pointer to the mapping, creating it on first call. More...
 

Static Public Attributes

static std::string className
 Name of the class. More...
 

Protected Member Functions

 MappingXYofZ (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields)
 
virtual GLOBAL_MAPPING_EXPORT void v_InitObject (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const TiXmlElement *pMapping) override
 
virtual GLOBAL_MAPPING_EXPORT void v_ContravarToCartesian (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray) override
 
virtual GLOBAL_MAPPING_EXPORT void v_CovarToCartesian (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray) override
 
virtual GLOBAL_MAPPING_EXPORT void v_ContravarFromCartesian (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray) override
 
virtual GLOBAL_MAPPING_EXPORT void v_CovarFromCartesian (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray) override
 
virtual GLOBAL_MAPPING_EXPORT void v_GetJacobian (Array< OneD, NekDouble > &outarray) override
 
virtual GLOBAL_MAPPING_EXPORT void v_DotGradJacobian (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, NekDouble > &outarray) override
 
virtual GLOBAL_MAPPING_EXPORT void v_GetMetricTensor (Array< OneD, Array< OneD, NekDouble >> &outarray) override
 
virtual GLOBAL_MAPPING_EXPORT void v_GetInvMetricTensor (Array< OneD, Array< OneD, NekDouble >> &outarray) override
 
virtual GLOBAL_MAPPING_EXPORT void v_ApplyChristoffelContravar (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray) override
 
virtual GLOBAL_MAPPING_EXPORT void v_ApplyChristoffelCovar (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray) override
 
virtual GLOBAL_MAPPING_EXPORT void v_UpdateGeomInfo () override
 
- Protected Member Functions inherited from Nektar::GlobalMapping::Mapping
GLOBAL_MAPPING_EXPORT Mapping (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields)
 Constructor. More...
 
GLOBAL_MAPPING_EXPORT void EvaluateFunction (Array< OneD, MultiRegions::ExpListSharedPtr > pFields, LibUtilities::SessionReaderSharedPtr pSession, std::string pFieldName, Array< OneD, NekDouble > &pArray, const std::string &pFunctionName, NekDouble pTime=NekDouble(0))
 
GLOBAL_MAPPING_EXPORT void EvaluateTimeFunction (LibUtilities::SessionReaderSharedPtr pSession, std::string pFieldName, Array< OneD, NekDouble > &pArray, const std::string &pFunctionName, NekDouble pTime=NekDouble(0))
 
virtual GLOBAL_MAPPING_EXPORT void v_GetCartesianCoordinates (Array< OneD, NekDouble > &out0, Array< OneD, NekDouble > &out1, Array< OneD, NekDouble > &out2)
 
virtual GLOBAL_MAPPING_EXPORT void v_GetCoordVelocity (Array< OneD, Array< OneD, NekDouble >> &outarray)
 
virtual GLOBAL_MAPPING_EXPORT void v_LowerIndex (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
 
virtual GLOBAL_MAPPING_EXPORT void v_RaiseIndex (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
 
virtual GLOBAL_MAPPING_EXPORT void v_Divergence (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, NekDouble > &outarray)
 
virtual GLOBAL_MAPPING_EXPORT void v_VelocityLaplacian (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray, const NekDouble alpha)
 
virtual GLOBAL_MAPPING_EXPORT void v_gradgradU (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
 
virtual GLOBAL_MAPPING_EXPORT void v_CurlCurlField (Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray, const bool generalized)
 
virtual GLOBAL_MAPPING_EXPORT void v_UpdateMapping (const NekDouble time, const Array< OneD, Array< OneD, NekDouble >> &coords=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble >> &coordsVel=NullNekDoubleArrayOfArray)
 
virtual GLOBAL_MAPPING_EXPORT void v_UpdateBCs (const NekDouble time)
 

Friends

class MemoryManager< MappingXYofZ >
 

Additional Inherited Members

- Public Member Functions inherited from Nektar::GlobalMapping::Mapping
virtual GLOBAL_MAPPING_EXPORT ~Mapping ()
 Destructor. More...
 
GLOBAL_MAPPING_EXPORT void InitObject (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const TiXmlElement *pMapping)
 Initialise the mapping object. More...
 
GLOBAL_MAPPING_EXPORT void ReplaceField (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields)
 Replace the Expansion List used by the mapping. More...
 
GLOBAL_MAPPING_EXPORT void Output (LibUtilities::FieldMetaDataMap &fieldMetaDataMap, const std::string &outname)
 Output function called when a chk or fld file is written. More...
 
GLOBAL_MAPPING_EXPORT void ContravarToCartesian (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
 Convert a contravariant vector to the Cartesian system. More...
 
GLOBAL_MAPPING_EXPORT void CovarToCartesian (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
 Convert a covariant vector to the Cartesian system. More...
 
GLOBAL_MAPPING_EXPORT void ContravarFromCartesian (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
 Convert a contravariant vector to the transformed system. More...
 
GLOBAL_MAPPING_EXPORT void CovarFromCartesian (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
 Convert a covariant vector to the transformed system. More...
 
GLOBAL_MAPPING_EXPORT void GetCartesianCoordinates (Array< OneD, NekDouble > &out0, Array< OneD, NekDouble > &out1, Array< OneD, NekDouble > &out2)
 Get the Cartesian coordinates in the field. More...
 
GLOBAL_MAPPING_EXPORT void GetJacobian (Array< OneD, NekDouble > &outarray)
 Get the Jacobian of the transformation. More...
 
GLOBAL_MAPPING_EXPORT void DotGradJacobian (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, NekDouble > &outarray)
 Calculate the dot product with the gradient of the Jacobian. More...
 
GLOBAL_MAPPING_EXPORT void GetMetricTensor (Array< OneD, Array< OneD, NekDouble >> &outarray)
 Get the metric tensor \(g_{ij}\). More...
 
GLOBAL_MAPPING_EXPORT void GetInvMetricTensor (Array< OneD, Array< OneD, NekDouble >> &outarray)
 Get the inverse of metric tensor \(g^{ij}\). More...
 
GLOBAL_MAPPING_EXPORT void LowerIndex (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
 Lower index of vector: \(v_{i} = g_{ij}*v^{j}\). More...
 
GLOBAL_MAPPING_EXPORT void RaiseIndex (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
 Raise index of vector: \(v^{i} = g^{ij}*v_{j}\). More...
 
GLOBAL_MAPPING_EXPORT void ApplyChristoffelContravar (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
 Apply the Christoffel symbols to a contravariant vector. More...
 
GLOBAL_MAPPING_EXPORT void ApplyChristoffelCovar (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
 Apply the Christoffel symbols to a covariant vector. More...
 
GLOBAL_MAPPING_EXPORT void GetCoordVelocity (Array< OneD, Array< OneD, NekDouble >> &outarray)
 Obtain the velocity of the coordinates. More...
 
GLOBAL_MAPPING_EXPORT void Divergence (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, NekDouble > &outarray)
 Calculate the generalised divergence operator. More...
 
GLOBAL_MAPPING_EXPORT void VelocityLaplacian (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray, const NekDouble alpha=0.0)
 Generalised (correction to the) velocity Laplacian operator. More...
 
GLOBAL_MAPPING_EXPORT void gradgradU (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
 Second order covariant derivatives of a contravariant vector. More...
 
GLOBAL_MAPPING_EXPORT void CurlCurlField (Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray, const bool generalized)
 CurlCurl calculated on the whole field. More...
 
GLOBAL_MAPPING_EXPORT bool IsTimeDependent ()
 Get flag defining if mapping is time-dependent. More...
 
GLOBAL_MAPPING_EXPORT void SetTimeDependent (const bool value)
 Set flag defining if mapping is time-dependent. More...
 
GLOBAL_MAPPING_EXPORT bool IsFromFunction ()
 Get flag defining if mapping is defined by a function. More...
 
GLOBAL_MAPPING_EXPORT void SetFromFunction (const bool value)
 Set flag defining if mapping is defined by a function. More...
 
GLOBAL_MAPPING_EXPORT bool HasConstantJacobian ()
 Get flag defining if mapping has constant Jacobian. More...
 
GLOBAL_MAPPING_EXPORT bool IsDefined ()
 Get flag determining if the mapping was defined or is trivial. More...
 
GLOBAL_MAPPING_EXPORT void UpdateBCs (const NekDouble time)
 Update the Dirichlet Boundary Conditions when using Mappings. More...
 
GLOBAL_MAPPING_EXPORT void UpdateMapping (const NekDouble time, const Array< OneD, Array< OneD, NekDouble >> &coords=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble >> &coordsVel=NullNekDoubleArrayOfArray)
 Update the Mapping with new coordinates. More...
 
GLOBAL_MAPPING_EXPORT void UpdateGeomInfo ()
 Recompute the metric terms of the Mapping. More...
 
- Protected Attributes inherited from Nektar::GlobalMapping::Mapping
LibUtilities::SessionReaderSharedPtr m_session
 Session reader. More...
 
LibUtilities::FieldIOSharedPtr m_fld
 
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
 
Array< OneD, Array< OneD, NekDouble > > m_coords
 Array with the Cartesian coordinates. More...
 
Array< OneD, Array< OneD, NekDouble > > m_coordsVel
 Array with the velocity of the coordinates. More...
 
Array< OneD, Array< OneD, NekDouble > > m_GeometricInfo
 Array with metric terms of the mapping. More...
 
int m_nConvectiveFields
 Number of velocity components. More...
 
std::string m_funcName
 Name of the function containing the coordinates. More...
 
std::string m_velFuncName
 Name of the function containing the velocity of the coordinates. More...
 
bool m_constantJacobian
 Flag defining if the Jacobian is constant. More...
 
bool m_timeDependent
 Flag defining if the Mapping is time-dependent. More...
 
bool m_fromFunction
 Flag defining if the Mapping is defined by a function. More...
 
Array< OneD, Array< OneD, NekDouble > > m_wk1
 
Array< OneD, Array< OneD, NekDouble > > m_wk2
 
Array< OneD, Array< OneD, NekDouble > > m_tmp
 
- Static Protected Attributes inherited from Nektar::GlobalMapping::Mapping
static MappingSharedPtr m_mappingPtr = MappingSharedPtr()
 
static bool m_init = false
 
static bool m_isDefined = false
 

Detailed Description

This class implements a constant-Jacobian mapping defined by

\[ \bar{x} = \bar{x}(x,z) = x + f(z) \]

\[ \bar{y} = \bar{y}(y,z) = y + g(z) \]

\[ \bar{z} = z \]

where \((\bar{x},\bar{y},\bar{z})\) are the Cartesian (physical) coordinates and \((x,y,z)\) are the transformed (computational) coordinates.

Definition at line 50 of file MappingXYofZ.h.

Constructor & Destructor Documentation

◆ MappingXYofZ()

Nektar::GlobalMapping::MappingXYofZ::MappingXYofZ ( const LibUtilities::SessionReaderSharedPtr &  pSession,
const Array< OneD, MultiRegions::ExpListSharedPtr > &  pFields 
)
protected

Definition at line 59 of file MappingXYofZ.cpp.

62  : Mapping(pSession, pFields)
63 {
64 }
Nektar::GlobalMapping::Mapping::Mapping
GLOBAL_MAPPING_EXPORT Mapping(const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields)
Constructor.
Definition: Mapping.cpp:59

Member Function Documentation

◆ create()

static GLOBAL_MAPPING_EXPORT MappingSharedPtr Nektar::GlobalMapping::MappingXYofZ::create ( const LibUtilities::SessionReaderSharedPtr &  pSession,
const Array< OneD, MultiRegions::ExpListSharedPtr > &  pFields,
const TiXmlElement *  pMapping 
)
inlinestatic

Creates an instance of this class.

Definition at line 57 of file MappingXYofZ.h.

61  {
62  MappingSharedPtr p =
63  MemoryManager<MappingXYofZ>::AllocateSharedPtr(pSession, pFields);
64  p->InitObject(pFields, pMapping);
65  return p;
66  }
Nektar::MemoryManager::AllocateSharedPtr
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:168
CellMLToNektar.cellml_metadata.p
p
Definition: cellml_metadata.py:350
Nektar::GlobalMapping::MappingSharedPtr
GLOBAL_MAPPING_EXPORT typedef std::shared_ptr< Mapping > MappingSharedPtr
A shared pointer to a Mapping object.
Definition: Mapping.h:50

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::GlobalMapping::MappingSharedPtr, and CellMLToNektar.cellml_metadata::p.

◆ v_ApplyChristoffelContravar()

void Nektar::GlobalMapping::MappingXYofZ::v_ApplyChristoffelContravar ( const Array< OneD, Array< OneD, NekDouble >> &  inarray,
Array< OneD, Array< OneD, NekDouble >> &  outarray 
)
overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 248 of file MappingXYofZ.cpp.

251 {
252  int physTot = m_fields[0]->GetTotPoints();
253  int nvel = m_nConvectiveFields;
254 
255  for (int i = 0; i < nvel; i++)
256  {
257  for (int j = 0; j < nvel; j++)
258  {
259  outarray[i * nvel + j] = Array<OneD, NekDouble>(physTot, 0.0);
260  }
261  }
262 
263  // Calculate non-zero terms
264 
265  // outarray(0,2) = U3 * fzz
266  Vmath::Vmul(physTot, m_GeometricInfo[1], 1, inarray[2], 1,
267  outarray[0 * nvel + 2], 1);
268 
269  // outarray(1,2) = U3 * gzz
270  Vmath::Vmul(physTot, m_GeometricInfo[4], 1, inarray[2], 1,
271  outarray[1 * nvel + 2], 1);
272 }
Nektar::GlobalMapping::Mapping::m_nConvectiveFields
int m_nConvectiveFields
Number of velocity components.
Definition: Mapping.h:414
Nektar::GlobalMapping::Mapping::m_GeometricInfo
Array< OneD, Array< OneD, NekDouble > > m_GeometricInfo
Array with metric terms of the mapping.
Definition: Mapping.h:412
Nektar::GlobalMapping::Mapping::m_fields
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
Definition: Mapping.h:406
Vmath::Vmul
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.
Definition: Vmath.cpp:209

References Nektar::GlobalMapping::Mapping::m_fields, Nektar::GlobalMapping::Mapping::m_GeometricInfo, Nektar::GlobalMapping::Mapping::m_nConvectiveFields, and Vmath::Vmul().

◆ v_ApplyChristoffelCovar()

void Nektar::GlobalMapping::MappingXYofZ::v_ApplyChristoffelCovar ( const Array< OneD, Array< OneD, NekDouble >> &  inarray,
Array< OneD, Array< OneD, NekDouble >> &  outarray 
)
overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 274 of file MappingXYofZ.cpp.

277 {
278  int physTot = m_fields[0]->GetTotPoints();
279  int nvel = m_nConvectiveFields;
280 
281  for (int i = 0; i < nvel; i++)
282  {
283  for (int j = 0; j < nvel; j++)
284  {
285  outarray[i * nvel + j] = Array<OneD, NekDouble>(physTot, 0.0);
286  }
287  }
288 
289  // Calculate non-zero terms
290 
291  // outarray(2,2) = U1 * fzz + U^2 * gzz
292  Vmath::Vmul(physTot, m_GeometricInfo[1], 1, inarray[0], 1,
293  outarray[2 * nvel + 2], 1);
294  Vmath::Vvtvp(physTot, m_GeometricInfo[4], 1, inarray[1], 1,
295  outarray[2 * nvel + 2], 1, outarray[2 * nvel + 2], 1);
296 }
Vmath::Vvtvp
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
Definition: Vmath.cpp:574

References Nektar::GlobalMapping::Mapping::m_fields, Nektar::GlobalMapping::Mapping::m_GeometricInfo, Nektar::GlobalMapping::Mapping::m_nConvectiveFields, Vmath::Vmul(), and Vmath::Vvtvp().

◆ v_ContravarFromCartesian()

void Nektar::GlobalMapping::MappingXYofZ::v_ContravarFromCartesian ( const Array< OneD, Array< OneD, NekDouble >> &  inarray,
Array< OneD, Array< OneD, NekDouble >> &  outarray 
)
overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 119 of file MappingXYofZ.cpp.

122 {
123  int physTot = m_fields[0]->GetTotPoints();
124  Array<OneD, NekDouble> wk(physTot, 0.0);
125 
126  // U1 = u1 - fz * u3
127  Vmath::Vmul(physTot, m_GeometricInfo[0], 1, inarray[2], 1, wk, 1);
128  Vmath::Vsub(physTot, inarray[0], 1, wk, 1, outarray[0], 1);
129 
130  // U2 = u2 - gz*u3
131  Vmath::Vmul(physTot, m_GeometricInfo[3], 1, inarray[2], 1, wk, 1);
132  Vmath::Vsub(physTot, inarray[1], 1, wk, 1, outarray[1], 1);
133 
134  // U3 = u3
135  Vmath::Vcopy(physTot, inarray[2], 1, outarray[2], 1);
136 }
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1255
Vmath::Vsub
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.
Definition: Vmath.cpp:419

References Nektar::GlobalMapping::Mapping::m_fields, Nektar::GlobalMapping::Mapping::m_GeometricInfo, Vmath::Vcopy(), Vmath::Vmul(), and Vmath::Vsub().

◆ v_ContravarToCartesian()

void Nektar::GlobalMapping::MappingXYofZ::v_ContravarToCartesian ( const Array< OneD, Array< OneD, NekDouble >> &  inarray,
Array< OneD, Array< OneD, NekDouble >> &  outarray 
)
overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 81 of file MappingXYofZ.cpp.

84 {
85  int physTot = m_fields[0]->GetTotPoints();
86 
87  // U1 = u1 + fz*u3
88  Vmath::Vvtvp(physTot, m_GeometricInfo[0], 1, inarray[2], 1, inarray[0], 1,
89  outarray[0], 1);
90 
91  // U2 = u2 + gz*u3
92  Vmath::Vvtvp(physTot, m_GeometricInfo[3], 1, inarray[2], 1, inarray[1], 1,
93  outarray[1], 1);
94 
95  // U3 = u3
96  Vmath::Vcopy(physTot, inarray[2], 1, outarray[2], 1);
97 }

References Nektar::GlobalMapping::Mapping::m_fields, Nektar::GlobalMapping::Mapping::m_GeometricInfo, Vmath::Vcopy(), and Vmath::Vvtvp().

◆ v_CovarFromCartesian()

void Nektar::GlobalMapping::MappingXYofZ::v_CovarFromCartesian ( const Array< OneD, Array< OneD, NekDouble >> &  inarray,
Array< OneD, Array< OneD, NekDouble >> &  outarray 
)
overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 138 of file MappingXYofZ.cpp.

141 {
142  int physTot = m_fields[0]->GetTotPoints();
143 
144  // U1 = u1
145  Vmath::Vcopy(physTot, inarray[0], 1, outarray[0], 1);
146 
147  // U2 = u2
148  Vmath::Vcopy(physTot, inarray[1], 1, outarray[1], 1);
149 
150  // U3 = u3 + fz*u1 + gz*u2
151  Vmath::Vmul(physTot, m_GeometricInfo[0], 1, inarray[0], 1, outarray[2], 1);
152  Vmath::Vvtvp(physTot, m_GeometricInfo[3], 1, inarray[1], 1, outarray[2], 1,
153  outarray[2], 1);
154  Vmath::Vadd(physTot, inarray[2], 1, outarray[2], 1, outarray[2], 1);
155 }
Vmath::Vadd
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.
Definition: Vmath.cpp:359

References Nektar::GlobalMapping::Mapping::m_fields, Nektar::GlobalMapping::Mapping::m_GeometricInfo, Vmath::Vadd(), Vmath::Vcopy(), Vmath::Vmul(), and Vmath::Vvtvp().

◆ v_CovarToCartesian()

void Nektar::GlobalMapping::MappingXYofZ::v_CovarToCartesian ( const Array< OneD, Array< OneD, NekDouble >> &  inarray,
Array< OneD, Array< OneD, NekDouble >> &  outarray 
)
overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 99 of file MappingXYofZ.cpp.

102 {
103  int physTot = m_fields[0]->GetTotPoints();
104  Array<OneD, NekDouble> wk(physTot, 0.0);
105 
106  // U1 = u1
107  Vmath::Vcopy(physTot, inarray[0], 1, outarray[0], 1);
108 
109  // U2 = u2
110  Vmath::Vcopy(physTot, inarray[1], 1, outarray[1], 1);
111 
112  // U3 = u3 - fz*u1 - gz*u2
113  Vmath::Vmul(physTot, m_GeometricInfo[0], 1, inarray[0], 1, wk, 1);
114  Vmath::Vsub(physTot, inarray[2], 1, wk, 1, outarray[2], 1);
115  Vmath::Vmul(physTot, m_GeometricInfo[3], 1, inarray[1], 1, wk, 1);
116  Vmath::Vsub(physTot, inarray[2], 1, wk, 1, outarray[2], 1);
117 }

References Nektar::GlobalMapping::Mapping::m_fields, Nektar::GlobalMapping::Mapping::m_GeometricInfo, Vmath::Vcopy(), Vmath::Vmul(), and Vmath::Vsub().

◆ v_DotGradJacobian()

void Nektar::GlobalMapping::MappingXYofZ::v_DotGradJacobian ( const Array< OneD, Array< OneD, NekDouble >> &  inarray,
Array< OneD, NekDouble > &  outarray 
)
overrideprotectedvirtual

Reimplemented from Nektar::GlobalMapping::Mapping.

Definition at line 163 of file MappingXYofZ.cpp.

166 {
167  boost::ignore_unused(inarray);
168 
169  int physTot = m_fields[0]->GetTotPoints();
170  Vmath::Zero(physTot, outarray, 1);
171 }
Vmath::Zero
void Zero(int n, T *x, const int incx)
Zero vector.
Definition: Vmath.cpp:492

References Nektar::GlobalMapping::Mapping::m_fields, and Vmath::Zero().

◆ v_GetInvMetricTensor()

void Nektar::GlobalMapping::MappingXYofZ::v_GetInvMetricTensor ( Array< OneD, Array< OneD, NekDouble >> &  outarray)
overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 205 of file MappingXYofZ.cpp.

207 {
208  int physTot = m_fields[0]->GetTotPoints();
209  int nvel = m_nConvectiveFields;
210  Array<OneD, NekDouble> wk(physTot, 0.0);
211 
212  for (int i = 0; i < nvel * nvel; i++)
213  {
214  outarray[i] = Array<OneD, NekDouble>(physTot, 0.0);
215  }
216  // Fill diagonal with 1.0
217  for (int i = 0; i < nvel; i++)
218  {
219  Vmath::Sadd(physTot, 1.0, outarray[i * nvel + i], 1,
220  outarray[i * nvel + i], 1);
221  }
222 
223  // G^{11} = 1+fz^2
224  Vmath::Vadd(physTot, outarray[0 * nvel + 0], 1, m_GeometricInfo[2], 1,
225  outarray[0 * nvel + 0], 1);
226 
227  // G^{22} = 1+gz^2
228  Vmath::Vadd(physTot, outarray[1 * nvel + 1], 1, m_GeometricInfo[5], 1,
229  outarray[1 * nvel + 1], 1);
230 
231  // G^{12} and G^{21} = fz*gz
232  Vmath::Vcopy(physTot, m_GeometricInfo[6], 1, outarray[0 * nvel + 1], 1);
233  Vmath::Vcopy(physTot, outarray[0 * nvel + 1], 1, outarray[1 * nvel + 0], 1);
234 
235  // G^{13} and G^{31} = -fz
236  Vmath::Vcopy(physTot, m_GeometricInfo[0], 1, wk, 1); // fz
237  Vmath::Neg(physTot, wk, 1);
238  Vmath::Vcopy(physTot, wk, 1, outarray[0 * nvel + 2], 1);
239  Vmath::Vcopy(physTot, wk, 1, outarray[2 * nvel + 0], 1);
240 
241  // G^{23} and G^{32} = -gz
242  Vmath::Vcopy(physTot, m_GeometricInfo[3], 1, wk, 1); // fz
243  Vmath::Neg(physTot, wk, 1);
244  Vmath::Vcopy(physTot, wk, 1, outarray[1 * nvel + 2], 1);
245  Vmath::Vcopy(physTot, wk, 1, outarray[2 * nvel + 1], 1);
246 }
Vmath::Neg
void Neg(int n, T *x, const int incx)
Negate x = -x.
Definition: Vmath.cpp:518
Vmath::Sadd
void Sadd(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Add scalar y = alpha + x.
Definition: Vmath.cpp:384

References Nektar::GlobalMapping::Mapping::m_fields, Nektar::GlobalMapping::Mapping::m_GeometricInfo, Nektar::GlobalMapping::Mapping::m_nConvectiveFields, Vmath::Neg(), Vmath::Sadd(), Vmath::Vadd(), and Vmath::Vcopy().

◆ v_GetJacobian()

void Nektar::GlobalMapping::MappingXYofZ::v_GetJacobian ( Array< OneD, NekDouble > &  outarray)
overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 157 of file MappingXYofZ.cpp.

158 {
159  int physTot = m_fields[0]->GetTotPoints();
160  Vmath::Fill(physTot, 1.0, outarray, 1);
161 }
Vmath::Fill
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
Definition: Vmath.cpp:45

References Vmath::Fill(), and Nektar::GlobalMapping::Mapping::m_fields.

◆ v_GetMetricTensor()

void Nektar::GlobalMapping::MappingXYofZ::v_GetMetricTensor ( Array< OneD, Array< OneD, NekDouble >> &  outarray)
overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 173 of file MappingXYofZ.cpp.

175 {
176  int physTot = m_fields[0]->GetTotPoints();
177  int nvel = m_nConvectiveFields;
178 
179  for (int i = 0; i < nvel * nvel; i++)
180  {
181  outarray[i] = Array<OneD, NekDouble>(physTot, 0.0);
182  }
183  // Fill diagonal with 1.0
184  for (int i = 0; i < nvel; i++)
185  {
186  Vmath::Sadd(physTot, 1.0, outarray[i * nvel + i], 1,
187  outarray[i * nvel + i], 1);
188  }
189 
190  // G_{13} and G_{31} = fz
191  Vmath::Vcopy(physTot, m_GeometricInfo[0], 1, outarray[0 * nvel + 2], 1);
192  Vmath::Vcopy(physTot, m_GeometricInfo[0], 1, outarray[2 * nvel + 0], 1);
193 
194  // G_{23} and G_{32} = gz
195  Vmath::Vcopy(physTot, m_GeometricInfo[3], 1, outarray[1 * nvel + 2], 1);
196  Vmath::Vcopy(physTot, m_GeometricInfo[3], 1, outarray[2 * nvel + 1], 1);
197 
198  // G^{33} = (1+fz^2 + gz^2)
199  Vmath::Vadd(physTot, m_GeometricInfo[2], 1, outarray[2 * nvel + 2], 1,
200  outarray[2 * nvel + 2], 1);
201  Vmath::Vadd(physTot, m_GeometricInfo[5], 1, outarray[2 * nvel + 2], 1,
202  outarray[2 * nvel + 2], 1);
203 }

References Nektar::GlobalMapping::Mapping::m_fields, Nektar::GlobalMapping::Mapping::m_GeometricInfo, Nektar::GlobalMapping::Mapping::m_nConvectiveFields, Vmath::Sadd(), Vmath::Vadd(), and Vmath::Vcopy().

◆ v_InitObject()

void Nektar::GlobalMapping::MappingXYofZ::v_InitObject ( const Array< OneD, MultiRegions::ExpListSharedPtr > &  pFields,
const TiXmlElement *  pMapping 
)
overrideprotectedvirtual

This function initialises the Mapping object. It computes the coordinates and velocity coordinates, initialises the workspace variables, and calls UpdateGeomInfo, which will perform the calculations specific for each type of Mapping.

Parameters
pFieldsExpList array used in the mapping
pMappingxml element describing the mapping

Reimplemented from Nektar::GlobalMapping::Mapping.

Definition at line 69 of file MappingXYofZ.cpp.

72 {
73  Mapping::v_InitObject(pFields, pMapping);
74 
75  m_constantJacobian = true;
76 
77  ASSERTL0(m_nConvectiveFields == 3,
78  "Mapping X = x + f(z), Y = y+g(z) needs 3 velocity components.");
79 }
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:215
Nektar::GlobalMapping::Mapping::v_InitObject
virtual GLOBAL_MAPPING_EXPORT void v_InitObject(const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const TiXmlElement *pMapping)
Definition: Mapping.cpp:101
Nektar::GlobalMapping::Mapping::m_constantJacobian
bool m_constantJacobian
Flag defining if the Jacobian is constant.
Definition: Mapping.h:423

References ASSERTL0, Nektar::GlobalMapping::Mapping::m_constantJacobian, Nektar::GlobalMapping::Mapping::m_nConvectiveFields, and Nektar::GlobalMapping::Mapping::v_InitObject().

◆ v_UpdateGeomInfo()

void Nektar::GlobalMapping::MappingXYofZ::v_UpdateGeomInfo ( )
overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 298 of file MappingXYofZ.cpp.

299 {
300  int phystot = m_fields[0]->GetTotPoints();
301  // Allocation of geometry memory
302  m_GeometricInfo = Array<OneD, Array<OneD, NekDouble>>(7);
303  for (int i = 0; i < m_GeometricInfo.size(); i++)
304  {
305  m_GeometricInfo[i] = Array<OneD, NekDouble>(phystot, 0.0);
306  }
307 
308  bool waveSpace = m_fields[0]->GetWaveSpace();
309  m_fields[0]->SetWaveSpace(false);
310 
311  // Calculate derivatives of x transformation --> m_GeometricInfo 0-1
312  m_fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[2], m_coords[0],
313  m_GeometricInfo[0]);
314  m_fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[2], m_GeometricInfo[0],
315  m_GeometricInfo[1]);
316  // m_GeometricInfo[2] = fz^2
317  Vmath::Vmul(phystot, m_GeometricInfo[0], 1, m_GeometricInfo[0], 1,
318  m_GeometricInfo[2], 1);
319 
320  // Calculate derivatives of transformation -> m_GeometricInfo 3-4
321  m_fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[2], m_coords[1],
322  m_GeometricInfo[3]);
323  m_fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[2], m_GeometricInfo[3],
324  m_GeometricInfo[4]);
325  // m_GeometricInfo[5] = gz^2
326  Vmath::Vmul(phystot, m_GeometricInfo[3], 1, m_GeometricInfo[3], 1,
327  m_GeometricInfo[5], 1);
328 
329  // m_GeometricInfo[6] = gz*fz
330  Vmath::Vmul(phystot, m_GeometricInfo[0], 1, m_GeometricInfo[3], 1,
331  m_GeometricInfo[6], 1);
332 
333  m_fields[0]->SetWaveSpace(waveSpace);
334 }
Nektar::GlobalMapping::Mapping::m_coords
Array< OneD, Array< OneD, NekDouble > > m_coords
Array with the Cartesian coordinates.
Definition: Mapping.h:408
Nektar::MultiRegions::DirCartesianMap
MultiRegions::Direction const DirCartesianMap[]
Definition: ExpList.h:91

References Nektar::MultiRegions::DirCartesianMap, Nektar::GlobalMapping::Mapping::m_coords, Nektar::GlobalMapping::Mapping::m_fields, Nektar::GlobalMapping::Mapping::m_GeometricInfo, and Vmath::Vmul().

Friends And Related Function Documentation

◆ MemoryManager< MappingXYofZ >

friend class MemoryManager< MappingXYofZ >
friend

Definition at line 1 of file MappingXYofZ.h.

Member Data Documentation

◆ className

std::string Nektar::GlobalMapping::MappingXYofZ::className
static
Initial value:
=
GetMappingFactory().RegisterCreatorFunction("XYofZ", MappingXYofZ::create,
"X = x + f(z), Y = y +g(z)")
Nektar::GlobalMapping::MappingXYofZ::create
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.
Definition: MappingXYofZ.h:57
Nektar::LibUtilities::NekFactory::RegisterCreatorFunction
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
Definition: NekFactory.hpp:198
Nektar::GlobalMapping::GetMappingFactory
MappingFactory & GetMappingFactory()
Declaration of the mapping factory singleton.
Definition: Mapping.cpp:53

Name of the class.

Definition at line 69 of file MappingXYofZ.h.

Generated on Sat Mar 25 2023 12:35:57 for Nektar++ by doxygen 1.9.1