45 T
computeEquilibrium(
int iPop, T rho,
const T u[DESCRIPTOR::d], T uSqr)
const override;
54 void computeApproxMomentum(T approxMomentum[DESCRIPTOR::d],
56 const T &rho,
const T u[DESCRIPTOR::d],
const T xi[DESCRIPTOR::d],
57 const std::vector<int> knownIndexes,
const std::vector<int> missingIndexes);
60 T computeError(
const T &rho,
const T u[DESCRIPTOR::d],
const T approxMomentum[DESCRIPTOR::d]);
62 void computeGradGradError(T gradGradError[DESCRIPTOR::d][DESCRIPTOR::d],
63 T gradError[DESCRIPTOR::d],
64 const T &rho,
const T u[DESCRIPTOR::d],
const T xi[DESCRIPTOR::d],
65 const T approxMomentum[DESCRIPTOR::d],
66 const std::vector<int> missingIndexes);
69 bool newtonRaphson(T xi[DESCRIPTOR::d],
70 const T gradError[DESCRIPTOR::d],
71 const T gradGradError[DESCRIPTOR::d][DESCRIPTOR::d]);
73 bool invert(
const T a[2][2],T b[2][2]);
75 bool invert(
const T a[3][3],T b[3][3]);