p_vspogt_elements.H

/* p_vspogt_elements.H */
/* $Id: p_vspogt_elements.H,v 2.6 2001/09/25 17:58:05 hoh Exp $ */
#ifndef __p_vspogt_elements_H
#define __p_vspogt_elements_H

#include <mix/IntAr.H>
#include <mix/IndexAr.H>
#include <mod3d/p_raytr_model.H>
#include <mod3d/p_raytr_elements.H>
#include <mod3d/rays.H>

// ´ÐÝÝëÕ ÔÛï ßàÞÔÞÛÖÕÝØï ÛãçÐ ÒÝãâàØ âÕÛÐstruct extend_data{
  double vray[3];
  double dTime;
  point pStart;
};

class time_depth_list_local{
protected:  
  Logic lPrepared;
  IndexAr ia;
  DoubleAr daTimes;

  DoubleAr daFi;
  DoubleAr daAz;

  DoubleAr daProp;
  DoubleAr daDiver;

  DoubleAr daRefRat;

  DoubleAr davpolx;
  DoubleAr davpoly;
  DoubleAr davpolz;

  aPoints paReflections;
  p_pl_splines paReflectingSplines;

  // ¿àÞØÝâÕàßÞÛØàÞÒÐÛØ âÞçÚã ßÞ áÞáÕÔïÜ
  int interpolate (int iPoint,point&pRes) const;

  // ¸ÝÞÓÔÐ ßàØåÞÔØâáï Ø ãÔÐÛïâì ÛØèÝØÕ íÛÕÜÕÝâë
  void remove(int i);
public:

  // ½Ð áÛãçÐÙ, ÕáÛØ áÞ×ÔÐÕâáï ÛÞÚÐÛìÝëÙ áßØáÞÚ, ÚÞâÞàëÙ ÝÕ ÑãÔÕÜ Ò ±´ ßØáÐâì
  time_depth_list_local(unsigned quant = DEFAULT_QUANT,
            unsigned volume = START_VOLUME);
  
  // ¿Þ ÒàÕÜÕÝØ - ÛØÝÕÙÝÐï ØÝâÕàßÞÛïæØï ßÐàÐÜÕâàÞÒ
  void interpolate(double dTime,double& prop,double& diver,double& refrat,
           point&pRefl,pl_spline& sReflect) const; 
  void interpolate(double dTime, double* vpol) const;
  double start_time()const{return daTimes(ia(0));};
  double finish_time()const{return daTimes(ia(_last_));};


  // ²àÕÜï áÕàñÔÚØ áÕÓÜÕÝâÐ á ßÕàÕÔÐÒÐÕÜëÜ ÝÞÜÕàÞÜ
  double t_of_new_border(int) const;

  // ÂØßÐ ÒÝÕáÛØ ÝÐçÐÛìÝëÕ ÞâàÐÖÕÝØï ØâÐÔÐÛÕÕ
  void prepared(){lPrepared=TRUE;};

  int n_of_segments() const; // ºÞÛØçÕáâÒÞ áÕÓÜÕÝâÞÒ Ò áßØáÚÕ

  // ÔÞÑÐÒØâì ÒàÕÜï-âÞçÚã-áßÛÐÙÝ
  void add_new_node(double, const point3d &,
            p_pl_intern_spline,ray_to_receiver &,Logic); 
  void add_new_node(double, const point3d &,
            p_pl_intern_spline,ray_to_receiver &,double); 

  void merge(const time_depth_list_local&); // ¿àØÑÐÒØâì Ú ÞÔÝÞÜã áßØáÚã ÔàãÓÞÙ

  // ¿Þ ÒàÕÜÕÝØ - ÝÞÜÕà ÞâàÕ×ÚÐ
  int t2i(double)const;

  // ÁÛØïÝØÕ ÔÒãå âØáßÐáßØáÚÞÒ ÑÕ× ßÞÒâÞàÞÒ
  void merge_unique(const time_depth_list_local&); 

  // ÂÞ ÖÕ, ÝÞ ÔÞÑÐÒÛïÕâ ÝÞÒëÕ íÛÕÜÕÝâë ØáÚÛîçØâÕÛìÝÞ ßÞ ÚàÐïÜ
  void merge_unique_monotone(const time_depth_list_local&); 

  // ÀÐ×ÝÞÓÞ àÞÔÐ çØáâÚØ  
  void clear();
  void pack();

  // ÀÐ×ÛØçÝëÕ ÞâÛÐÔÞçÝëÕ äØçÚØ
  void debug_output() const;

  friend class time_depth_lists;
};

class time_depth_list: public db_container,public time_depth_list_local{
protected:
  dbHandle dbhReceivers;
  dbHandle dbhTD;
  static int iLast;
  int N1,N2;
  int iTDLStoreRegime;

  // ÀÐÑÞâÐ á ±´
  int read_from_db (dbItem& dbiTD, dbItem& dbiRN);
  int write_to_db (dbItem& dbiTD, dbItem& dbiRN)const;
  void get_item (dbItem& dbiRes, dbItem& dbiRN) const;
public:
  // µáÛØ âÐÚØÜ ÞÑàÐ×ÞÜ áÞ×ÔÐÝ - ÜÞÖÕÜ ßØáÐâì Ò ±´
  time_depth_list(db_results_manager& d,char*sp, int id,int n1,int n2,
          int iTDLSR=PM_ADD,
          unsigned quant = DEFAULT_QUANT,unsigned volume = START_VOLUME);
  // ±´
  virtual int store() const;
  virtual int load();
  void prepare_pack_load (dbHandle&dbhR,dbHandle&dbhT);
  void done_loading();
  void pack_load(dbHandle&,dbHandle&);

  friend class time_depth_lists;
};

typedef time_depth_list* p_time_depth_list;
typedef ObjectsAr<p_time_depth_list> time_depth_lists_dynar;

class time_depth_lists:public time_depth_lists_dynar{
public:
  p_time_depth_list get(int,int)const;
};

class pl_vspogt_body:public pl_intern_body,public time_depth_list_local{
protected:
  pl_intern_borders                                       paActiveBorders;  
  DoubleAr daActualTimes;           // ²àÕÜÕÝÐ, ÝÐ ÚÞâÞàëå ÞâàÐÖÐÛØáì ^ ÓàÐÝØæë

  Logic interpolated(int) const;
  p_pl_fict_border get_fict_border(double)const; // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  extend_data *ed;
public:
  pl_vspogt_body(geophys_params&g);
  ~pl_vspogt_body(){paActiveBorders.clean();delete ed;};


  void prepare_for_vspogt(const pl_intern_borders&,rays_to_receiver* prtr,
              const DoubleAr& daRT);
  void get_initial_fict_borders(pl_fict_borders &)const;

  // ¿ÐàÐ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ÝÞÜÕàã ßÛÞåÞØÝâÕàßÞÛØàÞÒÐÝÝÞÓÞ ¶8() ÒàÕÜÕÝØ
  pl_fict_borders* get_fict_borders(int) const;
  // ½ÕÚÞâÞàÞÕ ÚÞÛØçÕáâÒÞ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ã×Ûã ib
  void add_n_borders(int iNNewPointsLeft,int ib,pl_fict_borders*,
             ray_parameters_dynar *)const;

  // ½ÞÜÕà íâÞÓÞáÐÜÞÓÞ ã×ÛÐ
  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

struct extend_data{
  double vray[3];
  double dTime;
  point pStart;
};

class time_depth_list_local{
protected:  
  Logic lPrepared;
  IndexAr ia;
  DoubleAr daTimes;

  DoubleAr daFi;
  DoubleAr daAz;

  DoubleAr daProp;
  DoubleAr daDiver;

  DoubleAr daRefRat;

  DoubleAr davpolx;
  DoubleAr davpoly;
  DoubleAr davpolz;

  aPoints paReflections;
  p_pl_splines paReflectingSplines;

  // ¿àÞØÝâÕàßÞÛØàÞÒÐÛØ âÞçÚã ßÞ áÞáÕÔïÜ  int interpolate (int iPoint,point&pRes) const;

  // ¸ÝÞÓÔÐ ßàØåÞÔØâáï Ø ãÔÐÛïâì ÛØèÝØÕ íÛÕÜÕÝâë
  void remove(int i);
public:

  // ½Ð áÛãçÐÙ, ÕáÛØ áÞ×ÔÐÕâáï ÛÞÚÐÛìÝëÙ áßØáÞÚ, ÚÞâÞàëÙ ÝÕ ÑãÔÕÜ Ò ±´ ßØáÐâì
  time_depth_list_local(unsigned quant = DEFAULT_QUANT,
            unsigned volume = START_VOLUME);
  
  // ¿Þ ÒàÕÜÕÝØ - ÛØÝÕÙÝÐï ØÝâÕàßÞÛïæØï ßÐàÐÜÕâàÞÒ
  void interpolate(double dTime,double& prop,double& diver,double& refrat,
           point&pRefl,pl_spline& sReflect) const; 
  void interpolate(double dTime, double* vpol) const;
  double start_time()const{return daTimes(ia(0));};
  double finish_time()const{return daTimes(ia(_last_));};


  // ²àÕÜï áÕàñÔÚØ áÕÓÜÕÝâÐ á ßÕàÕÔÐÒÐÕÜëÜ ÝÞÜÕàÞÜ
  double t_of_new_border(int) const;

  // ÂØßÐ ÒÝÕáÛØ ÝÐçÐÛìÝëÕ ÞâàÐÖÕÝØï ØâÐÔÐÛÕÕ
  void prepared(){lPrepared=TRUE;};

  int n_of_segments() const; // ºÞÛØçÕáâÒÞ áÕÓÜÕÝâÞÒ Ò áßØáÚÕ

  // ÔÞÑÐÒØâì ÒàÕÜï-âÞçÚã-áßÛÐÙÝ
  void add_new_node(double, const point3d &,
            p_pl_intern_spline,ray_to_receiver &,Logic); 
  void add_new_node(double, const point3d &,
            p_pl_intern_spline,ray_to_receiver &,double); 

  void merge(const time_depth_list_local&); // ¿àØÑÐÒØâì Ú ÞÔÝÞÜã áßØáÚã ÔàãÓÞÙ

  // ¿Þ ÒàÕÜÕÝØ - ÝÞÜÕà ÞâàÕ×ÚÐ
  int t2i(double)const;

  // ÁÛØïÝØÕ ÔÒãå âØáßÐáßØáÚÞÒ ÑÕ× ßÞÒâÞàÞÒ
  void merge_unique(const time_depth_list_local&); 

  // ÂÞ ÖÕ, ÝÞ ÔÞÑÐÒÛïÕâ ÝÞÒëÕ íÛÕÜÕÝâë ØáÚÛîçØâÕÛìÝÞ ßÞ ÚàÐïÜ
  void merge_unique_monotone(const time_depth_list_local&); 

  // ÀÐ×ÝÞÓÞ àÞÔÐ çØáâÚØ  
  void clear();
  void pack();

  // ÀÐ×ÛØçÝëÕ ÞâÛÐÔÞçÝëÕ äØçÚØ
  void debug_output() const;

  friend class time_depth_lists;
};

class time_depth_list: public db_container,public time_depth_list_local{
protected:
  dbHandle dbhReceivers;
  dbHandle dbhTD;
  static int iLast;
  int N1,N2;
  int iTDLStoreRegime;

  // ÀÐÑÞâÐ á ±´
  int read_from_db (dbItem& dbiTD, dbItem& dbiRN);
  int write_to_db (dbItem& dbiTD, dbItem& dbiRN)const;
  void get_item (dbItem& dbiRes, dbItem& dbiRN) const;
public:
  // µáÛØ âÐÚØÜ ÞÑàÐ×ÞÜ áÞ×ÔÐÝ - ÜÞÖÕÜ ßØáÐâì Ò ±´
  time_depth_list(db_results_manager& d,char*sp, int id,int n1,int n2,
          int iTDLSR=PM_ADD,
          unsigned quant = DEFAULT_QUANT,unsigned volume = START_VOLUME);
  // ±´
  virtual int store() const;
  virtual int load();
  void prepare_pack_load (dbHandle&dbhR,dbHandle&dbhT);
  void done_loading();
  void pack_load(dbHandle&,dbHandle&);

  friend class time_depth_lists;
};

typedef time_depth_list* p_time_depth_list;
typedef ObjectsAr<p_time_depth_list> time_depth_lists_dynar;

class time_depth_lists:public time_depth_lists_dynar{
public:
  p_time_depth_list get(int,int)const;
};

class pl_vspogt_body:public pl_intern_body,public time_depth_list_local{
protected:
  pl_intern_borders                                       paActiveBorders;  
  DoubleAr daActualTimes;           // ²àÕÜÕÝÐ, ÝÐ ÚÞâÞàëå ÞâàÐÖÐÛØáì ^ ÓàÐÝØæë

  Logic interpolated(int) const;
  p_pl_fict_border get_fict_border(double)const; // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  extend_data *ed;
public:
  pl_vspogt_body(geophys_params&g);
  ~pl_vspogt_body(){paActiveBorders.clean();delete ed;};


  void prepare_for_vspogt(const pl_intern_borders&,rays_to_receiver* prtr,
              const DoubleAr& daRT);
  void get_initial_fict_borders(pl_fict_borders &)const;

  // ¿ÐàÐ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ÝÞÜÕàã ßÛÞåÞØÝâÕàßÞÛØàÞÒÐÝÝÞÓÞ ¶8() ÒàÕÜÕÝØ
  pl_fict_borders* get_fict_borders(int) const;
  // ½ÕÚÞâÞàÞÕ ÚÞÛØçÕáâÒÞ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ã×Ûã ib
  void add_n_borders(int iNNewPointsLeft,int ib,pl_fict_borders*,
             ray_parameters_dynar *)const;

  // ½ÞÜÕà íâÞÓÞáÐÜÞÓÞ ã×ÛÐ
  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  int interpolate (int iPoint,point&pRes) const;

  // ¸ÝÞÓÔÐ ßàØåÞÔØâáï Ø ãÔÐÛïâì ÛØèÝØÕ íÛÕÜÕÝâë  void remove(int i);
public:

  // ½Ð áÛãçÐÙ, ÕáÛØ áÞ×ÔÐÕâáï ÛÞÚÐÛìÝëÙ áßØáÞÚ, ÚÞâÞàëÙ ÝÕ ÑãÔÕÜ Ò ±´ ßØáÐâì
  time_depth_list_local(unsigned quant = DEFAULT_QUANT,
            unsigned volume = START_VOLUME);
  
  // ¿Þ ÒàÕÜÕÝØ - ÛØÝÕÙÝÐï ØÝâÕàßÞÛïæØï ßÐàÐÜÕâàÞÒ
  void interpolate(double dTime,double& prop,double& diver,double& refrat,
           point&pRefl,pl_spline& sReflect) const; 
  void interpolate(double dTime, double* vpol) const;
  double start_time()const{return daTimes(ia(0));};
  double finish_time()const{return daTimes(ia(_last_));};


  // ²àÕÜï áÕàñÔÚØ áÕÓÜÕÝâÐ á ßÕàÕÔÐÒÐÕÜëÜ ÝÞÜÕàÞÜ
  double t_of_new_border(int) const;

  // ÂØßÐ ÒÝÕáÛØ ÝÐçÐÛìÝëÕ ÞâàÐÖÕÝØï ØâÐÔÐÛÕÕ
  void prepared(){lPrepared=TRUE;};

  int n_of_segments() const; // ºÞÛØçÕáâÒÞ áÕÓÜÕÝâÞÒ Ò áßØáÚÕ

  // ÔÞÑÐÒØâì ÒàÕÜï-âÞçÚã-áßÛÐÙÝ
  void add_new_node(double, const point3d &,
            p_pl_intern_spline,ray_to_receiver &,Logic); 
  void add_new_node(double, const point3d &,
            p_pl_intern_spline,ray_to_receiver &,double); 

  void merge(const time_depth_list_local&); // ¿àØÑÐÒØâì Ú ÞÔÝÞÜã áßØáÚã ÔàãÓÞÙ

  // ¿Þ ÒàÕÜÕÝØ - ÝÞÜÕà ÞâàÕ×ÚÐ
  int t2i(double)const;

  // ÁÛØïÝØÕ ÔÒãå âØáßÐáßØáÚÞÒ ÑÕ× ßÞÒâÞàÞÒ
  void merge_unique(const time_depth_list_local&); 

  // ÂÞ ÖÕ, ÝÞ ÔÞÑÐÒÛïÕâ ÝÞÒëÕ íÛÕÜÕÝâë ØáÚÛîçØâÕÛìÝÞ ßÞ ÚàÐïÜ
  void merge_unique_monotone(const time_depth_list_local&); 

  // ÀÐ×ÝÞÓÞ àÞÔÐ çØáâÚØ  
  void clear();
  void pack();

  // ÀÐ×ÛØçÝëÕ ÞâÛÐÔÞçÝëÕ äØçÚØ
  void debug_output() const;

  friend class time_depth_lists;
};

class time_depth_list: public db_container,public time_depth_list_local{
protected:
  dbHandle dbhReceivers;
  dbHandle dbhTD;
  static int iLast;
  int N1,N2;
  int iTDLStoreRegime;

  // ÀÐÑÞâÐ á ±´
  int read_from_db (dbItem& dbiTD, dbItem& dbiRN);
  int write_to_db (dbItem& dbiTD, dbItem& dbiRN)const;
  void get_item (dbItem& dbiRes, dbItem& dbiRN) const;
public:
  // µáÛØ âÐÚØÜ ÞÑàÐ×ÞÜ áÞ×ÔÐÝ - ÜÞÖÕÜ ßØáÐâì Ò ±´
  time_depth_list(db_results_manager& d,char*sp, int id,int n1,int n2,
          int iTDLSR=PM_ADD,
          unsigned quant = DEFAULT_QUANT,unsigned volume = START_VOLUME);
  // ±´
  virtual int store() const;
  virtual int load();
  void prepare_pack_load (dbHandle&dbhR,dbHandle&dbhT);
  void done_loading();
  void pack_load(dbHandle&,dbHandle&);

  friend class time_depth_lists;
};

typedef time_depth_list* p_time_depth_list;
typedef ObjectsAr<p_time_depth_list> time_depth_lists_dynar;

class time_depth_lists:public time_depth_lists_dynar{
public:
  p_time_depth_list get(int,int)const;
};

class pl_vspogt_body:public pl_intern_body,public time_depth_list_local{
protected:
  pl_intern_borders                                       paActiveBorders;  
  DoubleAr daActualTimes;           // ²àÕÜÕÝÐ, ÝÐ ÚÞâÞàëå ÞâàÐÖÐÛØáì ^ ÓàÐÝØæë

  Logic interpolated(int) const;
  p_pl_fict_border get_fict_border(double)const; // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  extend_data *ed;
public:
  pl_vspogt_body(geophys_params&g);
  ~pl_vspogt_body(){paActiveBorders.clean();delete ed;};


  void prepare_for_vspogt(const pl_intern_borders&,rays_to_receiver* prtr,
              const DoubleAr& daRT);
  void get_initial_fict_borders(pl_fict_borders &)const;

  // ¿ÐàÐ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ÝÞÜÕàã ßÛÞåÞØÝâÕàßÞÛØàÞÒÐÝÝÞÓÞ ¶8() ÒàÕÜÕÝØ
  pl_fict_borders* get_fict_borders(int) const;
  // ½ÕÚÞâÞàÞÕ ÚÞÛØçÕáâÒÞ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ã×Ûã ib
  void add_n_borders(int iNNewPointsLeft,int ib,pl_fict_borders*,
             ray_parameters_dynar *)const;

  // ½ÞÜÕà íâÞÓÞáÐÜÞÓÞ ã×ÛÐ
  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  void remove(int i);
public:

  // ½Ð áÛãçÐÙ, ÕáÛØ áÞ×ÔÐÕâáï ÛÞÚÐÛìÝëÙ áßØáÞÚ, ÚÞâÞàëÙ ÝÕ ÑãÔÕÜ Ò ±´ ßØáÐâì  time_depth_list_local(unsigned quant = DEFAULT_QUANT,
            unsigned volume = START_VOLUME);
  
  // ¿Þ ÒàÕÜÕÝØ - ÛØÝÕÙÝÐï ØÝâÕàßÞÛïæØï ßÐàÐÜÕâàÞÒ
  void interpolate(double dTime,double& prop,double& diver,double& refrat,
           point&pRefl,pl_spline& sReflect) const; 
  void interpolate(double dTime, double* vpol) const;
  double start_time()const{return daTimes(ia(0));};
  double finish_time()const{return daTimes(ia(_last_));};


  // ²àÕÜï áÕàñÔÚØ áÕÓÜÕÝâÐ á ßÕàÕÔÐÒÐÕÜëÜ ÝÞÜÕàÞÜ
  double t_of_new_border(int) const;

  // ÂØßÐ ÒÝÕáÛØ ÝÐçÐÛìÝëÕ ÞâàÐÖÕÝØï ØâÐÔÐÛÕÕ
  void prepared(){lPrepared=TRUE;};

  int n_of_segments() const; // ºÞÛØçÕáâÒÞ áÕÓÜÕÝâÞÒ Ò áßØáÚÕ

  // ÔÞÑÐÒØâì ÒàÕÜï-âÞçÚã-áßÛÐÙÝ
  void add_new_node(double, const point3d &,
            p_pl_intern_spline,ray_to_receiver &,Logic); 
  void add_new_node(double, const point3d &,
            p_pl_intern_spline,ray_to_receiver &,double); 

  void merge(const time_depth_list_local&); // ¿àØÑÐÒØâì Ú ÞÔÝÞÜã áßØáÚã ÔàãÓÞÙ

  // ¿Þ ÒàÕÜÕÝØ - ÝÞÜÕà ÞâàÕ×ÚÐ
  int t2i(double)const;

  // ÁÛØïÝØÕ ÔÒãå âØáßÐáßØáÚÞÒ ÑÕ× ßÞÒâÞàÞÒ
  void merge_unique(const time_depth_list_local&); 

  // ÂÞ ÖÕ, ÝÞ ÔÞÑÐÒÛïÕâ ÝÞÒëÕ íÛÕÜÕÝâë ØáÚÛîçØâÕÛìÝÞ ßÞ ÚàÐïÜ
  void merge_unique_monotone(const time_depth_list_local&); 

  // ÀÐ×ÝÞÓÞ àÞÔÐ çØáâÚØ  
  void clear();
  void pack();

  // ÀÐ×ÛØçÝëÕ ÞâÛÐÔÞçÝëÕ äØçÚØ
  void debug_output() const;

  friend class time_depth_lists;
};

class time_depth_list: public db_container,public time_depth_list_local{
protected:
  dbHandle dbhReceivers;
  dbHandle dbhTD;
  static int iLast;
  int N1,N2;
  int iTDLStoreRegime;

  // ÀÐÑÞâÐ á ±´
  int read_from_db (dbItem& dbiTD, dbItem& dbiRN);
  int write_to_db (dbItem& dbiTD, dbItem& dbiRN)const;
  void get_item (dbItem& dbiRes, dbItem& dbiRN) const;
public:
  // µáÛØ âÐÚØÜ ÞÑàÐ×ÞÜ áÞ×ÔÐÝ - ÜÞÖÕÜ ßØáÐâì Ò ±´
  time_depth_list(db_results_manager& d,char*sp, int id,int n1,int n2,
          int iTDLSR=PM_ADD,
          unsigned quant = DEFAULT_QUANT,unsigned volume = START_VOLUME);
  // ±´
  virtual int store() const;
  virtual int load();
  void prepare_pack_load (dbHandle&dbhR,dbHandle&dbhT);
  void done_loading();
  void pack_load(dbHandle&,dbHandle&);

  friend class time_depth_lists;
};

typedef time_depth_list* p_time_depth_list;
typedef ObjectsAr<p_time_depth_list> time_depth_lists_dynar;

class time_depth_lists:public time_depth_lists_dynar{
public:
  p_time_depth_list get(int,int)const;
};

class pl_vspogt_body:public pl_intern_body,public time_depth_list_local{
protected:
  pl_intern_borders                                       paActiveBorders;  
  DoubleAr daActualTimes;           // ²àÕÜÕÝÐ, ÝÐ ÚÞâÞàëå ÞâàÐÖÐÛØáì ^ ÓàÐÝØæë

  Logic interpolated(int) const;
  p_pl_fict_border get_fict_border(double)const; // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  extend_data *ed;
public:
  pl_vspogt_body(geophys_params&g);
  ~pl_vspogt_body(){paActiveBorders.clean();delete ed;};


  void prepare_for_vspogt(const pl_intern_borders&,rays_to_receiver* prtr,
              const DoubleAr& daRT);
  void get_initial_fict_borders(pl_fict_borders &)const;

  // ¿ÐàÐ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ÝÞÜÕàã ßÛÞåÞØÝâÕàßÞÛØàÞÒÐÝÝÞÓÞ ¶8() ÒàÕÜÕÝØ
  pl_fict_borders* get_fict_borders(int) const;
  // ½ÕÚÞâÞàÞÕ ÚÞÛØçÕáâÒÞ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ã×Ûã ib
  void add_n_borders(int iNNewPointsLeft,int ib,pl_fict_borders*,
             ray_parameters_dynar *)const;

  // ½ÞÜÕà íâÞÓÞáÐÜÞÓÞ ã×ÛÐ
  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  time_depth_list_local(unsigned quant = DEFAULT_QUANT,
            unsigned volume = START_VOLUME);
  
  // ¿Þ ÒàÕÜÕÝØ - ÛØÝÕÙÝÐï ØÝâÕàßÞÛïæØï ßÐàÐÜÕâàÞÒ
  void interpolate(double dTime,double& prop,double& diver,double& refrat,
           point&pRefl,pl_spline& sReflect) const; 
  void interpolate(double dTime, double* vpol) const;
  double start_time()const{return daTimes(ia(0));};
  double finish_time()const{return daTimes(ia(_last_));};


  // ²àÕÜï áÕàñÔÚØ áÕÓÜÕÝâÐ á ßÕàÕÔÐÒÐÕÜëÜ ÝÞÜÕàÞÜ
  double t_of_new_border(int) const;

  // ÂØßÐ ÒÝÕáÛØ ÝÐçÐÛìÝëÕ ÞâàÐÖÕÝØï ØâÐÔÐÛÕÕ
  void prepared(){lPrepared=TRUE;};

  int n_of_segments() const; // ºÞÛØçÕáâÒÞ áÕÓÜÕÝâÞÒ Ò áßØáÚÕ

  // ÔÞÑÐÒØâì ÒàÕÜï-âÞçÚã-áßÛÐÙÝ  void add_new_node(double, const point3d &,
            p_pl_intern_spline,ray_to_receiver &,Logic); 
  void add_new_node(double, const point3d &,
            p_pl_intern_spline,ray_to_receiver &,double); 

  void merge(const time_depth_list_local&); // ¿àØÑÐÒØâì Ú ÞÔÝÞÜã áßØáÚã ÔàãÓÞÙ

  // ¿Þ ÒàÕÜÕÝØ - ÝÞÜÕà ÞâàÕ×ÚÐ
  int t2i(double)const;

  // ÁÛØïÝØÕ ÔÒãå âØáßÐáßØáÚÞÒ ÑÕ× ßÞÒâÞàÞÒ
  void merge_unique(const time_depth_list_local&); 

  // ÂÞ ÖÕ, ÝÞ ÔÞÑÐÒÛïÕâ ÝÞÒëÕ íÛÕÜÕÝâë ØáÚÛîçØâÕÛìÝÞ ßÞ ÚàÐïÜ
  void merge_unique_monotone(const time_depth_list_local&); 

  // ÀÐ×ÝÞÓÞ àÞÔÐ çØáâÚØ  
  void clear();
  void pack();

  // ÀÐ×ÛØçÝëÕ ÞâÛÐÔÞçÝëÕ äØçÚØ
  void debug_output() const;

  friend class time_depth_lists;
};

class time_depth_list: public db_container,public time_depth_list_local{
protected:
  dbHandle dbhReceivers;
  dbHandle dbhTD;
  static int iLast;
  int N1,N2;
  int iTDLStoreRegime;

  // ÀÐÑÞâÐ á ±´
  int read_from_db (dbItem& dbiTD, dbItem& dbiRN);
  int write_to_db (dbItem& dbiTD, dbItem& dbiRN)const;
  void get_item (dbItem& dbiRes, dbItem& dbiRN) const;
public:
  // µáÛØ âÐÚØÜ ÞÑàÐ×ÞÜ áÞ×ÔÐÝ - ÜÞÖÕÜ ßØáÐâì Ò ±´
  time_depth_list(db_results_manager& d,char*sp, int id,int n1,int n2,
          int iTDLSR=PM_ADD,
          unsigned quant = DEFAULT_QUANT,unsigned volume = START_VOLUME);
  // ±´
  virtual int store() const;
  virtual int load();
  void prepare_pack_load (dbHandle&dbhR,dbHandle&dbhT);
  void done_loading();
  void pack_load(dbHandle&,dbHandle&);

  friend class time_depth_lists;
};

typedef time_depth_list* p_time_depth_list;
typedef ObjectsAr<p_time_depth_list> time_depth_lists_dynar;

class time_depth_lists:public time_depth_lists_dynar{
public:
  p_time_depth_list get(int,int)const;
};

class pl_vspogt_body:public pl_intern_body,public time_depth_list_local{
protected:
  pl_intern_borders                                       paActiveBorders;  
  DoubleAr daActualTimes;           // ²àÕÜÕÝÐ, ÝÐ ÚÞâÞàëå ÞâàÐÖÐÛØáì ^ ÓàÐÝØæë

  Logic interpolated(int) const;
  p_pl_fict_border get_fict_border(double)const; // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  extend_data *ed;
public:
  pl_vspogt_body(geophys_params&g);
  ~pl_vspogt_body(){paActiveBorders.clean();delete ed;};


  void prepare_for_vspogt(const pl_intern_borders&,rays_to_receiver* prtr,
              const DoubleAr& daRT);
  void get_initial_fict_borders(pl_fict_borders &)const;

  // ¿ÐàÐ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ÝÞÜÕàã ßÛÞåÞØÝâÕàßÞÛØàÞÒÐÝÝÞÓÞ ¶8() ÒàÕÜÕÝØ
  pl_fict_borders* get_fict_borders(int) const;
  // ½ÕÚÞâÞàÞÕ ÚÞÛØçÕáâÒÞ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ã×Ûã ib
  void add_n_borders(int iNNewPointsLeft,int ib,pl_fict_borders*,
             ray_parameters_dynar *)const;

  // ½ÞÜÕà íâÞÓÞáÐÜÞÓÞ ã×ÛÐ
  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  void add_new_node(double, const point3d &,
            p_pl_intern_spline,ray_to_receiver &,Logic); 
  void add_new_node(double, const point3d &,
            p_pl_intern_spline,ray_to_receiver &,double); 

  void merge(const time_depth_list_local&); // ¿àØÑÐÒØâì Ú ÞÔÝÞÜã áßØáÚã ÔàãÓÞÙ
  // ¿Þ ÒàÕÜÕÝØ - ÝÞÜÕà ÞâàÕ×ÚÐ
  int t2i(double)const;

  // ÁÛØïÝØÕ ÔÒãå âØáßÐáßØáÚÞÒ ÑÕ× ßÞÒâÞàÞÒ
  void merge_unique(const time_depth_list_local&); 

  // ÂÞ ÖÕ, ÝÞ ÔÞÑÐÒÛïÕâ ÝÞÒëÕ íÛÕÜÕÝâë ØáÚÛîçØâÕÛìÝÞ ßÞ ÚàÐïÜ
  void merge_unique_monotone(const time_depth_list_local&); 

  // ÀÐ×ÝÞÓÞ àÞÔÐ çØáâÚØ  
  void clear();
  void pack();

  // ÀÐ×ÛØçÝëÕ ÞâÛÐÔÞçÝëÕ äØçÚØ
  void debug_output() const;

  friend class time_depth_lists;
};

class time_depth_list: public db_container,public time_depth_list_local{
protected:
  dbHandle dbhReceivers;
  dbHandle dbhTD;
  static int iLast;
  int N1,N2;
  int iTDLStoreRegime;

  // ÀÐÑÞâÐ á ±´
  int read_from_db (dbItem& dbiTD, dbItem& dbiRN);
  int write_to_db (dbItem& dbiTD, dbItem& dbiRN)const;
  void get_item (dbItem& dbiRes, dbItem& dbiRN) const;
public:
  // µáÛØ âÐÚØÜ ÞÑàÐ×ÞÜ áÞ×ÔÐÝ - ÜÞÖÕÜ ßØáÐâì Ò ±´
  time_depth_list(db_results_manager& d,char*sp, int id,int n1,int n2,
          int iTDLSR=PM_ADD,
          unsigned quant = DEFAULT_QUANT,unsigned volume = START_VOLUME);
  // ±´
  virtual int store() const;
  virtual int load();
  void prepare_pack_load (dbHandle&dbhR,dbHandle&dbhT);
  void done_loading();
  void pack_load(dbHandle&,dbHandle&);

  friend class time_depth_lists;
};

typedef time_depth_list* p_time_depth_list;
typedef ObjectsAr<p_time_depth_list> time_depth_lists_dynar;

class time_depth_lists:public time_depth_lists_dynar{
public:
  p_time_depth_list get(int,int)const;
};

class pl_vspogt_body:public pl_intern_body,public time_depth_list_local{
protected:
  pl_intern_borders                                       paActiveBorders;  
  DoubleAr daActualTimes;           // ²àÕÜÕÝÐ, ÝÐ ÚÞâÞàëå ÞâàÐÖÐÛØáì ^ ÓàÐÝØæë

  Logic interpolated(int) const;
  p_pl_fict_border get_fict_border(double)const; // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  extend_data *ed;
public:
  pl_vspogt_body(geophys_params&g);
  ~pl_vspogt_body(){paActiveBorders.clean();delete ed;};


  void prepare_for_vspogt(const pl_intern_borders&,rays_to_receiver* prtr,
              const DoubleAr& daRT);
  void get_initial_fict_borders(pl_fict_borders &)const;

  // ¿ÐàÐ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ÝÞÜÕàã ßÛÞåÞØÝâÕàßÞÛØàÞÒÐÝÝÞÓÞ ¶8() ÒàÕÜÕÝØ
  pl_fict_borders* get_fict_borders(int) const;
  // ½ÕÚÞâÞàÞÕ ÚÞÛØçÕáâÒÞ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ã×Ûã ib
  void add_n_borders(int iNNewPointsLeft,int ib,pl_fict_borders*,
             ray_parameters_dynar *)const;

  // ½ÞÜÕà íâÞÓÞáÐÜÞÓÞ ã×ÛÐ
  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 


  // ¿Þ ÒàÕÜÕÝØ - ÝÞÜÕà ÞâàÕ×ÚÐ  int t2i(double)const;

  // ÁÛØïÝØÕ ÔÒãå âØáßÐáßØáÚÞÒ ÑÕ× ßÞÒâÞàÞÒ
  void merge_unique(const time_depth_list_local&); 

  // ÂÞ ÖÕ, ÝÞ ÔÞÑÐÒÛïÕâ ÝÞÒëÕ íÛÕÜÕÝâë ØáÚÛîçØâÕÛìÝÞ ßÞ ÚàÐïÜ
  void merge_unique_monotone(const time_depth_list_local&); 

  // ÀÐ×ÝÞÓÞ àÞÔÐ çØáâÚØ  
  void clear();
  void pack();

  // ÀÐ×ÛØçÝëÕ ÞâÛÐÔÞçÝëÕ äØçÚØ
  void debug_output() const;

  friend class time_depth_lists;
};

class time_depth_list: public db_container,public time_depth_list_local{
protected:
  dbHandle dbhReceivers;
  dbHandle dbhTD;
  static int iLast;
  int N1,N2;
  int iTDLStoreRegime;

  // ÀÐÑÞâÐ á ±´
  int read_from_db (dbItem& dbiTD, dbItem& dbiRN);
  int write_to_db (dbItem& dbiTD, dbItem& dbiRN)const;
  void get_item (dbItem& dbiRes, dbItem& dbiRN) const;
public:
  // µáÛØ âÐÚØÜ ÞÑàÐ×ÞÜ áÞ×ÔÐÝ - ÜÞÖÕÜ ßØáÐâì Ò ±´
  time_depth_list(db_results_manager& d,char*sp, int id,int n1,int n2,
          int iTDLSR=PM_ADD,
          unsigned quant = DEFAULT_QUANT,unsigned volume = START_VOLUME);
  // ±´
  virtual int store() const;
  virtual int load();
  void prepare_pack_load (dbHandle&dbhR,dbHandle&dbhT);
  void done_loading();
  void pack_load(dbHandle&,dbHandle&);

  friend class time_depth_lists;
};

typedef time_depth_list* p_time_depth_list;
typedef ObjectsAr<p_time_depth_list> time_depth_lists_dynar;

class time_depth_lists:public time_depth_lists_dynar{
public:
  p_time_depth_list get(int,int)const;
};

class pl_vspogt_body:public pl_intern_body,public time_depth_list_local{
protected:
  pl_intern_borders                                       paActiveBorders;  
  DoubleAr daActualTimes;           // ²àÕÜÕÝÐ, ÝÐ ÚÞâÞàëå ÞâàÐÖÐÛØáì ^ ÓàÐÝØæë

  Logic interpolated(int) const;
  p_pl_fict_border get_fict_border(double)const; // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  extend_data *ed;
public:
  pl_vspogt_body(geophys_params&g);
  ~pl_vspogt_body(){paActiveBorders.clean();delete ed;};


  void prepare_for_vspogt(const pl_intern_borders&,rays_to_receiver* prtr,
              const DoubleAr& daRT);
  void get_initial_fict_borders(pl_fict_borders &)const;

  // ¿ÐàÐ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ÝÞÜÕàã ßÛÞåÞØÝâÕàßÞÛØàÞÒÐÝÝÞÓÞ ¶8() ÒàÕÜÕÝØ
  pl_fict_borders* get_fict_borders(int) const;
  // ½ÕÚÞâÞàÞÕ ÚÞÛØçÕáâÒÞ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ã×Ûã ib
  void add_n_borders(int iNNewPointsLeft,int ib,pl_fict_borders*,
             ray_parameters_dynar *)const;

  // ½ÞÜÕà íâÞÓÞáÐÜÞÓÞ ã×ÛÐ
  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  int t2i(double)const;

  // ÁÛØïÝØÕ ÔÒãå âØáßÐáßØáÚÞÒ ÑÕ× ßÞÒâÞàÞÒ
  void merge_unique(const time_depth_list_local&); 

  // ÂÞ ÖÕ, ÝÞ ÔÞÑÐÒÛïÕâ ÝÞÒëÕ íÛÕÜÕÝâë ØáÚÛîçØâÕÛìÝÞ ßÞ ÚàÐïÜ  void merge_unique_monotone(const time_depth_list_local&); 

  // ÀÐ×ÝÞÓÞ àÞÔÐ çØáâÚØ  
  void clear();
  void pack();

  // ÀÐ×ÛØçÝëÕ ÞâÛÐÔÞçÝëÕ äØçÚØ
  void debug_output() const;

  friend class time_depth_lists;
};

class time_depth_list: public db_container,public time_depth_list_local{
protected:
  dbHandle dbhReceivers;
  dbHandle dbhTD;
  static int iLast;
  int N1,N2;
  int iTDLStoreRegime;

  // ÀÐÑÞâÐ á ±´
  int read_from_db (dbItem& dbiTD, dbItem& dbiRN);
  int write_to_db (dbItem& dbiTD, dbItem& dbiRN)const;
  void get_item (dbItem& dbiRes, dbItem& dbiRN) const;
public:
  // µáÛØ âÐÚØÜ ÞÑàÐ×ÞÜ áÞ×ÔÐÝ - ÜÞÖÕÜ ßØáÐâì Ò ±´
  time_depth_list(db_results_manager& d,char*sp, int id,int n1,int n2,
          int iTDLSR=PM_ADD,
          unsigned quant = DEFAULT_QUANT,unsigned volume = START_VOLUME);
  // ±´
  virtual int store() const;
  virtual int load();
  void prepare_pack_load (dbHandle&dbhR,dbHandle&dbhT);
  void done_loading();
  void pack_load(dbHandle&,dbHandle&);

  friend class time_depth_lists;
};

typedef time_depth_list* p_time_depth_list;
typedef ObjectsAr<p_time_depth_list> time_depth_lists_dynar;

class time_depth_lists:public time_depth_lists_dynar{
public:
  p_time_depth_list get(int,int)const;
};

class pl_vspogt_body:public pl_intern_body,public time_depth_list_local{
protected:
  pl_intern_borders                                       paActiveBorders;  
  DoubleAr daActualTimes;           // ²àÕÜÕÝÐ, ÝÐ ÚÞâÞàëå ÞâàÐÖÐÛØáì ^ ÓàÐÝØæë

  Logic interpolated(int) const;
  p_pl_fict_border get_fict_border(double)const; // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  extend_data *ed;
public:
  pl_vspogt_body(geophys_params&g);
  ~pl_vspogt_body(){paActiveBorders.clean();delete ed;};


  void prepare_for_vspogt(const pl_intern_borders&,rays_to_receiver* prtr,
              const DoubleAr& daRT);
  void get_initial_fict_borders(pl_fict_borders &)const;

  // ¿ÐàÐ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ÝÞÜÕàã ßÛÞåÞØÝâÕàßÞÛØàÞÒÐÝÝÞÓÞ ¶8() ÒàÕÜÕÝØ
  pl_fict_borders* get_fict_borders(int) const;
  // ½ÕÚÞâÞàÞÕ ÚÞÛØçÕáâÒÞ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ã×Ûã ib
  void add_n_borders(int iNNewPointsLeft,int ib,pl_fict_borders*,
             ray_parameters_dynar *)const;

  // ½ÞÜÕà íâÞÓÞáÐÜÞÓÞ ã×ÛÐ
  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  void merge_unique_monotone(const time_depth_list_local&); 

  // ÀÐ×ÝÞÓÞ àÞÔÐ çØáâÚØ  
  void clear();
  void pack();

  // ÀÐ×ÛØçÝëÕ ÞâÛÐÔÞçÝëÕ äØçÚØ
  void debug_output() const;

  friend class time_depth_lists;
};

class time_depth_list: public db_container,public time_depth_list_local{
protected:
  dbHandle dbhReceivers;
  dbHandle dbhTD;
  static int iLast;
  int N1,N2;
  int iTDLStoreRegime;

  // ÀÐÑÞâÐ á ±´  int read_from_db (dbItem& dbiTD, dbItem& dbiRN);
  int write_to_db (dbItem& dbiTD, dbItem& dbiRN)const;
  void get_item (dbItem& dbiRes, dbItem& dbiRN) const;
public:
  // µáÛØ âÐÚØÜ ÞÑàÐ×ÞÜ áÞ×ÔÐÝ - ÜÞÖÕÜ ßØáÐâì Ò ±´
  time_depth_list(db_results_manager& d,char*sp, int id,int n1,int n2,
          int iTDLSR=PM_ADD,
          unsigned quant = DEFAULT_QUANT,unsigned volume = START_VOLUME);
  // ±´
  virtual int store() const;
  virtual int load();
  void prepare_pack_load (dbHandle&dbhR,dbHandle&dbhT);
  void done_loading();
  void pack_load(dbHandle&,dbHandle&);

  friend class time_depth_lists;
};

typedef time_depth_list* p_time_depth_list;
typedef ObjectsAr<p_time_depth_list> time_depth_lists_dynar;

class time_depth_lists:public time_depth_lists_dynar{
public:
  p_time_depth_list get(int,int)const;
};

class pl_vspogt_body:public pl_intern_body,public time_depth_list_local{
protected:
  pl_intern_borders                                       paActiveBorders;  
  DoubleAr daActualTimes;           // ²àÕÜÕÝÐ, ÝÐ ÚÞâÞàëå ÞâàÐÖÐÛØáì ^ ÓàÐÝØæë

  Logic interpolated(int) const;
  p_pl_fict_border get_fict_border(double)const; // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  extend_data *ed;
public:
  pl_vspogt_body(geophys_params&g);
  ~pl_vspogt_body(){paActiveBorders.clean();delete ed;};


  void prepare_for_vspogt(const pl_intern_borders&,rays_to_receiver* prtr,
              const DoubleAr& daRT);
  void get_initial_fict_borders(pl_fict_borders &)const;

  // ¿ÐàÐ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ÝÞÜÕàã ßÛÞåÞØÝâÕàßÞÛØàÞÒÐÝÝÞÓÞ ¶8() ÒàÕÜÕÝØ
  pl_fict_borders* get_fict_borders(int) const;
  // ½ÕÚÞâÞàÞÕ ÚÞÛØçÕáâÒÞ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ã×Ûã ib
  void add_n_borders(int iNNewPointsLeft,int ib,pl_fict_borders*,
             ray_parameters_dynar *)const;

  // ½ÞÜÕà íâÞÓÞáÐÜÞÓÞ ã×ÛÐ
  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  int read_from_db (dbItem& dbiTD, dbItem& dbiRN);
  int write_to_db (dbItem& dbiTD, dbItem& dbiRN)const;
  void get_item (dbItem& dbiRes, dbItem& dbiRN) const;
public:
  // µáÛØ âÐÚØÜ ÞÑàÐ×ÞÜ áÞ×ÔÐÝ - ÜÞÖÕÜ ßØáÐâì Ò ±´
  time_depth_list(db_results_manager& d,char*sp, int id,int n1,int n2,
          int iTDLSR=PM_ADD,
          unsigned quant = DEFAULT_QUANT,unsigned volume = START_VOLUME);
  // ±´
  virtual int store() const;
  virtual int load();
  void prepare_pack_load (dbHandle&dbhR,dbHandle&dbhT);
  void done_loading();
  void pack_load(dbHandle&,dbHandle&);

  friend class time_depth_lists;
};

typedef time_depth_list* p_time_depth_list;
typedef ObjectsAr<p_time_depth_list> time_depth_lists_dynar;

class time_depth_lists:public time_depth_lists_dynar{
public:
  p_time_depth_list get(int,int)const;
};

class pl_vspogt_body:public pl_intern_body,public time_depth_list_local{
protected:
  pl_intern_borders                                       paActiveBorders;  
  DoubleAr daActualTimes;           // ²àÕÜÕÝÐ, ÝÐ ÚÞâÞàëå ÞâàÐÖÐÛØáì ^ ÓàÐÝØæë
  Logic interpolated(int) const;
  p_pl_fict_border get_fict_border(double)const; // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  extend_data *ed;
public:
  pl_vspogt_body(geophys_params&g);
  ~pl_vspogt_body(){paActiveBorders.clean();delete ed;};


  void prepare_for_vspogt(const pl_intern_borders&,rays_to_receiver* prtr,
              const DoubleAr& daRT);
  void get_initial_fict_borders(pl_fict_borders &)const;

  // ¿ÐàÐ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ÝÞÜÕàã ßÛÞåÞØÝâÕàßÞÛØàÞÒÐÝÝÞÓÞ ¶8() ÒàÕÜÕÝØ
  pl_fict_borders* get_fict_borders(int) const;
  // ½ÕÚÞâÞàÞÕ ÚÞÛØçÕáâÒÞ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ã×Ûã ib
  void add_n_borders(int iNNewPointsLeft,int ib,pl_fict_borders*,
             ray_parameters_dynar *)const;

  // ½ÞÜÕà íâÞÓÞáÐÜÞÓÞ ã×ÛÐ
  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 


  Logic interpolated(int) const;
  p_pl_fict_border get_fict_border(double)const; // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ  extend_data *ed;
public:
  pl_vspogt_body(geophys_params&g);
  ~pl_vspogt_body(){paActiveBorders.clean();delete ed;};


  void prepare_for_vspogt(const pl_intern_borders&,rays_to_receiver* prtr,
              const DoubleAr& daRT);
  void get_initial_fict_borders(pl_fict_borders &)const;

  // ¿ÐàÐ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ÝÞÜÕàã ßÛÞåÞØÝâÕàßÞÛØàÞÒÐÝÝÞÓÞ ¶8() ÒàÕÜÕÝØ
  pl_fict_borders* get_fict_borders(int) const;
  // ½ÕÚÞâÞàÞÕ ÚÞÛØçÕáâÒÞ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ã×Ûã ib
  void add_n_borders(int iNNewPointsLeft,int ib,pl_fict_borders*,
             ray_parameters_dynar *)const;

  // ½ÞÜÕà íâÞÓÞáÐÜÞÓÞ ã×ÛÐ
  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  extend_data *ed;
public:
  pl_vspogt_body(geophys_params&g);
  ~pl_vspogt_body(){paActiveBorders.clean();delete ed;};


  void prepare_for_vspogt(const pl_intern_borders&,rays_to_receiver* prtr,
              const DoubleAr& daRT);
  void get_initial_fict_borders(pl_fict_borders &)const;

  // ¿ÐàÐ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ÝÞÜÕàã ßÛÞåÞØÝâÕàßÞÛØàÞÒÐÝÝÞÓÞ ¶8() ÒàÕÜÕÝØ  pl_fict_borders* get_fict_borders(int) const;
  // ½ÕÚÞâÞàÞÕ ÚÞÛØçÕáâÒÞ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ã×Ûã ib
  void add_n_borders(int iNNewPointsLeft,int ib,pl_fict_borders*,
             ray_parameters_dynar *)const;

  // ½ÞÜÕà íâÞÓÞáÐÜÞÓÞ ã×ÛÐ
  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  pl_fict_borders* get_fict_borders(int) const;
  // ½ÕÚÞâÞàÞÕ ÚÞÛØçÕáâÒÞ äØÚâØÒÝëå ÓàÐÝØæ ßÞ ã×Ûã ib
  void add_n_borders(int iNNewPointsLeft,int ib,pl_fict_borders*,
             ray_parameters_dynar *)const;

  // ½ÞÜÕà íâÞÓÞáÐÜÞÓÞ ã×ÛÐ  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ
  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  int bad_interpolation();

  void clear();

  //¿ÞÛãçÕÝØÕ âàÞÙÚØ äØÚâØÒÝëå ÓàÐÝØæ ÝØÖÕ ÞâàÐÖÐîéÕÙ (áÛãçÐÙ ÓÛãÑÞçÐÙèÕÓÞ âÕÛÐ)
  pl_fict_borders*extend_initial(ray* prInitial,double dTLast,DoubleAr&daTimes,
                 p_pl_intern_border);

  // ¿àÞÔÞÛÖÕÝØÕ ÛãçÐ ÒÝãâàØ âÕÛÐ (ßÞ ÔÐÝÝëÜ, ãÖÕ ßÞÔÓÞâÞÒÛÕÝÝëÜ ^^^)
  Logic extend(double dWhere,point& pRes);

  // ÄØÚâØÒÝÐï ÓàÐÝØæÐ ßÞ ÒàÕÜÕÝØ  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ
  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  p_pl_fict_border get_fict_border(pl_fict_border&,point)const; 

  p_pl_fict_border get_fict_horizontal_border(point)const;
  
};

typedef pl_vspogt_body* p_pl_vspogt_body;
typedef ObjectsAr<p_pl_vspogt_body> pl_vspogt_bodies;

typedef ObjectsAr<IntAr*> aKindaBorders;
class mixer:public aSegments2D{
protected:
  int *iaFlags;
  aKindaBorders akbCurrent;

  void clear_data();

  // ¾ßàÕÔÕÛïÕâ, ÛÕÖØâ ÛØ âÞçÚÐ Ò áÞ×ÔÐÝÝÞÜ âÕÛÕ  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã
  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  Logic point_in_body(const point2D& pWhat,int iB);

  // ¾ßàÕÔÕÛïÕâ ÝÐßàÐÒÛÕÝØÕ ßÞ âÕÛã  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ
  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  void get_direction_body(const point2D& pWhat,int iB,double& A,double& B);

  // ¿ëâÐÕâáï áÞ×ÔÐâì ÕéÕ ÞÔÝÞ âÕÛÞ  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ
  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  IntAr* new_body (const point2D& pWhat);

  // ·ÐáãÒÐÕâ âÕÛÞ Ò áßØáÞÚ
  void add_body(IntAr*piaBody);

  inline int sign(const point2D&,segment2D&);
  int segment_other_body_line(int iB);
public:
  mixer(){iaFlags=NULL;};
  // ¿Þ âÞçÚÕ ÒëÔÐÕâ ÝÐßàÐÒÛÕÝØÕ áÜÕáØâÕÛï Ò ÝÕÙ  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ
  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  void get_mix_direction(const point2D& pWhere,double& A,double& B);

  // ¿ÕàÕßØáëÒÐÕÜ ÒáâÐÒÚã - ãÔÐÛÕÝØÕ  virtual void addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 

  virtual void  addh(const segment2D & val);
  virtual void  addl(const segment2D& val);

  virtual void  addh(const aSegments2D& pArray){Throw(ERROR__N_A);};
  virtual void  addl(const aSegments2D& pArray){Throw(ERROR__N_A);};

  virtual void  insert (unsigned i, const aSegments2D& val){Throw(ERROR__N_A);};
  // remove ßÞÚÐ çâÞ ÝÕ ×ÝÐî, ÚÐÚ ßÕàÕÔÕÛÐâì, ÔÐ Ø ÝÐÔÞ ÛØ
};
inline int mixer::sign (const point2D&p,segment2D&s){
  double dEps=0.0001;
  double A,B,C;
  s.get_coefs(A,B,C);
  double S=(A*p.X()+B*p.Y()+C);
  if(S>dEps)
    return 1;
  else if(S<-dEps)
    return -1;
  else
    return 0;
};

class fict_bord_hit_info:public hit_point_info{
protected:
  rays_shoot * pr;
  int iBody1;
  int iBody2;
  const pl_fict_border &fbReflect;
  point pRes; 
  double dTime;
  double dAcc;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  fict_bord_hit_info(rays_shoot*,int,int,const pl_fict_border &,double,
             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);
  point3d stopped_at()const;
  double time()const;
};

#endif /* p_vspogt_elements.H */ 
/*  Melnikov George Yu. 19-Mar-1999  */ 
/*  Melnikov George Yu. 23-Mar-1999  */ 
/*  Melnikov George Yu. 24-Mar-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */