Perform operation.
1092 int nPhys =
m_stdExp->GetTotPoints();
1094 Array<OneD, NekDouble> tmp0,tmp1,tmp2;
1095 Array<OneD, Array<OneD, NekDouble> > Diff(3);
1096 Array<OneD, Array<OneD, NekDouble> > out(3);
1097 out[0] = output0; out[1] = output1; out[2] = output2;
1099 for(
int i = 0; i <
m_dim; ++i)
1101 Diff[i] = wsp + i*ntot;
1113 1.0, &input[i*nPhys],
m_nquad0*m_nquad1,
1127 &Diff[1][i*nPhys+j*
m_nquad0*m_nquad1],
1133 Diff[1].get() + i*nPhys, 1,
1134 Diff[2].get() + i*nPhys, 1,
1135 Diff[2].get() + i*nPhys, 1);
1139 Diff[1].get() + i*nPhys, 1,
1140 Diff[1].get() + i*nPhys, 1);
1144 Diff[0].get() + i*nPhys, 1,
1145 Diff[1].get() + i*nPhys, 1,
1146 Diff[1].get() + i*nPhys, 1);
1150 Diff[0].get() + i*nPhys, 1,
1151 Diff[2].get() + i*nPhys, 1,
1152 Diff[2].get() + i*nPhys, 1);
1156 Diff[0].get() + i*nPhys, 1,
1157 Diff[0].get() + i*nPhys, 1);
1165 for(
int j = 1; j <
m_dim; ++j)
1168 Diff[j], 1, out[i], 1, out[i], 1);
Array< OneD, NekDouble > m_fac2
void Vvtvp(int n, const T *w, const int incw, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
vvtvp (vector times vector plus vector): z = w*x + y
StdRegions::StdExpansionSharedPtr m_stdExp
Array< OneD, NekDouble > m_fac1
Array< OneD, NekDouble > m_fac0
Array< TwoD, const NekDouble > m_derivFac
Array< OneD, NekDouble > m_fac3
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.