dflinearmobilityproduct.cpp

Go to the documentation of this file.

#include "dflinearmobilityproduct.h"
//#include "exception.h"
//#include "numericmethods.h"



DFLinearMobilityProduct::DFLinearMobilityProduct(double vW,double vO,double srw, double sro)
  :Function1D(1)
{

  _vw = vW;
  _vo = vO;
  _MaxSw=1.-sro;
  _Srw=srw;
  _c1=_MaxSw*_MaxSw*_vo;
  _c2=(1.0-_Srw)*(1.0-_Srw)*_vw;
  _prob.sampleFunction(*this,0,0,1,100000);
}


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

  // register double krw = x-_Srw;
  // register double kro = _MaxSw-x;
  
  assert(cmp==0);
  return (_MaxSw-x)/(_c2*_c1*((x-_Srw)/_c2+(_MaxSw-x)/_c1))-(x-_Srw)/(_c2*_c1*((x-_Srw)/_c2+(_MaxSw-x)/_c1))-(x-_Srw)*(_MaxSw-x)*(1/_c2-1/_c1)/(_c2*_c1*pow(((x-_Srw)/_c2+(_MaxSw-x)/_c1),2));
}

void DFLinearMobilityProduct::getMinMaxValues(double a, double b,double &min,double &max) const
{
  _prob.getMinMaxValues(a,b,min,max);
}


DFLinearMobilityProduct::~DFLinearMobilityProduct() 
{
 
}

---