90 int nquad =
m_base[0]->GetNumPoints();
92 if(
m_base[0]->Collocation())
99 #ifdef NEKTAR_USING_DIRECT_BLAS_CALLS
101 Blas::Dgemv(
'N',nquad,
m_base[0]->GetNumModes(),1.0, (
m_base[0]->GetBdata()).
get(),
102 nquad,&inarray[0],1,0.0,&outarray[0],1);
104 #else //NEKTAR_USING_DIRECT_BLAS_CALLS
111 #endif //NEKTAR_USING_DIRECT_BLAS_CALLS
120 if(
m_base[0]->Collocation())
143 if(
m_base[0]->Collocation())
166 ASSERTL0(
false,
"This type of FwdTrans is not defined for this shapex type");
169 fill(outarray.get(), outarray.get()+
m_ncoeffs, 0.0 );
190 Blas::Dgemv(
'N',nInteriorDofs,nInteriorDofs,1.0, &(matsys->GetPtr())[0],
191 nInteriorDofs,tmp1.get()+offset,1,0.0,outarray.get()+offset,1);
210 int nquad =
m_base[0]->GetNumPoints();
217 if(coll_check&&
m_base[0]->Collocation())
223 Blas::Dgemv(
'T',nquad,
m_ncoeffs,1.0,base.get(),nquad,
224 &tmp[0],1,0.0,outarray.get(),1);
239 ASSERTL1(dir >= 0 && dir < 1,
"input dir is out of range");
246 bool multiplybyweights)
267 (*Mat) = Imass*Iprod;
virtual LibUtilities::ShapeType v_DetShapeType() const
#define ASSERTL0(condition, msg)
MatrixType GetMatrixType() const
void MassMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
virtual void v_BwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
int GetNumPoints(const int dir) const
This function returns the number of quadrature points in the dir direction.
Principle Modified Functions .
virtual void v_IProductWRTBase(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Calculates the inner product of a given function f with the different modes of the expansion...
virtual void v_GetCoords(Array< OneD, NekDouble > &coords_0, Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2)
LibUtilities::NekManager< StdMatrixKey, DNekBlkMat, StdMatrixKey::opLess > m_stdStaticCondMatrixManager
boost::shared_ptr< DNekMat > DNekMatSharedPtr
DNekMatSharedPtr GetStdMatrix(const StdMatrixKey &mkey)
virtual void v_FwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Transform a given function from physical quadrature space to coefficient space.
virtual void v_BwdTrans_SumFac(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
virtual void v_FwdTrans_BndConstrained(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
The base class for all shapes.
int GetVertexMap(const int localVertexId, bool useCoeffPacking=false)
Principle Modified Functions .
DNekMatSharedPtr CreateGeneralMatrix(const StdMatrixKey &mkey)
this function generates the mass matrix
virtual void v_IProductWRTBase_SumFac(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true)
virtual DNekMatSharedPtr v_CreateStdMatrix(const StdMatrixKey &mkey)
LibUtilities::BasisType GetBasisType(const int dir) const
This function returns the type of basis used in the dir direction.
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.
virtual void v_IProductWRTDerivBase(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
#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 Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Describes the specification for a Basis.
virtual DNekMatSharedPtr v_GenMatrix(const StdMatrixKey &mkey)
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.