p_straight_model.H

/* p_straight_model.H */
/* $Id: p_straight_model.H,v 1.8 2005/07/18 07:16:08 nikita Exp $ */
#ifndef __p_straight_model_H
#define __p_straight_model_H


#include <mod3d/p_raytr_elements.H>
#include <mod3d/p_raytr_model.H>

#include "p_straight_elements.H"
#include "rays_straight.H"

#define PM_OUT_OF_DATE  -6
class pl_straight_model:public pl_raytr_model{
protected:  

  // ¸ÝäÞàÜÐæØï ÔÛï ÑÞÛÕÕ ÑëáâàÞÓÞ ßÞÛãçÕÝØï ßÕàÒëå ßàØÑÛØÖÕÝØÙ ßàØ ÝÐáâàÕÛÕ  // ¾ÑéÐï
  int iBFBLast;       // ½ÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ Ò ßÞáÛÕÔÝØÙ àÐ× ßÞßÐÛØ Ò ßàØÕÜÝØÚ

  // ´Ûï ßÐÔÐîéÕÙ ÒÞÛÝë...
  rays raStraightLast;
  crossings_ptrs caCrossingsStraight;

  // ´Ûï ÞâàÐÖÕÝÚØ...
  IntAr iaReflBordersLast;
  rays raReflsLast;
  crossings_ptrs caCrossingsRefl;

  // ÄÛÐÓØ ÝÐçÐÛìÝÞÙ ÝÐáâàÕÛïÝÝÞáâØ [8()]
  int iDWShotFlag;  // µáÛØ ßÞÛÞÖØâÕÛìÝëÙ - ÝÐáâàÕÛïÝë, ÕáÛØ ÝÕâ, âÞ ÝÕâ
  // ½ÞÛì - ÒáÕ ÝÐáâàÕÛïÝë, çØáÛÞ - ÞÔÝÐ ÓàÐÝØæÐ, Þâà. - ÝØ äØÓÐ ÝÕ ÝÐáâàÕÛïÝÞ
  int iUWShotFlag;  

  pl_straight_bodies abNewBodies;
  pl_straight_bodies abNewStraightBodies; // ÁßØáÞÚ ÝÞÒëå âÕÛ ÝÐ ÞÑàÐÑÞâÚã

  // ¿ÞÛãçØâì ÞÔÝÐÖÔë ÞâàÐÖÕÝÝëÙ Ûãç ØáâÞçÝØÚ-ßàØÕÜÝØÚ á ÞâàÐÖÕÝØÕÜ
  int get_reflected_ray(int,int,ray_parameters,double dAcc);
  
  void cover_possible_borders (int iReceiver,rays&,
                   IntAr&,double dAcc) const;
  void cover_one_border(int iReceiver,int iBR,ray*&rRes,double dAcc)const;
  // ÁÞ×ÔÐÝØÕ âàÐááë ßÞ ÝÞÜÕàã ØáâÞçÝØÚÐ Ø ßÐàÐÜÕâàÐÜ áÕÙáÜÞÓàÐÜÜë
  StorageID create_trace (int iReceiver, int iSize,double dDt,double *dComp[3],
              double* daSignal,int iSigSize,double dAcc,
              int iWFlags);

  // ÁÞ×ÔÐÕâ àÐ×ÛØçÝëÕ ÜÞÔÕÛìÝëÕ áÕÙáÜÞÓàÐÜÜë Ò ØáâÞçÝØÚÕ,
  // ÒÞ×ÒàÐéÐÕâ ÚÞÔ ÞèØÑÚØ
  StorageID create_iis(int iReciever, meta_iis & miRes,double dAcc,int iWFlags);

  // ÁÞ×ÔÐÝØÕ âàÐááë ßÞ ßÞáçØâÐÝÝÞÙ ÜÕâÐ ØØá Ø ßÐàÐÜÕâàÐÜ áÕÙáÜÞÓàÐÜÜë
  int form_trace (int iSize,double dDt,const meta_iis & miRes,
          double *dComp[3],double *daSignal,int iSigSize);

  // ²ëçØáÛïÕâ ÒàÕÜï ßÕàÒÞÓÞ ÒáâãßÛÕÝØï  Ò ïáÝÞÚÐÚÞÜ ßàØÕÜÝØÚÕ
  void calculate_first_break (int iN1,int iN2,const point& pReceiver,
                  double& dTime,double dAcc);

  // ´Ûï ØáßÞÛì×ÞÒÐÝØï Ò ÚÞÝáâàãÚâÞàÕ ßÞ ÑÐ×Õ ÔÐÝÝëå
  virtual p_pl_intern_border new_border(p_pl_intern_spline,border_flag);
  virtual p_pl_intern_body new_body(geophys_params&);

  void reflect_ray_from_border_to_border (ray*,int J,int K);
  void prepare_reflected_rays ();
  void push_body (int ib);
  int GetNumHitBorder(int j, const IntAr& iaBordersToExec, 
     const rays& arRaysToExec, const IntAr& iaCrossingsToExec); // fern  
  int GetNumClosestRay(const rays& arRaysToExec,const IntAr& iaCrossingsToExec,
               int iNewBorder);  // fern
  Logic extend_parallel (point pOnBorder, border_hit_info* pbhi);
  void extend_ray_to_border (ray* prToPass, int iNewBorder, int iBodyN,
                 crossing*cStart,
                 border_hit_info* pbhi=NULL);
  void pass_body (p_pl_straight_body pBCurrent);
  void init_shot_point (const point& pShotPoint,int ik);
  void init_straight_tracing_process (const point &pShotPoint,int ik);
  void init_reflected_tracing_process (int ik);
  rays* get_straight_rays(int iReceiver,double dAcc);

  void try_to_pass_big_angle();
public:
  pl_straight_model (const dbItem &dbiWhere);
  pl_straight_model (pl_model& am);
  pl_straight_model ();
  ~pl_straight_model ();

  // ²ÞÞÑçÕ âÞ íâÞ ÔÛï ÒÝãâàÕÝÝÕÓÞ ãßÞâàÕÑÛÕÝØï...
  Logic store_ray_and_node (int iBodN,ray* pr, int iNodeN);

  // ÁÞ×ÔÐÝØÕ âàÐááë ßÞ ÒÞÕÝÝëÜ ÝÞÜÕàÐÜ ØáâÞçÝØÚÐ Ø ßÐàÐÜÕâàÐÜ áÕÙáÜÞÓàÐÜÜë
  StorageID create_trace (int N1,int N2, int iSize,double dDt,double *dComp[3],
              double* daSignal,int iSigSize,double dAcc,
              int iWFlags);
  // ¼ÞÔÕÛØàÞÒÐÝØÕ ÑÕ× áÞ×ÔÐÝØï áÕÙáÜØçÕáÚÞÙ âàÐááë
  StorageID straight_modeling (int N1,int N2,double dAcc);

  // ¾ßàÕÔÕÛÕÝØÕ ßÕàÒÞÓÞ ÒáâãßÛÕÝØï á àÐ×ÝëÜØ ßàØÕÜÝØÚÐÜØ
  void  first_break_time (int iN1,int iN2,const point& pReceiver,
              double& dTime,double dAcc);
  void first_break_time(int iN1,int iN2,double& dTime);
  void uw_break_time(int iN1,int iN2,point&pR,int iB,double& dTime,double dA);
  void trace_reflected(int iUWKind=PM_P);
  void trace_reflected(int iBorder, int iUWKind);
  void trace_straight(const point &pShotPoint,int iDWKind);
  Logic register_receiver(int N1,int N2,point pRP);
  Logic straight_rays_shot()const;
  Logic reflected_rays_shot(int iBorder=0)const;

  virtual const char* class_id(){return "pl_straight_model";};

  virtual Logic isDerivedFrom(const char* class_name);
};

#endif /* p_straight_model.H */
/*  Melnikov George Yu. 31-Jan-2000  */ 

  // ¾ÑéÐï
  int iBFBLast;       // ½ÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ Ò ßÞáÛÕÔÝØÙ àÐ× ßÞßÐÛØ Ò ßàØÕÜÝØÚ

  // ´Ûï ßÐÔÐîéÕÙ ÒÞÛÝë...
  rays raStraightLast;
  crossings_ptrs caCrossingsStraight;

  // ´Ûï ÞâàÐÖÕÝÚØ...
  IntAr iaReflBordersLast;
  rays raReflsLast;
  crossings_ptrs caCrossingsRefl;

  // ÄÛÐÓØ ÝÐçÐÛìÝÞÙ ÝÐáâàÕÛïÝÝÞáâØ [8()]
  int iDWShotFlag;  // µáÛØ ßÞÛÞÖØâÕÛìÝëÙ - ÝÐáâàÕÛïÝë, ÕáÛØ ÝÕâ, âÞ ÝÕâ  // ½ÞÛì - ÒáÕ ÝÐáâàÕÛïÝë, çØáÛÞ - ÞÔÝÐ ÓàÐÝØæÐ, Þâà. - ÝØ äØÓÐ ÝÕ ÝÐáâàÕÛïÝÞ
  int iUWShotFlag;  

  pl_straight_bodies abNewBodies;
  pl_straight_bodies abNewStraightBodies; // ÁßØáÞÚ ÝÞÒëå âÕÛ ÝÐ ÞÑàÐÑÞâÚã

  // ¿ÞÛãçØâì ÞÔÝÐÖÔë ÞâàÐÖÕÝÝëÙ Ûãç ØáâÞçÝØÚ-ßàØÕÜÝØÚ á ÞâàÐÖÕÝØÕÜ
  int get_reflected_ray(int,int,ray_parameters,double dAcc);
  
  void cover_possible_borders (int iReceiver,rays&,
                   IntAr&,double dAcc) const;
  void cover_one_border(int iReceiver,int iBR,ray*&rRes,double dAcc)const;
  // ÁÞ×ÔÐÝØÕ âàÐááë ßÞ ÝÞÜÕàã ØáâÞçÝØÚÐ Ø ßÐàÐÜÕâàÐÜ áÕÙáÜÞÓàÐÜÜë
  StorageID create_trace (int iReceiver, int iSize,double dDt,double *dComp[3],
              double* daSignal,int iSigSize,double dAcc,
              int iWFlags);

  // ÁÞ×ÔÐÕâ àÐ×ÛØçÝëÕ ÜÞÔÕÛìÝëÕ áÕÙáÜÞÓàÐÜÜë Ò ØáâÞçÝØÚÕ,
  // ÒÞ×ÒàÐéÐÕâ ÚÞÔ ÞèØÑÚØ
  StorageID create_iis(int iReciever, meta_iis & miRes,double dAcc,int iWFlags);

  // ÁÞ×ÔÐÝØÕ âàÐááë ßÞ ßÞáçØâÐÝÝÞÙ ÜÕâÐ ØØá Ø ßÐàÐÜÕâàÐÜ áÕÙáÜÞÓàÐÜÜë
  int form_trace (int iSize,double dDt,const meta_iis & miRes,
          double *dComp[3],double *daSignal,int iSigSize);

  // ²ëçØáÛïÕâ ÒàÕÜï ßÕàÒÞÓÞ ÒáâãßÛÕÝØï  Ò ïáÝÞÚÐÚÞÜ ßàØÕÜÝØÚÕ
  void calculate_first_break (int iN1,int iN2,const point& pReceiver,
                  double& dTime,double dAcc);

  // ´Ûï ØáßÞÛì×ÞÒÐÝØï Ò ÚÞÝáâàãÚâÞàÕ ßÞ ÑÐ×Õ ÔÐÝÝëå
  virtual p_pl_intern_border new_border(p_pl_intern_spline,border_flag);
  virtual p_pl_intern_body new_body(geophys_params&);

  void reflect_ray_from_border_to_border (ray*,int J,int K);
  void prepare_reflected_rays ();
  void push_body (int ib);
  int GetNumHitBorder(int j, const IntAr& iaBordersToExec, 
     const rays& arRaysToExec, const IntAr& iaCrossingsToExec); // fern  
  int GetNumClosestRay(const rays& arRaysToExec,const IntAr& iaCrossingsToExec,
               int iNewBorder);  // fern
  Logic extend_parallel (point pOnBorder, border_hit_info* pbhi);
  void extend_ray_to_border (ray* prToPass, int iNewBorder, int iBodyN,
                 crossing*cStart,
                 border_hit_info* pbhi=NULL);
  void pass_body (p_pl_straight_body pBCurrent);
  void init_shot_point (const point& pShotPoint,int ik);
  void init_straight_tracing_process (const point &pShotPoint,int ik);
  void init_reflected_tracing_process (int ik);
  rays* get_straight_rays(int iReceiver,double dAcc);

  void try_to_pass_big_angle();
public:
  pl_straight_model (const dbItem &dbiWhere);
  pl_straight_model (pl_model& am);
  pl_straight_model ();
  ~pl_straight_model ();

  // ²ÞÞÑçÕ âÞ íâÞ ÔÛï ÒÝãâàÕÝÝÕÓÞ ãßÞâàÕÑÛÕÝØï...
  Logic store_ray_and_node (int iBodN,ray* pr, int iNodeN);

  // ÁÞ×ÔÐÝØÕ âàÐááë ßÞ ÒÞÕÝÝëÜ ÝÞÜÕàÐÜ ØáâÞçÝØÚÐ Ø ßÐàÐÜÕâàÐÜ áÕÙáÜÞÓàÐÜÜë
  StorageID create_trace (int N1,int N2, int iSize,double dDt,double *dComp[3],
              double* daSignal,int iSigSize,double dAcc,
              int iWFlags);
  // ¼ÞÔÕÛØàÞÒÐÝØÕ ÑÕ× áÞ×ÔÐÝØï áÕÙáÜØçÕáÚÞÙ âàÐááë
  StorageID straight_modeling (int N1,int N2,double dAcc);

  // ¾ßàÕÔÕÛÕÝØÕ ßÕàÒÞÓÞ ÒáâãßÛÕÝØï á àÐ×ÝëÜØ ßàØÕÜÝØÚÐÜØ
  void  first_break_time (int iN1,int iN2,const point& pReceiver,
              double& dTime,double dAcc);
  void first_break_time(int iN1,int iN2,double& dTime);
  void uw_break_time(int iN1,int iN2,point&pR,int iB,double& dTime,double dA);
  void trace_reflected(int iUWKind=PM_P);
  void trace_reflected(int iBorder, int iUWKind);
  void trace_straight(const point &pShotPoint,int iDWKind);
  Logic register_receiver(int N1,int N2,point pRP);
  Logic straight_rays_shot()const;
  Logic reflected_rays_shot(int iBorder=0)const;

  virtual const char* class_id(){return "pl_straight_model";};

  virtual Logic isDerivedFrom(const char* class_name);
};

#endif /* p_straight_model.H */
/*  Melnikov George Yu. 31-Jan-2000  */ 

  // ½ÞÛì - ÒáÕ ÝÐáâàÕÛïÝë, çØáÛÞ - ÞÔÝÐ ÓàÐÝØæÐ, Þâà. - ÝØ äØÓÐ ÝÕ ÝÐáâàÕÛïÝÞ
  int iUWShotFlag;  

  pl_straight_bodies abNewBodies;
  pl_straight_bodies abNewStraightBodies; // ÁßØáÞÚ ÝÞÒëå âÕÛ ÝÐ ÞÑàÐÑÞâÚã

  // ¿ÞÛãçØâì ÞÔÝÐÖÔë ÞâàÐÖÕÝÝëÙ Ûãç ØáâÞçÝØÚ-ßàØÕÜÝØÚ á ÞâàÐÖÕÝØÕÜ
  int get_reflected_ray(int,int,ray_parameters,double dAcc);
  
  void cover_possible_borders (int iReceiver,rays&,
                   IntAr&,double dAcc) const;
  void cover_one_border(int iReceiver,int iBR,ray*&rRes,double dAcc)const;
  // ÁÞ×ÔÐÝØÕ âàÐááë ßÞ ÝÞÜÕàã ØáâÞçÝØÚÐ Ø ßÐàÐÜÕâàÐÜ áÕÙáÜÞÓàÐÜÜë
  StorageID create_trace (int iReceiver, int iSize,double dDt,double *dComp[3],
              double* daSignal,int iSigSize,double dAcc,
              int iWFlags);

  // ÁÞ×ÔÐÕâ àÐ×ÛØçÝëÕ ÜÞÔÕÛìÝëÕ áÕÙáÜÞÓàÐÜÜë Ò ØáâÞçÝØÚÕ,
  // ÒÞ×ÒàÐéÐÕâ ÚÞÔ ÞèØÑÚØ
  StorageID create_iis(int iReciever, meta_iis & miRes,double dAcc,int iWFlags);

  // ÁÞ×ÔÐÝØÕ âàÐááë ßÞ ßÞáçØâÐÝÝÞÙ ÜÕâÐ ØØá Ø ßÐàÐÜÕâàÐÜ áÕÙáÜÞÓàÐÜÜë
  int form_trace (int iSize,double dDt,const meta_iis & miRes,
          double *dComp[3],double *daSignal,int iSigSize);

  // ²ëçØáÛïÕâ ÒàÕÜï ßÕàÒÞÓÞ ÒáâãßÛÕÝØï  Ò ïáÝÞÚÐÚÞÜ ßàØÕÜÝØÚÕ  void calculate_first_break (int iN1,int iN2,const point& pReceiver,
                  double& dTime,double dAcc);

  // ´Ûï ØáßÞÛì×ÞÒÐÝØï Ò ÚÞÝáâàãÚâÞàÕ ßÞ ÑÐ×Õ ÔÐÝÝëå
  virtual p_pl_intern_border new_border(p_pl_intern_spline,border_flag);
  virtual p_pl_intern_body new_body(geophys_params&);

  void reflect_ray_from_border_to_border (ray*,int J,int K);
  void prepare_reflected_rays ();
  void push_body (int ib);
  int GetNumHitBorder(int j, const IntAr& iaBordersToExec, 
     const rays& arRaysToExec, const IntAr& iaCrossingsToExec); // fern  
  int GetNumClosestRay(const rays& arRaysToExec,const IntAr& iaCrossingsToExec,
               int iNewBorder);  // fern
  Logic extend_parallel (point pOnBorder, border_hit_info* pbhi);
  void extend_ray_to_border (ray* prToPass, int iNewBorder, int iBodyN,
                 crossing*cStart,
                 border_hit_info* pbhi=NULL);
  void pass_body (p_pl_straight_body pBCurrent);
  void init_shot_point (const point& pShotPoint,int ik);
  void init_straight_tracing_process (const point &pShotPoint,int ik);
  void init_reflected_tracing_process (int ik);
  rays* get_straight_rays(int iReceiver,double dAcc);

  void try_to_pass_big_angle();
public:
  pl_straight_model (const dbItem &dbiWhere);
  pl_straight_model (pl_model& am);
  pl_straight_model ();
  ~pl_straight_model ();

  // ²ÞÞÑçÕ âÞ íâÞ ÔÛï ÒÝãâàÕÝÝÕÓÞ ãßÞâàÕÑÛÕÝØï...
  Logic store_ray_and_node (int iBodN,ray* pr, int iNodeN);

  // ÁÞ×ÔÐÝØÕ âàÐááë ßÞ ÒÞÕÝÝëÜ ÝÞÜÕàÐÜ ØáâÞçÝØÚÐ Ø ßÐàÐÜÕâàÐÜ áÕÙáÜÞÓàÐÜÜë
  StorageID create_trace (int N1,int N2, int iSize,double dDt,double *dComp[3],
              double* daSignal,int iSigSize,double dAcc,
              int iWFlags);
  // ¼ÞÔÕÛØàÞÒÐÝØÕ ÑÕ× áÞ×ÔÐÝØï áÕÙáÜØçÕáÚÞÙ âàÐááë
  StorageID straight_modeling (int N1,int N2,double dAcc);

  // ¾ßàÕÔÕÛÕÝØÕ ßÕàÒÞÓÞ ÒáâãßÛÕÝØï á àÐ×ÝëÜØ ßàØÕÜÝØÚÐÜØ
  void  first_break_time (int iN1,int iN2,const point& pReceiver,
              double& dTime,double dAcc);
  void first_break_time(int iN1,int iN2,double& dTime);
  void uw_break_time(int iN1,int iN2,point&pR,int iB,double& dTime,double dA);
  void trace_reflected(int iUWKind=PM_P);
  void trace_reflected(int iBorder, int iUWKind);
  void trace_straight(const point &pShotPoint,int iDWKind);
  Logic register_receiver(int N1,int N2,point pRP);
  Logic straight_rays_shot()const;
  Logic reflected_rays_shot(int iBorder=0)const;

  virtual const char* class_id(){return "pl_straight_model";};

  virtual Logic isDerivedFrom(const char* class_name);
};

#endif /* p_straight_model.H */
/*  Melnikov George Yu. 31-Jan-2000  */ 

  void calculate_first_break (int iN1,int iN2,const point& pReceiver,
                  double& dTime,double dAcc);

  // ´Ûï ØáßÞÛì×ÞÒÐÝØï Ò ÚÞÝáâàãÚâÞàÕ ßÞ ÑÐ×Õ ÔÐÝÝëå  virtual p_pl_intern_border new_border(p_pl_intern_spline,border_flag);
  virtual p_pl_intern_body new_body(geophys_params&);

  void reflect_ray_from_border_to_border (ray*,int J,int K);
  void prepare_reflected_rays ();
  void push_body (int ib);
  int GetNumHitBorder(int j, const IntAr& iaBordersToExec, 
     const rays& arRaysToExec, const IntAr& iaCrossingsToExec); // fern  
  int GetNumClosestRay(const rays& arRaysToExec,const IntAr& iaCrossingsToExec,
               int iNewBorder);  // fern
  Logic extend_parallel (point pOnBorder, border_hit_info* pbhi);
  void extend_ray_to_border (ray* prToPass, int iNewBorder, int iBodyN,
                 crossing*cStart,
                 border_hit_info* pbhi=NULL);
  void pass_body (p_pl_straight_body pBCurrent);
  void init_shot_point (const point& pShotPoint,int ik);
  void init_straight_tracing_process (const point &pShotPoint,int ik);
  void init_reflected_tracing_process (int ik);
  rays* get_straight_rays(int iReceiver,double dAcc);

  void try_to_pass_big_angle();
public:
  pl_straight_model (const dbItem &dbiWhere);
  pl_straight_model (pl_model& am);
  pl_straight_model ();
  ~pl_straight_model ();

  // ²ÞÞÑçÕ âÞ íâÞ ÔÛï ÒÝãâàÕÝÝÕÓÞ ãßÞâàÕÑÛÕÝØï...
  Logic store_ray_and_node (int iBodN,ray* pr, int iNodeN);

  // ÁÞ×ÔÐÝØÕ âàÐááë ßÞ ÒÞÕÝÝëÜ ÝÞÜÕàÐÜ ØáâÞçÝØÚÐ Ø ßÐàÐÜÕâàÐÜ áÕÙáÜÞÓàÐÜÜë
  StorageID create_trace (int N1,int N2, int iSize,double dDt,double *dComp[3],
              double* daSignal,int iSigSize,double dAcc,
              int iWFlags);
  // ¼ÞÔÕÛØàÞÒÐÝØÕ ÑÕ× áÞ×ÔÐÝØï áÕÙáÜØçÕáÚÞÙ âàÐááë
  StorageID straight_modeling (int N1,int N2,double dAcc);

  // ¾ßàÕÔÕÛÕÝØÕ ßÕàÒÞÓÞ ÒáâãßÛÕÝØï á àÐ×ÝëÜØ ßàØÕÜÝØÚÐÜØ
  void  first_break_time (int iN1,int iN2,const point& pReceiver,
              double& dTime,double dAcc);
  void first_break_time(int iN1,int iN2,double& dTime);
  void uw_break_time(int iN1,int iN2,point&pR,int iB,double& dTime,double dA);
  void trace_reflected(int iUWKind=PM_P);
  void trace_reflected(int iBorder, int iUWKind);
  void trace_straight(const point &pShotPoint,int iDWKind);
  Logic register_receiver(int N1,int N2,point pRP);
  Logic straight_rays_shot()const;
  Logic reflected_rays_shot(int iBorder=0)const;

  virtual const char* class_id(){return "pl_straight_model";};

  virtual Logic isDerivedFrom(const char* class_name);
};

#endif /* p_straight_model.H */
/*  Melnikov George Yu. 31-Jan-2000  */ 

  virtual p_pl_intern_border new_border(p_pl_intern_spline,border_flag);
  virtual p_pl_intern_body new_body(geophys_params&);

  void reflect_ray_from_border_to_border (ray*,int J,int K);
  void prepare_reflected_rays ();
  void push_body (int ib);
  int GetNumHitBorder(int j, const IntAr& iaBordersToExec, 
     const rays& arRaysToExec, const IntAr& iaCrossingsToExec); // fern  
  int GetNumClosestRay(const rays& arRaysToExec,const IntAr& iaCrossingsToExec,
               int iNewBorder);  // fern
  Logic extend_parallel (point pOnBorder, border_hit_info* pbhi);
  void extend_ray_to_border (ray* prToPass, int iNewBorder, int iBodyN,
                 crossing*cStart,
                 border_hit_info* pbhi=NULL);
  void pass_body (p_pl_straight_body pBCurrent);
  void init_shot_point (const point& pShotPoint,int ik);
  void init_straight_tracing_process (const point &pShotPoint,int ik);
  void init_reflected_tracing_process (int ik);
  rays* get_straight_rays(int iReceiver,double dAcc);

  void try_to_pass_big_angle();
public:
  pl_straight_model (const dbItem &dbiWhere);
  pl_straight_model (pl_model& am);
  pl_straight_model ();
  ~pl_straight_model ();

  // ²ÞÞÑçÕ âÞ íâÞ ÔÛï ÒÝãâàÕÝÝÕÓÞ ãßÞâàÕÑÛÕÝØï...
  Logic store_ray_and_node (int iBodN,ray* pr, int iNodeN);

  // ÁÞ×ÔÐÝØÕ âàÐááë ßÞ ÒÞÕÝÝëÜ ÝÞÜÕàÐÜ ØáâÞçÝØÚÐ Ø ßÐàÐÜÕâàÐÜ áÕÙáÜÞÓàÐÜÜë
  StorageID create_trace (int N1,int N2, int iSize,double dDt,double *dComp[3],
              double* daSignal,int iSigSize,double dAcc,
              int iWFlags);
  // ¼ÞÔÕÛØàÞÒÐÝØÕ ÑÕ× áÞ×ÔÐÝØï áÕÙáÜØçÕáÚÞÙ âàÐááë
  StorageID straight_modeling (int N1,int N2,double dAcc);

  // ¾ßàÕÔÕÛÕÝØÕ ßÕàÒÞÓÞ ÒáâãßÛÕÝØï á àÐ×ÝëÜØ ßàØÕÜÝØÚÐÜØ  void  first_break_time (int iN1,int iN2,const point& pReceiver,
              double& dTime,double dAcc);
  void first_break_time(int iN1,int iN2,double& dTime);
  void uw_break_time(int iN1,int iN2,point&pR,int iB,double& dTime,double dA);
  void trace_reflected(int iUWKind=PM_P);
  void trace_reflected(int iBorder, int iUWKind);
  void trace_straight(const point &pShotPoint,int iDWKind);
  Logic register_receiver(int N1,int N2,point pRP);
  Logic straight_rays_shot()const;
  Logic reflected_rays_shot(int iBorder=0)const;

  virtual const char* class_id(){return "pl_straight_model";};

  virtual Logic isDerivedFrom(const char* class_name);
};

#endif /* p_straight_model.H */
/*  Melnikov George Yu. 31-Jan-2000  */ 

  void  first_break_time (int iN1,int iN2,const point& pReceiver,
              double& dTime,double dAcc);
  void first_break_time(int iN1,int iN2,double& dTime);
  void uw_break_time(int iN1,int iN2,point&pR,int iB,double& dTime,double dA);
  void trace_reflected(int iUWKind=PM_P);
  void trace_reflected(int iBorder, int iUWKind);
  void trace_straight(const point &pShotPoint,int iDWKind);
  Logic register_receiver(int N1,int N2,point pRP);
  Logic straight_rays_shot()const;
  Logic reflected_rays_shot(int iBorder=0)const;

  virtual const char* class_id(){return "pl_straight_model";};

  virtual Logic isDerivedFrom(const char* class_name);
};

#endif /* p_straight_model.H */
/*  Melnikov George Yu. 31-Jan-2000  */