u3d_raytr_defines.hpp

/* u3d_raytr_defines.hpp */
/* $Id: u3d_raytr_defines.hpp,v 1.26 2008/04/25 11:37:42 guser1 Exp $ */

#ifndef __u3d_raytr_defines_hpp
#define __u3d_raytr_defines_hpp

#define U3D_RAY_ANGLE_CHORD_STEP 0.01
#define U3D_RAY_LINE_CHORD_STEP 50.0
#define U3D_MIN_RAY_LINE_STEP 0.5

#define U3D_SEISMIC_TIME_MULTIPLIER 1000.0

#define U3D_NHB_CONF_NUM 6


#define U3D_RAY_STEP_ERROR_MESSAGE


#include <complex>
#define dcomplex complex<double>

#include <u3d_model_defines.hpp> /* for U3D_WAVE_TYPE */


//#ifdef __cplusplus
//extern "C" {
//#endif /* __cplusplus */

const double PI_2 = M_PI*2.;

//#ifdef __cplusplus
//};
//#endif /* __cplusplus */



enum U3D_STEP_RESULT { U3D_STEP_DONE, U3D_OUT_OF_MODEL, U3D_BOUND_INTERSECTION, 
               U3D_WELL_INTERSECTION, U3D_STEP_ERROR, U3D_UPGOING_WAVE, 
               U3D_COMPLEX_CELL, U3D_NO_PARAMS }; 

//============================================================================
struct U3dCodeItem
{
  int        bNum;   
  U3D_WAVE_TYPE  wType;  
  bool       trFlag; 
  bool operator==(const U3dCodeItem& p) const {
    return ((bNum == p.bNum ) && (wType == p.wType) && (trFlag == p.trFlag) );
  }

};


typedef vector<U3dCodeItem> U3dCodeAr;




//============================================================================
//static 
bool analyseRayCode(U3dCodeAr code, int ciId, int &type);

void analyseAngles(double &f, double &az);
  
string printRayType(int ray_id);




struct U3dRayPoint {
  U3dPoint p;     
  double f;       
  double az;      
  double t;       
  bool operator==(const U3dRayPoint&) const;
  bool operator!=(const U3dRayPoint&) const;
};

typedef vector<U3dRayPoint> U3dRayPointsAr;



//============================================================================
struct U3dRayIntSt
{
  int          bNum1, bNum2;    
  int          boundNum;        
  U3dPoint     point;           
  double       coeff;            
  double       t;               
  double       f1;              
  double       az1;             
  double       f2;              
  double       az2;             
  double       iAngle;          
  double       iAz;             
  double       iCurv;          
};

typedef vector<U3dRayIntSt> U3dRayIntStAr;



//============================================================================
struct U3dRayWIntSt
{
  U3dPoint p;
  double   f, az;
  double   t;
  double   ampl;
};

typedef vector<U3dRayWIntSt> U3dRayWIntStAr;



//============================================================================
enum U3dReflCoefType {
  ZOEPPRITZ,
  AKI_RICHARDS,
  SHUEY
};

//===============================================================================
#define U3D_P      0
#define U3D_PPup   1
#define U3D_PSup   2
#define U3D_PSdown 3
//---------------------- exotic ---------------------------------------
#define U3D_S      4
#define U3D_SSup   5
#define U3D_SPup   6
#define U3D_SPdown 7


#endif  /* u3d_raytr_defines.hpp */