#include <MappingXofXZ.h>

Inheritance diagram for Nektar::GlobalMapping::MappingXofXZ:

Collaboration 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)

virtual GLOBAL_MAPPING_EXPORT void	v_ContravarToCartesian (const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray)

virtual GLOBAL_MAPPING_EXPORT void	v_CovarToCartesian (const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray)

virtual GLOBAL_MAPPING_EXPORT void	v_ContravarFromCartesian (const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray)

virtual GLOBAL_MAPPING_EXPORT void	v_CovarFromCartesian (const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray)

virtual GLOBAL_MAPPING_EXPORT void	v_GetJacobian (Array< OneD, NekDouble > &outarray)

virtual GLOBAL_MAPPING_EXPORT void	v_DotGradJacobian (const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, NekDouble > &outarray)

virtual GLOBAL_MAPPING_EXPORT void	v_GetMetricTensor (Array< OneD, Array< OneD, NekDouble > > &outarray)

virtual GLOBAL_MAPPING_EXPORT void	v_GetInvMetricTensor (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_ApplyChristoffelContravar (const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray)

virtual GLOBAL_MAPPING_EXPORT void	v_ApplyChristoffelCovar (const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray)

virtual GLOBAL_MAPPING_EXPORT void	v_UpdateGeomInfo ()

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 (const 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 (const 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 52 of file MappingXofXZ.h.

Constructor & Destructor Documentation

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

protected

Definition at line 58 of file MappingXofXZ.cpp.

     : Mapping(pSession, pFields)
 {
 }

Member Function Documentation

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 60 of file MappingXofXZ.h.

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

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

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

protectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 314 of file MappingXofXZ.cpp.

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

 {
     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); 
 
 }

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

protectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 347 of file MappingXofXZ.cpp.

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

 {
     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);  
 }

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

protectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 121 of file MappingXofXZ.cpp.

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

 {
     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);        
 }

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

protectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 81 of file MappingXofXZ.cpp.

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

 {
     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);
 }

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

protectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 143 of file MappingXofXZ.cpp.

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

 {
     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);
 }

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

protectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 100 of file MappingXofXZ.cpp.

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

 {
     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);
 }

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

protectedvirtual

Reimplemented from Nektar::GlobalMapping::Mapping.

Definition at line 169 of file MappingXofXZ.cpp.

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

 {
     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);
 }

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

protectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 215 of file MappingXofXZ.cpp.

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().

 {
     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);
 }

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

protectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 162 of file MappingXofXZ.cpp.

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

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

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

protectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 180 of file MappingXofXZ.cpp.

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().

 {
     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);
 }

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

protectedvirtual

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 69 of file MappingXofXZ.cpp.

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

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

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

protectedvirtual

Reimplemented from Nektar::GlobalMapping::Mapping.

Definition at line 252 of file MappingXofXZ.cpp.

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

 {
     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); 
 }

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

protectedvirtual

Reimplemented from Nektar::GlobalMapping::Mapping.

Definition at line 283 of file MappingXofXZ.cpp.

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

 {
     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); 
 
 
 }

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

protectedvirtual

Implements Nektar::GlobalMapping::Mapping.

Definition at line 379 of file MappingXofXZ.cpp.

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

 {
     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.num_elements(); 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);
 }

Friends And Related Function Documentation

friend class MemoryManager< MappingXofXZ >

friend

Definition at line 56 of file MappingXofXZ.h.

Member Data Documentation

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 73 of file MappingXofXZ.h.

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation