73 #ifdef NEKTAR_USING_DIRECT_BLAS_CALLS
75 if( inarray.data() == outarray.data())
79 Blas::Dgemv(
'N',nquad,nquad,1.0,&(D->GetPtr())[0],nquad,
80 &wsp[0],1,0.0,&outarray[0],1);
84 Blas::Dgemv(
'N',nquad,nquad,1.0,&(D->GetPtr())[0],nquad,
85 &inarray[0],1,0.0,&outarray[0],1);
88 #else //NEKTAR_USING_DIRECT_BLAS_CALLS
92 if(inarray.data() == outarray.data())
103 #endif //NEKTAR_USING_DIRECT_BLAS_CALLS
112 ASSERTL2(Lcoord[0] >= -1,
"Lcoord[0] < -1");
113 ASSERTL2(Lcoord[0] <= 1,
"Lcoord[0] > 1");
115 val =
Blas::Ddot(nquad, I->GetPtr(), 1, physvals, 1);
boost::shared_ptr< DNekMat > DNekMatSharedPtr
void PhysTensorDeriv(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
int GetTotPoints() const
This function returns the total number of quadrature points used in the element.
virtual NekDouble v_PhysEvaluate(const Array< OneD, const NekDouble > &coords, const Array< OneD, const NekDouble > &physvals)
The base class for all shapes.
virtual ~StdExpansion0D()
void CopyArray(const Array< OneD, ConstDataType > &source, Array< OneD, DataType > &dest)
T Ddot(int n, const Array< OneD, const T > &w, const int incw, const Array< OneD, const T > &x, const int incx, const Array< OneD, const int > &y, const int incy)
#define ASSERTL2(condition, msg)
Assert Level 2 – Debugging which is used FULLDEBUG compilation mode. This level assert is designed t...
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Describes the specification for a Basis.