Example #1
0
    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]);
        }
    }
Example #2
0
    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);

        }