p_straight_elements.H

/* p_straight_elements.H */
/* $Id: p_straight_elements.H,v 1.6 2001/05/07 17:12:10 hoh Exp $ */
#ifndef __p_straight_elements_H
#define __p_straight_elements_H

#include <mod3d/mod3d_raytrace.H>
#include "rays_straight.H"

class pl_straight_body:public pl_intern_body{
protected:
  // ¾ÑàÐÑÞâÐÝÝëÕ Ø ÕéÕ ÝÕ ÞÑàÐÑÞâÐÝÝëÕ ÛãçØ  rays arProcessedRays;
  rays arNewRays;

  // ÁÞÞâÒÕâáâÒãîéØÕ /\ ßÕàÕáÕçÕÝØï ÝÐ ÒåÞÔÕ Ò âÕÛÞ
  IntAr aiProcessedCrossings;
  IntAr aiNewCrossings;
    
  void addh_ray (ray *pRNew,int iCN);
public:
  pl_straight_body(geophys_params gp):pl_intern_body(gp){};
  pl_straight_body(any_borders &ab,geophys_params gp);
  ~pl_straight_body(){arProcessedRays.clean();arNewRays.clean();};
  void commit_pass ();
  const IntAr& get_all_crossings ()const;
  const rays& get_all_rays ()const;
  // ÍâÞâ ÜÐááØÒ, Ú áÞÖÐÛÕÝØî, ØÝÞÓÔÐ ÝãÖÝÞ ØáßàÐÒÛïâì...
  IntAr& get_new_crossings ();
  const rays &get_new_rays ()const;
  Logic store_ray (ray* pRNew,int iBorder);
  void clear_rays(){arNewRays.clean(); arProcessedRays.clean();};
  Logic reflecting_border(int iB)const;
};

typedef pl_straight_body* p_pl_straight_body;
typedef ObjectsAr<p_pl_straight_body> pl_straight_bodies;




// ¿àÞåÞÖÔÕÝØÕ ßÞ ßãâØ, ×ÐÔÐÝÝÞÜã èÐÑÛÞÝÞÜ
class border_hit_info:public hit_point_info{
protected:
  rays_straight * pr;
  IntAr iaTemplate;
  int iReflectBorder;
  int iHitBorder;
  int iBody;

  Logic verify_ray (ray*r)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  border_hit_info(rays_straight*,int,int,
          ray*prToPass,int,
          double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
          int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
          double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);  
  border_hit_info(rays_straight*,ray*,int,int,
          double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
          int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
          double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
  int get_ref_bord_num();
  int get_hit_bord_num();
};

// ¿àÞåÞÖÔÕÝØÕ ßÞ ßãâØ, ×ÐÔÐÝÝÞÜã èÐÑÛÞÝÞÜ, áâàÞÓÞ ÔÐÛìèÕ, ÝÕÖÕÛØ íâÞâ èÐÑÛÞÝ
class hit_through_border_info:public general_hit_info{
protected:
  rays_straight * pr;
  ray * prTemplate;
  int iReflectBorder;
  geophys_params gp;
  raytr_params rp;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);

  virtual void
  get_ray_parameters (const point3d &pWhere,const crossing&c,double &dFi,
              double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);

  hit_through_border_info(rays_straight *,ray *,
              double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
              int i3=PM_MAX_ITERATIONS_COUNT_D,
              double d2=PM_ALPHA_STOP_D/10.,
              int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D);
};

#endif /* p_straight_elements.H */

/*  Melnikov George Yu.  2-Feb-2000  */ 
/*          ...          */

  rays arProcessedRays;
  rays arNewRays;

  // ÁÞÞâÒÕâáâÒãîéØÕ /\ ßÕàÕáÕçÕÝØï ÝÐ ÒåÞÔÕ Ò âÕÛÞ  IntAr aiProcessedCrossings;
  IntAr aiNewCrossings;
    
  void addh_ray (ray *pRNew,int iCN);
public:
  pl_straight_body(geophys_params gp):pl_intern_body(gp){};
  pl_straight_body(any_borders &ab,geophys_params gp);
  ~pl_straight_body(){arProcessedRays.clean();arNewRays.clean();};
  void commit_pass ();
  const IntAr& get_all_crossings ()const;
  const rays& get_all_rays ()const;
  // ÍâÞâ ÜÐááØÒ, Ú áÞÖÐÛÕÝØî, ØÝÞÓÔÐ ÝãÖÝÞ ØáßàÐÒÛïâì...
  IntAr& get_new_crossings ();
  const rays &get_new_rays ()const;
  Logic store_ray (ray* pRNew,int iBorder);
  void clear_rays(){arNewRays.clean(); arProcessedRays.clean();};
  Logic reflecting_border(int iB)const;
};

typedef pl_straight_body* p_pl_straight_body;
typedef ObjectsAr<p_pl_straight_body> pl_straight_bodies;




// ¿àÞåÞÖÔÕÝØÕ ßÞ ßãâØ, ×ÐÔÐÝÝÞÜã èÐÑÛÞÝÞÜ
class border_hit_info:public hit_point_info{
protected:
  rays_straight * pr;
  IntAr iaTemplate;
  int iReflectBorder;
  int iHitBorder;
  int iBody;

  Logic verify_ray (ray*r)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  border_hit_info(rays_straight*,int,int,
          ray*prToPass,int,
          double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
          int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
          double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);  
  border_hit_info(rays_straight*,ray*,int,int,
          double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
          int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
          double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
  int get_ref_bord_num();
  int get_hit_bord_num();
};

// ¿àÞåÞÖÔÕÝØÕ ßÞ ßãâØ, ×ÐÔÐÝÝÞÜã èÐÑÛÞÝÞÜ, áâàÞÓÞ ÔÐÛìèÕ, ÝÕÖÕÛØ íâÞâ èÐÑÛÞÝ
class hit_through_border_info:public general_hit_info{
protected:
  rays_straight * pr;
  ray * prTemplate;
  int iReflectBorder;
  geophys_params gp;
  raytr_params rp;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);

  virtual void
  get_ray_parameters (const point3d &pWhere,const crossing&c,double &dFi,
              double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);

  hit_through_border_info(rays_straight *,ray *,
              double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
              int i3=PM_MAX_ITERATIONS_COUNT_D,
              double d2=PM_ALPHA_STOP_D/10.,
              int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D);
};

#endif /* p_straight_elements.H */

/*  Melnikov George Yu.  2-Feb-2000  */ 
/*          ...          */

  IntAr aiProcessedCrossings;
  IntAr aiNewCrossings;
    
  void addh_ray (ray *pRNew,int iCN);
public:
  pl_straight_body(geophys_params gp):pl_intern_body(gp){};
  pl_straight_body(any_borders &ab,geophys_params gp);
  ~pl_straight_body(){arProcessedRays.clean();arNewRays.clean();};
  void commit_pass ();
  const IntAr& get_all_crossings ()const;
  const rays& get_all_rays ()const;
  // ÍâÞâ ÜÐááØÒ, Ú áÞÖÐÛÕÝØî, ØÝÞÓÔÐ ÝãÖÝÞ ØáßàÐÒÛïâì...
  IntAr& get_new_crossings ();
  const rays &get_new_rays ()const;
  Logic store_ray (ray* pRNew,int iBorder);
  void clear_rays(){arNewRays.clean(); arProcessedRays.clean();};
  Logic reflecting_border(int iB)const;
};

typedef pl_straight_body* p_pl_straight_body;
typedef ObjectsAr<p_pl_straight_body> pl_straight_bodies;




// ¿àÞåÞÖÔÕÝØÕ ßÞ ßãâØ, ×ÐÔÐÝÝÞÜã èÐÑÛÞÝÞÜclass border_hit_info:public hit_point_info{
protected:
  rays_straight * pr;
  IntAr iaTemplate;
  int iReflectBorder;
  int iHitBorder;
  int iBody;

  Logic verify_ray (ray*r)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  border_hit_info(rays_straight*,int,int,
          ray*prToPass,int,
          double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
          int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
          double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);  
  border_hit_info(rays_straight*,ray*,int,int,
          double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
          int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
          double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
  int get_ref_bord_num();
  int get_hit_bord_num();
};

// ¿àÞåÞÖÔÕÝØÕ ßÞ ßãâØ, ×ÐÔÐÝÝÞÜã èÐÑÛÞÝÞÜ, áâàÞÓÞ ÔÐÛìèÕ, ÝÕÖÕÛØ íâÞâ èÐÑÛÞÝ
class hit_through_border_info:public general_hit_info{
protected:
  rays_straight * pr;
  ray * prTemplate;
  int iReflectBorder;
  geophys_params gp;
  raytr_params rp;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);

  virtual void
  get_ray_parameters (const point3d &pWhere,const crossing&c,double &dFi,
              double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);

  hit_through_border_info(rays_straight *,ray *,
              double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
              int i3=PM_MAX_ITERATIONS_COUNT_D,
              double d2=PM_ALPHA_STOP_D/10.,
              int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D);
};

#endif /* p_straight_elements.H */

/*  Melnikov George Yu.  2-Feb-2000  */ 
/*          ...          */

class border_hit_info:public hit_point_info{
protected:
  rays_straight * pr;
  IntAr iaTemplate;
  int iReflectBorder;
  int iHitBorder;
  int iBody;

  Logic verify_ray (ray*r)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  border_hit_info(rays_straight*,int,int,
          ray*prToPass,int,
          double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
          int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
          double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);  
  border_hit_info(rays_straight*,ray*,int,int,
          double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
          int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
          double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
  int get_ref_bord_num();
  int get_hit_bord_num();
};

// ¿àÞåÞÖÔÕÝØÕ ßÞ ßãâØ, ×ÐÔÐÝÝÞÜã èÐÑÛÞÝÞÜ, áâàÞÓÞ ÔÐÛìèÕ, ÝÕÖÕÛØ íâÞâ èÐÑÛÞÝclass hit_through_border_info:public general_hit_info{
protected:
  rays_straight * pr;
  ray * prTemplate;
  int iReflectBorder;
  geophys_params gp;
  raytr_params rp;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);

  virtual void
  get_ray_parameters (const point3d &pWhere,const crossing&c,double &dFi,
              double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);

  hit_through_border_info(rays_straight *,ray *,
              double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
              int i3=PM_MAX_ITERATIONS_COUNT_D,
              double d2=PM_ALPHA_STOP_D/10.,
              int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D);
};

#endif /* p_straight_elements.H */

/*  Melnikov George Yu.  2-Feb-2000  */ 
/*          ...          */

class hit_through_border_info:public general_hit_info{
protected:
  rays_straight * pr;
  ray * prTemplate;
  int iReflectBorder;
  geophys_params gp;
  raytr_params rp;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);

  virtual void
  get_ray_parameters (const point3d &pWhere,const crossing&c,double &dFi,
              double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);

  hit_through_border_info(rays_straight *,ray *,
              double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
              int i3=PM_MAX_ITERATIONS_COUNT_D,
              double d2=PM_ALPHA_STOP_D/10.,
              int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D);
};

#endif /* p_straight_elements.H */

/*  Melnikov George Yu.  2-Feb-2000  */ 
/*          ...          */