pcinvfunc.cpp

Go to the documentation of this file.

#include "pcinvfunc.h"

InvCapPress::InvCapPress(Function1D &P, double a, double b, double smax, unsigned N, double Tol)
:_P(P)
{
        _a=a;
        _b=b;
        assert(_a < _b);
        _N=N;
        _Tol=Tol;
        _smax=smax;
        _Pmax=_P(_smax);
        // build appropriate bisection inverse and its linear interpolate for p on [0,Pmax]
        PInv = new BSInverse(_P,_smax,_b,_Tol);
        PInvInt =  new FInterpolate((*PInv),0.0,_Pmax,_N);
        // build bisection inverse for p on (Pmax,infty)
        PInvL = new BSInverse(_P,_a,_smax,_Tol);
}

double InvCapPress::operator()(double P, unsigned cmp=0) const
{
        assert(isInDomain(P,cmp));
        // use interpolation if P <= Pmax
        if (P <= _Pmax) return (*PInvInt)(P,cmp);
        // use bisection of P > Pmax
        return (*PInvL)(P,cmp);
}

InvCapPress::~InvCapPress(){}

---