#include <MappingXofXZ.h>

Inheritance diagram for Nektar::GlobalMapping::MappingXofXZ:

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
	MappingXofXZ (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_LowerIndex (const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray) override

virtual GLOBAL_MAPPING_EXPORT void	v_RaiseIndex (const Array< OneD, Array< OneD, NekDouble >> &inarray, 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_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< MappingXofXZ >

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 mapping defined by a transformation of the type

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

\[ \bar{y} = y \]

\[ \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 51 of file MappingXofXZ.h.

Constructor & Destructor Documentation

◆ MappingXofXZ()

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

protected

Definition at line 57 of file MappingXofXZ.cpp.

     : Mapping(pSession, pFields)
 {
 }

Member Function Documentation

◆ create()

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

inlinestatic

Creates an instance of this class.

Definition at line 58 of file MappingXofXZ.h.

     {
         MappingSharedPtr p =
             MemoryManager<MappingXofXZ>::AllocateSharedPtr(pSession, pFields);
         p->InitObject(pFields, pMapping);
         return p;
     }

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

◆ v_ApplyChristoffelContravar()

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

overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 299 of file MappingXofXZ.cpp.

 {
     int physTot = m_fields[0]->GetTotPoints();
     int nvel    = m_nConvectiveFields;
     Array<OneD, NekDouble> wk(physTot, 0.0);
  
     for (int i = 0; i < nvel; i++)
     {
         for (int j = 0; j < nvel; j++)
         {
             outarray[i * nvel + j] = Array<OneD, NekDouble>(physTot, 0.0);
         }
     }
  
     // Calculate non-zero terms
  
     // outarray(0,0) = U1 * fxx/fx + U3 * fxz/fx
     Vmath::Vdiv(physTot, m_GeometricInfo[2], 1, m_GeometricInfo[0], 1, wk, 1);
     Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[0 * nvel + 0], 1);
     Vmath::Vdiv(physTot, m_GeometricInfo[3], 1, m_GeometricInfo[0], 1, wk, 1);
     Vmath::Vvtvp(physTot, wk, 1, inarray[2], 1, outarray[0 * nvel + 0], 1,
                  outarray[0 * nvel + 0], 1);
  
     // outarray(0,2) = U1 * fxz/fx + U3 * fzz/fx
     Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[0 * nvel + 2], 1);
     Vmath::Vdiv(physTot, m_GeometricInfo[4], 1, m_GeometricInfo[0], 1, wk, 1);
     Vmath::Vvtvp(physTot, wk, 1, inarray[2], 1, outarray[0 * nvel + 2], 1,
                  outarray[0 * nvel + 2], 1);
 }

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

◆ v_ApplyChristoffelCovar()

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

overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 331 of file MappingXofXZ.cpp.

 {
     int physTot = m_fields[0]->GetTotPoints();
     int nvel    = m_nConvectiveFields;
     Array<OneD, NekDouble> wk(physTot, 0.0);
  
     for (int i = 0; i < nvel; i++)
     {
         for (int j = 0; j < nvel; j++)
         {
             outarray[i * nvel + j] = Array<OneD, NekDouble>(physTot, 0.0);
         }
     }
  
     // Calculate non-zero terms
  
     // outarray(0,0) = U1 * fxx/fx
     Vmath::Vdiv(physTot, m_GeometricInfo[2], 1, m_GeometricInfo[0], 1, wk, 1);
     Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[0 * nvel + 0], 1);
  
     // outarray(0,2) = outarray(2,0) = U1 * fxz/fx
     Vmath::Vdiv(physTot, m_GeometricInfo[3], 1, m_GeometricInfo[0], 1, wk, 1);
     Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[0 * nvel + 2], 1);
     Vmath::Vcopy(physTot, outarray[0 * nvel + 2], 1, outarray[2 * nvel + 0], 1);
  
     // outarray(2,2) = U1 * fzz/fx
     Vmath::Vdiv(physTot, m_GeometricInfo[4], 1, m_GeometricInfo[0], 1, wk, 1);
     Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[2 * nvel + 2], 1);
 }

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

◆ v_ContravarFromCartesian()

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

overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 116 of file MappingXofXZ.cpp.

 {
     int physTot = m_fields[0]->GetTotPoints();
     Array<OneD, NekDouble> wk(physTot, 0.0);
  
     // U1 = u1/fx - fz/fx * u3
     Vmath::Vdiv(physTot, inarray[0], 1, m_GeometricInfo[0], 1, outarray[0], 1);
     Vmath::Vdiv(physTot, m_GeometricInfo[1], 1, m_GeometricInfo[0], 1, wk, 1);
     Vmath::Vmul(physTot, wk, 1, inarray[2], 1, wk, 1);
     Vmath::Vsub(physTot, outarray[0], 1, wk, 1, outarray[0], 1);
  
     // U2 = u2
     Vmath::Vcopy(physTot, inarray[1], 1, outarray[1], 1);
  
     // U3 = u3
     Vmath::Vcopy(physTot, inarray[2], 1, outarray[2], 1);
 }

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

◆ v_ContravarToCartesian()

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

overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 79 of file MappingXofXZ.cpp.

 {
     int physTot = m_fields[0]->GetTotPoints();
  
     // U1 = fx*u1 + fz*u3
     Vmath::Vmul(physTot, m_GeometricInfo[0], 1, inarray[0], 1, outarray[0], 1);
     Vmath::Vvtvp(physTot, m_GeometricInfo[1], 1, inarray[2], 1, outarray[0], 1,
                  outarray[0], 1);
  
     // U2 = u2
     Vmath::Vcopy(physTot, inarray[1], 1, outarray[1], 1);
  
     // U3 = u3
     Vmath::Vcopy(physTot, inarray[2], 1, outarray[2], 1);
 }

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

◆ v_CovarFromCartesian()

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

overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 136 of file MappingXofXZ.cpp.

 {
     int physTot = m_fields[0]->GetTotPoints();
  
     // U1 = u1*fx
     Vmath::Vmul(physTot, inarray[0], 1, m_GeometricInfo[0], 1, outarray[0], 1);
  
     // U2 = u2
     Vmath::Vcopy(physTot, inarray[1], 1, outarray[1], 1);
  
     // U3 = u3 + fz*u1
     Vmath::Vmul(physTot, m_GeometricInfo[1], 1, inarray[0], 1, outarray[2], 1);
     Vmath::Vadd(physTot, inarray[2], 1, outarray[2], 1, outarray[2], 1);
 }

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

◆ v_CovarToCartesian()

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

overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 97 of file MappingXofXZ.cpp.

 {
     int physTot = m_fields[0]->GetTotPoints();
     Array<OneD, NekDouble> wk(physTot, 0.0);
  
     // U1 = u1/fx
     Vmath::Vdiv(physTot, inarray[0], 1, m_GeometricInfo[0], 1, outarray[0], 1);
  
     // U2 = u2
     Vmath::Vcopy(physTot, inarray[1], 1, outarray[1], 1);
  
     // U3 = u3 - fz/fx*u1
     Vmath::Vdiv(physTot, m_GeometricInfo[1], 1, m_GeometricInfo[0], 1, wk, 1);
     Vmath::Vmul(physTot, wk, 1, inarray[0], 1, wk, 1);
     Vmath::Vsub(physTot, inarray[2], 1, wk, 1, outarray[2], 1);
 }

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

◆ v_DotGradJacobian()

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

overrideprotectedvirtual

Reimplemented from Nektar::GlobalMapping::Mapping.

Definition at line 159 of file MappingXofXZ.cpp.

 {
     int physTot = m_fields[0]->GetTotPoints();
  
     Vmath::Vmul(physTot, m_GeometricInfo[2], 1, inarray[0], 1, outarray, 1);
     Vmath::Vvtvp(physTot, m_GeometricInfo[3], 1, inarray[2], 1, outarray, 1,
                  outarray, 1);
 }

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

◆ v_GetInvMetricTensor()

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

overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 205 of file MappingXofXZ.cpp.

 {
     int physTot = m_fields[0]->GetTotPoints();
     int nvel    = m_nConvectiveFields;
     Array<OneD, NekDouble> wk(physTot, 0.0);
  
     for (int i = 0; i < nvel * nvel; i++)
     {
         outarray[i] = Array<OneD, NekDouble>(physTot, 0.0);
     }
     // Fill diagonal with 1.0
     for (int i = 0; i < nvel; i++)
     {
         Vmath::Sadd(physTot, 1.0, outarray[i + nvel * i], 1,
                     outarray[i + nvel * i], 1);
     }
  
     // G^{13} and G^{31} = -fz/fx
     Vmath::Vdiv(physTot, m_GeometricInfo[1], 1, m_GeometricInfo[0], 1, wk,
                 1); // fz/fx
     Vmath::Neg(physTot, wk, 1);
     Vmath::Vcopy(physTot, wk, 1, outarray[0 * nvel + 2], 1);
     Vmath::Vcopy(physTot, wk, 1, outarray[2 * nvel + 0], 1);
  
     // G^{11} = (1+fz^2)/(fx^2)
     Vmath::Vmul(physTot, m_GeometricInfo[1], 1, m_GeometricInfo[1], 1, wk,
                 1); // fz^2
     Vmath::Vadd(physTot, wk, 1, outarray[0 * nvel + 0], 1,
                 outarray[0 * nvel + 0], 1);
  
     Vmath::Vmul(physTot, m_GeometricInfo[0], 1, m_GeometricInfo[0], 1, wk,
                 1); // fx^2
     Vmath::Vdiv(physTot, outarray[0 * nvel + 0], 1, wk, 1,
                 outarray[0 * nvel + 0], 1);
 }

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

◆ v_GetJacobian()

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

overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 153 of file MappingXofXZ.cpp.

 {
     int physTot = m_fields[0]->GetTotPoints();
     Vmath::Vcopy(physTot, m_GeometricInfo[0], 1, outarray, 1);
 }

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

◆ v_GetMetricTensor()

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

overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 170 of file MappingXofXZ.cpp.

 {
     int physTot = m_fields[0]->GetTotPoints();
     int nvel    = m_nConvectiveFields;
     Array<OneD, NekDouble> wk(physTot, 0.0);
  
     for (int i = 0; i < nvel * nvel; i++)
     {
         outarray[i] = Array<OneD, NekDouble>(physTot, 0.0);
     }
     // Fill G^{22} and G^{33} with 1.0
     for (int i = 1; i < nvel; i++)
     {
         Vmath::Sadd(physTot, 1.0, outarray[i + nvel * i], 1,
                     outarray[i + nvel * i], 1);
     }
  
     // G_{13} and G_{31} = fz*fx
     Vmath::Vmul(physTot, m_GeometricInfo[1], 1, m_GeometricInfo[0], 1, wk,
                 1); // fz*fx
     Vmath::Vcopy(physTot, wk, 1, outarray[0 * nvel + 2], 1);
     Vmath::Vcopy(physTot, wk, 1, outarray[2 * nvel + 0], 1);
  
     // G^{11} = (fx^2)
     Vmath::Vmul(physTot, m_GeometricInfo[0], 1, m_GeometricInfo[0], 1,
                 outarray[0 * nvel + 0], 1);
  
     // G^{33} = (1+fz^2)
     Vmath::Vmul(physTot, m_GeometricInfo[1], 1, m_GeometricInfo[1], 1, wk,
                 1); // fz^2
     Vmath::Vadd(physTot, wk, 1, outarray[2 * nvel + 2], 1,
                 outarray[2 * nvel + 2], 1);
 }

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

◆ v_InitObject()

void Nektar::GlobalMapping::MappingXofXZ::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

pFields	ExpList array used in the mapping
pMapping	xml element describing the mapping

Reimplemented from Nektar::GlobalMapping::Mapping.

Definition at line 67 of file MappingXofXZ.cpp.

 {
     Mapping::v_InitObject(pFields, pMapping);
  
     m_constantJacobian = false;
  
     ASSERTL0(m_nConvectiveFields == 3,
              "Mapping X = X(x,z) needs 3 velocity components.");
 }

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

◆ v_LowerIndex()

void Nektar::GlobalMapping::MappingXofXZ::v_LowerIndex	(	const Array< OneD, Array< OneD, NekDouble >> &	inarray,
		Array< OneD, Array< OneD, NekDouble >> &	outarray
	)

overrideprotectedvirtual

Reimplemented from Nektar::GlobalMapping::Mapping.

Definition at line 242 of file MappingXofXZ.cpp.

 {
     int physTot = m_fields[0]->GetTotPoints();
     Array<OneD, NekDouble> wk(physTot, 0.0);
  
     // out[0] = in[0]*fx^2 + in[2] * fz*fx
     Vmath::Vmul(physTot, m_GeometricInfo[1], 1, m_GeometricInfo[0], 1, wk,
                 1);                                             // fz*fx
     Vmath::Vmul(physTot, wk, 1, inarray[2], 1, outarray[0], 1); // in[2]*fz*fx
     Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[2], 1); // in[0]*fz*fx
  
     Vmath::Vmul(physTot, m_GeometricInfo[0], 1, m_GeometricInfo[0], 1, wk,
                 1);                                    // fx^2
     Vmath::Vmul(physTot, wk, 1, inarray[0], 1, wk, 1); // in[0]*fx^2
  
     Vmath::Vadd(physTot, outarray[0], 1, wk, 1, outarray[0], 1);
  
     // out[1] = in[1]
     Vmath::Vcopy(physTot, inarray[1], 1, outarray[1], 1);
  
     // out[2] = fx*fz*in[0] + (1+fz^2)*in[2]
     Vmath::Vmul(physTot, m_GeometricInfo[1], 1, m_GeometricInfo[1], 1, wk,
                 1);                                    // fz^2
     Vmath::Sadd(physTot, 1.0, wk, 1, wk, 1);           // 1+fz^2
     Vmath::Vmul(physTot, wk, 1, inarray[2], 1, wk, 1); // (1+fz^2)*in[2]
  
     Vmath::Vadd(physTot, wk, 1, outarray[2], 1, outarray[2], 1);
 }

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

◆ v_RaiseIndex()

void Nektar::GlobalMapping::MappingXofXZ::v_RaiseIndex	(	const Array< OneD, Array< OneD, NekDouble >> &	inarray,
		Array< OneD, Array< OneD, NekDouble >> &	outarray
	)

overrideprotectedvirtual

Reimplemented from Nektar::GlobalMapping::Mapping.

Definition at line 273 of file MappingXofXZ.cpp.

 {
     int physTot = m_fields[0]->GetTotPoints();
     Array<OneD, NekDouble> wk(physTot, 0.0);
     Array<OneD, NekDouble> wk_2(physTot, 0.0);
  
     // out[2] = in[2] - in[0] * fz/fx
     Vmath::Vdiv(physTot, m_GeometricInfo[1], 1, m_GeometricInfo[0], 1, wk, 1);
     Vmath::Vmul(physTot, wk, 1, inarray[0], 1, outarray[2], 1);
     Vmath::Vsub(physTot, inarray[2], 1, outarray[2], 1, outarray[2], 1);
  
     // out[0] = in[0]*(1+fz^2)/(fx^2) - in[2] * fz/fx
     Vmath::Vmul(physTot, wk, 1, inarray[2], 1, outarray[0], 1);
     Vmath::Vmul(physTot, m_GeometricInfo[1], 1, m_GeometricInfo[1], 1, wk, 1);
     Vmath::Sadd(physTot, 1.0, wk, 1, wk, 1);
     Vmath::Vmul(physTot, m_GeometricInfo[0], 1, m_GeometricInfo[0], 1, wk_2, 1);
     Vmath::Vdiv(physTot, wk, 1, wk_2, 1, wk, 1);
     Vmath::Vmul(physTot, wk, 1, inarray[0], 1, wk, 1);
     Vmath::Vsub(physTot, wk, 1, outarray[0], 1, outarray[0], 1);
  
     // out[1] = in[1]
     Vmath::Vcopy(physTot, inarray[1], 1, outarray[1], 1);
 }

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

◆ v_UpdateGeomInfo()

void Nektar::GlobalMapping::MappingXofXZ::v_UpdateGeomInfo ( )

overrideprotectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 363 of file MappingXofXZ.cpp.

 {
     int phystot = m_fields[0]->GetTotPoints();
     // Allocation of geometry memory
     m_GeometricInfo = Array<OneD, Array<OneD, NekDouble>>(5);
     for (int i = 0; i < m_GeometricInfo.size(); i++)
     {
         m_GeometricInfo[i] = Array<OneD, NekDouble>(phystot, 0.0);
     }
  
     bool waveSpace = m_fields[0]->GetWaveSpace();
     m_fields[0]->SetWaveSpace(false);
  
     // Calculate derivatives of transformation
     m_fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[0], m_coords[0],
                            m_GeometricInfo[0]); // f_x
     m_fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[2], m_coords[0],
                            m_GeometricInfo[1]); // f_z
  
     m_fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[0], m_GeometricInfo[0],
                            m_GeometricInfo[2]); // f_xx
     m_fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[2], m_GeometricInfo[0],
                            m_GeometricInfo[3]); // f_xz
     m_fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[2], m_GeometricInfo[1],
                            m_GeometricInfo[4]); // f_zz
  
     m_fields[0]->SetWaveSpace(waveSpace);
 }

References Nektar::MultiRegions::DirCartesianMap, Nektar::GlobalMapping::Mapping::m_coords, Nektar::GlobalMapping::Mapping::m_fields, and Nektar::GlobalMapping::Mapping::m_GeometricInfo.

Friends And Related Function Documentation

◆ MemoryManager< MappingXofXZ >

friend class MemoryManager< MappingXofXZ >

friend

Definition at line 1 of file MappingXofXZ.h.

Member Data Documentation

◆ className

std::string Nektar::GlobalMapping::MappingXofXZ::className

static

Initial value:

=
    GetMappingFactory().RegisterCreatorFunction("XofXZ", MappingXofXZ::create,
                                                "X = X(x,z)")

Name of the class.

Definition at line 70 of file MappingXofXZ.h.

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MappingXofXZ()

Member Function Documentation

◆ create()

◆ v_ApplyChristoffelContravar()

◆ v_ApplyChristoffelCovar()

◆ v_ContravarFromCartesian()

◆ v_ContravarToCartesian()

◆ v_CovarFromCartesian()

◆ v_CovarToCartesian()

◆ v_DotGradJacobian()

◆ v_GetInvMetricTensor()

◆ v_GetJacobian()

◆ v_GetMetricTensor()

◆ v_InitObject()

◆ v_LowerIndex()

◆ v_RaiseIndex()

◆ v_UpdateGeomInfo()

Friends And Related Function Documentation

◆ MemoryManager< MappingXofXZ >

Member Data Documentation

◆ className