ray_tracing.hpp

// ray_tracing.hpp
// $Id: ray_tracing.hpp,v 1.11 2008/07/30 14:29:22 guser1 Exp $
//------------------------------------------------------------------------------
#ifndef __RayTracingH
#define __RayTracingH

#include "vel2dmod_raytr_global.hpp"
#include "rt_storage.hpp"


//---------------------------------------------------------------------
void OverturnD(double *&M, int N);
void OverturnI(int *&M, int N);

//class RayInfo;

//==============================================================================

class RayTracing
{
public:
  
  RayTracing(Vel2dModel *mod2d);
  ~RayTracing();

  RayInfo*     RayShoot(double f, double startX, double startZ, bool waySave, CodeAr rCode , bool calcSD,
            int_ar *bNums=NULL);

  void         AddObSysPtr(ObservSys*);  

  void         SetDL(double dl);
  double       GetDL();

  
  Vel2dModel *Mod2d;

  ObservSys    *ObsSys;

  void         CalcIntParams2(double x, double z, double f, int rtr,
                  double &x0,double &z0, double &x1, double &z1);

  static bool Snellius(double vin, double teta0, double vout, int rtr,
               double iAngle, double &teta1);

  
protected:

  STEP_RESULT  RayStep(double &x, double &z, double &f, double &v, double df, double &dt,
               double &resultAngle, double &resultCurvity, int &iBoundNum,
               WAVE_TYPE vPS, bool prevInt, bool noIntersect=false);
  
  bool         CalcRayAmpl(RayInfo *Ray);    
  
  double       CalcRayCoeff(double x, double z, double x0, double z0, double x1, double z1, int bn1,
                int bn2, int boundN, WAVE_TYPE wt1, WAVE_TYPE wt2, double iAngle);

  void         CalcIntParams(double x, double z, double x0, double z0, int rtr,
                 double &x1, double &z1);
  void         CalcIntParams(double x, double z, double &x0, double &z0, int rtr,
                 double f, double &x1, double &z1);

  //  /** calculation of the ray points before and after bound intersection (rt_calc.cpp) */
//    void         CalcIntParams2(double x, double z, double f, int rtr,
//                    double &x0,double &z0, double &x1, double &z1);

//    /** calculation of the ray reflection/refraction angle on the bound (ray_tracing.cpp) */
//    bool         Snellius(double vin, double teta0, double vout, int rtr,
//              double iAngle, double &teta1);

 private:

  
  double       Time(double x0, double z0, double x1, double z1, double v0, double v1);
  

  int          Chorda1(double x0, double z0, double k, double alpha, double F, double teta0,
               double v0, double &x1, double &z1, double &teta1);
  
  bool         StepOnRay(double x1, double z1, double &x2, double &z2, double f1, double &f2,
             double df, WAVE_TYPE vPS, bool prevInt);

  void         ModLimInt(double x1, double z1, double x2, double z2, double &res_x, double &res_z);
  
  bool         CalcRaySD(RayInfo *Ray);        
  bool         CalcRaySD2(RayInfo *Ray);       
  bool         CalcRaySD3(RayInfo *Ray);
  
  double      gridSize;    
  bool        isObSys;     
  double      DL;          
  double      iX[50], iZ[50], iAngles[50], iCurvity[50];
  int         ibNum[50], ibrd[50];
};

#endif