u3d_raytracing.hpp

/* u3d_raytracing.hpp */
/* $Id: u3d_raytracing.hpp,v 1.3 2008/06/02 12:02:11 sasha Exp $ */

#ifndef __u3d_raytracing_hpp
#define __u3d_raytracing_hpp

#include <cmath>
#include <rbm/u3d_model.hpp>

#include "u3d_ray.hpp"
#include "u3d_raytr_defines.hpp"

//============================================================================
class U3dRayTracing
{ 
 public:

  U3dRayTracing(U3dModel *mod3d);                                         
  virtual ~U3dRayTracing();

  virtual U3dRay* rayShoot(U3dPoint startP, double f, double az, 
               U3dCodeAr rCode, /*IntAr *bNums=NULL,*/ 
               bool waySave=false);
  bool         snellius(double vin, double teta0, double az, 
            double vout, int rtr, double iAngle, 
            double iAz, double &teta1, double &az1);

  U3dModel* getModel() const { return m_model; };
  
  
 protected:

  U3dModel   *m_model;      
  U3dBox     *m_modelBox;   
  U3D_STEP_RESULT  RayStep(U3dRayPoint &rp, double df, double &dt, 
               double &iAngle, double &iAz, 
               double &iCurvity, int &iBoundNum, 
               U3D_WAVE_TYPE vPS, bool prevInt, 
               bool noIntersect=false);
  
  U3D_STEP_RESULT  rayStep(U3dRayPoint &rp, double df, double &dt, 
               double &iAngle, double &iAz, 
               double &iCurvity, int &iBoundNum, int &bn1, int &bn2,
               U3D_WAVE_TYPE vPS, bool prevInt, 
               bool noIntersect=false);
  
  bool         intersectCellTriangles(const U3dGridCell* cell, 
                      U3dPoint p1, 
                      U3dPoint p2, 
                      U3dPoint &ip, int &iBoundNum, double &iAngle, double &iAz, 
                      int debug_level=0);
  
  bool         processNeighbourCells(const U3dGridCell* cell1, const U3dGridCell* cell2, 
                     U3dPoint p1, U3dPoint p2, 
                     U3dPoint &ip, int &iBoundNum, double &iAngle, double &iAz, 
                     int debug_level=0);
  
  double       getTime(U3dPoint p0, U3dPoint p1, double v0);
  
  bool         stepOnRay(U3dRayPoint p1, U3dRayPoint &p2,
             double df, U3D_WAVE_TYPE vPS, bool prevInt, 
             double &vel);
  
  void         stepOnChord(double k, double fK, double azK, double df, double vel, 
               U3dRayPoint rp0, U3dRayPoint &rp1);
  
  void         CalcIntParams(U3dRayPoint p, U3dPoint p0, int rtr, U3dPoint &p1);
  
  void         calcIntNeighPoints(U3dPoint p_int, double f, double az, U3dPoint &p0, U3dPoint &p1);
  
  void         getEndPoint(U3dPoint p0, U3dPoint &p1);
  
  double     dl;             
   
private:
  
  double     BIGGEST_RADIUS; 
 
};


#endif /* u3d_raytracing.hpp */