47 namespace LibUtilities
67 if(fpoints0 == tpoints0)
97 if(fpoints0 == tpoints0)
111 from, 1, 0.0, to, 1);
134 Interp2D(fpoints0,fpoints1,from.data(),tpoints0,tpoints1,to.data());
145 if((fpoints0 == tpoints0)&&(fpoints1 == tpoints1))
160 if(fpoints1 == tpoints1)
167 Blas::Dgemm(
'N',
'T', fnp0, tnp1, fnp1, 1.0, from, fnp0,
168 I1->GetPtr().get(), tnp1, 0.0, wsp.get(), fnp0);
171 if(fpoints0 == tpoints0)
178 Blas::Dgemm(
'N',
'N', tnp0, tnp1, fnp0, 1.0, I0->GetPtr().get(),
179 tnp0, wsp.get(), fnp0, 0.0, to, tnp0);
211 Interp3D(fpoints0,fpoints1,fpoints2,from.data(),
212 tpoints0,tpoints1,tpoints2,to.data());
241 Blas::Dgemm(
'N',
'N', tnp0, fnp1*fnp2, fnp0, 1.0, I0->GetPtr().get(),
242 tnp0, from, fnp0, 0.0, wsp2.get(), tnp0);
244 for(i = 0; i < fnp2; i++)
246 Blas::Dgemm(
'N',
'T', tnp0, tnp1, fnp1, 1.0, wsp2.get()+i*tnp0*fnp1,
247 tnp0, I1->GetPtr().get(), tnp1, 0.0, wsp1.get()+i*tnp0*tnp1, tnp0);
250 Blas::Dgemm(
'N',
'T', tnp0*tnp1, tnp2, fnp2, 1.0, wsp1.get(),
251 tnp0*tnp1, I2->GetPtr().get(), tnp2, 0.0, to, tnp0*tnp1);
boost::shared_ptr< DNekMat > DNekMatSharedPtr
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...
PointsManagerT & PointsManager(void)
Defines a specification for a set of points.
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...
PointsKey GetPointsKey() const
Return distribution of points.
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 ...
unsigned int GetNumPoints() const
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Describes the specification for a Basis.