p_raytr_elements.H

/* p_raytr_elements.H */
/* $Id: p_raytr_elements.H,v 2.5 2001/10/09 14:54:42 hoh Exp $ */
#ifndef __p_raytr_elements_H
#define __p_raytr_elements_H
/*
   ÍÛÕÜÕÝâë âàÕåÜÕàÝÞÙ ÜÞÔÕÛØ á ßÛÞáÚØÜØ ÓàÐÝØæÐÜØ, ÞÑÕáßÕçØÒÐîéØÕ âàÐááØàÞÒÚã
   ÛãçÕÙ.

   crossings etc -- ßÕàÕáÕçÕÝØï ÛãçÕÙ á ÓàÐÝØæÐÜØ.

   regular_set -- àÕÓãÛïàÝëÕ áÕâÚØ âÞçÚÕ ÝÐ ÓàÐÝØæÐå, ÞÑÕáßÕçØÒÐîéØÕ ßÞÛÝÞÕ
   ßÞÚàëâØÕ ÓàÐÝØæ ÝÐÑÞàÞÜ ÛãçÕÙ.

   pl_intern_raytr_âØß -- ÕéÕ ÞÔÝÞ ßÞàÞÖÔÕÝØÕ Þâ any_âØß, ÔÛï ßÞÔÔÕàÖÚØ ÝÐáâàÕÛÐ   ÛãçÕÙ áÞÑáâÒÕÝÝÞ ÜÞÔÕÛìî

   meta_iis -- áâàãÚâãàÐ, áÞÔÕàÖÐéÐï ÔÐÝÝëÕ, ßÞ ÚÞâÞàëÜ ÑãÔÕâ äÞàÜØàÞÒÐâìáï
   ØÔÕÐÛìÝÐï ØÜßãÛìáÝÐï áÕÙáÜÞÓàÐÜÜÐ
 */

#include <mix/wenv.h>
#include <mix/ObjectsAr.H>

#include <mod3d/any_elements.H>
#include <mod3d/p_elements.H>

#include "rays.H"
#include "p_raytr_hits.H"

#define PM_TRANSMISSION 0
#define PM_REFLECTION 1

// ºÛÐááë, ãßÞÜØÝÐîéØÕáï Ò äÐÙÛÕ
class pl_intern_raytr_border;
class rays_to_receiver;

// ¾ßØáÐÝØÕ ßàØÕÜÝØÚÐ
class receiver_info{
protected:  
  int iN1;
  int iN2;
  point pRP;
public: 
  receiver_info(int n1,int n2,point p){iN1=n1;iN2=n2;pRP=p;};
  point rp()const{return pRP;};
  int n1(){return iN1;};
  int n2(){return iN2;};  
  Logic same_numbers(int n1,int n2)const{
    if((n1==iN1)&&(n2==iN2))return TRUE;
    return FALSE;
  };
};
#define AnyType receiver_info
#define AnyAr   receivers
#include <mix/AnyAr.H>

// ³ÕÞäØ×ØçÕáÚÐï ØÝäÞàÜÐæØï Þ ßÕàÕáÕçÕÝØØ ÛãçÐ á ÓàÐÝØæÕÙ (Þâ ¸ÓÞàï)
class cross_info{
public:
  // ²àÕÜï ßÕàÕáÕçÕÝØï
  double time;
  cross_info(){};
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï-ßàÕÛÞÜÛÕÝØï ÝÕ ×ÐÝÞáïâáï áàÐ×ã, Øå ÜÞÖÕâ Ø ÝÕ Ñëâì
  cross_info(double t,
         double dDx1,double dDy1,double dDz1,
         double dDx2,double dDy2,double dDz2,
         double dF  ,double dA,
         //  double dkdfi[3],double dkdt[3], double ksi[3], double vp[3],
         const double *dkdfi,const double *dkdt,const double *ksi,
         const double *vp,double *daDpdteta0,double *daDpdfi0,int ib,
         int w_t);
  void add_rc_coefs(double a[2],double b[2]);
  void add_rc_vectors(double[3],double[3]);
  void add_raypol(double da[3]);
  // ¿àÞØ×ÒÞÔÝëÕ ßÞ ãÓÛÐÜ
  double dDxDf,dDyDf,dDzDf;
  double dDxDaz,dDyDaz,dDzDaz;

  // ÁÞÑáâÒÕÝÝÞ ãÓÛë
  double dFi,dAz;

  // ÁÛãÖÕÑÝÐï ÔÛï ¸ÓÞàï
  double dksidfi[3];
  double dksidteta[3];
  double ksi[3];
  int iBody; // ½ÞÜÕà âÕÛÐ, ßÞ ÚÞâÞàÞÜã Ò íâÞ ßÕàÕáÕçÕÝØÕ Ûãç ßàØèÕÛ
  int wave_type; // âØß ÒÞÛÝë ÝÐ ÔÐÝÝÞÜ áÕÓÜÕÝâÕ ÛãçÐ

  // ¸ÝäÞàÜÐæØï, ÝÕÞÑåÞÔØÜÐï áÞÑáâÒÕÝÝÞ ÔÛï áÕÙáÜÞÓàÐÜÜ

  double vpol[3];
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï Ø ßàÕÛÞÜÛÕÝØï, ÚÞÜßÛÕÚáÝëÕ
  double dCRef[2],dCTrans[2];
  // »ãçÕÒëÕ ßÐàÐÜÕâàë ßàÞåÞÖÔÕÝØï Ø ßàÕÛÞÜÛÕÝØï ÝÐ ÓàÐÝØæÕ, ßÞáÛÕ refkin
  double vref[3],vtrans[3];
  // ²ÞÕÝÝÞÕ çØáÕÛÚÞ ßÞ ßÞÒÞÔã ÞâÝÞèÕÝØï ÚÞíääØæØÕÝâÞÒ ÞâàÐÖÕÝØï ÔÛï ¸ÓÞàï
  double dRefRat;

  double daDpdteta0[3];
  double daDpdfi0[3];

  // ²ÕÚâÞà ÔÛï ßÞÛïàØ×ÐæØØ ßÞßÕàÕçÝÞÙ ÒÞÛÝë ßÞâàÕÑÝëÙ...
  double daRayPol[3];
  int operator =(const cross_info&);
};


// ¾ßØáÐÝØÕ ÞÔÝÞÓÞ ßÕàÕáÕçÕÝØï ÛãçÐ á ÓàÐÝØæÕÙ

// ¾àÓÐÝØ×ÐæØï ÝÐÑÞàÐ ßÕàÕáÕçÕÝØÙ ÝÐ ßÞÒÕàåÝÞáâØ

// ÁÛãÖÕÑÝëÙ ÚÛÐáá, ÝãÖÕÝ ÝÕ ÒáÕÜ
class crossing:public point3d{
protected:
  int iBorder; // ½ÞÜÕà ÓàÐÝØæë, ÝÐ ÚÞâÞàÞÙ ßàÞØ×ÞèÛÞ ßÕàÕáÕçÕÝØÕ
  cross_info cr;
  ray_parameters rp;
  int iKind; // 0 - ßàÕÛÞÜÛÕÝØÕ, 1 - ÞâàÐÖÕÝØÕ
public:
  crossing(const ray_parameters &rp, cross_info cr_,int ib,
       double x,double y, double z,int ik);
  crossing();
  cross_info crossing_info()const {return cr;};
  ray_parameters ray_pars(){return rp;};
  int border()const{return iBorder;};
  void add_rc_coefs( double c1[2], double c2[2]){cr.add_rc_coefs(c1,c2);};
  void add_rc_vectors(double v1[3],double v2[3]){cr.add_rc_vectors(v1,v2);};
  void add_raypol(double da[3]){cr.add_raypol(da);};
  int reflected() const {return iKind;};
  int operator =(const crossing&);
  void set_refl_type(int ik);
  void set_refrat(double rr);
  int kind() {return iKind;};
  int get_wave_type() {return cr.wave_type;};
  ray_parameters rp_as_sp()const;
};

typedef crossing* p_crossing;

#define AnyType crossing
#define AnyAr   crossings_ptrs
#include <mix/AnyAr.H>

struct meta_iis_element{
  double dTime;
  double dValueX[2];
  double dValueY[2];
  double dValueZ[2];
};

#define AnyType meta_iis_element
#define AnyAr   meta_iis
#include <mix/AnyAr.H>



/*
  ±Õ×ÚÛÐááÞÒëÕ äãÝÚæØØ 
 */
void general_get_ray_parameters (const point3d &pWhere,
                 const crossing &cClosest,
                 double &dFi,double &dAz,
                 double &dDFi,double& dDAz,double &dRo,
                 double dfl=0.8,double dal=0.8);

double compare_reflected(const crossing&,const crossing&);
double compare_transmitted(const crossing&,const crossing&);
void get_ray_angles (double& dFi,double& dAz,const point&pStart,
             const point& pCurr,const geophys_params&gp);
#endif /* p_raytr_elements.H */

/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 31-Jan-2000  */ 
/*          ...          */

   ÛãçÕÙ áÞÑáâÒÕÝÝÞ ÜÞÔÕÛìî
   meta_iis -- áâàãÚâãàÐ, áÞÔÕàÖÐéÐï ÔÐÝÝëÕ, ßÞ ÚÞâÞàëÜ ÑãÔÕâ äÞàÜØàÞÒÐâìáï
   ØÔÕÐÛìÝÐï ØÜßãÛìáÝÐï áÕÙáÜÞÓàÐÜÜÐ
 */

#include <mix/wenv.h>
#include <mix/ObjectsAr.H>

#include <mod3d/any_elements.H>
#include <mod3d/p_elements.H>

#include "rays.H"
#include "p_raytr_hits.H"

#define PM_TRANSMISSION 0
#define PM_REFLECTION 1

// ºÛÐááë, ãßÞÜØÝÐîéØÕáï Ò äÐÙÛÕ
class pl_intern_raytr_border;
class rays_to_receiver;

// ¾ßØáÐÝØÕ ßàØÕÜÝØÚÐ
class receiver_info{
protected:  
  int iN1;
  int iN2;
  point pRP;
public: 
  receiver_info(int n1,int n2,point p){iN1=n1;iN2=n2;pRP=p;};
  point rp()const{return pRP;};
  int n1(){return iN1;};
  int n2(){return iN2;};  
  Logic same_numbers(int n1,int n2)const{
    if((n1==iN1)&&(n2==iN2))return TRUE;
    return FALSE;
  };
};
#define AnyType receiver_info
#define AnyAr   receivers
#include <mix/AnyAr.H>

// ³ÕÞäØ×ØçÕáÚÐï ØÝäÞàÜÐæØï Þ ßÕàÕáÕçÕÝØØ ÛãçÐ á ÓàÐÝØæÕÙ (Þâ ¸ÓÞàï)
class cross_info{
public:
  // ²àÕÜï ßÕàÕáÕçÕÝØï
  double time;
  cross_info(){};
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï-ßàÕÛÞÜÛÕÝØï ÝÕ ×ÐÝÞáïâáï áàÐ×ã, Øå ÜÞÖÕâ Ø ÝÕ Ñëâì
  cross_info(double t,
         double dDx1,double dDy1,double dDz1,
         double dDx2,double dDy2,double dDz2,
         double dF  ,double dA,
         //  double dkdfi[3],double dkdt[3], double ksi[3], double vp[3],
         const double *dkdfi,const double *dkdt,const double *ksi,
         const double *vp,double *daDpdteta0,double *daDpdfi0,int ib,
         int w_t);
  void add_rc_coefs(double a[2],double b[2]);
  void add_rc_vectors(double[3],double[3]);
  void add_raypol(double da[3]);
  // ¿àÞØ×ÒÞÔÝëÕ ßÞ ãÓÛÐÜ
  double dDxDf,dDyDf,dDzDf;
  double dDxDaz,dDyDaz,dDzDaz;

  // ÁÞÑáâÒÕÝÝÞ ãÓÛë
  double dFi,dAz;

  // ÁÛãÖÕÑÝÐï ÔÛï ¸ÓÞàï
  double dksidfi[3];
  double dksidteta[3];
  double ksi[3];
  int iBody; // ½ÞÜÕà âÕÛÐ, ßÞ ÚÞâÞàÞÜã Ò íâÞ ßÕàÕáÕçÕÝØÕ Ûãç ßàØèÕÛ
  int wave_type; // âØß ÒÞÛÝë ÝÐ ÔÐÝÝÞÜ áÕÓÜÕÝâÕ ÛãçÐ

  // ¸ÝäÞàÜÐæØï, ÝÕÞÑåÞÔØÜÐï áÞÑáâÒÕÝÝÞ ÔÛï áÕÙáÜÞÓàÐÜÜ

  double vpol[3];
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï Ø ßàÕÛÞÜÛÕÝØï, ÚÞÜßÛÕÚáÝëÕ
  double dCRef[2],dCTrans[2];
  // »ãçÕÒëÕ ßÐàÐÜÕâàë ßàÞåÞÖÔÕÝØï Ø ßàÕÛÞÜÛÕÝØï ÝÐ ÓàÐÝØæÕ, ßÞáÛÕ refkin
  double vref[3],vtrans[3];
  // ²ÞÕÝÝÞÕ çØáÕÛÚÞ ßÞ ßÞÒÞÔã ÞâÝÞèÕÝØï ÚÞíääØæØÕÝâÞÒ ÞâàÐÖÕÝØï ÔÛï ¸ÓÞàï
  double dRefRat;

  double daDpdteta0[3];
  double daDpdfi0[3];

  // ²ÕÚâÞà ÔÛï ßÞÛïàØ×ÐæØØ ßÞßÕàÕçÝÞÙ ÒÞÛÝë ßÞâàÕÑÝëÙ...
  double daRayPol[3];
  int operator =(const cross_info&);
};


// ¾ßØáÐÝØÕ ÞÔÝÞÓÞ ßÕàÕáÕçÕÝØï ÛãçÐ á ÓàÐÝØæÕÙ

// ¾àÓÐÝØ×ÐæØï ÝÐÑÞàÐ ßÕàÕáÕçÕÝØÙ ÝÐ ßÞÒÕàåÝÞáâØ

// ÁÛãÖÕÑÝëÙ ÚÛÐáá, ÝãÖÕÝ ÝÕ ÒáÕÜ
class crossing:public point3d{
protected:
  int iBorder; // ½ÞÜÕà ÓàÐÝØæë, ÝÐ ÚÞâÞàÞÙ ßàÞØ×ÞèÛÞ ßÕàÕáÕçÕÝØÕ
  cross_info cr;
  ray_parameters rp;
  int iKind; // 0 - ßàÕÛÞÜÛÕÝØÕ, 1 - ÞâàÐÖÕÝØÕ
public:
  crossing(const ray_parameters &rp, cross_info cr_,int ib,
       double x,double y, double z,int ik);
  crossing();
  cross_info crossing_info()const {return cr;};
  ray_parameters ray_pars(){return rp;};
  int border()const{return iBorder;};
  void add_rc_coefs( double c1[2], double c2[2]){cr.add_rc_coefs(c1,c2);};
  void add_rc_vectors(double v1[3],double v2[3]){cr.add_rc_vectors(v1,v2);};
  void add_raypol(double da[3]){cr.add_raypol(da);};
  int reflected() const {return iKind;};
  int operator =(const crossing&);
  void set_refl_type(int ik);
  void set_refrat(double rr);
  int kind() {return iKind;};
  int get_wave_type() {return cr.wave_type;};
  ray_parameters rp_as_sp()const;
};

typedef crossing* p_crossing;

#define AnyType crossing
#define AnyAr   crossings_ptrs
#include <mix/AnyAr.H>

struct meta_iis_element{
  double dTime;
  double dValueX[2];
  double dValueY[2];
  double dValueZ[2];
};

#define AnyType meta_iis_element
#define AnyAr   meta_iis
#include <mix/AnyAr.H>



/*
  ±Õ×ÚÛÐááÞÒëÕ äãÝÚæØØ 
 */
void general_get_ray_parameters (const point3d &pWhere,
                 const crossing &cClosest,
                 double &dFi,double &dAz,
                 double &dDFi,double& dDAz,double &dRo,
                 double dfl=0.8,double dal=0.8);

double compare_reflected(const crossing&,const crossing&);
double compare_transmitted(const crossing&,const crossing&);
void get_ray_angles (double& dFi,double& dAz,const point&pStart,
             const point& pCurr,const geophys_params&gp);
#endif /* p_raytr_elements.H */

/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 31-Jan-2000  */ 
/*          ...          */


   meta_iis -- áâàãÚâãàÐ, áÞÔÕàÖÐéÐï ÔÐÝÝëÕ, ßÞ ÚÞâÞàëÜ ÑãÔÕâ äÞàÜØàÞÒÐâìáï
   ØÔÕÐÛìÝÐï ØÜßãÛìáÝÐï áÕÙáÜÞÓàÐÜÜÐ
 */

#include <mix/wenv.h>
#include <mix/ObjectsAr.H>

#include <mod3d/any_elements.H>
#include <mod3d/p_elements.H>

#include "rays.H"
#include "p_raytr_hits.H"

#define PM_TRANSMISSION 0
#define PM_REFLECTION 1

// ºÛÐááë, ãßÞÜØÝÐîéØÕáï Ò äÐÙÛÕ
class pl_intern_raytr_border;
class rays_to_receiver;

// ¾ßØáÐÝØÕ ßàØÕÜÝØÚÐclass receiver_info{
protected:  
  int iN1;
  int iN2;
  point pRP;
public: 
  receiver_info(int n1,int n2,point p){iN1=n1;iN2=n2;pRP=p;};
  point rp()const{return pRP;};
  int n1(){return iN1;};
  int n2(){return iN2;};  
  Logic same_numbers(int n1,int n2)const{
    if((n1==iN1)&&(n2==iN2))return TRUE;
    return FALSE;
  };
};
#define AnyType receiver_info
#define AnyAr   receivers
#include <mix/AnyAr.H>

// ³ÕÞäØ×ØçÕáÚÐï ØÝäÞàÜÐæØï Þ ßÕàÕáÕçÕÝØØ ÛãçÐ á ÓàÐÝØæÕÙ (Þâ ¸ÓÞàï)
class cross_info{
public:
  // ²àÕÜï ßÕàÕáÕçÕÝØï
  double time;
  cross_info(){};
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï-ßàÕÛÞÜÛÕÝØï ÝÕ ×ÐÝÞáïâáï áàÐ×ã, Øå ÜÞÖÕâ Ø ÝÕ Ñëâì
  cross_info(double t,
         double dDx1,double dDy1,double dDz1,
         double dDx2,double dDy2,double dDz2,
         double dF  ,double dA,
         //  double dkdfi[3],double dkdt[3], double ksi[3], double vp[3],
         const double *dkdfi,const double *dkdt,const double *ksi,
         const double *vp,double *daDpdteta0,double *daDpdfi0,int ib,
         int w_t);
  void add_rc_coefs(double a[2],double b[2]);
  void add_rc_vectors(double[3],double[3]);
  void add_raypol(double da[3]);
  // ¿àÞØ×ÒÞÔÝëÕ ßÞ ãÓÛÐÜ
  double dDxDf,dDyDf,dDzDf;
  double dDxDaz,dDyDaz,dDzDaz;

  // ÁÞÑáâÒÕÝÝÞ ãÓÛë
  double dFi,dAz;

  // ÁÛãÖÕÑÝÐï ÔÛï ¸ÓÞàï
  double dksidfi[3];
  double dksidteta[3];
  double ksi[3];
  int iBody; // ½ÞÜÕà âÕÛÐ, ßÞ ÚÞâÞàÞÜã Ò íâÞ ßÕàÕáÕçÕÝØÕ Ûãç ßàØèÕÛ
  int wave_type; // âØß ÒÞÛÝë ÝÐ ÔÐÝÝÞÜ áÕÓÜÕÝâÕ ÛãçÐ

  // ¸ÝäÞàÜÐæØï, ÝÕÞÑåÞÔØÜÐï áÞÑáâÒÕÝÝÞ ÔÛï áÕÙáÜÞÓàÐÜÜ

  double vpol[3];
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï Ø ßàÕÛÞÜÛÕÝØï, ÚÞÜßÛÕÚáÝëÕ
  double dCRef[2],dCTrans[2];
  // »ãçÕÒëÕ ßÐàÐÜÕâàë ßàÞåÞÖÔÕÝØï Ø ßàÕÛÞÜÛÕÝØï ÝÐ ÓàÐÝØæÕ, ßÞáÛÕ refkin
  double vref[3],vtrans[3];
  // ²ÞÕÝÝÞÕ çØáÕÛÚÞ ßÞ ßÞÒÞÔã ÞâÝÞèÕÝØï ÚÞíääØæØÕÝâÞÒ ÞâàÐÖÕÝØï ÔÛï ¸ÓÞàï
  double dRefRat;

  double daDpdteta0[3];
  double daDpdfi0[3];

  // ²ÕÚâÞà ÔÛï ßÞÛïàØ×ÐæØØ ßÞßÕàÕçÝÞÙ ÒÞÛÝë ßÞâàÕÑÝëÙ...
  double daRayPol[3];
  int operator =(const cross_info&);
};


// ¾ßØáÐÝØÕ ÞÔÝÞÓÞ ßÕàÕáÕçÕÝØï ÛãçÐ á ÓàÐÝØæÕÙ

// ¾àÓÐÝØ×ÐæØï ÝÐÑÞàÐ ßÕàÕáÕçÕÝØÙ ÝÐ ßÞÒÕàåÝÞáâØ

// ÁÛãÖÕÑÝëÙ ÚÛÐáá, ÝãÖÕÝ ÝÕ ÒáÕÜ
class crossing:public point3d{
protected:
  int iBorder; // ½ÞÜÕà ÓàÐÝØæë, ÝÐ ÚÞâÞàÞÙ ßàÞØ×ÞèÛÞ ßÕàÕáÕçÕÝØÕ
  cross_info cr;
  ray_parameters rp;
  int iKind; // 0 - ßàÕÛÞÜÛÕÝØÕ, 1 - ÞâàÐÖÕÝØÕ
public:
  crossing(const ray_parameters &rp, cross_info cr_,int ib,
       double x,double y, double z,int ik);
  crossing();
  cross_info crossing_info()const {return cr;};
  ray_parameters ray_pars(){return rp;};
  int border()const{return iBorder;};
  void add_rc_coefs( double c1[2], double c2[2]){cr.add_rc_coefs(c1,c2);};
  void add_rc_vectors(double v1[3],double v2[3]){cr.add_rc_vectors(v1,v2);};
  void add_raypol(double da[3]){cr.add_raypol(da);};
  int reflected() const {return iKind;};
  int operator =(const crossing&);
  void set_refl_type(int ik);
  void set_refrat(double rr);
  int kind() {return iKind;};
  int get_wave_type() {return cr.wave_type;};
  ray_parameters rp_as_sp()const;
};

typedef crossing* p_crossing;

#define AnyType crossing
#define AnyAr   crossings_ptrs
#include <mix/AnyAr.H>

struct meta_iis_element{
  double dTime;
  double dValueX[2];
  double dValueY[2];
  double dValueZ[2];
};

#define AnyType meta_iis_element
#define AnyAr   meta_iis
#include <mix/AnyAr.H>



/*
  ±Õ×ÚÛÐááÞÒëÕ äãÝÚæØØ 
 */
void general_get_ray_parameters (const point3d &pWhere,
                 const crossing &cClosest,
                 double &dFi,double &dAz,
                 double &dDFi,double& dDAz,double &dRo,
                 double dfl=0.8,double dal=0.8);

double compare_reflected(const crossing&,const crossing&);
double compare_transmitted(const crossing&,const crossing&);
void get_ray_angles (double& dFi,double& dAz,const point&pStart,
             const point& pCurr,const geophys_params&gp);
#endif /* p_raytr_elements.H */

/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 31-Jan-2000  */ 
/*          ...          */

class receiver_info{
protected:  
  int iN1;
  int iN2;
  point pRP;
public: 
  receiver_info(int n1,int n2,point p){iN1=n1;iN2=n2;pRP=p;};
  point rp()const{return pRP;};
  int n1(){return iN1;};
  int n2(){return iN2;};  
  Logic same_numbers(int n1,int n2)const{
    if((n1==iN1)&&(n2==iN2))return TRUE;
    return FALSE;
  };
};
#define AnyType receiver_info
#define AnyAr   receivers
#include <mix/AnyAr.H>

// ³ÕÞäØ×ØçÕáÚÐï ØÝäÞàÜÐæØï Þ ßÕàÕáÕçÕÝØØ ÛãçÐ á ÓàÐÝØæÕÙ (Þâ ¸ÓÞàï)class cross_info{
public:
  // ²àÕÜï ßÕàÕáÕçÕÝØï
  double time;
  cross_info(){};
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï-ßàÕÛÞÜÛÕÝØï ÝÕ ×ÐÝÞáïâáï áàÐ×ã, Øå ÜÞÖÕâ Ø ÝÕ Ñëâì
  cross_info(double t,
         double dDx1,double dDy1,double dDz1,
         double dDx2,double dDy2,double dDz2,
         double dF  ,double dA,
         //  double dkdfi[3],double dkdt[3], double ksi[3], double vp[3],
         const double *dkdfi,const double *dkdt,const double *ksi,
         const double *vp,double *daDpdteta0,double *daDpdfi0,int ib,
         int w_t);
  void add_rc_coefs(double a[2],double b[2]);
  void add_rc_vectors(double[3],double[3]);
  void add_raypol(double da[3]);
  // ¿àÞØ×ÒÞÔÝëÕ ßÞ ãÓÛÐÜ
  double dDxDf,dDyDf,dDzDf;
  double dDxDaz,dDyDaz,dDzDaz;

  // ÁÞÑáâÒÕÝÝÞ ãÓÛë
  double dFi,dAz;

  // ÁÛãÖÕÑÝÐï ÔÛï ¸ÓÞàï
  double dksidfi[3];
  double dksidteta[3];
  double ksi[3];
  int iBody; // ½ÞÜÕà âÕÛÐ, ßÞ ÚÞâÞàÞÜã Ò íâÞ ßÕàÕáÕçÕÝØÕ Ûãç ßàØèÕÛ
  int wave_type; // âØß ÒÞÛÝë ÝÐ ÔÐÝÝÞÜ áÕÓÜÕÝâÕ ÛãçÐ

  // ¸ÝäÞàÜÐæØï, ÝÕÞÑåÞÔØÜÐï áÞÑáâÒÕÝÝÞ ÔÛï áÕÙáÜÞÓàÐÜÜ

  double vpol[3];
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï Ø ßàÕÛÞÜÛÕÝØï, ÚÞÜßÛÕÚáÝëÕ
  double dCRef[2],dCTrans[2];
  // »ãçÕÒëÕ ßÐàÐÜÕâàë ßàÞåÞÖÔÕÝØï Ø ßàÕÛÞÜÛÕÝØï ÝÐ ÓàÐÝØæÕ, ßÞáÛÕ refkin
  double vref[3],vtrans[3];
  // ²ÞÕÝÝÞÕ çØáÕÛÚÞ ßÞ ßÞÒÞÔã ÞâÝÞèÕÝØï ÚÞíääØæØÕÝâÞÒ ÞâàÐÖÕÝØï ÔÛï ¸ÓÞàï
  double dRefRat;

  double daDpdteta0[3];
  double daDpdfi0[3];

  // ²ÕÚâÞà ÔÛï ßÞÛïàØ×ÐæØØ ßÞßÕàÕçÝÞÙ ÒÞÛÝë ßÞâàÕÑÝëÙ...
  double daRayPol[3];
  int operator =(const cross_info&);
};


// ¾ßØáÐÝØÕ ÞÔÝÞÓÞ ßÕàÕáÕçÕÝØï ÛãçÐ á ÓàÐÝØæÕÙ

// ¾àÓÐÝØ×ÐæØï ÝÐÑÞàÐ ßÕàÕáÕçÕÝØÙ ÝÐ ßÞÒÕàåÝÞáâØ

// ÁÛãÖÕÑÝëÙ ÚÛÐáá, ÝãÖÕÝ ÝÕ ÒáÕÜ
class crossing:public point3d{
protected:
  int iBorder; // ½ÞÜÕà ÓàÐÝØæë, ÝÐ ÚÞâÞàÞÙ ßàÞØ×ÞèÛÞ ßÕàÕáÕçÕÝØÕ
  cross_info cr;
  ray_parameters rp;
  int iKind; // 0 - ßàÕÛÞÜÛÕÝØÕ, 1 - ÞâàÐÖÕÝØÕ
public:
  crossing(const ray_parameters &rp, cross_info cr_,int ib,
       double x,double y, double z,int ik);
  crossing();
  cross_info crossing_info()const {return cr;};
  ray_parameters ray_pars(){return rp;};
  int border()const{return iBorder;};
  void add_rc_coefs( double c1[2], double c2[2]){cr.add_rc_coefs(c1,c2);};
  void add_rc_vectors(double v1[3],double v2[3]){cr.add_rc_vectors(v1,v2);};
  void add_raypol(double da[3]){cr.add_raypol(da);};
  int reflected() const {return iKind;};
  int operator =(const crossing&);
  void set_refl_type(int ik);
  void set_refrat(double rr);
  int kind() {return iKind;};
  int get_wave_type() {return cr.wave_type;};
  ray_parameters rp_as_sp()const;
};

typedef crossing* p_crossing;

#define AnyType crossing
#define AnyAr   crossings_ptrs
#include <mix/AnyAr.H>

struct meta_iis_element{
  double dTime;
  double dValueX[2];
  double dValueY[2];
  double dValueZ[2];
};

#define AnyType meta_iis_element
#define AnyAr   meta_iis
#include <mix/AnyAr.H>



/*
  ±Õ×ÚÛÐááÞÒëÕ äãÝÚæØØ 
 */
void general_get_ray_parameters (const point3d &pWhere,
                 const crossing &cClosest,
                 double &dFi,double &dAz,
                 double &dDFi,double& dDAz,double &dRo,
                 double dfl=0.8,double dal=0.8);

double compare_reflected(const crossing&,const crossing&);
double compare_transmitted(const crossing&,const crossing&);
void get_ray_angles (double& dFi,double& dAz,const point&pStart,
             const point& pCurr,const geophys_params&gp);
#endif /* p_raytr_elements.H */

/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 31-Jan-2000  */ 
/*          ...          */

class cross_info{
public:
  // ²àÕÜï ßÕàÕáÕçÕÝØï
  double time;
  cross_info(){};
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï-ßàÕÛÞÜÛÕÝØï ÝÕ ×ÐÝÞáïâáï áàÐ×ã, Øå ÜÞÖÕâ Ø ÝÕ Ñëâì  cross_info(double t,
         double dDx1,double dDy1,double dDz1,
         double dDx2,double dDy2,double dDz2,
         double dF  ,double dA,
         //  double dkdfi[3],double dkdt[3], double ksi[3], double vp[3],
         const double *dkdfi,const double *dkdt,const double *ksi,
         const double *vp,double *daDpdteta0,double *daDpdfi0,int ib,
         int w_t);
  void add_rc_coefs(double a[2],double b[2]);
  void add_rc_vectors(double[3],double[3]);
  void add_raypol(double da[3]);
  // ¿àÞØ×ÒÞÔÝëÕ ßÞ ãÓÛÐÜ
  double dDxDf,dDyDf,dDzDf;
  double dDxDaz,dDyDaz,dDzDaz;

  // ÁÞÑáâÒÕÝÝÞ ãÓÛë
  double dFi,dAz;

  // ÁÛãÖÕÑÝÐï ÔÛï ¸ÓÞàï
  double dksidfi[3];
  double dksidteta[3];
  double ksi[3];
  int iBody; // ½ÞÜÕà âÕÛÐ, ßÞ ÚÞâÞàÞÜã Ò íâÞ ßÕàÕáÕçÕÝØÕ Ûãç ßàØèÕÛ
  int wave_type; // âØß ÒÞÛÝë ÝÐ ÔÐÝÝÞÜ áÕÓÜÕÝâÕ ÛãçÐ

  // ¸ÝäÞàÜÐæØï, ÝÕÞÑåÞÔØÜÐï áÞÑáâÒÕÝÝÞ ÔÛï áÕÙáÜÞÓàÐÜÜ

  double vpol[3];
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï Ø ßàÕÛÞÜÛÕÝØï, ÚÞÜßÛÕÚáÝëÕ
  double dCRef[2],dCTrans[2];
  // »ãçÕÒëÕ ßÐàÐÜÕâàë ßàÞåÞÖÔÕÝØï Ø ßàÕÛÞÜÛÕÝØï ÝÐ ÓàÐÝØæÕ, ßÞáÛÕ refkin
  double vref[3],vtrans[3];
  // ²ÞÕÝÝÞÕ çØáÕÛÚÞ ßÞ ßÞÒÞÔã ÞâÝÞèÕÝØï ÚÞíääØæØÕÝâÞÒ ÞâàÐÖÕÝØï ÔÛï ¸ÓÞàï
  double dRefRat;

  double daDpdteta0[3];
  double daDpdfi0[3];

  // ²ÕÚâÞà ÔÛï ßÞÛïàØ×ÐæØØ ßÞßÕàÕçÝÞÙ ÒÞÛÝë ßÞâàÕÑÝëÙ...
  double daRayPol[3];
  int operator =(const cross_info&);
};


// ¾ßØáÐÝØÕ ÞÔÝÞÓÞ ßÕàÕáÕçÕÝØï ÛãçÐ á ÓàÐÝØæÕÙ

// ¾àÓÐÝØ×ÐæØï ÝÐÑÞàÐ ßÕàÕáÕçÕÝØÙ ÝÐ ßÞÒÕàåÝÞáâØ

// ÁÛãÖÕÑÝëÙ ÚÛÐáá, ÝãÖÕÝ ÝÕ ÒáÕÜ
class crossing:public point3d{
protected:
  int iBorder; // ½ÞÜÕà ÓàÐÝØæë, ÝÐ ÚÞâÞàÞÙ ßàÞØ×ÞèÛÞ ßÕàÕáÕçÕÝØÕ
  cross_info cr;
  ray_parameters rp;
  int iKind; // 0 - ßàÕÛÞÜÛÕÝØÕ, 1 - ÞâàÐÖÕÝØÕ
public:
  crossing(const ray_parameters &rp, cross_info cr_,int ib,
       double x,double y, double z,int ik);
  crossing();
  cross_info crossing_info()const {return cr;};
  ray_parameters ray_pars(){return rp;};
  int border()const{return iBorder;};
  void add_rc_coefs( double c1[2], double c2[2]){cr.add_rc_coefs(c1,c2);};
  void add_rc_vectors(double v1[3],double v2[3]){cr.add_rc_vectors(v1,v2);};
  void add_raypol(double da[3]){cr.add_raypol(da);};
  int reflected() const {return iKind;};
  int operator =(const crossing&);
  void set_refl_type(int ik);
  void set_refrat(double rr);
  int kind() {return iKind;};
  int get_wave_type() {return cr.wave_type;};
  ray_parameters rp_as_sp()const;
};

typedef crossing* p_crossing;

#define AnyType crossing
#define AnyAr   crossings_ptrs
#include <mix/AnyAr.H>

struct meta_iis_element{
  double dTime;
  double dValueX[2];
  double dValueY[2];
  double dValueZ[2];
};

#define AnyType meta_iis_element
#define AnyAr   meta_iis
#include <mix/AnyAr.H>



/*
  ±Õ×ÚÛÐááÞÒëÕ äãÝÚæØØ 
 */
void general_get_ray_parameters (const point3d &pWhere,
                 const crossing &cClosest,
                 double &dFi,double &dAz,
                 double &dDFi,double& dDAz,double &dRo,
                 double dfl=0.8,double dal=0.8);

double compare_reflected(const crossing&,const crossing&);
double compare_transmitted(const crossing&,const crossing&);
void get_ray_angles (double& dFi,double& dAz,const point&pStart,
             const point& pCurr,const geophys_params&gp);
#endif /* p_raytr_elements.H */

/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 31-Jan-2000  */ 
/*          ...          */

  cross_info(double t,
         double dDx1,double dDy1,double dDz1,
         double dDx2,double dDy2,double dDz2,
         double dF  ,double dA,
         //  double dkdfi[3],double dkdt[3], double ksi[3], double vp[3],
         const double *dkdfi,const double *dkdt,const double *ksi,
         const double *vp,double *daDpdteta0,double *daDpdfi0,int ib,
         int w_t);
  void add_rc_coefs(double a[2],double b[2]);
  void add_rc_vectors(double[3],double[3]);
  void add_raypol(double da[3]);
  // ¿àÞØ×ÒÞÔÝëÕ ßÞ ãÓÛÐÜ
  double dDxDf,dDyDf,dDzDf;
  double dDxDaz,dDyDaz,dDzDaz;

  // ÁÞÑáâÒÕÝÝÞ ãÓÛë  double dFi,dAz;

  // ÁÛãÖÕÑÝÐï ÔÛï ¸ÓÞàï
  double dksidfi[3];
  double dksidteta[3];
  double ksi[3];
  int iBody; // ½ÞÜÕà âÕÛÐ, ßÞ ÚÞâÞàÞÜã Ò íâÞ ßÕàÕáÕçÕÝØÕ Ûãç ßàØèÕÛ
  int wave_type; // âØß ÒÞÛÝë ÝÐ ÔÐÝÝÞÜ áÕÓÜÕÝâÕ ÛãçÐ

  // ¸ÝäÞàÜÐæØï, ÝÕÞÑåÞÔØÜÐï áÞÑáâÒÕÝÝÞ ÔÛï áÕÙáÜÞÓàÐÜÜ

  double vpol[3];
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï Ø ßàÕÛÞÜÛÕÝØï, ÚÞÜßÛÕÚáÝëÕ
  double dCRef[2],dCTrans[2];
  // »ãçÕÒëÕ ßÐàÐÜÕâàë ßàÞåÞÖÔÕÝØï Ø ßàÕÛÞÜÛÕÝØï ÝÐ ÓàÐÝØæÕ, ßÞáÛÕ refkin
  double vref[3],vtrans[3];
  // ²ÞÕÝÝÞÕ çØáÕÛÚÞ ßÞ ßÞÒÞÔã ÞâÝÞèÕÝØï ÚÞíääØæØÕÝâÞÒ ÞâàÐÖÕÝØï ÔÛï ¸ÓÞàï
  double dRefRat;

  double daDpdteta0[3];
  double daDpdfi0[3];

  // ²ÕÚâÞà ÔÛï ßÞÛïàØ×ÐæØØ ßÞßÕàÕçÝÞÙ ÒÞÛÝë ßÞâàÕÑÝëÙ...
  double daRayPol[3];
  int operator =(const cross_info&);
};


// ¾ßØáÐÝØÕ ÞÔÝÞÓÞ ßÕàÕáÕçÕÝØï ÛãçÐ á ÓàÐÝØæÕÙ

// ¾àÓÐÝØ×ÐæØï ÝÐÑÞàÐ ßÕàÕáÕçÕÝØÙ ÝÐ ßÞÒÕàåÝÞáâØ

// ÁÛãÖÕÑÝëÙ ÚÛÐáá, ÝãÖÕÝ ÝÕ ÒáÕÜ
class crossing:public point3d{
protected:
  int iBorder; // ½ÞÜÕà ÓàÐÝØæë, ÝÐ ÚÞâÞàÞÙ ßàÞØ×ÞèÛÞ ßÕàÕáÕçÕÝØÕ
  cross_info cr;
  ray_parameters rp;
  int iKind; // 0 - ßàÕÛÞÜÛÕÝØÕ, 1 - ÞâàÐÖÕÝØÕ
public:
  crossing(const ray_parameters &rp, cross_info cr_,int ib,
       double x,double y, double z,int ik);
  crossing();
  cross_info crossing_info()const {return cr;};
  ray_parameters ray_pars(){return rp;};
  int border()const{return iBorder;};
  void add_rc_coefs( double c1[2], double c2[2]){cr.add_rc_coefs(c1,c2);};
  void add_rc_vectors(double v1[3],double v2[3]){cr.add_rc_vectors(v1,v2);};
  void add_raypol(double da[3]){cr.add_raypol(da);};
  int reflected() const {return iKind;};
  int operator =(const crossing&);
  void set_refl_type(int ik);
  void set_refrat(double rr);
  int kind() {return iKind;};
  int get_wave_type() {return cr.wave_type;};
  ray_parameters rp_as_sp()const;
};

typedef crossing* p_crossing;

#define AnyType crossing
#define AnyAr   crossings_ptrs
#include <mix/AnyAr.H>

struct meta_iis_element{
  double dTime;
  double dValueX[2];
  double dValueY[2];
  double dValueZ[2];
};

#define AnyType meta_iis_element
#define AnyAr   meta_iis
#include <mix/AnyAr.H>



/*
  ±Õ×ÚÛÐááÞÒëÕ äãÝÚæØØ 
 */
void general_get_ray_parameters (const point3d &pWhere,
                 const crossing &cClosest,
                 double &dFi,double &dAz,
                 double &dDFi,double& dDAz,double &dRo,
                 double dfl=0.8,double dal=0.8);

double compare_reflected(const crossing&,const crossing&);
double compare_transmitted(const crossing&,const crossing&);
void get_ray_angles (double& dFi,double& dAz,const point&pStart,
             const point& pCurr,const geophys_params&gp);
#endif /* p_raytr_elements.H */

/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 31-Jan-2000  */ 
/*          ...          */

  double dFi,dAz;

  // ÁÛãÖÕÑÝÐï ÔÛï ¸ÓÞàï  double dksidfi[3];
  double dksidteta[3];
  double ksi[3];
  int iBody; // ½ÞÜÕà âÕÛÐ, ßÞ ÚÞâÞàÞÜã Ò íâÞ ßÕàÕáÕçÕÝØÕ Ûãç ßàØèÕÛ
  int wave_type; // âØß ÒÞÛÝë ÝÐ ÔÐÝÝÞÜ áÕÓÜÕÝâÕ ÛãçÐ

  // ¸ÝäÞàÜÐæØï, ÝÕÞÑåÞÔØÜÐï áÞÑáâÒÕÝÝÞ ÔÛï áÕÙáÜÞÓàÐÜÜ

  double vpol[3];
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï Ø ßàÕÛÞÜÛÕÝØï, ÚÞÜßÛÕÚáÝëÕ
  double dCRef[2],dCTrans[2];
  // »ãçÕÒëÕ ßÐàÐÜÕâàë ßàÞåÞÖÔÕÝØï Ø ßàÕÛÞÜÛÕÝØï ÝÐ ÓàÐÝØæÕ, ßÞáÛÕ refkin
  double vref[3],vtrans[3];
  // ²ÞÕÝÝÞÕ çØáÕÛÚÞ ßÞ ßÞÒÞÔã ÞâÝÞèÕÝØï ÚÞíääØæØÕÝâÞÒ ÞâàÐÖÕÝØï ÔÛï ¸ÓÞàï
  double dRefRat;

  double daDpdteta0[3];
  double daDpdfi0[3];

  // ²ÕÚâÞà ÔÛï ßÞÛïàØ×ÐæØØ ßÞßÕàÕçÝÞÙ ÒÞÛÝë ßÞâàÕÑÝëÙ...
  double daRayPol[3];
  int operator =(const cross_info&);
};


// ¾ßØáÐÝØÕ ÞÔÝÞÓÞ ßÕàÕáÕçÕÝØï ÛãçÐ á ÓàÐÝØæÕÙ

// ¾àÓÐÝØ×ÐæØï ÝÐÑÞàÐ ßÕàÕáÕçÕÝØÙ ÝÐ ßÞÒÕàåÝÞáâØ

// ÁÛãÖÕÑÝëÙ ÚÛÐáá, ÝãÖÕÝ ÝÕ ÒáÕÜ
class crossing:public point3d{
protected:
  int iBorder; // ½ÞÜÕà ÓàÐÝØæë, ÝÐ ÚÞâÞàÞÙ ßàÞØ×ÞèÛÞ ßÕàÕáÕçÕÝØÕ
  cross_info cr;
  ray_parameters rp;
  int iKind; // 0 - ßàÕÛÞÜÛÕÝØÕ, 1 - ÞâàÐÖÕÝØÕ
public:
  crossing(const ray_parameters &rp, cross_info cr_,int ib,
       double x,double y, double z,int ik);
  crossing();
  cross_info crossing_info()const {return cr;};
  ray_parameters ray_pars(){return rp;};
  int border()const{return iBorder;};
  void add_rc_coefs( double c1[2], double c2[2]){cr.add_rc_coefs(c1,c2);};
  void add_rc_vectors(double v1[3],double v2[3]){cr.add_rc_vectors(v1,v2);};
  void add_raypol(double da[3]){cr.add_raypol(da);};
  int reflected() const {return iKind;};
  int operator =(const crossing&);
  void set_refl_type(int ik);
  void set_refrat(double rr);
  int kind() {return iKind;};
  int get_wave_type() {return cr.wave_type;};
  ray_parameters rp_as_sp()const;
};

typedef crossing* p_crossing;

#define AnyType crossing
#define AnyAr   crossings_ptrs
#include <mix/AnyAr.H>

struct meta_iis_element{
  double dTime;
  double dValueX[2];
  double dValueY[2];
  double dValueZ[2];
};

#define AnyType meta_iis_element
#define AnyAr   meta_iis
#include <mix/AnyAr.H>



/*
  ±Õ×ÚÛÐááÞÒëÕ äãÝÚæØØ 
 */
void general_get_ray_parameters (const point3d &pWhere,
                 const crossing &cClosest,
                 double &dFi,double &dAz,
                 double &dDFi,double& dDAz,double &dRo,
                 double dfl=0.8,double dal=0.8);

double compare_reflected(const crossing&,const crossing&);
double compare_transmitted(const crossing&,const crossing&);
void get_ray_angles (double& dFi,double& dAz,const point&pStart,
             const point& pCurr,const geophys_params&gp);
#endif /* p_raytr_elements.H */

/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 31-Jan-2000  */ 
/*          ...          */

  double dksidfi[3];
  double dksidteta[3];
  double ksi[3];
  int iBody; // ½ÞÜÕà âÕÛÐ, ßÞ ÚÞâÞàÞÜã Ò íâÞ ßÕàÕáÕçÕÝØÕ Ûãç ßàØèÕÛ
  int wave_type; // âØß ÒÞÛÝë ÝÐ ÔÐÝÝÞÜ áÕÓÜÕÝâÕ ÛãçÐ
  // ¸ÝäÞàÜÐæØï, ÝÕÞÑåÞÔØÜÐï áÞÑáâÒÕÝÝÞ ÔÛï áÕÙáÜÞÓàÐÜÜ

  double vpol[3];
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï Ø ßàÕÛÞÜÛÕÝØï, ÚÞÜßÛÕÚáÝëÕ
  double dCRef[2],dCTrans[2];
  // »ãçÕÒëÕ ßÐàÐÜÕâàë ßàÞåÞÖÔÕÝØï Ø ßàÕÛÞÜÛÕÝØï ÝÐ ÓàÐÝØæÕ, ßÞáÛÕ refkin
  double vref[3],vtrans[3];
  // ²ÞÕÝÝÞÕ çØáÕÛÚÞ ßÞ ßÞÒÞÔã ÞâÝÞèÕÝØï ÚÞíääØæØÕÝâÞÒ ÞâàÐÖÕÝØï ÔÛï ¸ÓÞàï
  double dRefRat;

  double daDpdteta0[3];
  double daDpdfi0[3];

  // ²ÕÚâÞà ÔÛï ßÞÛïàØ×ÐæØØ ßÞßÕàÕçÝÞÙ ÒÞÛÝë ßÞâàÕÑÝëÙ...
  double daRayPol[3];
  int operator =(const cross_info&);
};


// ¾ßØáÐÝØÕ ÞÔÝÞÓÞ ßÕàÕáÕçÕÝØï ÛãçÐ á ÓàÐÝØæÕÙ

// ¾àÓÐÝØ×ÐæØï ÝÐÑÞàÐ ßÕàÕáÕçÕÝØÙ ÝÐ ßÞÒÕàåÝÞáâØ

// ÁÛãÖÕÑÝëÙ ÚÛÐáá, ÝãÖÕÝ ÝÕ ÒáÕÜ
class crossing:public point3d{
protected:
  int iBorder; // ½ÞÜÕà ÓàÐÝØæë, ÝÐ ÚÞâÞàÞÙ ßàÞØ×ÞèÛÞ ßÕàÕáÕçÕÝØÕ
  cross_info cr;
  ray_parameters rp;
  int iKind; // 0 - ßàÕÛÞÜÛÕÝØÕ, 1 - ÞâàÐÖÕÝØÕ
public:
  crossing(const ray_parameters &rp, cross_info cr_,int ib,
       double x,double y, double z,int ik);
  crossing();
  cross_info crossing_info()const {return cr;};
  ray_parameters ray_pars(){return rp;};
  int border()const{return iBorder;};
  void add_rc_coefs( double c1[2], double c2[2]){cr.add_rc_coefs(c1,c2);};
  void add_rc_vectors(double v1[3],double v2[3]){cr.add_rc_vectors(v1,v2);};
  void add_raypol(double da[3]){cr.add_raypol(da);};
  int reflected() const {return iKind;};
  int operator =(const crossing&);
  void set_refl_type(int ik);
  void set_refrat(double rr);
  int kind() {return iKind;};
  int get_wave_type() {return cr.wave_type;};
  ray_parameters rp_as_sp()const;
};

typedef crossing* p_crossing;

#define AnyType crossing
#define AnyAr   crossings_ptrs
#include <mix/AnyAr.H>

struct meta_iis_element{
  double dTime;
  double dValueX[2];
  double dValueY[2];
  double dValueZ[2];
};

#define AnyType meta_iis_element
#define AnyAr   meta_iis
#include <mix/AnyAr.H>



/*
  ±Õ×ÚÛÐááÞÒëÕ äãÝÚæØØ 
 */
void general_get_ray_parameters (const point3d &pWhere,
                 const crossing &cClosest,
                 double &dFi,double &dAz,
                 double &dDFi,double& dDAz,double &dRo,
                 double dfl=0.8,double dal=0.8);

double compare_reflected(const crossing&,const crossing&);
double compare_transmitted(const crossing&,const crossing&);
void get_ray_angles (double& dFi,double& dAz,const point&pStart,
             const point& pCurr,const geophys_params&gp);
#endif /* p_raytr_elements.H */

/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 31-Jan-2000  */ 
/*          ...          */


  // ¸ÝäÞàÜÐæØï, ÝÕÞÑåÞÔØÜÐï áÞÑáâÒÕÝÝÞ ÔÛï áÕÙáÜÞÓàÐÜÜ
  double vpol[3];
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï Ø ßàÕÛÞÜÛÕÝØï, ÚÞÜßÛÕÚáÝëÕ
  double dCRef[2],dCTrans[2];
  // »ãçÕÒëÕ ßÐàÐÜÕâàë ßàÞåÞÖÔÕÝØï Ø ßàÕÛÞÜÛÕÝØï ÝÐ ÓàÐÝØæÕ, ßÞáÛÕ refkin
  double vref[3],vtrans[3];
  // ²ÞÕÝÝÞÕ çØáÕÛÚÞ ßÞ ßÞÒÞÔã ÞâÝÞèÕÝØï ÚÞíääØæØÕÝâÞÒ ÞâàÐÖÕÝØï ÔÛï ¸ÓÞàï
  double dRefRat;

  double daDpdteta0[3];
  double daDpdfi0[3];

  // ²ÕÚâÞà ÔÛï ßÞÛïàØ×ÐæØØ ßÞßÕàÕçÝÞÙ ÒÞÛÝë ßÞâàÕÑÝëÙ...
  double daRayPol[3];
  int operator =(const cross_info&);
};


// ¾ßØáÐÝØÕ ÞÔÝÞÓÞ ßÕàÕáÕçÕÝØï ÛãçÐ á ÓàÐÝØæÕÙ

// ¾àÓÐÝØ×ÐæØï ÝÐÑÞàÐ ßÕàÕáÕçÕÝØÙ ÝÐ ßÞÒÕàåÝÞáâØ

// ÁÛãÖÕÑÝëÙ ÚÛÐáá, ÝãÖÕÝ ÝÕ ÒáÕÜ
class crossing:public point3d{
protected:
  int iBorder; // ½ÞÜÕà ÓàÐÝØæë, ÝÐ ÚÞâÞàÞÙ ßàÞØ×ÞèÛÞ ßÕàÕáÕçÕÝØÕ
  cross_info cr;
  ray_parameters rp;
  int iKind; // 0 - ßàÕÛÞÜÛÕÝØÕ, 1 - ÞâàÐÖÕÝØÕ
public:
  crossing(const ray_parameters &rp, cross_info cr_,int ib,
       double x,double y, double z,int ik);
  crossing();
  cross_info crossing_info()const {return cr;};
  ray_parameters ray_pars(){return rp;};
  int border()const{return iBorder;};
  void add_rc_coefs( double c1[2], double c2[2]){cr.add_rc_coefs(c1,c2);};
  void add_rc_vectors(double v1[3],double v2[3]){cr.add_rc_vectors(v1,v2);};
  void add_raypol(double da[3]){cr.add_raypol(da);};
  int reflected() const {return iKind;};
  int operator =(const crossing&);
  void set_refl_type(int ik);
  void set_refrat(double rr);
  int kind() {return iKind;};
  int get_wave_type() {return cr.wave_type;};
  ray_parameters rp_as_sp()const;
};

typedef crossing* p_crossing;

#define AnyType crossing
#define AnyAr   crossings_ptrs
#include <mix/AnyAr.H>

struct meta_iis_element{
  double dTime;
  double dValueX[2];
  double dValueY[2];
  double dValueZ[2];
};

#define AnyType meta_iis_element
#define AnyAr   meta_iis
#include <mix/AnyAr.H>



/*
  ±Õ×ÚÛÐááÞÒëÕ äãÝÚæØØ 
 */
void general_get_ray_parameters (const point3d &pWhere,
                 const crossing &cClosest,
                 double &dFi,double &dAz,
                 double &dDFi,double& dDAz,double &dRo,
                 double dfl=0.8,double dal=0.8);

double compare_reflected(const crossing&,const crossing&);
double compare_transmitted(const crossing&,const crossing&);
void get_ray_angles (double& dFi,double& dAz,const point&pStart,
             const point& pCurr,const geophys_params&gp);
#endif /* p_raytr_elements.H */

/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 31-Jan-2000  */ 
/*          ...          */


  double vpol[3];
  // ºÞíääØæØÕÝâë ÞâàÐÖÕÝØï Ø ßàÕÛÞÜÛÕÝØï, ÚÞÜßÛÕÚáÝëÕ  double dCRef[2],dCTrans[2];
  // »ãçÕÒëÕ ßÐàÐÜÕâàë ßàÞåÞÖÔÕÝØï Ø ßàÕÛÞÜÛÕÝØï ÝÐ ÓàÐÝØæÕ, ßÞáÛÕ refkin
  double vref[3],vtrans[3];
  // ²ÞÕÝÝÞÕ çØáÕÛÚÞ ßÞ ßÞÒÞÔã ÞâÝÞèÕÝØï ÚÞíääØæØÕÝâÞÒ ÞâàÐÖÕÝØï ÔÛï ¸ÓÞàï
  double dRefRat;

  double daDpdteta0[3];
  double daDpdfi0[3];

  // ²ÕÚâÞà ÔÛï ßÞÛïàØ×ÐæØØ ßÞßÕàÕçÝÞÙ ÒÞÛÝë ßÞâàÕÑÝëÙ...
  double daRayPol[3];
  int operator =(const cross_info&);
};


// ¾ßØáÐÝØÕ ÞÔÝÞÓÞ ßÕàÕáÕçÕÝØï ÛãçÐ á ÓàÐÝØæÕÙ

// ¾àÓÐÝØ×ÐæØï ÝÐÑÞàÐ ßÕàÕáÕçÕÝØÙ ÝÐ ßÞÒÕàåÝÞáâØ

// ÁÛãÖÕÑÝëÙ ÚÛÐáá, ÝãÖÕÝ ÝÕ ÒáÕÜ
class crossing:public point3d{
protected:
  int iBorder; // ½ÞÜÕà ÓàÐÝØæë, ÝÐ ÚÞâÞàÞÙ ßàÞØ×ÞèÛÞ ßÕàÕáÕçÕÝØÕ
  cross_info cr;
  ray_parameters rp;
  int iKind; // 0 - ßàÕÛÞÜÛÕÝØÕ, 1 - ÞâàÐÖÕÝØÕ
public:
  crossing(const ray_parameters &rp, cross_info cr_,int ib,
       double x,double y, double z,int ik);
  crossing();
  cross_info crossing_info()const {return cr;};
  ray_parameters ray_pars(){return rp;};
  int border()const{return iBorder;};
  void add_rc_coefs( double c1[2], double c2[2]){cr.add_rc_coefs(c1,c2);};
  void add_rc_vectors(double v1[3],double v2[3]){cr.add_rc_vectors(v1,v2);};
  void add_raypol(double da[3]){cr.add_raypol(da);};
  int reflected() const {return iKind;};
  int operator =(const crossing&);
  void set_refl_type(int ik);
  void set_refrat(double rr);
  int kind() {return iKind;};
  int get_wave_type() {return cr.wave_type;};
  ray_parameters rp_as_sp()const;
};

typedef crossing* p_crossing;

#define AnyType crossing
#define AnyAr   crossings_ptrs
#include <mix/AnyAr.H>

struct meta_iis_element{
  double dTime;
  double dValueX[2];
  double dValueY[2];
  double dValueZ[2];
};

#define AnyType meta_iis_element
#define AnyAr   meta_iis
#include <mix/AnyAr.H>



/*
  ±Õ×ÚÛÐááÞÒëÕ äãÝÚæØØ 
 */
void general_get_ray_parameters (const point3d &pWhere,
                 const crossing &cClosest,
                 double &dFi,double &dAz,
                 double &dDFi,double& dDAz,double &dRo,
                 double dfl=0.8,double dal=0.8);

double compare_reflected(const crossing&,const crossing&);
double compare_transmitted(const crossing&,const crossing&);
void get_ray_angles (double& dFi,double& dAz,const point&pStart,
             const point& pCurr,const geophys_params&gp);
#endif /* p_raytr_elements.H */

/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 31-Jan-2000  */ 
/*          ...          */

  double dCRef[2],dCTrans[2];
  // »ãçÕÒëÕ ßÐàÐÜÕâàë ßàÞåÞÖÔÕÝØï Ø ßàÕÛÞÜÛÕÝØï ÝÐ ÓàÐÝØæÕ, ßÞáÛÕ refkin
  double vref[3],vtrans[3];
  // ²ÞÕÝÝÞÕ çØáÕÛÚÞ ßÞ ßÞÒÞÔã ÞâÝÞèÕÝØï ÚÞíääØæØÕÝâÞÒ ÞâàÐÖÕÝØï ÔÛï ¸ÓÞàï  double dRefRat;

  double daDpdteta0[3];
  double daDpdfi0[3];

  // ²ÕÚâÞà ÔÛï ßÞÛïàØ×ÐæØØ ßÞßÕàÕçÝÞÙ ÒÞÛÝë ßÞâàÕÑÝëÙ...
  double daRayPol[3];
  int operator =(const cross_info&);
};


// ¾ßØáÐÝØÕ ÞÔÝÞÓÞ ßÕàÕáÕçÕÝØï ÛãçÐ á ÓàÐÝØæÕÙ

// ¾àÓÐÝØ×ÐæØï ÝÐÑÞàÐ ßÕàÕáÕçÕÝØÙ ÝÐ ßÞÒÕàåÝÞáâØ

// ÁÛãÖÕÑÝëÙ ÚÛÐáá, ÝãÖÕÝ ÝÕ ÒáÕÜ
class crossing:public point3d{
protected:
  int iBorder; // ½ÞÜÕà ÓàÐÝØæë, ÝÐ ÚÞâÞàÞÙ ßàÞØ×ÞèÛÞ ßÕàÕáÕçÕÝØÕ
  cross_info cr;
  ray_parameters rp;
  int iKind; // 0 - ßàÕÛÞÜÛÕÝØÕ, 1 - ÞâàÐÖÕÝØÕ
public:
  crossing(const ray_parameters &rp, cross_info cr_,int ib,
       double x,double y, double z,int ik);
  crossing();
  cross_info crossing_info()const {return cr;};
  ray_parameters ray_pars(){return rp;};
  int border()const{return iBorder;};
  void add_rc_coefs( double c1[2], double c2[2]){cr.add_rc_coefs(c1,c2);};
  void add_rc_vectors(double v1[3],double v2[3]){cr.add_rc_vectors(v1,v2);};
  void add_raypol(double da[3]){cr.add_raypol(da);};
  int reflected() const {return iKind;};
  int operator =(const crossing&);
  void set_refl_type(int ik);
  void set_refrat(double rr);
  int kind() {return iKind;};
  int get_wave_type() {return cr.wave_type;};
  ray_parameters rp_as_sp()const;
};

typedef crossing* p_crossing;

#define AnyType crossing
#define AnyAr   crossings_ptrs
#include <mix/AnyAr.H>

struct meta_iis_element{
  double dTime;
  double dValueX[2];
  double dValueY[2];
  double dValueZ[2];
};

#define AnyType meta_iis_element
#define AnyAr   meta_iis
#include <mix/AnyAr.H>



/*
  ±Õ×ÚÛÐááÞÒëÕ äãÝÚæØØ 
 */
void general_get_ray_parameters (const point3d &pWhere,
                 const crossing &cClosest,
                 double &dFi,double &dAz,
                 double &dDFi,double& dDAz,double &dRo,
                 double dfl=0.8,double dal=0.8);

double compare_reflected(const crossing&,const crossing&);
double compare_transmitted(const crossing&,const crossing&);
void get_ray_angles (double& dFi,double& dAz,const point&pStart,
             const point& pCurr,const geophys_params&gp);
#endif /* p_raytr_elements.H */

/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 31-Jan-2000  */ 
/*          ...          */

  double dRefRat;

  double daDpdteta0[3];
  double daDpdfi0[3];

  // ²ÕÚâÞà ÔÛï ßÞÛïàØ×ÐæØØ ßÞßÕàÕçÝÞÙ ÒÞÛÝë ßÞâàÕÑÝëÙ...
  double daRayPol[3];
  int operator =(const cross_info&);
};


// ¾ßØáÐÝØÕ ÞÔÝÞÓÞ ßÕàÕáÕçÕÝØï ÛãçÐ á ÓàÐÝØæÕÙ

// ¾àÓÐÝØ×ÐæØï ÝÐÑÞàÐ ßÕàÕáÕçÕÝØÙ ÝÐ ßÞÒÕàåÝÞáâØ
// ÁÛãÖÕÑÝëÙ ÚÛÐáá, ÝãÖÕÝ ÝÕ ÒáÕÜ
class crossing:public point3d{
protected:
  int iBorder; // ½ÞÜÕà ÓàÐÝØæë, ÝÐ ÚÞâÞàÞÙ ßàÞØ×ÞèÛÞ ßÕàÕáÕçÕÝØÕ
  cross_info cr;
  ray_parameters rp;
  int iKind; // 0 - ßàÕÛÞÜÛÕÝØÕ, 1 - ÞâàÐÖÕÝØÕ
public:
  crossing(const ray_parameters &rp, cross_info cr_,int ib,
       double x,double y, double z,int ik);
  crossing();
  cross_info crossing_info()const {return cr;};
  ray_parameters ray_pars(){return rp;};
  int border()const{return iBorder;};
  void add_rc_coefs( double c1[2], double c2[2]){cr.add_rc_coefs(c1,c2);};
  void add_rc_vectors(double v1[3],double v2[3]){cr.add_rc_vectors(v1,v2);};
  void add_raypol(double da[3]){cr.add_raypol(da);};
  int reflected() const {return iKind;};
  int operator =(const crossing&);
  void set_refl_type(int ik);
  void set_refrat(double rr);
  int kind() {return iKind;};
  int get_wave_type() {return cr.wave_type;};
  ray_parameters rp_as_sp()const;
};

typedef crossing* p_crossing;

#define AnyType crossing
#define AnyAr   crossings_ptrs
#include <mix/AnyAr.H>

struct meta_iis_element{
  double dTime;
  double dValueX[2];
  double dValueY[2];
  double dValueZ[2];
};

#define AnyType meta_iis_element
#define AnyAr   meta_iis
#include <mix/AnyAr.H>



/*
  ±Õ×ÚÛÐááÞÒëÕ äãÝÚæØØ 
 */
void general_get_ray_parameters (const point3d &pWhere,
                 const crossing &cClosest,
                 double &dFi,double &dAz,
                 double &dDFi,double& dDAz,double &dRo,
                 double dfl=0.8,double dal=0.8);

double compare_reflected(const crossing&,const crossing&);
double compare_transmitted(const crossing&,const crossing&);
void get_ray_angles (double& dFi,double& dAz,const point&pStart,
             const point& pCurr,const geophys_params&gp);
#endif /* p_raytr_elements.H */

/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 31-Jan-2000  */ 
/*          ...          */


// ÁÛãÖÕÑÝëÙ ÚÛÐáá, ÝãÖÕÝ ÝÕ ÒáÕÜ
class crossing:public point3d{
protected:
  int iBorder; // ½ÞÜÕà ÓàÐÝØæë, ÝÐ ÚÞâÞàÞÙ ßàÞØ×ÞèÛÞ ßÕàÕáÕçÕÝØÕ
  cross_info cr;
  ray_parameters rp;
  int iKind; // 0 - ßàÕÛÞÜÛÕÝØÕ, 1 - ÞâàÐÖÕÝØÕ
public:
  crossing(const ray_parameters &rp, cross_info cr_,int ib,
       double x,double y, double z,int ik);
  crossing();
  cross_info crossing_info()const {return cr;};
  ray_parameters ray_pars(){return rp;};
  int border()const{return iBorder;};
  void add_rc_coefs( double c1[2], double c2[2]){cr.add_rc_coefs(c1,c2);};
  void add_rc_vectors(double v1[3],double v2[3]){cr.add_rc_vectors(v1,v2);};
  void add_raypol(double da[3]){cr.add_raypol(da);};
  int reflected() const {return iKind;};
  int operator =(const crossing&);
  void set_refl_type(int ik);
  void set_refrat(double rr);
  int kind() {return iKind;};
  int get_wave_type() {return cr.wave_type;};
  ray_parameters rp_as_sp()const;
};

typedef crossing* p_crossing;

#define AnyType crossing
#define AnyAr   crossings_ptrs
#include <mix/AnyAr.H>

struct meta_iis_element{
  double dTime;
  double dValueX[2];
  double dValueY[2];
  double dValueZ[2];
};

#define AnyType meta_iis_element
#define AnyAr   meta_iis
#include <mix/AnyAr.H>



/*
  ±Õ×ÚÛÐááÞÒëÕ äãÝÚæØØ 
 */
void general_get_ray_parameters (const point3d &pWhere,
                 const crossing &cClosest,
                 double &dFi,double &dAz,
                 double &dDFi,double& dDAz,double &dRo,
                 double dfl=0.8,double dal=0.8);

double compare_reflected(const crossing&,const crossing&);
double compare_transmitted(const crossing&,const crossing&);
void get_ray_angles (double& dFi,double& dAz,const point&pStart,
             const point& pCurr,const geophys_params&gp);
#endif /* p_raytr_elements.H */

/*  Melnikov George Y.   2-Nov-1998  */
/*          ...          */
/*  Melnikov George Yu. 31-Jan-2000  */ 
/*          ...          */