pcinvfunc.h

Go to the documentation of this file.

#ifndef _InvCapPress
#define _InvCapPress

#include <assert.h>
#include "fbsinverse.h"
#include "finterpolate.h"
#include "sfunctions.h"
#include "opointer.h"

class InvCapPress : public Function1D
{
private:
        Function1D &_P;           // 1D function that maps [a,b] , 0 <= a < b <=1 , onto (0,infty)
                                  // Note: Pc(a) may be +infty
                                  //       Pc(b)=0
                                  // Pc must be decreasing on (a,b]
        double _a, _b;            // [a,b] is the domain of f
        double _Pmax;             // bisection will be used for p on (Pmax,infty)
                                  // interpolation will be used for p on [0,Pmax]
        double _smax;             // Pc(smax)=Pmax
        unsigned _N;              // number of subintervals used for the interpolation
        double _Tol;              // tolerance used in the bisection method

        OPointer<BSInverse> PInv, PInvL;  // Pointers to bisection inverse functions
                                          // PInv - bisection inverse function to be used for p on (Pmax,infty)
                                          // PInvL - bisection inverse function to be interpolated and then
                                          //         used for p on [0,Pmax]

        OPointer<FInterpolate> PInvInt;   // Pointer to linear interpolate of the bisection inverse function
                                          // PInvInt - Interpolate to be used for p on [0,Pmax]
        

public:
        InvCapPress(Function1D &P, double a, double b, double smax, unsigned N, double Tol);

        virtual double operator()(double P, unsigned cmp) const;

        virtual ~InvCapPress();

};

#endif

---