u3d_rtmod_defines.hpp

/* u3d_rtmod_defines.hpp */
/* $Id: u3d_rtmod_defines.hpp,v 1.1 2008/05/13 12:19:18 sad Exp $ */
#ifndef __u3d_rtmod_defines_hpp
#define __u3d_rtmod_defines_hpp

#include <mth/u3d_point.hpp>
#include <rbm/u3d_model_defines.hpp>
#include <u3d_raytr_defines.hpp>

#include <map>

/*                                                            DEFAULT this line must be straight */

//#define WF_DEBUG                 /* DEFAULT - OFF */                       /* | */
//#define WH_DEBUG                 /* DEFAULT - OFF */                       /* | */
//#define WF_DEBUG_SP_NUM 0        /* DEFAULT - OFF */                       /* | */
//#define U3D_DYN_DEBUG            /* DEFAULT - OFF */                       /* | */
//#define U3D_WF_FOR_MIGR          /* DEFAULT - OFF */                       /* | */

//#define U3D_SS_CODE_DEBUG          /* DEFAULT - ON */                        /* | */ /* DOES NOT WORK YET - assetion from broadcast stage appears */

#define U3D_SHEAR_ANGLES_RECALC    /* DEFAULT - ON */                        /* | */

#define U3D_NEW_RAY_STEP           /* DEFAULT - ON */                        /* | */
#define U3D_EXCLUDE_PERIMETER      /* DEFAULT - ?? */                        /* | */

#ifndef U3D_WF_FOR_MIGR
#define CALC_DYN                   /* DEFAULT - ON */                        /* | */
#define CALC_GEOM_DIV              /* DEFAULT - ON */                        /* | */
#endif

//#define SAVE_HODS                  /* DEFAULT - ON */                        /* | */





#ifdef WF_DEBUG
#define U3D_RTMOD_BOUND_0     1
#define U3D_RTMOD_BOUND_NUM   2
#endif






//---------------------------------------------------------------------
#define U3D_RTMOD_STEP_LEN 10.

#define U3D_RTMOD_F_EPS 10.                     /* For initial (rough) selection of incident cells (during two surface sources comparison) */

#define U3D_RTMOD_MIN_ANGLE_STEP 1.e-6          /* Minimal angle step (during direct wave upper part calculation) */

#define U3D_RTMOD_INIT_ANGLE_STEP 1.e-3         /* Initial angle step value (during direct wave upper part calculation) */
  
#define U3D_RTMOD_MAX_DIST 50.                  /* Maximal distance between the last point of the dyn.ray and receiver 
                           (to detect the evident error) */

#define U3D_RTMOD_LENGTH_EPS 5.                 /* Final precision of receivers shooting (probably must be equal to the SS cell size) */

#define U3D_RTMOD_SS_SHADOW_CELLS_PART 0.05     /* Allowed part of the shadowed cells on surface source */

#define U3D_FACE_ACTIVE_ZONES_NUMBER 10         /* Number of zones on surface in current direction (for surface source calculation) */

//===============================================================================
struct SSPQCell {
  double st_f;       
  double st_az;      
  double t;          
  double f;          
  double az;         
  double n_f;        
  double n_az;       
};


//===============================================================================
enum U3dWaveletType {
  RICKER,
  MIN_PHASE,
  VSP_DATA
};


//===============================================================================
struct U3dWaveInfoCell {
  double st_f;      
  double st_az;     
  double t;         
  double f;         
  double az;        
  int recId;        
  double ampl;      
};
typedef std::vector<U3dWaveInfoCell> U3dWaveInfo;



typedef std::vector<U3dPoint> U3dPoints;

//===============================================================================
struct U3dWaveHodoCell {
  U3dWaveInfoCell wi_cell;    
  U3dCodeAr code;             
  U3dPoints rp_ar;            
};
typedef std::vector<U3dWaveHodoCell> U3dWaveDynHodo;






//===============================================================================
struct U3dWaveHodItem {
  int recId;              
  U3dPoint recp;          
  double t;               
  U3dPoints ray_path;     
  U3dWaveHodItem() {}
  
  U3dWaveHodItem(int rec_id, U3dPoint p, double time, std::vector<U3dPoint> *rp = NULL) :
    recId(rec_id),
    recp(p),
    t(time)
  { if (rp)  ray_path = *rp; }
  
  //U3dWaveHodItem(const U3dWaveHodItem &item) :
  //recId(item.recId),
  //recp(item.recp),
  //t(item.t),
  //ray_path(item.ray_path)
  //{}
  
  ~U3dWaveHodItem() {}
  
};
typedef std::map<int,U3dWaveHodItem> U3dWaveHod;

//---------------------------------------------------------------------
struct U3dComplexWaveHod {
  U3dWaveHod hod;             
  U3dCodeAr code;             
};
typedef std::vector<U3dComplexWaveHod> U3dComplexWaveHodsAr;

//---------------------------------------------------------------------
struct U3dSPHods {
  int spId;
  U3dPoint sp;
  U3D_WAVE_TYPE dir_w_type;
  U3dWaveHod direct_hod;                          
  U3dComplexWaveHodsAr single_refl_hods;          
  U3dComplexWaveHodsAr single_refl_conv_hods;     
  U3dComplexWaveHodsAr single_trans_conv_hods;    
};
typedef std::vector<U3dSPHods> U3dSPHodsAr;



struct U3dRayPathItem {
  int spId;       /* Shot point identifier */
  int recId;      /* Receiver point identifier */
  int ray_type;   /* Type of ray path (see u3d_raytr_defines.hpp) */
  int boundId;    /* Incident bound identifier */
  U3dPoint rp;    /* Ray path point */
  int ray_pathId; /* Ray path identifier (unique among all ray paths) */
  int pointId;    /* Sequential index of the point in the ray path. */
};
typedef std::vector<U3dRayPathItem> U3dRayPath;


//---------------------------------------------------------------------
#endif /* u3d_migr_defines.hpp */