fchaventmobility.cpp

Go to the documentation of this file.

/*
 * fchaventmobility.cpp
 *
 *  Created on: Sep 15, 2008
 *      Author: marcosm
 */

#include "fchaventmobility.h"
FChaventMobility::FChaventMobility(double vw,double vo,double Srw,double Sro)
  :Function1D(1)
{
  m_M=vw/vo;

  MaxSo=1-Srw;
  MaxSw=1-Sro;
  MinSw=Srw;
  MinSo=Sro;


  //1_Srw2=1_Srw*1_Srw;
  //1_Sro=1-Sro;

  _c = m_M*MaxSo*MaxSo/(MaxSw*MaxSw);

  
  

}

FChaventMobility::~FChaventMobility()
{
}


double FChaventMobility::operator()(double x,unsigned cmp) const
{
  assert(cmp==0);
  if (x >= MaxSw)
     return 1.0;
  else if (x <= MinSw)
     return 0.0;
 double d1=MaxSw - x;
 double d2=x-MinSw;
 x=d1/d2;
 x*=x;
 return 1.0/(_c*x+1.0);
}





void FChaventMobility::getMinMaxValues(double a, double b,double &min,double &max) const 
{
  assert(a>=0.0);
  assert(b<=1.0);
  assert(a<=b);
  min=(*this)(a);
  max=(*this)(b);

}




 void FChaventMobility::setParameters(const VecDouble &v) 
{
  assert(v.size() ==2);
  m_M=v(0)/v(1);
}

---