fbucleyleverettgravitymob.cpp

Go to the documentation of this file.

#include "fbucleyleverettgravitymob.h"
#include "exception.h"
#include "numericmethods.h"



FBucleyLeverettGravityMob::FBucleyLeverettGravityMob(double vW,double vO,double srw, double sro) 
{

  _vw = vW;
  _vo = vO;
  _MaxSw=1.-sro;
  _Srw=srw;
  _c1=_MaxSw*_MaxSw*_vo;
  _c2=(1.0-_Srw)*(1.0-_Srw)*_vw;
  
  //Coefficients for the polynomium in the expression of the first derivative
  //of the function
  //m_a = m_vo + m_vw;
  //m_b = -(m_vo + 4*m_vw);
  //m_c = 6*m_vw;
  //m_d = -4*m_vw;
          
  m_scanInf = new ProbeFunction1D(*this,0,0,1,10000);
  
}


 FBucleyLeverettGravityMob::FBucleyLeverettGravityMob(FBucleyLeverettGravityMob& f) 
{
  *this=f;
  m_scanInf = new ProbeFunction1D(*this,0,0,1,10000);
}


double FBucleyLeverettGravityMob::operator()(double x,unsigned cmp) const
{
   if (x >= _MaxSw)
     x=_MaxSw;
 else if (x < _Srw)
     x=_Srw;

  register double kw = x-_Srw;
  register double ko = x-_MaxSw;
  register double kw2=kw*kw;
  register double ko2=ko*ko;
  return kw2*ko2/(_c1*kw2 + _c2*ko2);
}



FBucleyLeverettGravityMob::~FBucleyLeverettGravityMob() 
{
  delete m_scanInf;
}

void FBucleyLeverettGravityMob::getMinMaxValues(double a, double b,double &min,double &max) const
{
  assert(a<=b);
  NumericMethods::adjustBounds(a,0,1);
  NumericMethods::adjustBounds(b,0,1);
  m_scanInf->getMinMaxValues(a,b,min,max);
}

---