48 StdQuadExp::StdQuadExp()
57 StdExpansion (Ba.GetNumModes()*Bb.GetNumModes(),2,Ba,Bb),
123 ASSERTL1(
false,
"input dir is out of range");
155 if(
m_base[0]->Collocation() &&
m_base[1]->Collocation())
158 inarray, 1, outarray, 1);
171 m_base[1]->GetNumModes());
175 inarray,outarray,wsp,
true,
true);
194 bool doCheckCollDir0,
195 bool doCheckCollDir1)
197 int nquad0 =
m_base[0]->GetNumPoints();
198 int nquad1 =
m_base[1]->GetNumPoints();
199 int nmodes0 =
m_base[0]->GetNumModes();
200 int nmodes1 =
m_base[1]->GetNumModes();
202 bool colldir0 = doCheckCollDir0?(
m_base[0]->Collocation()):
false;
203 bool colldir1 = doCheckCollDir1?(
m_base[1]->Collocation()):
false;
205 if(colldir0 && colldir1)
211 Blas::Dgemm(
'N',
'T', nquad0, nquad1,nmodes1, 1.0, &inarray[0], nquad0,
212 base1.get(), nquad1, 0.0, &outarray[0], nquad0);
216 Blas::Dgemm(
'N',
'N', nquad0,nmodes1,nmodes0,1.0, base0.get(),
217 nquad0, &inarray[0], nmodes0,0.0,&outarray[0], nquad0);
221 ASSERTL1(wsp.num_elements()>=nquad0*nmodes1,
"Workspace size is not sufficient");
225 Blas::Dgemm(
'N',
'N', nquad0,nmodes1,nmodes0,1.0, base0.get(),
226 nquad0, &inarray[0], nmodes0,0.0,&wsp[0], nquad0);
227 Blas::Dgemm(
'N',
'T', nquad0, nquad1,nmodes1, 1.0, &wsp[0], nquad0,
228 base1.get(), nquad1, 0.0, &outarray[0], nquad0);
236 if((
m_base[0]->Collocation())&&(
m_base[1]->Collocation()))
260 if((
m_base[0]->Collocation())&&(
m_base[1]->Collocation()))
267 int npoints[2] = {
m_base[0]->GetNumPoints(),
268 m_base[1]->GetNumPoints()};
269 int nmodes[2] = {
m_base[0]->GetNumModes(),
270 m_base[1]->GetNumModes()};
272 fill(outarray.get(), outarray.get()+
m_ncoeffs, 0.0 );
276 for(i = 0; i < 4; i++)
282 for(i = 0; i < npoints[0]; i++)
284 physEdge[0][i] = inarray[i];
285 physEdge[2][i] = inarray[npoints[0]*npoints[1]-1-i];
288 for(i = 0; i < npoints[1]; i++)
290 physEdge[1][i] = inarray[npoints[0]-1+i*npoints[0]];
291 physEdge[3][i] = inarray[(npoints[1]-1)*npoints[0]-i*npoints[0]];
301 for(i = 0; i < 4; i++)
303 segexp[i%2]->FwdTrans_BndConstrained(physEdge[i],coeffEdge[i]);
306 for(j=0; j < nmodes[i%2]; j++)
309 outarray[ mapArray[j] ] = sign * coeffEdge[i][j];
328 int nInteriorDofs =
m_ncoeffs - nBoundaryDofs;
336 for(i = 0; i < nInteriorDofs; i++)
338 rhs[i] = tmp1[ mapArray[i] ];
341 Blas::Dgemv(
'N',nInteriorDofs,nInteriorDofs,1.0, &(matsys->GetPtr())[0],
342 nInteriorDofs,rhs.get(),1,0.0,result.get(),1);
344 for(i = 0; i < nInteriorDofs; i++)
346 outarray[ mapArray[i] ] = result[i];
385 if(
m_base[0]->Collocation() &&
m_base[1]->Collocation())
398 bool multiplybyweights)
400 int nquad0 =
m_base[0]->GetNumPoints();
401 int nquad1 =
m_base[1]->GetNumPoints();
402 int order0 =
m_base[0]->GetNumModes();
404 if(multiplybyweights)
413 tmp,outarray,wsp,
true,
true);
420 inarray,outarray,wsp,
true,
true);
432 Blas::Dgemv(
'N',
m_ncoeffs,nq,1.0,iprodmat->GetPtr().get(),
433 m_ncoeffs, inarray.get(), 1, 0.0, outarray.get(), 1);
447 ASSERTL0((dir==0)||(dir==1),
"input dir is out of range");
449 int nquad0 =
m_base[0]->GetNumPoints();
450 int nquad1 =
m_base[1]->GetNumPoints();
451 int nqtot = nquad0*nquad1;
452 int order0 =
m_base[0]->GetNumModes();
464 tmp,outarray,wsp,
true,
false);
470 tmp,outarray,wsp,
false,
true);
478 ASSERTL0((dir==0)||(dir==1),
"input dir is out of range");
495 Blas::Dgemv(
'N',
m_ncoeffs,nq,1.0,iprodmat->GetPtr().get(),
496 m_ncoeffs, inarray.get(), 1, 0.0, outarray.get(), 1);
515 bool doCheckCollDir0,
516 bool doCheckCollDir1)
518 int nquad0 =
m_base[0]->GetNumPoints();
519 int nquad1 =
m_base[1]->GetNumPoints();
520 int nmodes0 =
m_base[0]->GetNumModes();
521 int nmodes1 =
m_base[1]->GetNumModes();
523 bool colldir0 = doCheckCollDir0?(
m_base[0]->Collocation()):
false;
524 bool colldir1 = doCheckCollDir1?(
m_base[1]->Collocation()):
false;
526 if(colldir0 && colldir1)
532 Blas::Dgemm(
'N',
'N',nmodes0,nmodes1, nquad1,1.0, inarray.get(),
533 nmodes0, base1.get(), nquad1, 0.0,outarray.get(),nmodes0);
537 Blas::Dgemm(
'T',
'N',nmodes0,nquad1,nquad0,1.0,base0.get(),
538 nquad0,inarray.get(),nquad0,0.0,outarray.get(),nmodes0);
542 ASSERTL1(wsp.num_elements()>=nquad1*nmodes0,
"Workspace size is not sufficient");
545 Blas::Dgemm(
'T',
'N',nmodes0,nquad1,nquad0,1.0,base0.get(),
546 nquad0,inarray.get(),nquad0,0.0,wsp.get(),nmodes0);
549 for(
int i = 0; i < nmodes0; ++i)
551 for(
int j = 0; j < nquad1; ++j)
554 base0.get()+i*nquad0,1,
555 inarray.get()+j*nquad0,1);
559 Blas::Dgemm(
'N',
'N',nmodes0,nmodes1, nquad1,1.0, wsp.get(),
560 nmodes0, base1.get(), nquad1, 0.0,outarray.get(),nmodes0);
586 int nquad0 =
m_base[0]->GetNumPoints();
587 int nquad1 =
m_base[1]->GetNumPoints();
590 int btmp0 =
m_base[0]->GetNumModes();
591 int mode0 = mode%btmp0;
592 int mode1 = mode/btmp0;
595 ASSERTL2(mode1 == (
int)floor((1.0*mode)/btmp0),
596 "Integer Truncation not Equiv to Floor");
599 "calling argument mode is larger than total expansion order");
601 for(i = 0; i < nquad1; ++i)
604 1, &outarray[0]+i*nquad0,1);
607 for(i = 0; i < nquad0; ++i)
610 &outarray[0]+i,nquad0,&outarray[0]+i,nquad0);
630 ASSERTL2((i >= 0)&&(i <= 3),
"edge id is out of range");
632 if((i == 0)||(i == 2))
644 ASSERTL2((i >= 0)&&(i <= 3),
"edge id is out of range");
646 if((i == 0)||(i == 2))
658 ASSERTL2((i >= 0)&&(i <= 3),
"edge id is out of range");
660 if((i == 0)||(i == 2))
672 ASSERTL2((edge >= 0)&&(edge <= 3),
"edge id is out of range");
674 if((edge == 0)||(edge == 2))
686 ASSERTL2((i >= 0)&&(i <= 3),
"edge id is out of range");
688 if((i == 0)||(i == 2))
710 "BasisType is not a boundary interior form");
714 "BasisType is not a boundary interior form");
724 "BasisType is not a boundary interior form");
728 "BasisType is not a boundary interior form");
734 const std::vector<unsigned int> &nummodes,
737 int nmodes = nummodes[modes_offset]*nummodes[modes_offset+1];
745 bool returnval =
false;
770 for(i = 0; i < nq1; ++i)
772 Blas::Dcopy(nq0,z0.get(), 1,&coords_0[0] + i*nq0,1);
785 int nummodes0, nummodes1;
786 int value1 = 0, value2 = 0;
792 nummodes0 =
m_base[0]->GetNumModes();
793 nummodes1 =
m_base[1]->GetNumModes();
805 value1 = 2*nummodes0;
808 ASSERTL0(0,
"Mapping array is not defined for this expansion");
812 for(i = 0; i < value1; i++)
822 value2 = value1+nummodes0-1;
828 ASSERTL0(0,
"Mapping array is not defined for this expansion");
832 for(i = 0; i < nummodes1-2; i++)
834 outarray[cnt++]=value1+i*nummodes0;
835 outarray[cnt++]=value2+i*nummodes0;
841 for(i = nummodes0*(nummodes1-1);i <
GetNcoeffs(); i++)
852 int nummodes0, nummodes1;
859 nummodes0 =
m_base[0]->GetNumModes();
860 nummodes1 =
m_base[1]->GetNumModes();
868 startvalue = nummodes0;
871 startvalue = 2*nummodes0;
874 ASSERTL0(0,
"Mapping array is not defined for this expansion");
887 ASSERTL0(0,
"Mapping array is not defined for this expansion");
891 for(i = 0; i < nummodes1-2; i++)
893 for(j = 0; j < nummodes0-2; j++)
895 outarray[cnt++]=startvalue+j;
897 startvalue+=nummodes0;
905 if(useCoeffPacking ==
true)
907 switch(localVertexId)
918 localDOF =
m_base[0]->GetNumModes()-1;
930 localDOF =
m_base[0]->GetNumModes() * (
m_base[1]->GetNumModes()-1);
934 localDOF =
m_base[0]->GetNumModes();
942 localDOF =
m_base[0]->GetNumModes()*
m_base[1]->GetNumModes()-1;
946 localDOF =
m_base[0]->GetNumModes()+1;
951 ASSERTL0(
false,
"eid must be between 0 and 3");
958 switch(localVertexId)
969 localDOF =
m_base[0]->GetNumModes()-1;
981 localDOF =
m_base[0]->GetNumModes()*
m_base[1]->GetNumModes()-1;
985 localDOF =
m_base[0]->GetNumModes()+1;
993 localDOF =
m_base[0]->GetNumModes() * (
m_base[1]->GetNumModes()-1);
997 localDOF =
m_base[0]->GetNumModes();
1002 ASSERTL0(
false,
"eid must be between 0 and 3");
1015 const int nummodes0 =
m_base[0]->GetNumModes();
1016 const int nummodes1 =
m_base[1]->GetNumModes();
1020 if(maparray.num_elements() != nEdgeIntCoeffs)
1025 if(signarray.num_elements() != nEdgeIntCoeffs)
1031 fill( signarray.get() , signarray.get()+nEdgeIntCoeffs, 1 );
1040 for(i = 0; i < nEdgeIntCoeffs; i++)
1047 for(i = 1; i < nEdgeIntCoeffs; i+=2)
1056 for(i = 0; i < nEdgeIntCoeffs; i++)
1058 maparray[i] = (i+2)*nummodes0 + 1;
1063 for(i = 1; i < nEdgeIntCoeffs; i+=2)
1072 for(i = 0; i < nEdgeIntCoeffs; i++)
1074 maparray[i] = nummodes0+i+2;
1079 for(i = 1; i < nEdgeIntCoeffs; i+=2)
1088 for(i = 0; i < nEdgeIntCoeffs; i++)
1090 maparray[i] = (i+2)*nummodes0;
1095 for(i = 1; i < nEdgeIntCoeffs; i+=2)
1103 ASSERTL0(
false,
"eid must be between 0 and 3");
1113 for(i = 0; i < nEdgeIntCoeffs; i++)
1121 for(i = 0; i < nEdgeIntCoeffs; i++)
1123 maparray[i] = (i+2)*nummodes0 - 1;
1129 for(i = 0; i < nEdgeIntCoeffs; i++)
1131 maparray[i] = nummodes0*nummodes1 - 2 - i;
1137 for(i = 0; i < nEdgeIntCoeffs; i++)
1139 maparray[i] = nummodes0*(nummodes1-2-i);
1144 ASSERTL0(
false,
"eid must be between 0 and 3");
1149 reverse( maparray.get() , maparray.get()+nEdgeIntCoeffs );
1154 ASSERTL0(
false,
"Mapping not defined for this type of basis");
1168 int order0 =
m_base[0]->GetNumModes();
1169 int order1 =
m_base[1]->GetNumModes();
1182 bool checkForZeroedModes =
false;
1187 else if(P != numModes)
1189 checkForZeroedModes =
true;
1194 if (maparray.num_elements() != P)
1199 if(signarray.num_elements() != P)
1205 fill(signarray.get(), signarray.get()+P, 1);
1214 for (i = 0; i < P; i++)
1221 swap(maparray[0], maparray[1]);
1223 for(i = 3; i < P; i+=2)
1232 for (i = 0; i < P; i++)
1234 maparray[i] = i*order0 + 1;
1239 swap(maparray[0], maparray[1]);
1241 for(i = 3; i < P; i+=2)
1250 for (i = 0; i < P; i++)
1252 maparray[i] = order0+i;
1257 swap(maparray[0], maparray[1]);
1259 for (i = 3; i < P; i+=2)
1268 for (i = 0; i < P; i++)
1270 maparray[i] = i*order0;
1275 swap(maparray[0], maparray[1]);
1277 for (i = 3; i < P; i+=2)
1285 ASSERTL0(
false,
"eid must be between 0 and 3");
1296 for (i = 0; i < P; i++)
1304 for (i = 0; i < P; i++)
1306 maparray[i] = (i+1)*order0 - 1;
1312 for (i = 0; i < P; i++)
1314 maparray[i] = order0*order1 - 1 - i;
1320 for (i = 0; i < P; i++)
1322 maparray[i] = order0*(order1-1-i);
1327 ASSERTL0(
false,
"eid must be between 0 and 3");
1332 reverse(maparray.get(), maparray.get()+P);
1337 ASSERTL0(
false,
"Mapping not defined for this type of basis");
1340 if (checkForZeroedModes)
1346 for (
int j = numModes; j < P; j++)
1349 maparray[j] = maparray[0];
1354 ASSERTL0(
false,
"Different trace space edge dimension "
1355 "and element edge dimension not possible "
1356 "for GLL-Lagrange bases");
1380 int nq0 =
m_base[0]->GetNumPoints();
1381 int nq1 =
m_base[1]->GetNumPoints();
1405 for(
int i = 0; i < nq; ++i)
1407 for(
int j = 0; j < nq; ++j,++cnt)
1410 coords[cnt][0] = -1.0 + 2*j/(
NekDouble)(nq-1);
1411 coords[cnt][1] = -1.0 + 2*i/(
NekDouble)(nq-1);
1415 for(
int i = 0; i < neq; ++i)
1419 I[0] =
m_base[0]->GetI(coll);
1420 I[1] =
m_base[1]->GetI(coll+1);
1423 for (
int j = 0; j < nq1; ++j)
1429 Mat->GetRawPtr()+j*nq0*neq+i,neq);
1442 for(i = 0; i < order1; ++i)
1444 (*Mat)(order0*i+1,i*order0+1) = 1.0;
1450 for(i = 0; i < order0; ++i)
1452 (*Mat)(order0+i ,order0+i) = 1.0;
1465 (*Mat) = Imass*Iprod;
1473 int dir = (edge + 1) % 2;
1474 int nCoeffs =
m_base[dir]->GetNumModes();
1482 coords[0] = (edge == 0 || edge == 3) ? -1.0 : 1.0;
1490 Vmath::Vcopy(nCoeffs, m_Ix->GetPtr(), 1, Mat->GetPtr(), 1);
1520 if(inarray.get() == outarray.get())
1526 m_ncoeffs, tmp.get(), 1, 0.0, outarray.get(), 1);
1531 m_ncoeffs, inarray.get(), 1, 0.0, outarray.get(), 1);
1540 int qa =
m_base[0]->GetNumPoints();
1541 int qb =
m_base[1]->GetNumPoints();
1542 int nmodes_a =
m_base[0]->GetNumModes();
1543 int nmodes_b =
m_base[1]->GetNumModes();
1544 int nmodes = min(nmodes_a,nmodes_b);
1573 OrthoExp.
FwdTrans(tmp,orthocoeffs);
1575 for(
int j = 0; j < nmodes_a; ++j)
1577 for(
int k = 0; k < nmodes_b; ++k)
1580 orthocoeffs[j*nmodes_b+k] *=
1582 pow(j/(nmodes_a-1)+k/(nmodes_b-1),0.5*nmodes));
1592 OrthoExp.
BwdTrans(orthocoeffs,tmp);
1602 min(nmodes_a,nmodes_b));
1605 OrthoExp.
FwdTrans(array,orthocoeffs);
1611 for(j = 0; j < nmodes_a; ++j)
1613 for(k = 0; k < nmodes_b; ++k)
1617 orthocoeffs[j*nmodes_b+k] *=
1618 (SvvDiffCoeff*exp(-(j+k-nmodes)*(j+k-nmodes)/
1624 orthocoeffs[j*nmodes_b+k] *= 0.0;
1631 OrthoExp.
BwdTrans(orthocoeffs,array);
1640 int n_coeffs = inarray.num_elements();
1648 int nmodes0 =
m_base[0]->GetNumModes();
1649 int nmodes1 =
m_base[1]->GetNumModes();
1650 int numMax = nmodes0;
1666 b0, b1, coeff_tmp, bortho0, bortho1, coeff);
1671 for (
int i = 0; i < numMin+1; ++i)
1675 tmp2 = coeff_tmp+cnt,1);
1681 bortho0, bortho1, coeff_tmp, b0, b1, outarray);
1730 int nquad0 =
m_base[0]->GetNumPoints();
1731 int nquad1 =
m_base[1]->GetNumPoints();
1737 for(i = 0; i < nquad1; ++i)
1740 w0.get(),1,outarray.get()+i*nquad0,1);
1743 for(i = 0; i < nquad0; ++i)
1745 Vmath::Vmul(nquad1,outarray.get()+i,nquad0,w1.get(),1,
1746 outarray.get()+i,nquad0);
1754 int np1 =
m_base[0]->GetNumPoints();
1755 int np2 =
m_base[1]->GetNumPoints();
1756 int np = max(np1,np2);
1763 for(
int i = 0; i < np-1; ++i)
1766 for(
int j = 0; j < np-1; ++j)
1768 conn[cnt++] = row +j;
1769 conn[cnt++] = row +j+1;
1770 conn[cnt++] = rowp1 +j;
1772 conn[cnt++] = rowp1 +j+1;
1773 conn[cnt++] = rowp1 +j;
1774 conn[cnt++] = row +j+1;
virtual DNekMatSharedPtr v_GenMatrix(const StdMatrixKey &mkey)
LibUtilities::ShapeType DetShapeType() const
This function returns the shape of the expansion domain.
NekDouble GetConstFactor(const ConstFactorType &factor) const
DNekMatSharedPtr GenMatrix(const StdMatrixKey &mkey)
#define ASSERTL0(condition, msg)
const LibUtilities::BasisSharedPtr & GetBasis(int dir) const
This function gets the shared point to basis in the dir direction.
const ConstFactorMap & GetConstFactors() const
virtual void v_IProductWRTDerivBase_SumFac(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
MatrixType GetMatrixType() const
void MassMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
static Array< OneD, NekDouble > NullNekDouble1DArray
void BwdTrans_SumFacKernel(const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0=true, bool doCheckCollDir1=true)
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
int GetBasisNumModes(const int dir) const
This function returns the number of expansion modes in the dir direction.
bool HasVarCoeff(const StdRegions::VarCoeffType &coeff) const
#define sign(a, b)
return the sign(b)*a
virtual void v_PhysDeriv(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)
Calculate the derivative of the physical points.
void IProductWRTDerivBase_SumFac(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
void Vsqrt(int n, const T *x, const int incx, T *y, const int incy)
sqrt y = sqrt(x)
void MultiplyByQuadratureMetric(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
virtual DNekMatSharedPtr v_CreateStdMatrix(const StdMatrixKey &mkey)
void IProductWRTBase(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
this function calculates the inner product of a given function f with the different modes of the expa...
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
virtual void v_ReduceOrderCoeffs(int numMin, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
void LaplacianMatrixOp_MatFree(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
const Array< OneD, const NekDouble > & GetVarCoeff(const StdRegions::VarCoeffType &coeff) const
int GetNumPoints(const int dir) const
This function returns the number of quadrature points in the dir direction.
Principle Modified Functions .
virtual bool v_IsBoundaryInteriorExpansion()
Lagrange Polynomials using the Gauss points .
void LocCoordToLocCollapsed(const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta)
Convert local cartesian coordinate xi into local collapsed coordinates eta.
std::map< ConstFactorType, NekDouble > ConstFactorMap
NekDouble Integral(const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &w0, const Array< OneD, const NekDouble > &w1)
virtual void v_GetInteriorMap(Array< OneD, unsigned int > &outarray)
virtual void v_SVVLaplacianFilter(Array< OneD, NekDouble > &array, const StdMatrixKey &mkey)
virtual void v_LaplacianMatrixOp_MatFree(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
LibUtilities::NekManager< StdMatrixKey, DNekBlkMat, StdMatrixKey::opLess > m_stdStaticCondMatrixManager
int GetEdgeNcoeffs(const int i) const
This function returns the number of expansion coefficients belonging to the i-th edge.
boost::shared_ptr< StdSegExp > StdSegExpSharedPtr
boost::shared_ptr< DNekMat > DNekMatSharedPtr
virtual void v_StdPhysDeriv(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)
DNekMatSharedPtr GetStdMatrix(const StdMatrixKey &mkey)
void PhysTensorDeriv(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray_d0, Array< OneD, NekDouble > &outarray_d1)
Calculate the 2D derivative in the local tensor/collapsed coordinate at the physical points...
1D Gauss-Gauss-Legendre quadrature points
virtual void v_FwdTrans_BndConstrained(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
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.
void v_GeneralMatrixOp_MatOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
bool ConstFactorExists(const ConstFactorType &factor) const
BasisManagerT & BasisManager(void)
virtual int v_NumDGBndryCoeffs() const
int GetTotPoints() const
This function returns the total number of quadrature points used in the element.
int NumBndryCoeffs(void) const
virtual void v_GetBoundaryMap(Array< OneD, unsigned int > &outarray)
LibUtilities::BasisType GetEdgeBasisType(const int i) const
This function returns the type of expansion basis on the i-th edge.
void v_MultiplyByStdQuadratureMetric(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
virtual void v_IProductWRTBase_SumFacKernel(const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1)
The base class for all shapes.
virtual const LibUtilities::BasisKey v_DetEdgeBasisKey(const int i) const
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*y.
virtual void v_HelmholtzMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
virtual void v_IProductWRTBase_MatOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
virtual void v_BwdTrans_SumFac(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
virtual int v_DetCartesianDirOfEdge(const int edge)
void WeakDerivMatrixOp_MatFree(const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
void GetInteriorMap(Array< OneD, unsigned int > &outarray)
virtual void v_GetEdgeInteriorMap(const int eid, const Orientation edgeOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray)
DNekMatSharedPtr CreateGeneralMatrix(const StdMatrixKey &mkey)
this function generates the mass matrix
Principle Orthogonal Functions .
virtual void v_IProductWRTDerivBase(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
void v_GetEdgeToElementMap(const int eid, const Orientation edgeOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, int P=-1)
virtual void v_LocCoordToLocCollapsed(const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta)
Defines a specification for a set of points.
virtual void v_WeakDerivMatrixOp(const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
virtual void v_GetCoords(Array< OneD, NekDouble > &coords_0, Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2)
virtual int v_GetEdgeNcoeffs(const int i) const
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)
LibUtilities::BasisType GetBasisType(const int dir) const
This function returns the type of basis used in the dir direction.
virtual void v_BwdTrans_SumFacKernel(const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1)
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_MatOp(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
virtual int v_GetNedges() const
virtual int v_CalcNumberOfCoefficients(const std::vector< unsigned int > &nummodes, int &modes_offset)
void GetEdgeToElementMap(const int eid, const Orientation edgeOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, int P=-1)
virtual LibUtilities::BasisType v_GetEdgeBasisType(const int i) const
virtual int v_NumBndryCoeffs() const
int getNumberOfCoefficients(int Na, int Nb)
#define ASSERTL2(condition, msg)
Assert Level 2 – Debugging which is used FULLDEBUG compilation mode. This level assert is designed t...
virtual void v_IProductWRTBase(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Calculate the inner product of inarray with respect to the basis B=base0*base1 and put into outarray...
void MassMatrixOp_MatFree(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
virtual int v_GetVertexMap(int localVertexId, bool useCoeffPacking=false)
void BwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
This function performs the Backward transformation from coefficient space to physical space...
virtual void v_IProductWRTBase_SumFac(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true)
int GetNcoeffs(void) const
This function returns the total number of coefficients used in the expansion.
LibUtilities::NekManager< StdMatrixKey, DNekMat, StdMatrixKey::opLess > m_stdMatrixManager
NekDouble v_Integral(const Array< OneD, const NekDouble > &inarray)
Integrates the specified function over the domain.
virtual void v_BwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
virtual void v_FillMode(const int mode, Array< OneD, NekDouble > &array)
Fill outarray with mode mode of expansion.
void IProductWRTBase_SumFacKernel(const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0=true, bool doCheckCollDir1=true)
virtual void v_MassMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
boost::shared_ptr< Basis > BasisSharedPtr
virtual int v_GetEdgeNumPoints(const int i) const
void Zero(int n, T *x, const int incx)
Zero vector.
LibUtilities::PointsType GetPointsType(const int dir) const
This function returns the type of quadrature points used in the dir direction.
#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 InterpCoeff2D(const BasisKey &fbasis0, const BasisKey &fbasis1, const Array< OneD, const NekDouble > &from, const BasisKey &tbasis0, const BasisKey &tbasis1, Array< OneD, NekDouble > &to)
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Describes the specification for a Basis.
virtual LibUtilities::ShapeType v_DetShapeType() const
virtual int v_GetNverts() const
1D Gauss-Lobatto-Legendre quadrature points
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.
virtual void v_GetSimplexEquiSpacedConnectivity(Array< OneD, int > &conn, bool standard=true)
virtual void v_HelmholtzMatrixOp_MatFree(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
virtual void v_LaplacianMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
void FwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
This function performs the Forward transformation from physical space to coefficient space...