u_Sh_rd_calc.hpp

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

#ifndef __u_sh_rd_calc_hpp
#define __u_sh_rd_calc_hpp

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



double getA(double B, double sigma);
  
double getR0(double vp0, double vp1,  double den0, double den1);

double getB(double vp0, double vp1, double R0);

double getPoissonsRatio(double vp, double vs);



class URDCalcShuey {

public:
  
  URDCalcShuey();
  
  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 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