void SinglePvtDeadSpline::mu(const int n, const double* p, const double* /*z*/, double* output_mu) const { // #pragma omp parallel for for (int i = 0; i < n; ++i) { output_mu[i] = viscosity_(p[i]); } }
void SinglePvtDead::mu(const int n, const double* p, const double* /*r*/, double* output_mu, double* output_dmudp, double* output_dmudr) const { // #pragma omp parallel for for (int i = 0; i < n; ++i) { output_mu[i] = viscosity_(p[i]); output_dmudp[i] = viscosity_.derivative(p[i]); } std::fill(output_dmudr, output_dmudr + n, 0.0); }