41 namespace GlobalMapping
71 const TiXmlElement *pMapping)
76 if (
m_session->DefinesElement(
"Nektar/Mapping"))
82 int phystot = pFields[0]->GetTotPoints();
88 for (
int i = 0; i < 3; i++)
110 int physTot =
m_fields[0]->GetTotPoints();
122 int physTot =
m_fields[0]->GetTotPoints();
134 int physTot =
m_fields[0]->GetTotPoints();
146 int physTot =
m_fields[0]->GetTotPoints();
157 int physTot =
m_fields[0]->GetTotPoints();
165 int physTot =
m_fields[0]->GetTotPoints();
173 int physTot =
m_fields[0]->GetTotPoints();
176 for (
int i=0; i<nvel*nvel; i++)
181 for (
int i=0; i<nvel; i++)
184 outarray[i+nvel*i], 1);
191 int physTot =
m_fields[0]->GetTotPoints();
194 for (
int i=0; i<nvel*nvel; i++)
199 for (
int i=0; i<nvel; i++)
202 outarray[i+nvel*i], 1);
210 int physTot =
m_fields[0]->GetTotPoints();
213 for (
int i=0; i<nvel; i++)
218 for (
int i=0; i<nvel; i++)
228 int physTot =
m_fields[0]->GetTotPoints();
231 for (
int i=0; i<nvel; i++)
236 for (
int i=0; i<nvel; i++)
246 int physTot =
m_fields[0]->GetTotPoints();
249 for (
int i = 0; i< nvel; i++)
251 for (
int j = 0; j< nvel; j++)
262 int physTot =
m_fields[0]->GetTotPoints();
265 for (
int i = 0; i< nvel; i++)
267 for (
int j = 0; j< nvel; j++)
Array< OneD, Array< OneD, NekDouble > > m_coords
Array with the Cartesian coordinates.
virtual GLOBAL_MAPPING_EXPORT void v_RaiseIndex(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray)
MappingFactory & GetMappingFactory()
Declaration of the mapping factory singleton.
virtual GLOBAL_MAPPING_EXPORT void v_ApplyChristoffelContravar(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray)
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
virtual GLOBAL_MAPPING_EXPORT void v_DotGradJacobian(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, NekDouble > &outarray)
virtual GLOBAL_MAPPING_EXPORT void v_CovarFromCartesian(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray)
MappingTranslation(const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields)
boost::shared_ptr< SessionReader > SessionReaderSharedPtr
int m_nConvectiveFields
Number of velocity components.
virtual GLOBAL_MAPPING_EXPORT void v_CovarToCartesian(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray)
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
LibUtilities::SessionReaderSharedPtr m_session
Session reader.
virtual GLOBAL_MAPPING_EXPORT void v_GetInvMetricTensor(Array< OneD, Array< OneD, NekDouble > > &outarray)
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.
bool m_timeDependent
Flag defining if the Mapping is time-dependent.
virtual GLOBAL_MAPPING_EXPORT void v_UpdateGeomInfo()
Array< OneD, Array< OneD, NekDouble > > m_coordsVel
Array with the velocity of the coordinates.
virtual GLOBAL_MAPPING_EXPORT void v_InitObject(const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const TiXmlElement *pMapping)
bool m_constantJacobian
Flag defining if the Jacobian is constant.
virtual GLOBAL_MAPPING_EXPORT void v_GetJacobian(Array< OneD, NekDouble > &outarray)
Array< OneD, Array< OneD, NekDouble > > m_tmp
void Sadd(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Add vector y = alpha + x.
virtual GLOBAL_MAPPING_EXPORT void v_ContravarToCartesian(const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray)
Base class for mapping to be applied to the coordinate system.
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_LowerIndex(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)
void Zero(int n, T *x, const int incx)
Zero vector.
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
virtual GLOBAL_MAPPING_EXPORT void v_GetMetricTensor(Array< OneD, Array< OneD, NekDouble > > &outarray)
virtual GLOBAL_MAPPING_EXPORT void v_InitObject(const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const TiXmlElement *pMapping)
static std::string className
Name of the class.
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, tDescription pDesc="")
Register a class with the factory.