u_AR_rd_calc.hpp

/* u__AR_rd_calc.hpp */
/* $Id: u_AR_rd_calc.hpp,v 1.1 2008/06/04 12:31:38 sasha Exp $ */

#ifndef __u_ar_rd_calc_hpp
#define __u_ar_rd_calc_hpp

#include "u_raytr_defines.hpp" /* FOR RTR_TYPES */

class URDCalcAkiRichards {

public:
  
  URDCalcAkiRichards();
  
  static double getCoeff(int rtr_type, double f0, double f1,
             double vp0, double vs0, double den0,
             double vp1, double vs1, double den1);
  
  static double getRDir(double f, 
            double v0, double den0,
            double v1, double den1);
  
  static double getRPPup(double f0, 
             double vp0, double vs0, double den0,
             double vp1, double vs1, double den1);
  
  static double getRPSup(double f0, double f1, 
             double vp0, double vs0, double den0,
             double vp1, double vs1, double den1);
  
  static double getRPSdown(double f0, double f1, 
               double vp0, double vs0, double den0,
               double vp1, double vs1, double den1);
  
  static double getRPS(double f0, double f1, 
               double vp0, double vs0, double den0,
               double vp1, double vs1, double den1,
               bool up = true);
  
  static double getRSSup(double f1,
             double vp0, double vs0, double den0,
             double vp1, double vs1, double den1);
  
  static double getRSPup(double f0, double f1,
             double vp0, double vs0, double den0,
             double vp1, double vs1, double den1);
  
  static double getRSPdown(double f0, double f1,
               double vp0, double vs0, double den0,
               double vp1, double vs1, double den1);
  
};



#endif