60 boost::ignore_unused(time);
63 int nVariables = physarray.size();
65 const Array<OneD, const int> &traceBndMap =
m_fields[0]->GetTraceBndMap();
69 int e, id1, id2, nBCEdgePts, eMax;
73 for (e = 0; e < eMax; ++e)
87 Vmath::Zero(nBCEdgePts, &Fwd[nVariables - 1][id2], 1);
91 Array<OneD, NekDouble> tmp(nBCEdgePts, 0.0);
101 Vmath::Smul(nBCEdgePts, -2.0, &tmp[0], 1, &tmp[0], 1);
107 &Fwd[1 + i][id2], 1, &Fwd[1 + i][id2], 1);
111 for (i = 0; i < nVariables; ++i)
116 ->UpdatePhys())[id1],
int m_spacedim
Space dimension.
Array< OneD, Array< OneD, NekDouble > > m_traceNormals
Trace normals.
int m_bcRegion
Id of the boundary region.
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
Array of fields.
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
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*x.
void Zero(int n, T *x, const int incx)
Zero vector.
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)