46 namespace LocalRegions
50 m_metricinfo(m_geom->GetGeomFactors())
59 int nDim =
m_base.num_elements();
60 string type =
"regular";
67 err << nDim <<
"D " << type <<
" Jacobian not positive "
68 <<
"(element ID = " <<
m_geom->GetGlobalID() <<
") "
69 <<
"(first vertex ID = " <<
m_geom->GetVid(0) <<
")";
76 m_metricinfo(pSrc.m_metricinfo)
107 const std::vector<unsigned int > &nummodes,
108 const int nmodes_offset,
110 std::vector<LibUtilities::BasisType> &fromType)
117 const boost::shared_ptr<Expansion> &EdgeExp,
127 const boost::shared_ptr<Expansion> &EdgeExp,
136 const boost::shared_ptr<Expansion> &FaceExp,
150 v_DGDeriv(dir, inarray, EdgeExp, coeffs, outarray);
234 const int expDim =
m_base.num_elements();
241 for (
int i = 0; i < expDim; ++i)
243 CBasis[i] =
m_geom->GetBasis(i);
244 nqGeom *= CBasis[i]->GetNumPoints();
245 doCopy = doCopy &&
m_base[i]->GetBasisKey().SamePoints(
246 CBasis[i]->GetBasisKey());
255 for (
int i = 0; i <
m_geom->GetCoordim(); ++i)
257 m_geom->GetXmap()->BwdTrans(
m_geom->GetCoeffs(i), tmp[i]);
262 for (
int i = 0; i <
m_geom->GetCoordim(); ++i)
265 m_geom->GetXmap()->BwdTrans(
m_geom->GetCoeffs(i), tmpGeom);
272 CBasis[0]->GetPointsKey(), &tmpGeom[0],
273 m_base[0]->GetPointsKey(), &tmp[i][0]);
279 CBasis[0]->GetPointsKey(),
280 CBasis[1]->GetPointsKey(),
282 m_base[0]->GetPointsKey(),
283 m_base[1]->GetPointsKey(),
290 CBasis[0]->GetPointsKey(),
291 CBasis[1]->GetPointsKey(),
292 CBasis[2]->GetPointsKey(),
294 m_base[0]->GetPointsKey(),
295 m_base[1]->GetPointsKey(),
296 m_base[2]->GetPointsKey(),
322 const std::vector<unsigned int > &nummodes,
323 const int nmodes_offset,
325 std::vector<LibUtilities::BasisType> &fromType)
332 const boost::shared_ptr<Expansion> &EdgeExp,
342 const boost::shared_ptr<Expansion> &EdgeExp,
351 const boost::shared_ptr<Expansion> &FaceExp,
void ComputeLaplacianMetric()
const LibUtilities::PointsKeyVector GetPointsKeys() const
LibUtilities::ShapeType DetShapeType() const
This function returns the shape of the expansion domain.
#define NEKERROR(type, msg)
Assert Level 0 – Fundamental assert which is used whether in FULLDEBUG, DEBUG or OPT compilation mod...
std::vector< PointsKey > PointsKeyVector
DNekMatSharedPtr BuildTransformationMatrix(const DNekScalMatSharedPtr &r_bnd, const StdRegions::MatrixType matrixType)
void ExtractDataToCoeffs(const NekDouble *data, const std::vector< unsigned int > &nummodes, const int nmodes_offset, NekDouble *coeffs, std::vector< LibUtilities::BasisType > &fromType)
void MultiplyByStdQuadratureMetric(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
virtual void v_MultiplyByQuadratureMetric(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
static DNekScalMatSharedPtr NullDNekScalMatSharedPtr
DNekMatSharedPtr BuildVertexMatrix(const DNekScalMatSharedPtr &r_bnd)
std::map< ConstFactorType, NekDouble > ConstFactorMap
SpatialDomains::GeometrySharedPtr m_geom
boost::shared_ptr< DNekMat > DNekMatSharedPtr
boost::shared_ptr< DNekScalMat > DNekScalMatSharedPtr
virtual DNekMatSharedPtr v_BuildTransformationMatrix(const DNekScalMatSharedPtr &r_bnd, const StdRegions::MatrixType matrixType)
virtual DNekMatSharedPtr v_BuildVertexMatrix(const DNekScalMatSharedPtr &r_bnd)
virtual void v_ExtractDataToCoeffs(const NekDouble *data, const std::vector< unsigned int > &nummodes, const int nmodes_offset, NekDouble *coeffs, std::vector< LibUtilities::BasisType > &fromType)
static DNekMatSharedPtr NullDNekMatSharedPtr
int GetTotPoints() const
This function returns the total number of quadrature points used in the element.
void Interp2D(const BasisKey &fbasis0, const BasisKey &fbasis1, const Array< OneD, const NekDouble > &from, const BasisKey &tbasis0, const BasisKey &tbasis1, Array< OneD, NekDouble > &to)
this function interpolates a 2D function evaluated at the quadrature points of the 2D basis...
void AddFaceNormBoundaryInt(const int face, const boost::shared_ptr< Expansion > &FaceExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)
Expansion(SpatialDomains::GeometrySharedPtr pGeom)
virtual void v_ComputeLaplacianMetric()
std::map< StdRegions::VarCoeffType, Array< OneD, NekDouble > > VarCoeffMap
virtual void v_GetCoords(Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2, Array< OneD, NekDouble > &coords_3)
virtual DNekScalMatSharedPtr v_GetLocMatrix(const LocalRegions::MatrixKey &mkey)
void DGDeriv(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, Array< OneD, NekDouble > > &coeffs, Array< OneD, NekDouble > &outarray)
virtual const SpatialDomains::GeomFactorsSharedPtr & v_GetMetricInfo() const
virtual void v_AddFaceNormBoundaryInt(const int face, const boost::shared_ptr< Expansion > &FaceExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)
DNekScalMatSharedPtr GetLocMatrix(const LocalRegions::MatrixKey &mkey)
void Interp3D(const BasisKey &fbasis0, const BasisKey &fbasis1, const BasisKey &fbasis2, const Array< OneD, const NekDouble > &from, const BasisKey &tbasis0, const BasisKey &tbasis1, const BasisKey &tbasis2, Array< OneD, NekDouble > &to)
this function interpolates a 3D function evaluated at the quadrature points of the 3D basis...
SpatialDomains::GeometrySharedPtr GetGeom() const
boost::shared_ptr< GeomFactors > GeomFactorsSharedPtr
Pointer to a GeomFactors object.
virtual void v_DGDeriv(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, Array< OneD, NekDouble > > &coeffs, Array< OneD, NekDouble > &outarray)
Geometry is straight-sided with constant geometric factors.
void Interp1D(const BasisKey &fbasis0, const Array< OneD, const NekDouble > &from, const BasisKey &tbasis0, Array< OneD, NekDouble > &to)
this function interpolates a 1D function evaluated at the quadrature points of the basis fbasis0 to ...
void ComputeQuadratureMetric()
GeomType
Indicates the type of element geometry.
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Array< OneD, LibUtilities::BasisSharedPtr > m_base
void AddEdgeNormBoundaryInt(const int edge, const boost::shared_ptr< Expansion > &EdgeExp, const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray)
Geometry is curved or has non-constant factors.
boost::shared_ptr< Geometry > GeometrySharedPtr
virtual void v_AddEdgeNormBoundaryInt(const int edge, const boost::shared_ptr< Expansion > &EdgeExp, const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray)
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.