fmobfromkpermlib.h

Go to the documentation of this file.

#include "sfunctions.h"
#include "numericmethods.h"



using NumericMethods::adjustBounds;

class FFracFluxFromRelK : public Function1D
{
 private:

 protected:
  Function1D &_fk1,&_fk2;
  double _sr1,_sr1Max;
  double _M;
 public:
  virtual double operator()(double x,unsigned cmp=0) const
  {
    /*
    if (x < _sr1) 
      return 0.0;
    if (x >= _sr1Max)
      return 1;
    */
    double fk1=_fk1(x);
    return fk1/(fk1 + _M*_fk2(x));
  }

 FFracFluxFromRelK(Function1D &fk1,Function1D &fk2,double sr1,double sr2,double v1,double v2)
   :_fk1(fk1),_fk2(fk2),_sr1(sr1),_sr1Max(1.0-sr2)
  {
    
    _M=v1/v2;
  }


  virtual ~FFracFluxFromRelK(){}

};



class FMobTFromRelK : public Function1D
{
 private:

 protected:
  Function1D &_fk1,&_fk2;
  VecDouble _viscs;
  double _sr1,_sr1Max;
 public:
  virtual double operator()(double x,unsigned cmp=0) const
  {
    /*if (x<=_sr1) //adjustBounds(x,_sr1,_sr1Max);
      return  _fk2(x)/_viscs(1);
    if (x>=_sr1Max)
    return _fk1(x)/_viscs(0);*/
    return _fk1(x)/_viscs(0) + _fk2(x)/_viscs(1);
  }

 FMobTFromRelK(Function1D &fk1,Function1D &fk2,double sr1,double sr2,double v1,double v2)
   :_fk1(fk1),_fk2(fk2),_sr1(sr1),_sr1Max(1.0-sr2)
  {
    _viscs.set(v1,v2);
  }
  virtual ~FMobTFromRelK(){}
};


class FFracGravFromRelK : public Function1D
{
 private:

 protected:
  Function1D &_fk1,&_fk2;
  VecDouble _viscs;
  double _sr1,_sr1Max;
 public:
  virtual double operator()(double x,unsigned cmp=0) const
  {
    /*
    if (x < _sr1 || x > _sr1Max)
    return 0.0;*/
    return _fk1(x)*_fk2(x)/(_viscs(1)*_fk1(x) + _viscs(0)*_fk2(x));
  }

 FFracGravFromRelK(Function1D &fk1,Function1D &fk2,double sr1,double sr2,double v1,double v2)
   :_fk1(fk1),_fk2(fk2),_sr1(sr1),_sr1Max(1.0-sr2)
  {
    _viscs.set(v1,v2);
  }


  virtual ~FFracGravFromRelK(){}

};



class DFFracGravFromRelK : public Function1D
{
 private:

 protected:
  Function1D &_fk1,&_fk2,&_dfk1,&_dfk2;
  double _sr1,_sr1Max;
  VecDouble _viscs;
 public:
  virtual double operator()(double x,unsigned cmp=0) const
  {
    /*
    if ( x < _sr1 || x > _sr1Max)
      return 0.0;
    */
    double f1=_fk1(x);
    double f2=_fk2(x);
    double df1=_dfk1(x);
    double df2=_dfk2(x);
    double denom=_viscs(1)*f1 + _viscs(0)*f2;
    
    return ((f1*df2+df1*f2)*denom - f1*f2*(_viscs(1)*df1 + _viscs(0)*df2))/(denom*denom); 
  }

 DFFracGravFromRelK(Function1D &fk1,Function1D &fk2,Function1D &dfk1,Function1D &dfk2,double sr1,double sr2,double v1,double v2)
   :_fk1(fk1),_fk2(fk2),_dfk1(dfk1),_dfk2(dfk2),_sr1(sr1),_sr1Max(1.0-sr2)
  {
    _viscs.set(v1,v2);
  }


  virtual ~DFFracGravFromRelK(){}

};







class DFFracFluxFromRelK : public Function1D
{
 private:

 protected:
  Function1D &_fk1,&_fk2,&_dfk1,&_dfk2;
  double _M,_sr1,_sr1Max;
 public:
  virtual double operator()(double x,unsigned cmp=0) const
  {
    /*
    if ( x < _sr1 || x > _sr1Max)
      return 0.0;
    */
    double fk1=_fk1(x);
    double dfk1=_dfk1(x);
    double dd = fk1 + _M*_fk2(x);
    return (dfk1*dd - fk1*(dfk1 + _M*_dfk2(x)))/(dd*dd);
  }

 DFFracFluxFromRelK(Function1D &fk1,Function1D &fk2,Function1D &dfk1,Function1D &dfk2,double sr1,double sr2,double v1,double v2)
   :_fk1(fk1),_fk2(fk2),_dfk1(dfk1),_dfk2(dfk2),_sr1(sr1),_sr1Max(1.0-sr2)
  {
    _M=v1/v2;
  }


  virtual ~DFFracFluxFromRelK(){}

};

---