rays.H

/* rays.H */
/* $Id: rays.H,v 2.12 2005/07/18 07:16:06 nikita Exp $ */
#ifndef __rays_H
#define __rays_H

/*
   ºÛÐááë, ÞßØáÐÝÝëÕ Ò ×ÐÓÞÛÞÒÚÕ:
   ray -- Ûãç, ÒëåÞÔïéØÙ Ø× ØáâÞçÝØÚÐ;

   ray_parameters -- ßÐàÐÜÕâàë ÛãçÐ, ßÞ ÚÞâÞàëÜ, ÒÝÕ ×ÐÒØáØÜÞáâØ Þâ âÕÚãéÕÓÞ
   áÞáâÞïÝØï rays_set ÜÞÖÝÞ ÒÞááâÐÝÞÒØâì ÝÐåÞÖÔÕÝØÕ ÛãçÐ Ò ßàïÜÞãÓÞÛìÝÞÙ áÕâÚÕ.

   rays_set -- ßàïÜÞãÓÞÛìÝÐï áÕâÚÐ ÛãçÕÙ, ÚÞâÞàÐï ßÞ ÜÕàÕ ÝÐÔÞÑÝÞáâØ   àÐ×àÕ×ÐÕâáï, Ø Ò ÚÞâÞàÞÙ åàÐÝïâáï ÒáÕ ÛãçØ, ÒëèÕÔèØÕ Ø× ØáâÞçÝØÚÐ.
 */


#include <ldb/dbtools.h>
#include <mix/DoubleAr.H>

#include <mix/ObjectsAr.H>
#include <mod3d/mod3d_basic.H>

#define PM_FB 1
#define PM_FULL 0

#define PM_CREATE 1
#define PM_ADD 2

// ºÛÐááë, ØáßÞÛì×ãÕÜëÕ Ò ×ÐÓÞÛÞÒÚÕ
class crossings;
class crossing;
class crossings_ptrs;
class pl_raytr_model;
class pl_fict_border;
class ray_specifics;

// start of insertion by nikita
typedef struct rayPathInfo {
  aPoints* points;
  int nReflectBorder;
  int receiver;
};

#define AssocAr rayPathAr

typedef struct {
  unsigned  key;//ray identifier in database RAYPATHINFO table
  rayPathInfo   value;//ray path points and receiver, reflection border binding
}  rayPath;
#define AssocItemType      rayPath
#define AssocKeyField      key
#define AssocKeyType       unsigned
#define AssocKeysCmp(key1,key2)    COMPAR(key1,key2)
#include <mix/AssocAr.H>

// end of insertion by nikita 



class ray_parameters{
protected:
  // ² àÐÔØÐÝÐå
  double dFi,dAz;
public:
  void set_pars(double,double);
  ray_parameters(){;};
  ray_parameters(double,double);
friend int operator==(const ray_parameters &rp1,
              const ray_parameters &rp2);
  double fi()const{return dFi;};
  double az()const{return dAz;};
};

#define AnyType ray_parameters
#define AnyAr ray_parameters_dynar
#include <mix/AnyAr.H>

#define R_BORDER    1
#define R_NODE      2
class ray{
protected:
  pl_raytr_model *pm;
  int iStored; // ¾ÑÞ×ÝÐçÐÕâ, Òë×ëÒÐÛØ ÛØ ÔÛï ÛãçÐ ×ÐßØáì ÕÓÞ ßÕàÕáÕçÕÝØÙ
  ray_parameters rp; // ¾âÚãÔÐ ÒëáâàÕÛØÛØ
  // ÁßØáÞÚ ãÚÐ×ÐâÕÛÕÙ ÝÐ ßÕàÕáÕçÕÝØï
  crossings_ptrs *cr;

  // ÁÞ×ÔÐÕâáï âÞÛìÚÞ Ø× rays_shoot
  ray(const ray_parameters &RP,crossings_ptrs *CR,pl_raytr_model*p);

  // ¿àÞÒÕàïÕâ, ßàÞåÞÔØâ ÛØ Ûãç çÕàÕ× ÓàÐÝØæã á ãÚÐ×ÐÝÝëÜ ÝÞÜÕàÞÜ
  Logic crosses(int i);

  // ²Þ×ÒàÐéÐÕâ ÑÛØÖÐÙèÕÕ ßÕàÕáÕçÕÝØÕ
  crossing closest_crossing(const point3d &pWhat)const;
  friend class rays_shoot;
  friend class rays_set;
public:  
  ray(ray*,int iTo,int iKind=R_BORDER);
  // ·ÐÓÛãèÕçÚÐ ÝÐ áÛãçÐÙ ßÞßÐÔÐÝØï Ò ÞÔÝÞÜ âÕÛÕ. ºàØÒÞ, ÚÞÝÕçÝÞ :(
  ray(ray_parameters _rp,pl_raytr_model *_pm);
  ~ray();
  // ²Þ×ÒàÐéÐîâ ÚÐÚãî-ÝÐÔÞ ØÝäÞàÜÐæØî Þ ÛãçÕ
  void add_RCoefs(  DoubleAr&) const;
  void add_Ro(DoubleAr&,const pl_raytr_model * const) const;
  void add_Ksi(DoubleAr&) const;
  void add_BorderNormals(DoubleAr&,const pl_raytr_model* const,
             double A=0,double B=0,double C=0) const;
  void add_last_dxdaz(DoubleAr&) const;
  void add_last_dxdfi(DoubleAr&) const;
  void add_last_vpol (DoubleAr&) const;
  double time() const;
  double time_without_last() const;
  crossing last_crossing() const;
  crossing next_to_last_crossing()const;
  int n_of_nodes() const;
  ray_parameters ray_params() const {return rp;};
  point reflection_point()const;
  double get_refrat() const;
  int reflection_number()const;
  Logic fictive_reflection()const;
  const crossings_ptrs * crossings()const{return cr;};
  void merge(const ray*);
  const pl_raytr_model *father ()const{return pm;};
};

typedef ray* p_ray;
typedef ObjectsAr<p_ray> rays;

// ³àïÔãâ ßÕàÕÔÕÛÚØ, ßÞáÚÞÛìÚã Ûãç ÜÞÖÕâ "ÛÞÜÐâìáï" - ßàÞåÞÔØâì çÕàÕ× àÕÑàÐ...
struct ray_to_receiver{
  void load_garries(ray*rTmp,pl_raytr_model*p,double A=0,double B=0,double C=0);
  void get_cross_points(aPoints* cr_points); //nikita 17.05.05
  static unsigned globalRayCounter;
  unsigned rayId;
#if 0
  // ÇâÞ-âÞ âÐÚÞÕ ÔÞÛÖÝÞ ØÜÕâì ÜÕáâÞ Ñëâì...
  PointsAr apBreaks;
  ray_parameters_dynar arbBreakAngles;
#endif
  point crossPoints[250];//nikita 17.05.05
  int nCrossPoints;//nikita 17.05.05
  double dFi;
  double dAz;
  double dTime;
  point pReflection;
  int iReflectionBorder;  

  // ´ÐÝÝëÕ ÔÛï àÐáçÕâÐ ÐÜßÛØâãÔë
  double diver,prop;
  double dRefRat;
  double vpol[3];

  ray_to_receiver(){/*crossPoints = NULL;*/ nCrossPoints=0;};
  ray_to_receiver(const ray_to_receiver&);//inserted by nikita 26.05.05
  ray_to_receiver(int ir,ray*,pl_raytr_model *);
  //  ~ray_to_receiver(){if (crossPoints!=NULL) free(crossPoints);};
  // ´Ûï ÒÒÞÔÐ Ø× ÑÐ×ë ÔÐÝÝëå
  ray_to_receiver(double,double,double,double,double,
          const point&,int,double, double vp[3]);

  // ½Ð áÛãçÐÙ äØÚâØÒÝÞÙ ÓàÐÝØæë
  ray_to_receiver (ray*rTmp,pl_raytr_model*, double A,double B,double C);
};
#define AnyType   ray_to_receiver
#define AnyAr     ray_to_receiver_ar
#include <mix/AnyAr.H>

class rays_to_receiver:public ray_to_receiver_ar,public db_container{
protected:
  int iHRStoreRegime;
  point3d pReceiverPoint;
  Logic lReceiver;
  int N1,N2;
  void get_item (dbItem& dbiRes, dbItem& dbiRN, dbItem* dbiRP=NULL)const;
  int read_from_db (dbItem& dbiRays,dbItem& dbiRN, dbItem* dbiRP=NULL);
  int write_to_db (dbItem& dbiRays,dbItem& dbiRN, const dbItem* dbiRP=NULL)const;
  rayPathAr* ray_path_ar;
public:  
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,
           const point3d& prp,int n1,int n2,int iHRR=PM_ADD):
  db_container(dbm,sP,id){
    iHRStoreRegime=iHRR;
    lReceiver=TRUE;
    N1=n1;N2=n2;
    pReceiverPoint=prp;
    ray_path_ar=NULL;
  };
  
  // ÇâÞÑë áçØâÐâì Ø× ÑÐ×ë ÔÐÝÝëå
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,int n1,int n2,
           int iHRR=PM_ADD,void (*idle_func)()=NULL);

  const rayPathAr* get_ray_paths(){return ray_path_ar;};

  ~rays_to_receiver();

  point3d receiver()const{return pReceiverPoint;};
  int n1() const {return N1;};
  int n2() const {return N2;};
  virtual int store() const;
  virtual int load();
};


class rays_to_receivers:public ObjectsAr<rays_to_receiver*>{
public:
  rays_to_receivers(){weFormatMsg(EL_DEBUG,"Creating HitRays...");};
  /* if ~rays_to_receivers is absent then g++ 2.7.2.1 generates:
     rays.H:338: Internal compiler error.
     rays.H:338: Please submit a full bug report to `bug-g++@prep.ai.mit.edu'.
   */
  ~rays_to_receivers(){
    for(register int i=0;i<count();i++){
      delete fetch(i);
    };
    clean();
  };
  rays_to_receiver *get(int N1,int N2)const;
};

class rays_shoot{
protected:
  point3d pShotPoint;        // ºÞÞàÔØÝÐâë ØáâÞçÝØÚÐ
  pl_raytr_model *pmFather;  // ¿ÐßÐ-ÜÞÔÕÛì    
  int iShotType;             // ÂØß âÞçÚØ - ØáâÞçÝØÚ/ßàØÕÜÝØÚ
  
  int N1,N2;                 // ²ÞÕÝÝëÕ ÝÞÜÕàÐ ÔÛï ØÔÕÝâØäØÚÐæØØ ²ÞÛÞÔÕÙ
  
  // ¿ÞÛãçØâì Ûãç (Ø ÒáÕ ÕÓÞ ßÕàÕáÕçÕÝØï)
  ray* shoot_out(const ray_parameters &rp,int iKind=PM_FULL)const; 
  ray* a_shot (const double& dFi,const double& dAz,ray_specifics& rsSpec,
           ray*r=NULL)const;
  int iUWKind;
  int iDWKind;
public:
  rays_shoot(const point3d& p,int,int,int,pl_raytr_model *);
  // ´Ûï ßÞÛãçÕÝØï áÕÙáÜÞÓàÐÜÜ
  ray *raw_ray(const ray_parameters& rp) const;
  ray *raw_ray(double df,double da) const;

  // ½ÐÑÞà àÐ×ÛØçÝëå äãÝÚæØÙ ßÞßÐÔÐÝØï Ò âÞçÚã
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point(ray_parameters, const point3d&,
             int iBNumber,Logic&,double& lE,IntAr *iat,ray*r=NULL) const;
  ray* raw_ray_point(double, double,
             int iBNumber, const point3d&, Logic&,double& lE,IntAr *iat,
             ray*r=NULL) const;


  // Á ÞâàÐÖÕÝØÕÜ
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point_reflected(ray_parameters, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;
  ray* raw_ray_point_reflected(double, double, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;

  // ¾âàÐÖÕÝØÕ Þâ äØÚâØÒÝÞÙ ÓàÐÝØæë Ò âÕÛÕ á Ø×ÒÕáâÝëÜ ÝÞÜÕàÞÜ

  ray* raw_ray_fict_bord(double, double, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  ray* raw_ray_fict_bord(ray_parameters, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  point sp()const{return pShotPoint;};  
  int n1()const{return N1;};
  int n2()const{return N2;};
  void set_uw_kind(int ik){iUWKind=ik;};
  void set_dw_kind(int ik){iDWKind=ik;};
  int uw_kind(){return iUWKind;};
  int dw_kind(){return iDWKind;};
};

// »îÜßÕÝë...
Logic equal_paths_to_border (ray* r1, ray*r2, int iB);
Logic equal_paths_to_node (ray* r1, ray*r2, int iN);


#endif /* rays.H */
/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 11-Jan-1999  */
/*          ...          */

   àÐ×àÕ×ÐÕâáï, Ø Ò ÚÞâÞàÞÙ åàÐÝïâáï ÒáÕ ÛãçØ, ÒëèÕÔèØÕ Ø× ØáâÞçÝØÚÐ.
 */


#include <ldb/dbtools.h>
#include <mix/DoubleAr.H>

#include <mix/ObjectsAr.H>
#include <mod3d/mod3d_basic.H>

#define PM_FB 1
#define PM_FULL 0

#define PM_CREATE 1
#define PM_ADD 2

// ºÛÐááë, ØáßÞÛì×ãÕÜëÕ Ò ×ÐÓÞÛÞÒÚÕclass crossings;
class crossing;
class crossings_ptrs;
class pl_raytr_model;
class pl_fict_border;
class ray_specifics;

// start of insertion by nikita
typedef struct rayPathInfo {
  aPoints* points;
  int nReflectBorder;
  int receiver;
};

#define AssocAr rayPathAr

typedef struct {
  unsigned  key;//ray identifier in database RAYPATHINFO table
  rayPathInfo   value;//ray path points and receiver, reflection border binding
}  rayPath;
#define AssocItemType      rayPath
#define AssocKeyField      key
#define AssocKeyType       unsigned
#define AssocKeysCmp(key1,key2)    COMPAR(key1,key2)
#include <mix/AssocAr.H>

// end of insertion by nikita 



class ray_parameters{
protected:
  // ² àÐÔØÐÝÐå
  double dFi,dAz;
public:
  void set_pars(double,double);
  ray_parameters(){;};
  ray_parameters(double,double);
friend int operator==(const ray_parameters &rp1,
              const ray_parameters &rp2);
  double fi()const{return dFi;};
  double az()const{return dAz;};
};

#define AnyType ray_parameters
#define AnyAr ray_parameters_dynar
#include <mix/AnyAr.H>

#define R_BORDER    1
#define R_NODE      2
class ray{
protected:
  pl_raytr_model *pm;
  int iStored; // ¾ÑÞ×ÝÐçÐÕâ, Òë×ëÒÐÛØ ÛØ ÔÛï ÛãçÐ ×ÐßØáì ÕÓÞ ßÕàÕáÕçÕÝØÙ
  ray_parameters rp; // ¾âÚãÔÐ ÒëáâàÕÛØÛØ
  // ÁßØáÞÚ ãÚÐ×ÐâÕÛÕÙ ÝÐ ßÕàÕáÕçÕÝØï
  crossings_ptrs *cr;

  // ÁÞ×ÔÐÕâáï âÞÛìÚÞ Ø× rays_shoot
  ray(const ray_parameters &RP,crossings_ptrs *CR,pl_raytr_model*p);

  // ¿àÞÒÕàïÕâ, ßàÞåÞÔØâ ÛØ Ûãç çÕàÕ× ÓàÐÝØæã á ãÚÐ×ÐÝÝëÜ ÝÞÜÕàÞÜ
  Logic crosses(int i);

  // ²Þ×ÒàÐéÐÕâ ÑÛØÖÐÙèÕÕ ßÕàÕáÕçÕÝØÕ
  crossing closest_crossing(const point3d &pWhat)const;
  friend class rays_shoot;
  friend class rays_set;
public:  
  ray(ray*,int iTo,int iKind=R_BORDER);
  // ·ÐÓÛãèÕçÚÐ ÝÐ áÛãçÐÙ ßÞßÐÔÐÝØï Ò ÞÔÝÞÜ âÕÛÕ. ºàØÒÞ, ÚÞÝÕçÝÞ :(
  ray(ray_parameters _rp,pl_raytr_model *_pm);
  ~ray();
  // ²Þ×ÒàÐéÐîâ ÚÐÚãî-ÝÐÔÞ ØÝäÞàÜÐæØî Þ ÛãçÕ
  void add_RCoefs(  DoubleAr&) const;
  void add_Ro(DoubleAr&,const pl_raytr_model * const) const;
  void add_Ksi(DoubleAr&) const;
  void add_BorderNormals(DoubleAr&,const pl_raytr_model* const,
             double A=0,double B=0,double C=0) const;
  void add_last_dxdaz(DoubleAr&) const;
  void add_last_dxdfi(DoubleAr&) const;
  void add_last_vpol (DoubleAr&) const;
  double time() const;
  double time_without_last() const;
  crossing last_crossing() const;
  crossing next_to_last_crossing()const;
  int n_of_nodes() const;
  ray_parameters ray_params() const {return rp;};
  point reflection_point()const;
  double get_refrat() const;
  int reflection_number()const;
  Logic fictive_reflection()const;
  const crossings_ptrs * crossings()const{return cr;};
  void merge(const ray*);
  const pl_raytr_model *father ()const{return pm;};
};

typedef ray* p_ray;
typedef ObjectsAr<p_ray> rays;

// ³àïÔãâ ßÕàÕÔÕÛÚØ, ßÞáÚÞÛìÚã Ûãç ÜÞÖÕâ "ÛÞÜÐâìáï" - ßàÞåÞÔØâì çÕàÕ× àÕÑàÐ...
struct ray_to_receiver{
  void load_garries(ray*rTmp,pl_raytr_model*p,double A=0,double B=0,double C=0);
  void get_cross_points(aPoints* cr_points); //nikita 17.05.05
  static unsigned globalRayCounter;
  unsigned rayId;
#if 0
  // ÇâÞ-âÞ âÐÚÞÕ ÔÞÛÖÝÞ ØÜÕâì ÜÕáâÞ Ñëâì...
  PointsAr apBreaks;
  ray_parameters_dynar arbBreakAngles;
#endif
  point crossPoints[250];//nikita 17.05.05
  int nCrossPoints;//nikita 17.05.05
  double dFi;
  double dAz;
  double dTime;
  point pReflection;
  int iReflectionBorder;  

  // ´ÐÝÝëÕ ÔÛï àÐáçÕâÐ ÐÜßÛØâãÔë
  double diver,prop;
  double dRefRat;
  double vpol[3];

  ray_to_receiver(){/*crossPoints = NULL;*/ nCrossPoints=0;};
  ray_to_receiver(const ray_to_receiver&);//inserted by nikita 26.05.05
  ray_to_receiver(int ir,ray*,pl_raytr_model *);
  //  ~ray_to_receiver(){if (crossPoints!=NULL) free(crossPoints);};
  // ´Ûï ÒÒÞÔÐ Ø× ÑÐ×ë ÔÐÝÝëå
  ray_to_receiver(double,double,double,double,double,
          const point&,int,double, double vp[3]);

  // ½Ð áÛãçÐÙ äØÚâØÒÝÞÙ ÓàÐÝØæë
  ray_to_receiver (ray*rTmp,pl_raytr_model*, double A,double B,double C);
};
#define AnyType   ray_to_receiver
#define AnyAr     ray_to_receiver_ar
#include <mix/AnyAr.H>

class rays_to_receiver:public ray_to_receiver_ar,public db_container{
protected:
  int iHRStoreRegime;
  point3d pReceiverPoint;
  Logic lReceiver;
  int N1,N2;
  void get_item (dbItem& dbiRes, dbItem& dbiRN, dbItem* dbiRP=NULL)const;
  int read_from_db (dbItem& dbiRays,dbItem& dbiRN, dbItem* dbiRP=NULL);
  int write_to_db (dbItem& dbiRays,dbItem& dbiRN, const dbItem* dbiRP=NULL)const;
  rayPathAr* ray_path_ar;
public:  
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,
           const point3d& prp,int n1,int n2,int iHRR=PM_ADD):
  db_container(dbm,sP,id){
    iHRStoreRegime=iHRR;
    lReceiver=TRUE;
    N1=n1;N2=n2;
    pReceiverPoint=prp;
    ray_path_ar=NULL;
  };
  
  // ÇâÞÑë áçØâÐâì Ø× ÑÐ×ë ÔÐÝÝëå
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,int n1,int n2,
           int iHRR=PM_ADD,void (*idle_func)()=NULL);

  const rayPathAr* get_ray_paths(){return ray_path_ar;};

  ~rays_to_receiver();

  point3d receiver()const{return pReceiverPoint;};
  int n1() const {return N1;};
  int n2() const {return N2;};
  virtual int store() const;
  virtual int load();
};


class rays_to_receivers:public ObjectsAr<rays_to_receiver*>{
public:
  rays_to_receivers(){weFormatMsg(EL_DEBUG,"Creating HitRays...");};
  /* if ~rays_to_receivers is absent then g++ 2.7.2.1 generates:
     rays.H:338: Internal compiler error.
     rays.H:338: Please submit a full bug report to `bug-g++@prep.ai.mit.edu'.
   */
  ~rays_to_receivers(){
    for(register int i=0;i<count();i++){
      delete fetch(i);
    };
    clean();
  };
  rays_to_receiver *get(int N1,int N2)const;
};

class rays_shoot{
protected:
  point3d pShotPoint;        // ºÞÞàÔØÝÐâë ØáâÞçÝØÚÐ
  pl_raytr_model *pmFather;  // ¿ÐßÐ-ÜÞÔÕÛì    
  int iShotType;             // ÂØß âÞçÚØ - ØáâÞçÝØÚ/ßàØÕÜÝØÚ
  
  int N1,N2;                 // ²ÞÕÝÝëÕ ÝÞÜÕàÐ ÔÛï ØÔÕÝâØäØÚÐæØØ ²ÞÛÞÔÕÙ
  
  // ¿ÞÛãçØâì Ûãç (Ø ÒáÕ ÕÓÞ ßÕàÕáÕçÕÝØï)
  ray* shoot_out(const ray_parameters &rp,int iKind=PM_FULL)const; 
  ray* a_shot (const double& dFi,const double& dAz,ray_specifics& rsSpec,
           ray*r=NULL)const;
  int iUWKind;
  int iDWKind;
public:
  rays_shoot(const point3d& p,int,int,int,pl_raytr_model *);
  // ´Ûï ßÞÛãçÕÝØï áÕÙáÜÞÓàÐÜÜ
  ray *raw_ray(const ray_parameters& rp) const;
  ray *raw_ray(double df,double da) const;

  // ½ÐÑÞà àÐ×ÛØçÝëå äãÝÚæØÙ ßÞßÐÔÐÝØï Ò âÞçÚã
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point(ray_parameters, const point3d&,
             int iBNumber,Logic&,double& lE,IntAr *iat,ray*r=NULL) const;
  ray* raw_ray_point(double, double,
             int iBNumber, const point3d&, Logic&,double& lE,IntAr *iat,
             ray*r=NULL) const;


  // Á ÞâàÐÖÕÝØÕÜ
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point_reflected(ray_parameters, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;
  ray* raw_ray_point_reflected(double, double, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;

  // ¾âàÐÖÕÝØÕ Þâ äØÚâØÒÝÞÙ ÓàÐÝØæë Ò âÕÛÕ á Ø×ÒÕáâÝëÜ ÝÞÜÕàÞÜ

  ray* raw_ray_fict_bord(double, double, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  ray* raw_ray_fict_bord(ray_parameters, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  point sp()const{return pShotPoint;};  
  int n1()const{return N1;};
  int n2()const{return N2;};
  void set_uw_kind(int ik){iUWKind=ik;};
  void set_dw_kind(int ik){iDWKind=ik;};
  int uw_kind(){return iUWKind;};
  int dw_kind(){return iDWKind;};
};

// »îÜßÕÝë...
Logic equal_paths_to_border (ray* r1, ray*r2, int iB);
Logic equal_paths_to_node (ray* r1, ray*r2, int iN);


#endif /* rays.H */
/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 11-Jan-1999  */
/*          ...          */

class crossings;
class crossing;
class crossings_ptrs;
class pl_raytr_model;
class pl_fict_border;
class ray_specifics;

// start of insertion by nikita
typedef struct rayPathInfo {
  aPoints* points;
  int nReflectBorder;
  int receiver;
};

#define AssocAr rayPathAr

typedef struct {
  unsigned  key;//ray identifier in database RAYPATHINFO table
  rayPathInfo   value;//ray path points and receiver, reflection border binding
}  rayPath;
#define AssocItemType      rayPath
#define AssocKeyField      key
#define AssocKeyType       unsigned
#define AssocKeysCmp(key1,key2)    COMPAR(key1,key2)
#include <mix/AssocAr.H>

// end of insertion by nikita 



class ray_parameters{
protected:
  // ² àÐÔØÐÝÐå  double dFi,dAz;
public:
  void set_pars(double,double);
  ray_parameters(){;};
  ray_parameters(double,double);
friend int operator==(const ray_parameters &rp1,
              const ray_parameters &rp2);
  double fi()const{return dFi;};
  double az()const{return dAz;};
};

#define AnyType ray_parameters
#define AnyAr ray_parameters_dynar
#include <mix/AnyAr.H>

#define R_BORDER    1
#define R_NODE      2
class ray{
protected:
  pl_raytr_model *pm;
  int iStored; // ¾ÑÞ×ÝÐçÐÕâ, Òë×ëÒÐÛØ ÛØ ÔÛï ÛãçÐ ×ÐßØáì ÕÓÞ ßÕàÕáÕçÕÝØÙ
  ray_parameters rp; // ¾âÚãÔÐ ÒëáâàÕÛØÛØ
  // ÁßØáÞÚ ãÚÐ×ÐâÕÛÕÙ ÝÐ ßÕàÕáÕçÕÝØï
  crossings_ptrs *cr;

  // ÁÞ×ÔÐÕâáï âÞÛìÚÞ Ø× rays_shoot
  ray(const ray_parameters &RP,crossings_ptrs *CR,pl_raytr_model*p);

  // ¿àÞÒÕàïÕâ, ßàÞåÞÔØâ ÛØ Ûãç çÕàÕ× ÓàÐÝØæã á ãÚÐ×ÐÝÝëÜ ÝÞÜÕàÞÜ
  Logic crosses(int i);

  // ²Þ×ÒàÐéÐÕâ ÑÛØÖÐÙèÕÕ ßÕàÕáÕçÕÝØÕ
  crossing closest_crossing(const point3d &pWhat)const;
  friend class rays_shoot;
  friend class rays_set;
public:  
  ray(ray*,int iTo,int iKind=R_BORDER);
  // ·ÐÓÛãèÕçÚÐ ÝÐ áÛãçÐÙ ßÞßÐÔÐÝØï Ò ÞÔÝÞÜ âÕÛÕ. ºàØÒÞ, ÚÞÝÕçÝÞ :(
  ray(ray_parameters _rp,pl_raytr_model *_pm);
  ~ray();
  // ²Þ×ÒàÐéÐîâ ÚÐÚãî-ÝÐÔÞ ØÝäÞàÜÐæØî Þ ÛãçÕ
  void add_RCoefs(  DoubleAr&) const;
  void add_Ro(DoubleAr&,const pl_raytr_model * const) const;
  void add_Ksi(DoubleAr&) const;
  void add_BorderNormals(DoubleAr&,const pl_raytr_model* const,
             double A=0,double B=0,double C=0) const;
  void add_last_dxdaz(DoubleAr&) const;
  void add_last_dxdfi(DoubleAr&) const;
  void add_last_vpol (DoubleAr&) const;
  double time() const;
  double time_without_last() const;
  crossing last_crossing() const;
  crossing next_to_last_crossing()const;
  int n_of_nodes() const;
  ray_parameters ray_params() const {return rp;};
  point reflection_point()const;
  double get_refrat() const;
  int reflection_number()const;
  Logic fictive_reflection()const;
  const crossings_ptrs * crossings()const{return cr;};
  void merge(const ray*);
  const pl_raytr_model *father ()const{return pm;};
};

typedef ray* p_ray;
typedef ObjectsAr<p_ray> rays;

// ³àïÔãâ ßÕàÕÔÕÛÚØ, ßÞáÚÞÛìÚã Ûãç ÜÞÖÕâ "ÛÞÜÐâìáï" - ßàÞåÞÔØâì çÕàÕ× àÕÑàÐ...
struct ray_to_receiver{
  void load_garries(ray*rTmp,pl_raytr_model*p,double A=0,double B=0,double C=0);
  void get_cross_points(aPoints* cr_points); //nikita 17.05.05
  static unsigned globalRayCounter;
  unsigned rayId;
#if 0
  // ÇâÞ-âÞ âÐÚÞÕ ÔÞÛÖÝÞ ØÜÕâì ÜÕáâÞ Ñëâì...
  PointsAr apBreaks;
  ray_parameters_dynar arbBreakAngles;
#endif
  point crossPoints[250];//nikita 17.05.05
  int nCrossPoints;//nikita 17.05.05
  double dFi;
  double dAz;
  double dTime;
  point pReflection;
  int iReflectionBorder;  

  // ´ÐÝÝëÕ ÔÛï àÐáçÕâÐ ÐÜßÛØâãÔë
  double diver,prop;
  double dRefRat;
  double vpol[3];

  ray_to_receiver(){/*crossPoints = NULL;*/ nCrossPoints=0;};
  ray_to_receiver(const ray_to_receiver&);//inserted by nikita 26.05.05
  ray_to_receiver(int ir,ray*,pl_raytr_model *);
  //  ~ray_to_receiver(){if (crossPoints!=NULL) free(crossPoints);};
  // ´Ûï ÒÒÞÔÐ Ø× ÑÐ×ë ÔÐÝÝëå
  ray_to_receiver(double,double,double,double,double,
          const point&,int,double, double vp[3]);

  // ½Ð áÛãçÐÙ äØÚâØÒÝÞÙ ÓàÐÝØæë
  ray_to_receiver (ray*rTmp,pl_raytr_model*, double A,double B,double C);
};
#define AnyType   ray_to_receiver
#define AnyAr     ray_to_receiver_ar
#include <mix/AnyAr.H>

class rays_to_receiver:public ray_to_receiver_ar,public db_container{
protected:
  int iHRStoreRegime;
  point3d pReceiverPoint;
  Logic lReceiver;
  int N1,N2;
  void get_item (dbItem& dbiRes, dbItem& dbiRN, dbItem* dbiRP=NULL)const;
  int read_from_db (dbItem& dbiRays,dbItem& dbiRN, dbItem* dbiRP=NULL);
  int write_to_db (dbItem& dbiRays,dbItem& dbiRN, const dbItem* dbiRP=NULL)const;
  rayPathAr* ray_path_ar;
public:  
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,
           const point3d& prp,int n1,int n2,int iHRR=PM_ADD):
  db_container(dbm,sP,id){
    iHRStoreRegime=iHRR;
    lReceiver=TRUE;
    N1=n1;N2=n2;
    pReceiverPoint=prp;
    ray_path_ar=NULL;
  };
  
  // ÇâÞÑë áçØâÐâì Ø× ÑÐ×ë ÔÐÝÝëå
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,int n1,int n2,
           int iHRR=PM_ADD,void (*idle_func)()=NULL);

  const rayPathAr* get_ray_paths(){return ray_path_ar;};

  ~rays_to_receiver();

  point3d receiver()const{return pReceiverPoint;};
  int n1() const {return N1;};
  int n2() const {return N2;};
  virtual int store() const;
  virtual int load();
};


class rays_to_receivers:public ObjectsAr<rays_to_receiver*>{
public:
  rays_to_receivers(){weFormatMsg(EL_DEBUG,"Creating HitRays...");};
  /* if ~rays_to_receivers is absent then g++ 2.7.2.1 generates:
     rays.H:338: Internal compiler error.
     rays.H:338: Please submit a full bug report to `bug-g++@prep.ai.mit.edu'.
   */
  ~rays_to_receivers(){
    for(register int i=0;i<count();i++){
      delete fetch(i);
    };
    clean();
  };
  rays_to_receiver *get(int N1,int N2)const;
};

class rays_shoot{
protected:
  point3d pShotPoint;        // ºÞÞàÔØÝÐâë ØáâÞçÝØÚÐ
  pl_raytr_model *pmFather;  // ¿ÐßÐ-ÜÞÔÕÛì    
  int iShotType;             // ÂØß âÞçÚØ - ØáâÞçÝØÚ/ßàØÕÜÝØÚ
  
  int N1,N2;                 // ²ÞÕÝÝëÕ ÝÞÜÕàÐ ÔÛï ØÔÕÝâØäØÚÐæØØ ²ÞÛÞÔÕÙ
  
  // ¿ÞÛãçØâì Ûãç (Ø ÒáÕ ÕÓÞ ßÕàÕáÕçÕÝØï)
  ray* shoot_out(const ray_parameters &rp,int iKind=PM_FULL)const; 
  ray* a_shot (const double& dFi,const double& dAz,ray_specifics& rsSpec,
           ray*r=NULL)const;
  int iUWKind;
  int iDWKind;
public:
  rays_shoot(const point3d& p,int,int,int,pl_raytr_model *);
  // ´Ûï ßÞÛãçÕÝØï áÕÙáÜÞÓàÐÜÜ
  ray *raw_ray(const ray_parameters& rp) const;
  ray *raw_ray(double df,double da) const;

  // ½ÐÑÞà àÐ×ÛØçÝëå äãÝÚæØÙ ßÞßÐÔÐÝØï Ò âÞçÚã
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point(ray_parameters, const point3d&,
             int iBNumber,Logic&,double& lE,IntAr *iat,ray*r=NULL) const;
  ray* raw_ray_point(double, double,
             int iBNumber, const point3d&, Logic&,double& lE,IntAr *iat,
             ray*r=NULL) const;


  // Á ÞâàÐÖÕÝØÕÜ
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point_reflected(ray_parameters, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;
  ray* raw_ray_point_reflected(double, double, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;

  // ¾âàÐÖÕÝØÕ Þâ äØÚâØÒÝÞÙ ÓàÐÝØæë Ò âÕÛÕ á Ø×ÒÕáâÝëÜ ÝÞÜÕàÞÜ

  ray* raw_ray_fict_bord(double, double, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  ray* raw_ray_fict_bord(ray_parameters, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  point sp()const{return pShotPoint;};  
  int n1()const{return N1;};
  int n2()const{return N2;};
  void set_uw_kind(int ik){iUWKind=ik;};
  void set_dw_kind(int ik){iDWKind=ik;};
  int uw_kind(){return iUWKind;};
  int dw_kind(){return iDWKind;};
};

// »îÜßÕÝë...
Logic equal_paths_to_border (ray* r1, ray*r2, int iB);
Logic equal_paths_to_node (ray* r1, ray*r2, int iN);


#endif /* rays.H */
/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 11-Jan-1999  */
/*          ...          */

  double dFi,dAz;
public:
  void set_pars(double,double);
  ray_parameters(){;};
  ray_parameters(double,double);
friend int operator==(const ray_parameters &rp1,
              const ray_parameters &rp2);
  double fi()const{return dFi;};
  double az()const{return dAz;};
};

#define AnyType ray_parameters
#define AnyAr ray_parameters_dynar
#include <mix/AnyAr.H>

#define R_BORDER    1
#define R_NODE      2
class ray{
protected:
  pl_raytr_model *pm;
  int iStored; // ¾ÑÞ×ÝÐçÐÕâ, Òë×ëÒÐÛØ ÛØ ÔÛï ÛãçÐ ×ÐßØáì ÕÓÞ ßÕàÕáÕçÕÝØÙ
  ray_parameters rp; // ¾âÚãÔÐ ÒëáâàÕÛØÛØ  // ÁßØáÞÚ ãÚÐ×ÐâÕÛÕÙ ÝÐ ßÕàÕáÕçÕÝØï
  crossings_ptrs *cr;

  // ÁÞ×ÔÐÕâáï âÞÛìÚÞ Ø× rays_shoot
  ray(const ray_parameters &RP,crossings_ptrs *CR,pl_raytr_model*p);

  // ¿àÞÒÕàïÕâ, ßàÞåÞÔØâ ÛØ Ûãç çÕàÕ× ÓàÐÝØæã á ãÚÐ×ÐÝÝëÜ ÝÞÜÕàÞÜ
  Logic crosses(int i);

  // ²Þ×ÒàÐéÐÕâ ÑÛØÖÐÙèÕÕ ßÕàÕáÕçÕÝØÕ
  crossing closest_crossing(const point3d &pWhat)const;
  friend class rays_shoot;
  friend class rays_set;
public:  
  ray(ray*,int iTo,int iKind=R_BORDER);
  // ·ÐÓÛãèÕçÚÐ ÝÐ áÛãçÐÙ ßÞßÐÔÐÝØï Ò ÞÔÝÞÜ âÕÛÕ. ºàØÒÞ, ÚÞÝÕçÝÞ :(
  ray(ray_parameters _rp,pl_raytr_model *_pm);
  ~ray();
  // ²Þ×ÒàÐéÐîâ ÚÐÚãî-ÝÐÔÞ ØÝäÞàÜÐæØî Þ ÛãçÕ
  void add_RCoefs(  DoubleAr&) const;
  void add_Ro(DoubleAr&,const pl_raytr_model * const) const;
  void add_Ksi(DoubleAr&) const;
  void add_BorderNormals(DoubleAr&,const pl_raytr_model* const,
             double A=0,double B=0,double C=0) const;
  void add_last_dxdaz(DoubleAr&) const;
  void add_last_dxdfi(DoubleAr&) const;
  void add_last_vpol (DoubleAr&) const;
  double time() const;
  double time_without_last() const;
  crossing last_crossing() const;
  crossing next_to_last_crossing()const;
  int n_of_nodes() const;
  ray_parameters ray_params() const {return rp;};
  point reflection_point()const;
  double get_refrat() const;
  int reflection_number()const;
  Logic fictive_reflection()const;
  const crossings_ptrs * crossings()const{return cr;};
  void merge(const ray*);
  const pl_raytr_model *father ()const{return pm;};
};

typedef ray* p_ray;
typedef ObjectsAr<p_ray> rays;

// ³àïÔãâ ßÕàÕÔÕÛÚØ, ßÞáÚÞÛìÚã Ûãç ÜÞÖÕâ "ÛÞÜÐâìáï" - ßàÞåÞÔØâì çÕàÕ× àÕÑàÐ...
struct ray_to_receiver{
  void load_garries(ray*rTmp,pl_raytr_model*p,double A=0,double B=0,double C=0);
  void get_cross_points(aPoints* cr_points); //nikita 17.05.05
  static unsigned globalRayCounter;
  unsigned rayId;
#if 0
  // ÇâÞ-âÞ âÐÚÞÕ ÔÞÛÖÝÞ ØÜÕâì ÜÕáâÞ Ñëâì...
  PointsAr apBreaks;
  ray_parameters_dynar arbBreakAngles;
#endif
  point crossPoints[250];//nikita 17.05.05
  int nCrossPoints;//nikita 17.05.05
  double dFi;
  double dAz;
  double dTime;
  point pReflection;
  int iReflectionBorder;  

  // ´ÐÝÝëÕ ÔÛï àÐáçÕâÐ ÐÜßÛØâãÔë
  double diver,prop;
  double dRefRat;
  double vpol[3];

  ray_to_receiver(){/*crossPoints = NULL;*/ nCrossPoints=0;};
  ray_to_receiver(const ray_to_receiver&);//inserted by nikita 26.05.05
  ray_to_receiver(int ir,ray*,pl_raytr_model *);
  //  ~ray_to_receiver(){if (crossPoints!=NULL) free(crossPoints);};
  // ´Ûï ÒÒÞÔÐ Ø× ÑÐ×ë ÔÐÝÝëå
  ray_to_receiver(double,double,double,double,double,
          const point&,int,double, double vp[3]);

  // ½Ð áÛãçÐÙ äØÚâØÒÝÞÙ ÓàÐÝØæë
  ray_to_receiver (ray*rTmp,pl_raytr_model*, double A,double B,double C);
};
#define AnyType   ray_to_receiver
#define AnyAr     ray_to_receiver_ar
#include <mix/AnyAr.H>

class rays_to_receiver:public ray_to_receiver_ar,public db_container{
protected:
  int iHRStoreRegime;
  point3d pReceiverPoint;
  Logic lReceiver;
  int N1,N2;
  void get_item (dbItem& dbiRes, dbItem& dbiRN, dbItem* dbiRP=NULL)const;
  int read_from_db (dbItem& dbiRays,dbItem& dbiRN, dbItem* dbiRP=NULL);
  int write_to_db (dbItem& dbiRays,dbItem& dbiRN, const dbItem* dbiRP=NULL)const;
  rayPathAr* ray_path_ar;
public:  
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,
           const point3d& prp,int n1,int n2,int iHRR=PM_ADD):
  db_container(dbm,sP,id){
    iHRStoreRegime=iHRR;
    lReceiver=TRUE;
    N1=n1;N2=n2;
    pReceiverPoint=prp;
    ray_path_ar=NULL;
  };
  
  // ÇâÞÑë áçØâÐâì Ø× ÑÐ×ë ÔÐÝÝëå
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,int n1,int n2,
           int iHRR=PM_ADD,void (*idle_func)()=NULL);

  const rayPathAr* get_ray_paths(){return ray_path_ar;};

  ~rays_to_receiver();

  point3d receiver()const{return pReceiverPoint;};
  int n1() const {return N1;};
  int n2() const {return N2;};
  virtual int store() const;
  virtual int load();
};


class rays_to_receivers:public ObjectsAr<rays_to_receiver*>{
public:
  rays_to_receivers(){weFormatMsg(EL_DEBUG,"Creating HitRays...");};
  /* if ~rays_to_receivers is absent then g++ 2.7.2.1 generates:
     rays.H:338: Internal compiler error.
     rays.H:338: Please submit a full bug report to `bug-g++@prep.ai.mit.edu'.
   */
  ~rays_to_receivers(){
    for(register int i=0;i<count();i++){
      delete fetch(i);
    };
    clean();
  };
  rays_to_receiver *get(int N1,int N2)const;
};

class rays_shoot{
protected:
  point3d pShotPoint;        // ºÞÞàÔØÝÐâë ØáâÞçÝØÚÐ
  pl_raytr_model *pmFather;  // ¿ÐßÐ-ÜÞÔÕÛì    
  int iShotType;             // ÂØß âÞçÚØ - ØáâÞçÝØÚ/ßàØÕÜÝØÚ
  
  int N1,N2;                 // ²ÞÕÝÝëÕ ÝÞÜÕàÐ ÔÛï ØÔÕÝâØäØÚÐæØØ ²ÞÛÞÔÕÙ
  
  // ¿ÞÛãçØâì Ûãç (Ø ÒáÕ ÕÓÞ ßÕàÕáÕçÕÝØï)
  ray* shoot_out(const ray_parameters &rp,int iKind=PM_FULL)const; 
  ray* a_shot (const double& dFi,const double& dAz,ray_specifics& rsSpec,
           ray*r=NULL)const;
  int iUWKind;
  int iDWKind;
public:
  rays_shoot(const point3d& p,int,int,int,pl_raytr_model *);
  // ´Ûï ßÞÛãçÕÝØï áÕÙáÜÞÓàÐÜÜ
  ray *raw_ray(const ray_parameters& rp) const;
  ray *raw_ray(double df,double da) const;

  // ½ÐÑÞà àÐ×ÛØçÝëå äãÝÚæØÙ ßÞßÐÔÐÝØï Ò âÞçÚã
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point(ray_parameters, const point3d&,
             int iBNumber,Logic&,double& lE,IntAr *iat,ray*r=NULL) const;
  ray* raw_ray_point(double, double,
             int iBNumber, const point3d&, Logic&,double& lE,IntAr *iat,
             ray*r=NULL) const;


  // Á ÞâàÐÖÕÝØÕÜ
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point_reflected(ray_parameters, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;
  ray* raw_ray_point_reflected(double, double, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;

  // ¾âàÐÖÕÝØÕ Þâ äØÚâØÒÝÞÙ ÓàÐÝØæë Ò âÕÛÕ á Ø×ÒÕáâÝëÜ ÝÞÜÕàÞÜ

  ray* raw_ray_fict_bord(double, double, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  ray* raw_ray_fict_bord(ray_parameters, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  point sp()const{return pShotPoint;};  
  int n1()const{return N1;};
  int n2()const{return N2;};
  void set_uw_kind(int ik){iUWKind=ik;};
  void set_dw_kind(int ik){iDWKind=ik;};
  int uw_kind(){return iUWKind;};
  int dw_kind(){return iDWKind;};
};

// »îÜßÕÝë...
Logic equal_paths_to_border (ray* r1, ray*r2, int iB);
Logic equal_paths_to_node (ray* r1, ray*r2, int iN);


#endif /* rays.H */
/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 11-Jan-1999  */
/*          ...          */

  // ÁßØáÞÚ ãÚÐ×ÐâÕÛÕÙ ÝÐ ßÕàÕáÕçÕÝØï
  crossings_ptrs *cr;

  // ÁÞ×ÔÐÕâáï âÞÛìÚÞ Ø× rays_shoot
  ray(const ray_parameters &RP,crossings_ptrs *CR,pl_raytr_model*p);

  // ¿àÞÒÕàïÕâ, ßàÞåÞÔØâ ÛØ Ûãç çÕàÕ× ÓàÐÝØæã á ãÚÐ×ÐÝÝëÜ ÝÞÜÕàÞÜ
  Logic crosses(int i);

  // ²Þ×ÒàÐéÐÕâ ÑÛØÖÐÙèÕÕ ßÕàÕáÕçÕÝØÕ
  crossing closest_crossing(const point3d &pWhat)const;
  friend class rays_shoot;
  friend class rays_set;
public:  
  ray(ray*,int iTo,int iKind=R_BORDER);
  // ·ÐÓÛãèÕçÚÐ ÝÐ áÛãçÐÙ ßÞßÐÔÐÝØï Ò ÞÔÝÞÜ âÕÛÕ. ºàØÒÞ, ÚÞÝÕçÝÞ :(
  ray(ray_parameters _rp,pl_raytr_model *_pm);
  ~ray();
  // ²Þ×ÒàÐéÐîâ ÚÐÚãî-ÝÐÔÞ ØÝäÞàÜÐæØî Þ ÛãçÕ  void add_RCoefs(  DoubleAr&) const;
  void add_Ro(DoubleAr&,const pl_raytr_model * const) const;
  void add_Ksi(DoubleAr&) const;
  void add_BorderNormals(DoubleAr&,const pl_raytr_model* const,
             double A=0,double B=0,double C=0) const;
  void add_last_dxdaz(DoubleAr&) const;
  void add_last_dxdfi(DoubleAr&) const;
  void add_last_vpol (DoubleAr&) const;
  double time() const;
  double time_without_last() const;
  crossing last_crossing() const;
  crossing next_to_last_crossing()const;
  int n_of_nodes() const;
  ray_parameters ray_params() const {return rp;};
  point reflection_point()const;
  double get_refrat() const;
  int reflection_number()const;
  Logic fictive_reflection()const;
  const crossings_ptrs * crossings()const{return cr;};
  void merge(const ray*);
  const pl_raytr_model *father ()const{return pm;};
};

typedef ray* p_ray;
typedef ObjectsAr<p_ray> rays;

// ³àïÔãâ ßÕàÕÔÕÛÚØ, ßÞáÚÞÛìÚã Ûãç ÜÞÖÕâ "ÛÞÜÐâìáï" - ßàÞåÞÔØâì çÕàÕ× àÕÑàÐ...
struct ray_to_receiver{
  void load_garries(ray*rTmp,pl_raytr_model*p,double A=0,double B=0,double C=0);
  void get_cross_points(aPoints* cr_points); //nikita 17.05.05
  static unsigned globalRayCounter;
  unsigned rayId;
#if 0
  // ÇâÞ-âÞ âÐÚÞÕ ÔÞÛÖÝÞ ØÜÕâì ÜÕáâÞ Ñëâì...
  PointsAr apBreaks;
  ray_parameters_dynar arbBreakAngles;
#endif
  point crossPoints[250];//nikita 17.05.05
  int nCrossPoints;//nikita 17.05.05
  double dFi;
  double dAz;
  double dTime;
  point pReflection;
  int iReflectionBorder;  

  // ´ÐÝÝëÕ ÔÛï àÐáçÕâÐ ÐÜßÛØâãÔë
  double diver,prop;
  double dRefRat;
  double vpol[3];

  ray_to_receiver(){/*crossPoints = NULL;*/ nCrossPoints=0;};
  ray_to_receiver(const ray_to_receiver&);//inserted by nikita 26.05.05
  ray_to_receiver(int ir,ray*,pl_raytr_model *);
  //  ~ray_to_receiver(){if (crossPoints!=NULL) free(crossPoints);};
  // ´Ûï ÒÒÞÔÐ Ø× ÑÐ×ë ÔÐÝÝëå
  ray_to_receiver(double,double,double,double,double,
          const point&,int,double, double vp[3]);

  // ½Ð áÛãçÐÙ äØÚâØÒÝÞÙ ÓàÐÝØæë
  ray_to_receiver (ray*rTmp,pl_raytr_model*, double A,double B,double C);
};
#define AnyType   ray_to_receiver
#define AnyAr     ray_to_receiver_ar
#include <mix/AnyAr.H>

class rays_to_receiver:public ray_to_receiver_ar,public db_container{
protected:
  int iHRStoreRegime;
  point3d pReceiverPoint;
  Logic lReceiver;
  int N1,N2;
  void get_item (dbItem& dbiRes, dbItem& dbiRN, dbItem* dbiRP=NULL)const;
  int read_from_db (dbItem& dbiRays,dbItem& dbiRN, dbItem* dbiRP=NULL);
  int write_to_db (dbItem& dbiRays,dbItem& dbiRN, const dbItem* dbiRP=NULL)const;
  rayPathAr* ray_path_ar;
public:  
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,
           const point3d& prp,int n1,int n2,int iHRR=PM_ADD):
  db_container(dbm,sP,id){
    iHRStoreRegime=iHRR;
    lReceiver=TRUE;
    N1=n1;N2=n2;
    pReceiverPoint=prp;
    ray_path_ar=NULL;
  };
  
  // ÇâÞÑë áçØâÐâì Ø× ÑÐ×ë ÔÐÝÝëå
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,int n1,int n2,
           int iHRR=PM_ADD,void (*idle_func)()=NULL);

  const rayPathAr* get_ray_paths(){return ray_path_ar;};

  ~rays_to_receiver();

  point3d receiver()const{return pReceiverPoint;};
  int n1() const {return N1;};
  int n2() const {return N2;};
  virtual int store() const;
  virtual int load();
};


class rays_to_receivers:public ObjectsAr<rays_to_receiver*>{
public:
  rays_to_receivers(){weFormatMsg(EL_DEBUG,"Creating HitRays...");};
  /* if ~rays_to_receivers is absent then g++ 2.7.2.1 generates:
     rays.H:338: Internal compiler error.
     rays.H:338: Please submit a full bug report to `bug-g++@prep.ai.mit.edu'.
   */
  ~rays_to_receivers(){
    for(register int i=0;i<count();i++){
      delete fetch(i);
    };
    clean();
  };
  rays_to_receiver *get(int N1,int N2)const;
};

class rays_shoot{
protected:
  point3d pShotPoint;        // ºÞÞàÔØÝÐâë ØáâÞçÝØÚÐ
  pl_raytr_model *pmFather;  // ¿ÐßÐ-ÜÞÔÕÛì    
  int iShotType;             // ÂØß âÞçÚØ - ØáâÞçÝØÚ/ßàØÕÜÝØÚ
  
  int N1,N2;                 // ²ÞÕÝÝëÕ ÝÞÜÕàÐ ÔÛï ØÔÕÝâØäØÚÐæØØ ²ÞÛÞÔÕÙ
  
  // ¿ÞÛãçØâì Ûãç (Ø ÒáÕ ÕÓÞ ßÕàÕáÕçÕÝØï)
  ray* shoot_out(const ray_parameters &rp,int iKind=PM_FULL)const; 
  ray* a_shot (const double& dFi,const double& dAz,ray_specifics& rsSpec,
           ray*r=NULL)const;
  int iUWKind;
  int iDWKind;
public:
  rays_shoot(const point3d& p,int,int,int,pl_raytr_model *);
  // ´Ûï ßÞÛãçÕÝØï áÕÙáÜÞÓàÐÜÜ
  ray *raw_ray(const ray_parameters& rp) const;
  ray *raw_ray(double df,double da) const;

  // ½ÐÑÞà àÐ×ÛØçÝëå äãÝÚæØÙ ßÞßÐÔÐÝØï Ò âÞçÚã
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point(ray_parameters, const point3d&,
             int iBNumber,Logic&,double& lE,IntAr *iat,ray*r=NULL) const;
  ray* raw_ray_point(double, double,
             int iBNumber, const point3d&, Logic&,double& lE,IntAr *iat,
             ray*r=NULL) const;


  // Á ÞâàÐÖÕÝØÕÜ
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point_reflected(ray_parameters, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;
  ray* raw_ray_point_reflected(double, double, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;

  // ¾âàÐÖÕÝØÕ Þâ äØÚâØÒÝÞÙ ÓàÐÝØæë Ò âÕÛÕ á Ø×ÒÕáâÝëÜ ÝÞÜÕàÞÜ

  ray* raw_ray_fict_bord(double, double, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  ray* raw_ray_fict_bord(ray_parameters, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  point sp()const{return pShotPoint;};  
  int n1()const{return N1;};
  int n2()const{return N2;};
  void set_uw_kind(int ik){iUWKind=ik;};
  void set_dw_kind(int ik){iDWKind=ik;};
  int uw_kind(){return iUWKind;};
  int dw_kind(){return iDWKind;};
};

// »îÜßÕÝë...
Logic equal_paths_to_border (ray* r1, ray*r2, int iB);
Logic equal_paths_to_node (ray* r1, ray*r2, int iN);


#endif /* rays.H */
/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 11-Jan-1999  */
/*          ...          */

  void add_RCoefs(  DoubleAr&) const;
  void add_Ro(DoubleAr&,const pl_raytr_model * const) const;
  void add_Ksi(DoubleAr&) const;
  void add_BorderNormals(DoubleAr&,const pl_raytr_model* const,
             double A=0,double B=0,double C=0) const;
  void add_last_dxdaz(DoubleAr&) const;
  void add_last_dxdfi(DoubleAr&) const;
  void add_last_vpol (DoubleAr&) const;
  double time() const;
  double time_without_last() const;
  crossing last_crossing() const;
  crossing next_to_last_crossing()const;
  int n_of_nodes() const;
  ray_parameters ray_params() const {return rp;};
  point reflection_point()const;
  double get_refrat() const;
  int reflection_number()const;
  Logic fictive_reflection()const;
  const crossings_ptrs * crossings()const{return cr;};
  void merge(const ray*);
  const pl_raytr_model *father ()const{return pm;};
};

typedef ray* p_ray;
typedef ObjectsAr<p_ray> rays;

// ³àïÔãâ ßÕàÕÔÕÛÚØ, ßÞáÚÞÛìÚã Ûãç ÜÞÖÕâ "ÛÞÜÐâìáï" - ßàÞåÞÔØâì çÕàÕ× àÕÑàÐ...
struct ray_to_receiver{
  void load_garries(ray*rTmp,pl_raytr_model*p,double A=0,double B=0,double C=0);
  void get_cross_points(aPoints* cr_points); //nikita 17.05.05
  static unsigned globalRayCounter;
  unsigned rayId;
#if 0
  // ÇâÞ-âÞ âÐÚÞÕ ÔÞÛÖÝÞ ØÜÕâì ÜÕáâÞ Ñëâì...
  PointsAr apBreaks;
  ray_parameters_dynar arbBreakAngles;
#endif
  point crossPoints[250];//nikita 17.05.05
  int nCrossPoints;//nikita 17.05.05
  double dFi;
  double dAz;
  double dTime;
  point pReflection;
  int iReflectionBorder;  

  // ´ÐÝÝëÕ ÔÛï àÐáçÕâÐ ÐÜßÛØâãÔë
  double diver,prop;
  double dRefRat;
  double vpol[3];

  ray_to_receiver(){/*crossPoints = NULL;*/ nCrossPoints=0;};
  ray_to_receiver(const ray_to_receiver&);//inserted by nikita 26.05.05
  ray_to_receiver(int ir,ray*,pl_raytr_model *);
  //  ~ray_to_receiver(){if (crossPoints!=NULL) free(crossPoints);};
  // ´Ûï ÒÒÞÔÐ Ø× ÑÐ×ë ÔÐÝÝëå  ray_to_receiver(double,double,double,double,double,
          const point&,int,double, double vp[3]);

  // ½Ð áÛãçÐÙ äØÚâØÒÝÞÙ ÓàÐÝØæë
  ray_to_receiver (ray*rTmp,pl_raytr_model*, double A,double B,double C);
};
#define AnyType   ray_to_receiver
#define AnyAr     ray_to_receiver_ar
#include <mix/AnyAr.H>

class rays_to_receiver:public ray_to_receiver_ar,public db_container{
protected:
  int iHRStoreRegime;
  point3d pReceiverPoint;
  Logic lReceiver;
  int N1,N2;
  void get_item (dbItem& dbiRes, dbItem& dbiRN, dbItem* dbiRP=NULL)const;
  int read_from_db (dbItem& dbiRays,dbItem& dbiRN, dbItem* dbiRP=NULL);
  int write_to_db (dbItem& dbiRays,dbItem& dbiRN, const dbItem* dbiRP=NULL)const;
  rayPathAr* ray_path_ar;
public:  
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,
           const point3d& prp,int n1,int n2,int iHRR=PM_ADD):
  db_container(dbm,sP,id){
    iHRStoreRegime=iHRR;
    lReceiver=TRUE;
    N1=n1;N2=n2;
    pReceiverPoint=prp;
    ray_path_ar=NULL;
  };
  
  // ÇâÞÑë áçØâÐâì Ø× ÑÐ×ë ÔÐÝÝëå
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,int n1,int n2,
           int iHRR=PM_ADD,void (*idle_func)()=NULL);

  const rayPathAr* get_ray_paths(){return ray_path_ar;};

  ~rays_to_receiver();

  point3d receiver()const{return pReceiverPoint;};
  int n1() const {return N1;};
  int n2() const {return N2;};
  virtual int store() const;
  virtual int load();
};


class rays_to_receivers:public ObjectsAr<rays_to_receiver*>{
public:
  rays_to_receivers(){weFormatMsg(EL_DEBUG,"Creating HitRays...");};
  /* if ~rays_to_receivers is absent then g++ 2.7.2.1 generates:
     rays.H:338: Internal compiler error.
     rays.H:338: Please submit a full bug report to `bug-g++@prep.ai.mit.edu'.
   */
  ~rays_to_receivers(){
    for(register int i=0;i<count();i++){
      delete fetch(i);
    };
    clean();
  };
  rays_to_receiver *get(int N1,int N2)const;
};

class rays_shoot{
protected:
  point3d pShotPoint;        // ºÞÞàÔØÝÐâë ØáâÞçÝØÚÐ
  pl_raytr_model *pmFather;  // ¿ÐßÐ-ÜÞÔÕÛì    
  int iShotType;             // ÂØß âÞçÚØ - ØáâÞçÝØÚ/ßàØÕÜÝØÚ
  
  int N1,N2;                 // ²ÞÕÝÝëÕ ÝÞÜÕàÐ ÔÛï ØÔÕÝâØäØÚÐæØØ ²ÞÛÞÔÕÙ
  
  // ¿ÞÛãçØâì Ûãç (Ø ÒáÕ ÕÓÞ ßÕàÕáÕçÕÝØï)
  ray* shoot_out(const ray_parameters &rp,int iKind=PM_FULL)const; 
  ray* a_shot (const double& dFi,const double& dAz,ray_specifics& rsSpec,
           ray*r=NULL)const;
  int iUWKind;
  int iDWKind;
public:
  rays_shoot(const point3d& p,int,int,int,pl_raytr_model *);
  // ´Ûï ßÞÛãçÕÝØï áÕÙáÜÞÓàÐÜÜ
  ray *raw_ray(const ray_parameters& rp) const;
  ray *raw_ray(double df,double da) const;

  // ½ÐÑÞà àÐ×ÛØçÝëå äãÝÚæØÙ ßÞßÐÔÐÝØï Ò âÞçÚã
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point(ray_parameters, const point3d&,
             int iBNumber,Logic&,double& lE,IntAr *iat,ray*r=NULL) const;
  ray* raw_ray_point(double, double,
             int iBNumber, const point3d&, Logic&,double& lE,IntAr *iat,
             ray*r=NULL) const;


  // Á ÞâàÐÖÕÝØÕÜ
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point_reflected(ray_parameters, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;
  ray* raw_ray_point_reflected(double, double, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;

  // ¾âàÐÖÕÝØÕ Þâ äØÚâØÒÝÞÙ ÓàÐÝØæë Ò âÕÛÕ á Ø×ÒÕáâÝëÜ ÝÞÜÕàÞÜ

  ray* raw_ray_fict_bord(double, double, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  ray* raw_ray_fict_bord(ray_parameters, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  point sp()const{return pShotPoint;};  
  int n1()const{return N1;};
  int n2()const{return N2;};
  void set_uw_kind(int ik){iUWKind=ik;};
  void set_dw_kind(int ik){iDWKind=ik;};
  int uw_kind(){return iUWKind;};
  int dw_kind(){return iDWKind;};
};

// »îÜßÕÝë...
Logic equal_paths_to_border (ray* r1, ray*r2, int iB);
Logic equal_paths_to_node (ray* r1, ray*r2, int iN);


#endif /* rays.H */
/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 11-Jan-1999  */
/*          ...          */

  ray_to_receiver(double,double,double,double,double,
          const point&,int,double, double vp[3]);

  // ½Ð áÛãçÐÙ äØÚâØÒÝÞÙ ÓàÐÝØæë  ray_to_receiver (ray*rTmp,pl_raytr_model*, double A,double B,double C);
};
#define AnyType   ray_to_receiver
#define AnyAr     ray_to_receiver_ar
#include <mix/AnyAr.H>

class rays_to_receiver:public ray_to_receiver_ar,public db_container{
protected:
  int iHRStoreRegime;
  point3d pReceiverPoint;
  Logic lReceiver;
  int N1,N2;
  void get_item (dbItem& dbiRes, dbItem& dbiRN, dbItem* dbiRP=NULL)const;
  int read_from_db (dbItem& dbiRays,dbItem& dbiRN, dbItem* dbiRP=NULL);
  int write_to_db (dbItem& dbiRays,dbItem& dbiRN, const dbItem* dbiRP=NULL)const;
  rayPathAr* ray_path_ar;
public:  
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,
           const point3d& prp,int n1,int n2,int iHRR=PM_ADD):
  db_container(dbm,sP,id){
    iHRStoreRegime=iHRR;
    lReceiver=TRUE;
    N1=n1;N2=n2;
    pReceiverPoint=prp;
    ray_path_ar=NULL;
  };
  
  // ÇâÞÑë áçØâÐâì Ø× ÑÐ×ë ÔÐÝÝëå
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,int n1,int n2,
           int iHRR=PM_ADD,void (*idle_func)()=NULL);

  const rayPathAr* get_ray_paths(){return ray_path_ar;};

  ~rays_to_receiver();

  point3d receiver()const{return pReceiverPoint;};
  int n1() const {return N1;};
  int n2() const {return N2;};
  virtual int store() const;
  virtual int load();
};


class rays_to_receivers:public ObjectsAr<rays_to_receiver*>{
public:
  rays_to_receivers(){weFormatMsg(EL_DEBUG,"Creating HitRays...");};
  /* if ~rays_to_receivers is absent then g++ 2.7.2.1 generates:
     rays.H:338: Internal compiler error.
     rays.H:338: Please submit a full bug report to `bug-g++@prep.ai.mit.edu'.
   */
  ~rays_to_receivers(){
    for(register int i=0;i<count();i++){
      delete fetch(i);
    };
    clean();
  };
  rays_to_receiver *get(int N1,int N2)const;
};

class rays_shoot{
protected:
  point3d pShotPoint;        // ºÞÞàÔØÝÐâë ØáâÞçÝØÚÐ
  pl_raytr_model *pmFather;  // ¿ÐßÐ-ÜÞÔÕÛì    
  int iShotType;             // ÂØß âÞçÚØ - ØáâÞçÝØÚ/ßàØÕÜÝØÚ
  
  int N1,N2;                 // ²ÞÕÝÝëÕ ÝÞÜÕàÐ ÔÛï ØÔÕÝâØäØÚÐæØØ ²ÞÛÞÔÕÙ
  
  // ¿ÞÛãçØâì Ûãç (Ø ÒáÕ ÕÓÞ ßÕàÕáÕçÕÝØï)
  ray* shoot_out(const ray_parameters &rp,int iKind=PM_FULL)const; 
  ray* a_shot (const double& dFi,const double& dAz,ray_specifics& rsSpec,
           ray*r=NULL)const;
  int iUWKind;
  int iDWKind;
public:
  rays_shoot(const point3d& p,int,int,int,pl_raytr_model *);
  // ´Ûï ßÞÛãçÕÝØï áÕÙáÜÞÓàÐÜÜ
  ray *raw_ray(const ray_parameters& rp) const;
  ray *raw_ray(double df,double da) const;

  // ½ÐÑÞà àÐ×ÛØçÝëå äãÝÚæØÙ ßÞßÐÔÐÝØï Ò âÞçÚã
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point(ray_parameters, const point3d&,
             int iBNumber,Logic&,double& lE,IntAr *iat,ray*r=NULL) const;
  ray* raw_ray_point(double, double,
             int iBNumber, const point3d&, Logic&,double& lE,IntAr *iat,
             ray*r=NULL) const;


  // Á ÞâàÐÖÕÝØÕÜ
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point_reflected(ray_parameters, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;
  ray* raw_ray_point_reflected(double, double, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;

  // ¾âàÐÖÕÝØÕ Þâ äØÚâØÒÝÞÙ ÓàÐÝØæë Ò âÕÛÕ á Ø×ÒÕáâÝëÜ ÝÞÜÕàÞÜ

  ray* raw_ray_fict_bord(double, double, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  ray* raw_ray_fict_bord(ray_parameters, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  point sp()const{return pShotPoint;};  
  int n1()const{return N1;};
  int n2()const{return N2;};
  void set_uw_kind(int ik){iUWKind=ik;};
  void set_dw_kind(int ik){iDWKind=ik;};
  int uw_kind(){return iUWKind;};
  int dw_kind(){return iDWKind;};
};

// »îÜßÕÝë...
Logic equal_paths_to_border (ray* r1, ray*r2, int iB);
Logic equal_paths_to_node (ray* r1, ray*r2, int iN);


#endif /* rays.H */
/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 11-Jan-1999  */
/*          ...          */

  ray_to_receiver (ray*rTmp,pl_raytr_model*, double A,double B,double C);
};
#define AnyType   ray_to_receiver
#define AnyAr     ray_to_receiver_ar
#include <mix/AnyAr.H>

class rays_to_receiver:public ray_to_receiver_ar,public db_container{
protected:
  int iHRStoreRegime;
  point3d pReceiverPoint;
  Logic lReceiver;
  int N1,N2;
  void get_item (dbItem& dbiRes, dbItem& dbiRN, dbItem* dbiRP=NULL)const;
  int read_from_db (dbItem& dbiRays,dbItem& dbiRN, dbItem* dbiRP=NULL);
  int write_to_db (dbItem& dbiRays,dbItem& dbiRN, const dbItem* dbiRP=NULL)const;
  rayPathAr* ray_path_ar;
public:  
  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,
           const point3d& prp,int n1,int n2,int iHRR=PM_ADD):
  db_container(dbm,sP,id){
    iHRStoreRegime=iHRR;
    lReceiver=TRUE;
    N1=n1;N2=n2;
    pReceiverPoint=prp;
    ray_path_ar=NULL;
  };
  
  // ÇâÞÑë áçØâÐâì Ø× ÑÐ×ë ÔÐÝÝëå  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,int n1,int n2,
           int iHRR=PM_ADD,void (*idle_func)()=NULL);

  const rayPathAr* get_ray_paths(){return ray_path_ar;};

  ~rays_to_receiver();

  point3d receiver()const{return pReceiverPoint;};
  int n1() const {return N1;};
  int n2() const {return N2;};
  virtual int store() const;
  virtual int load();
};


class rays_to_receivers:public ObjectsAr<rays_to_receiver*>{
public:
  rays_to_receivers(){weFormatMsg(EL_DEBUG,"Creating HitRays...");};
  /* if ~rays_to_receivers is absent then g++ 2.7.2.1 generates:
     rays.H:338: Internal compiler error.
     rays.H:338: Please submit a full bug report to `bug-g++@prep.ai.mit.edu'.
   */
  ~rays_to_receivers(){
    for(register int i=0;i<count();i++){
      delete fetch(i);
    };
    clean();
  };
  rays_to_receiver *get(int N1,int N2)const;
};

class rays_shoot{
protected:
  point3d pShotPoint;        // ºÞÞàÔØÝÐâë ØáâÞçÝØÚÐ
  pl_raytr_model *pmFather;  // ¿ÐßÐ-ÜÞÔÕÛì    
  int iShotType;             // ÂØß âÞçÚØ - ØáâÞçÝØÚ/ßàØÕÜÝØÚ
  
  int N1,N2;                 // ²ÞÕÝÝëÕ ÝÞÜÕàÐ ÔÛï ØÔÕÝâØäØÚÐæØØ ²ÞÛÞÔÕÙ
  
  // ¿ÞÛãçØâì Ûãç (Ø ÒáÕ ÕÓÞ ßÕàÕáÕçÕÝØï)
  ray* shoot_out(const ray_parameters &rp,int iKind=PM_FULL)const; 
  ray* a_shot (const double& dFi,const double& dAz,ray_specifics& rsSpec,
           ray*r=NULL)const;
  int iUWKind;
  int iDWKind;
public:
  rays_shoot(const point3d& p,int,int,int,pl_raytr_model *);
  // ´Ûï ßÞÛãçÕÝØï áÕÙáÜÞÓàÐÜÜ
  ray *raw_ray(const ray_parameters& rp) const;
  ray *raw_ray(double df,double da) const;

  // ½ÐÑÞà àÐ×ÛØçÝëå äãÝÚæØÙ ßÞßÐÔÐÝØï Ò âÞçÚã
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point(ray_parameters, const point3d&,
             int iBNumber,Logic&,double& lE,IntAr *iat,ray*r=NULL) const;
  ray* raw_ray_point(double, double,
             int iBNumber, const point3d&, Logic&,double& lE,IntAr *iat,
             ray*r=NULL) const;


  // Á ÞâàÐÖÕÝØÕÜ
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point_reflected(ray_parameters, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;
  ray* raw_ray_point_reflected(double, double, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;

  // ¾âàÐÖÕÝØÕ Þâ äØÚâØÒÝÞÙ ÓàÐÝØæë Ò âÕÛÕ á Ø×ÒÕáâÝëÜ ÝÞÜÕàÞÜ

  ray* raw_ray_fict_bord(double, double, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  ray* raw_ray_fict_bord(ray_parameters, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  point sp()const{return pShotPoint;};  
  int n1()const{return N1;};
  int n2()const{return N2;};
  void set_uw_kind(int ik){iUWKind=ik;};
  void set_dw_kind(int ik){iDWKind=ik;};
  int uw_kind(){return iUWKind;};
  int dw_kind(){return iDWKind;};
};

// »îÜßÕÝë...
Logic equal_paths_to_border (ray* r1, ray*r2, int iB);
Logic equal_paths_to_node (ray* r1, ray*r2, int iN);


#endif /* rays.H */
/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 11-Jan-1999  */
/*          ...          */

  rays_to_receiver(db_results_manager& dbm,const char* sP,int id,int n1,int n2,
           int iHRR=PM_ADD,void (*idle_func)()=NULL);

  const rayPathAr* get_ray_paths(){return ray_path_ar;};

  ~rays_to_receiver();

  point3d receiver()const{return pReceiverPoint;};
  int n1() const {return N1;};
  int n2() const {return N2;};
  virtual int store() const;
  virtual int load();
};


class rays_to_receivers:public ObjectsAr<rays_to_receiver*>{
public:
  rays_to_receivers(){weFormatMsg(EL_DEBUG,"Creating HitRays...");};
  /* if ~rays_to_receivers is absent then g++ 2.7.2.1 generates:
     rays.H:338: Internal compiler error.
     rays.H:338: Please submit a full bug report to `bug-g++@prep.ai.mit.edu'.
   */
  ~rays_to_receivers(){
    for(register int i=0;i<count();i++){
      delete fetch(i);
    };
    clean();
  };
  rays_to_receiver *get(int N1,int N2)const;
};

class rays_shoot{
protected:
  point3d pShotPoint;        // ºÞÞàÔØÝÐâë ØáâÞçÝØÚÐ
  pl_raytr_model *pmFather;  // ¿ÐßÐ-ÜÞÔÕÛì    
  int iShotType;             // ÂØß âÞçÚØ - ØáâÞçÝØÚ/ßàØÕÜÝØÚ
  
  int N1,N2;                 // ²ÞÕÝÝëÕ ÝÞÜÕàÐ ÔÛï ØÔÕÝâØäØÚÐæØØ ²ÞÛÞÔÕÙ  
  // ¿ÞÛãçØâì Ûãç (Ø ÒáÕ ÕÓÞ ßÕàÕáÕçÕÝØï)
  ray* shoot_out(const ray_parameters &rp,int iKind=PM_FULL)const; 
  ray* a_shot (const double& dFi,const double& dAz,ray_specifics& rsSpec,
           ray*r=NULL)const;
  int iUWKind;
  int iDWKind;
public:
  rays_shoot(const point3d& p,int,int,int,pl_raytr_model *);
  // ´Ûï ßÞÛãçÕÝØï áÕÙáÜÞÓàÐÜÜ
  ray *raw_ray(const ray_parameters& rp) const;
  ray *raw_ray(double df,double da) const;

  // ½ÐÑÞà àÐ×ÛØçÝëå äãÝÚæØÙ ßÞßÐÔÐÝØï Ò âÞçÚã
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point(ray_parameters, const point3d&,
             int iBNumber,Logic&,double& lE,IntAr *iat,ray*r=NULL) const;
  ray* raw_ray_point(double, double,
             int iBNumber, const point3d&, Logic&,double& lE,IntAr *iat,
             ray*r=NULL) const;


  // Á ÞâàÐÖÕÝØÕÜ
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point_reflected(ray_parameters, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;
  ray* raw_ray_point_reflected(double, double, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;

  // ¾âàÐÖÕÝØÕ Þâ äØÚâØÒÝÞÙ ÓàÐÝØæë Ò âÕÛÕ á Ø×ÒÕáâÝëÜ ÝÞÜÕàÞÜ

  ray* raw_ray_fict_bord(double, double, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  ray* raw_ray_fict_bord(ray_parameters, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  point sp()const{return pShotPoint;};  
  int n1()const{return N1;};
  int n2()const{return N2;};
  void set_uw_kind(int ik){iUWKind=ik;};
  void set_dw_kind(int ik){iDWKind=ik;};
  int uw_kind(){return iUWKind;};
  int dw_kind(){return iDWKind;};
};

// »îÜßÕÝë...
Logic equal_paths_to_border (ray* r1, ray*r2, int iB);
Logic equal_paths_to_node (ray* r1, ray*r2, int iN);


#endif /* rays.H */
/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 11-Jan-1999  */
/*          ...          */

  
  // ¿ÞÛãçØâì Ûãç (Ø ÒáÕ ÕÓÞ ßÕàÕáÕçÕÝØï)
  ray* shoot_out(const ray_parameters &rp,int iKind=PM_FULL)const; 
  ray* a_shot (const double& dFi,const double& dAz,ray_specifics& rsSpec,
           ray*r=NULL)const;
  int iUWKind;
  int iDWKind;
public:
  rays_shoot(const point3d& p,int,int,int,pl_raytr_model *);
  // ´Ûï ßÞÛãçÕÝØï áÕÙáÜÞÓàÐÜÜ  ray *raw_ray(const ray_parameters& rp) const;
  ray *raw_ray(double df,double da) const;

  // ½ÐÑÞà àÐ×ÛØçÝëå äãÝÚæØÙ ßÞßÐÔÐÝØï Ò âÞçÚã
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point(ray_parameters, const point3d&,
             int iBNumber,Logic&,double& lE,IntAr *iat,ray*r=NULL) const;
  ray* raw_ray_point(double, double,
             int iBNumber, const point3d&, Logic&,double& lE,IntAr *iat,
             ray*r=NULL) const;


  // Á ÞâàÐÖÕÝØÕÜ
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point_reflected(ray_parameters, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;
  ray* raw_ray_point_reflected(double, double, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;

  // ¾âàÐÖÕÝØÕ Þâ äØÚâØÒÝÞÙ ÓàÐÝØæë Ò âÕÛÕ á Ø×ÒÕáâÝëÜ ÝÞÜÕàÞÜ

  ray* raw_ray_fict_bord(double, double, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  ray* raw_ray_fict_bord(ray_parameters, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  point sp()const{return pShotPoint;};  
  int n1()const{return N1;};
  int n2()const{return N2;};
  void set_uw_kind(int ik){iUWKind=ik;};
  void set_dw_kind(int ik){iDWKind=ik;};
  int uw_kind(){return iUWKind;};
  int dw_kind(){return iDWKind;};
};

// »îÜßÕÝë...
Logic equal_paths_to_border (ray* r1, ray*r2, int iB);
Logic equal_paths_to_node (ray* r1, ray*r2, int iN);


#endif /* rays.H */
/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 11-Jan-1999  */
/*          ...          */

  ray *raw_ray(const ray_parameters& rp) const;
  ray *raw_ray(double df,double da) const;

  // ½ÐÑÞà àÐ×ÛØçÝëå äãÝÚæØÙ ßÞßÐÔÐÝØï Ò âÞçÚã
  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point(ray_parameters, const point3d&,
             int iBNumber,Logic&,double& lE,IntAr *iat,ray*r=NULL) const;
  ray* raw_ray_point(double, double,
             int iBNumber, const point3d&, Logic&,double& lE,IntAr *iat,
             ray*r=NULL) const;


  // Á ÞâàÐÖÕÝØÕÜ  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point_reflected(ray_parameters, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;
  ray* raw_ray_point_reflected(double, double, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;

  // ¾âàÐÖÕÝØÕ Þâ äØÚâØÒÝÞÙ ÓàÐÝØæë Ò âÕÛÕ á Ø×ÒÕáâÝëÜ ÝÞÜÕàÞÜ

  ray* raw_ray_fict_bord(double, double, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  ray* raw_ray_fict_bord(ray_parameters, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  point sp()const{return pShotPoint;};  
  int n1()const{return N1;};
  int n2()const{return N2;};
  void set_uw_kind(int ik){iUWKind=ik;};
  void set_dw_kind(int ik){iDWKind=ik;};
  int uw_kind(){return iUWKind;};
  int dw_kind(){return iDWKind;};
};

// »îÜßÕÝë...
Logic equal_paths_to_border (ray* r1, ray*r2, int iB);
Logic equal_paths_to_node (ray* r1, ray*r2, int iN);


#endif /* rays.H */
/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 11-Jan-1999  */
/*          ...          */

  // ¿ëâÐÕÜáï ßÞßÐáâì á âÞçÝÞáâìî lE
  // iBNumber - ÝÞÜÕà âÕÛÐ, Ò ÚÞâÞàÞÜ ÛÕÖØâ âÞçÚÐ - ßÕàÕÔÐÕâáï ÔÛï ÞßâØÜØ×ÐæØØ
  ray* raw_ray_point_reflected(ray_parameters, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;
  ray* raw_ray_point_reflected(double, double, const point3d&,int iBNumber,
                   Logic&,int,double& lE,IntAr *iat,
                   Logic lRefRat=FALSE,ray*r=NULL) const;

  // ¾âàÐÖÕÝØÕ Þâ äØÚâØÒÝÞÙ ÓàÐÝØæë Ò âÕÛÕ á Ø×ÒÕáâÝëÜ ÝÞÜÕàÞÜ

  ray* raw_ray_fict_bord(double, double, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  ray* raw_ray_fict_bord(ray_parameters, const point3d&,int iBInNumber,
             int iBReflNumber,const pl_fict_border& fbReflect,
             Logic&,double& lE,point3d&,ray*r=NULL) const;
  point sp()const{return pShotPoint;};  
  int n1()const{return N1;};
  int n2()const{return N2;};
  void set_uw_kind(int ik){iUWKind=ik;};
  void set_dw_kind(int ik){iDWKind=ik;};
  int uw_kind(){return iUWKind;};
  int dw_kind(){return iDWKind;};
};

// »îÜßÕÝë...
Logic equal_paths_to_border (ray* r1, ray*r2, int iB);
Logic equal_paths_to_node (ray* r1, ray*r2, int iN);


#endif /* rays.H */
/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 11-Jan-1999  */
/*          ...          */