physics.hpp

#ifndef PHYSICS_HPP
#define PHYSICS_HPP
 
class Model{
public:
  
  virtual double s(double)=0;
  virtual double ds(double)=0;
  virtual double kr(double)=0;
  virtual double dkr(double)=0;
};
 
 
class ModelBC{
public:
  double s(double);
  double ds(double);
  double kr(double);
  double dkr(double);
};
 
 
class Physics{
public:
 
  double g;
  double rho;
  
  double phi;
  double k0;
 
  double nivrivsat;
 
 
  double (*s)(double);
  double (*ds)(double);
  void (*s_ds)(double,double&,double&);
  
  double (*kr)(double);
  double (*dkr)(double);
  void (*kr_dkr)(double,double&,double&);
  
 
  //---------------------
  // Models descriptions 
  //---------------------
 
  double model_datas[6];
 
  //------------------------
  // Brooks and Corey model
  //------------------------
  
  double s_BC(double P);
  double ds_BC(double P);
  void (s_ds_BC)(double P,double& v,double& dv);
 
  double kr_BC(double P);
  double dkr_BC(double P);
  void (kr_dkr_BC)(double P,double& v,double& dv);
 
};
#endif