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
114 ASSERTL2(Lcoord[0] >= -1,
"Lcoord[0] < -1");
115 ASSERTL2(Lcoord[0] <= 1,
"Lcoord[0] > 1");
117 val =
Blas::Ddot(nquad, I->GetPtr(), 1, physvals, 1);
135 std::map<int, NormalVector>::const_iterator x;
138 "vertex normal not computed.");
#define ASSERTL0(condition, msg)
std::map< int, NormalVector > m_vertexNormals
virtual ~StdExpansion1D()
boost::shared_ptr< DNekMat > DNekMatSharedPtr
int GetTotPoints() const
This function returns the total number of quadrature points used in the element.
const NormalVector & v_GetVertexNormal(const int vertex) const
The base class for all shapes.
void ComputeVertexNormal(const int vertex)
const NormalVector & v_GetSurfaceNormal(const int id) const
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)
void PhysTensorDeriv(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Evaluate the derivative at the physical quadrature points given by inarray and return in outarray...
#define ASSERTL2(condition, msg)
Assert Level 2 – Debugging which is used FULLDEBUG compilation mode. This level assert is designed t...
virtual void v_SetUpPhysNormals(const int vertex)
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Describes the specification for a Basis.
virtual NekDouble v_PhysEvaluate(const Array< OneD, const NekDouble > &coords, const Array< OneD, const NekDouble > &physvals)