46 namespace LibUtilities
66 if(fpoints0 == tpoints0)
74 GP0 =
PointsManager()[tpoints0]->GetGalerkinProjection(fpoints0);
97 if(fpoints0 == tpoints0)
107 ->GetGalerkinProjection(fpoints0);
111 from, 1, 0.0, to, 1);
152 if(fpoints1 == tpoints1)
158 GP1 =
PointsManager()[tpoints1]->GetGalerkinProjection(fpoints1);
159 Blas::Dgemm(
'N',
'T', fnp0, tnp1, fnp1, 1.0, from, fnp0,
160 GP1->GetPtr().get(), tnp1, 0.0, wsp.get(), fnp0);
163 if(fpoints0 == tpoints0)
169 GP0 =
PointsManager()[tpoints0]->GetGalerkinProjection(fpoints0);
170 Blas::Dgemm(
'N',
'N', tnp0, tnp1, fnp0, 1.0,
172 tnp0, wsp.get(), fnp0, 0.0, to, tnp0);
207 tpoints0,tpoints1,tpoints2,to.data());
231 GP2 =
PointsManager()[tpoints2]->GetGalerkinProjection(fpoints2);
232 Blas::Dgemm(
'N',
'T', fnp0*fnp1, tnp2, fnp2, 1.0, from, fnp0*fnp1,
233 GP2->GetPtr().get(), tnp2, 0.0, wsp2.get(), fnp0*fnp1);
235 GP1 =
PointsManager()[tpoints1]->GetGalerkinProjection(fpoints1);
236 for(
int i = 0; i < tnp2; i++)
238 Blas::Dgemm(
'N',
'T', fnp0, tnp1, fnp1, 1.0,
239 wsp2.get()+i*fnp0*fnp1,
240 fnp0, GP1->GetPtr().get(),tnp1, 0.0,
241 wsp1.get()+i*fnp0*tnp1,
245 GP0 =
PointsManager()[tpoints0]->GetGalerkinProjection(fpoints0);
246 Blas::Dgemm(
'N',
'N', tnp0, tnp1*tnp2, fnp0, 1.0,
247 GP0->GetPtr().get(), tnp0, wsp1.get(), fnp0, 0.0,
void PhysGalerkinProject1D(const BasisKey &fbasis0, const Array< OneD, const NekDouble > &from, const BasisKey &tbasis0, Array< OneD, NekDouble > &to)
boost::shared_ptr< DNekMat > DNekMatSharedPtr
PointsManagerT & PointsManager(void)
void PhysGalerkinProject3D(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)
Defines a specification for a set of points.
PointsKey GetPointsKey() const
Return distribution of points.
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.
void PhysGalerkinProject2D(const BasisKey &fbasis0, const BasisKey &fbasis1, const Array< OneD, const NekDouble > &from, const BasisKey &tbasis0, const BasisKey &tbasis1, Array< OneD, NekDouble > &to)