any_elements.H

/* any_elements.H */
/* $Id: any_elements.H,v 2.21 2001/10/09 14:54:38 hoh Exp $ */
#ifndef __any_elements_H
#define __any_elements_H

#include <mix/IntAr.H>
#include <mix/DoubleAr.H>
#include <mix/ObjectsAr.H>
#include <mix/geometry.H>

#include "misc.H"

// ²ØÔ ÓàÐÝØæë: ÞÑëçÝÐï, ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï, áÒÞÑÞÔÝÐï ßÞÒÕàåÝÞáâìtypedef enum{NORMAL=0,MODEL=1,TOP=2, TEMP=3} border_flag;

// ²Þ×ÜÞÖÝëÕ âØßë ÒÞÛÝ - ÝØÚÐÚÐï ÛØÑÞ ÚÐÚÐï-âÞ ;)
#define PM_NOKIND   0
#define PM_P        1
#define PM_SV       2

// ¸ÝÚÐßáãÛïæØï ÓÕÞäØ×ØçÕáÚØå ßÐàÐÜÕâàÞÒ
typedef struct{
  double fVP;
  double fVS;
  double fDensity;
  double fAn_Azimuth;
  double fAn_Dip;
  double fAn_KP;
  double fAn_KS;
} geophys_params;

// ¸ÝÚÐßáãÛïæØï ÔÐÝÝëå ÔÛï ÝÐáâàÕÛÐ ÛãçÕÙ (ÕéÕ ÔÞÑÐÒØÜ ÒáïçÕáÚÞÙ àÐ×ÝÞáâØ...)
typedef struct{
  double dAniz[3];
  double dVP_Aniz;
  double dVS_Aniz;
  double*cl;
  double*cg;
  double*rt;
} raytr_params;

// Data incapsulation to store parallel model parameters
class parallel_model_params{
public:

  DoubleAr pVP;
  DoubleAr pVS;
  DoubleAr pKP;
  DoubleAr pKS;
  DoubleAr xTop;
  DoubleAr yTop;
  DoubleAr zTop;
  double normal[3];
  double *e,*d,*f,*h,*v,*ztop,*plr,*segm,*aux,*tgpsic;

  parallel_model_params ()
    :e(NULL),d(NULL),f(NULL),h(NULL),v(NULL),ztop(NULL),
     plr(NULL),segm(NULL),aux(NULL),tgpsic(NULL)
  {
    clean();
  };

  ~parallel_model_params (){
    if(NULL!=e) free(e);
    if(NULL!=d) free(d);
    if(NULL!=f) free(f);
    if(NULL!=h) free(h);
    if(NULL!=v) free(v);
    if(NULL!=ztop) free(ztop);
    if(NULL!=plr) free(plr);
    if(NULL!=segm) free(segm);
    if(NULL!=aux) free(aux);
    if(NULL!=tgpsic) free(tgpsic);
  };

  void clean (){
    pVP.clean(); pVS.clean(); pKP.clean(); pKS.clean();
    xTop.clean(); yTop.clean(); zTop.clean();
  };

  void allocate_mem_buf (int count){
    e=(double*) calloc(2*(count+1), sizeof(double));
    d=(double*) calloc(2*(count+1), sizeof(double));
    f=(double*) calloc(2*(count+1), sizeof(double));
    h=(double*) calloc(2*(count+1), sizeof(double));
    v=(double*) calloc(2*(count+1), sizeof(double));
    ztop=(double*) calloc(2*(count+1), sizeof(double));
    plr=(double*) calloc(2*(count+1), sizeof(double));
    segm=(double*) calloc(2*(count+1), sizeof(double));
    aux=(double*) calloc(2*(count+1), sizeof(double));
    tgpsic=(double*) calloc(2*(count+1), sizeof(double));    
  };
};

/****** ÍÛÕÜÕÝâë ÒÝãâàÕÝÝÕÓÞ ßàÕÔáâÐÒÛÕÝØï ÜÞÔÕÛØ ******/

// ¿àÞØ×ÒÞÛìÝÐï ÓÛÐÔÚÐï ßÞÒÕàåÝÞáâì - áßÛÐÙÝ
class any_spline{
protected:
  any_spline(){;};
};

typedef any_spline* p_any_spline;
typedef ObjectsAr<p_any_spline> any_splines;

// ¿àÞØ×ÒÞÛìÝÐï ÓàÐÝØæÐ àÞÒÝÞ ÜÕÖÔã ÔÒãÜï âÕÛÐÜØ
  
class any_border{
protected:
  border_flag bf;
  int iBodies[2];
  p_any_spline psSpline;
  any_border(any_border &bCopy);
  any_border(int ib[2],const p_any_spline ps,border_flag bf);
  any_border(const p_any_spline ps,border_flag bf);
public:
  int get_body(int ib) const;
  p_any_spline  spline()const{return psSpline;} ;
  void set_bodies(int i1,int i2){iBodies[0]=i1;iBodies[1]=i2;};
  border_flag kind() {return bf;};
  int another_body(int iB){
    if(iBodies[0]==iB)return iBodies[1];else
      if (iBodies[1]==iB)return iBodies[0];
    return -1;
  };
  ~any_border();
};

typedef any_border* p_any_border;
typedef ObjectsAr<p_any_border> any_borders;


// ÂÕÛÞ ÚÐÚ ÝÐÑÞà áÞÑáâÒÕÝÝëå ÓàÐÝØæ
class any_body{
protected:
  int iNumber;
  any_borders aBorders;
  IntAr aBordersNums;
  geophys_params gp;
  raytr_params rp;
  any_body(geophys_params _gp);
  any_body(any_body& bCopy);
  any_body(any_borders &ab,geophys_params gp);  
  void swap_borders(int,int);
  void gp2rp();
public:
  Logic numerated_borders()const{
    if(aBordersNums.count()>0)
      return TRUE;
    return FALSE;
  };
  void clear_borders();
  void set_number(int i){iNumber=i;};
  int n_of_borders() const {return aBorders.count();};
  void add_border(p_any_border pb){aBorders.addh(pb);};
  ~any_body();
  inline p_any_border get_border(int i)const ;
  geophys_params get_geophys() const {return gp;};
  void set_geophys(geophys_params GP) {gp=GP; gp2rp();};
  raytr_params get_raytr()const{return rp;};
  Logic contains_border(p_any_border pb)const;

  // ¾çÕÝÝÞ ÞßÐáÝÐï äãÝÚæØï, ÞÔÝÐÚÞ...
  void fix_numbers(const IntAr&aBorderNumbers);
  const IntAr& borders_numbers()const{return aBordersNums;};
  // ²ëÔÐÕâ áßØáÞÚ ÝÞÜÕàÞÒ âÕÛ-áÞáÕÔÕÙ (ÚàÞÜÕ -1, ÕááÝÞ)
  IntAr neighbors(int iNumber) const;

  // ·ÐÜÕÝØÛØ ÞÔÝã ÓàÐÝØæã ÝÐ ÔàãÓãî, ÕáÛØ çâÞ - ÚØÝãÛØ ØáÚÛîçÕÝØÕ
  void change_border(p_any_border,p_any_border); 

  // ¾âáÞàâØàÞÒÐâì ÓàÐÝØæë
  void sort_borders();
  int number()const{return iNumber;};

  void set_vp(double vp){gp.fVP=vp;gp2rp();};
  void set_kp(double kp){gp.fAn_KP=kp;gp2rp();};
  void set_vs(double vs){gp.fVS=vs;gp2rp();};
  void set_ks(double ks){gp.fAn_KS=ks;gp2rp();};  

  const any_borders& borders() const {return aBorders;};
};

// ¾çÕÝì ãÖ çÐáâÞ ÞÝÐ Òë×ëÒÐÕâáï...
inline p_any_border any_body::get_border (int i)const{
  // ¿àÞÒÕàÚØ ÝÐäØÓ, ßãáâì ObjectAr ßàÞÒÕàïÕâ...
#if 0
  if((i>=aBorders.count())||(i<0))
    return NULL;
#endif
  return aBorders(i);
};


typedef any_body* p_any_body;
typedef ObjectsAr<p_any_body> any_bodies;

// ÍÛÕÜÕÝâë ÒÝÕèÝÕÓÞ ßàÕÔáâÐÒÛÕÝØï ÜÞÔÕÛØ

// ÁßÛÐÙÝ
class any_ex_spline{
protected:
  any_ex_spline(){;};
};

// ³àÐÝØæÐ
class any_ex_border{
protected:
  border_flag bf;
  int iSpline;
  int iBodies[2];
  any_ex_border(int ib[2], int is, border_flag _bf);
  any_ex_border(){iSpline=-1; bf=NORMAL; iBodies[0]=iBodies[1]=-5;};
public:
  int get_body(int ib) ;
  int spline()const {return iSpline;};
  border_flag kind() {return bf;};
};

// ÂÕÛÞ
class any_ex_body{
protected:
  int iNumber;
  IntAr aBorders;
  geophys_params gp;
  raytr_params rp;
  any_ex_body(IntAr & ab, geophys_params const & gp, int in);
  any_ex_body(){;};
public:
  
  virtual int n_of_borders()const {return aBorders.count();};
  int make_body(const IntAr& ab, const geophys_params & gp,
        const raytr_params & rp,int in);
  geophys_params get_geophys()const{return gp;};
  raytr_params get_raytr()const{return rp;};
  int get_border(int i);
  int number()const{return iNumber;};
};

#endif /* any_elements.H */

/*  Melnikov George Y.  29-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 23-Nov-1999  */ 
/*          ...          */ 

typedef enum{NORMAL=0,MODEL=1,TOP=2, TEMP=3} border_flag;

// ²Þ×ÜÞÖÝëÕ âØßë ÒÞÛÝ - ÝØÚÐÚÐï ÛØÑÞ ÚÐÚÐï-âÞ ;)
#define PM_NOKIND   0
#define PM_P        1
#define PM_SV       2

// ¸ÝÚÐßáãÛïæØï ÓÕÞäØ×ØçÕáÚØå ßÐàÐÜÕâàÞÒ
typedef struct{
  double fVP;
  double fVS;
  double fDensity;
  double fAn_Azimuth;
  double fAn_Dip;
  double fAn_KP;
  double fAn_KS;
} geophys_params;

// ¸ÝÚÐßáãÛïæØï ÔÐÝÝëå ÔÛï ÝÐáâàÕÛÐ ÛãçÕÙ (ÕéÕ ÔÞÑÐÒØÜ ÒáïçÕáÚÞÙ àÐ×ÝÞáâØ...)
typedef struct{
  double dAniz[3];
  double dVP_Aniz;
  double dVS_Aniz;
  double*cl;
  double*cg;
  double*rt;
} raytr_params;

// Data incapsulation to store parallel model parameters
class parallel_model_params{
public:

  DoubleAr pVP;
  DoubleAr pVS;
  DoubleAr pKP;
  DoubleAr pKS;
  DoubleAr xTop;
  DoubleAr yTop;
  DoubleAr zTop;
  double normal[3];
  double *e,*d,*f,*h,*v,*ztop,*plr,*segm,*aux,*tgpsic;

  parallel_model_params ()
    :e(NULL),d(NULL),f(NULL),h(NULL),v(NULL),ztop(NULL),
     plr(NULL),segm(NULL),aux(NULL),tgpsic(NULL)
  {
    clean();
  };

  ~parallel_model_params (){
    if(NULL!=e) free(e);
    if(NULL!=d) free(d);
    if(NULL!=f) free(f);
    if(NULL!=h) free(h);
    if(NULL!=v) free(v);
    if(NULL!=ztop) free(ztop);
    if(NULL!=plr) free(plr);
    if(NULL!=segm) free(segm);
    if(NULL!=aux) free(aux);
    if(NULL!=tgpsic) free(tgpsic);
  };

  void clean (){
    pVP.clean(); pVS.clean(); pKP.clean(); pKS.clean();
    xTop.clean(); yTop.clean(); zTop.clean();
  };

  void allocate_mem_buf (int count){
    e=(double*) calloc(2*(count+1), sizeof(double));
    d=(double*) calloc(2*(count+1), sizeof(double));
    f=(double*) calloc(2*(count+1), sizeof(double));
    h=(double*) calloc(2*(count+1), sizeof(double));
    v=(double*) calloc(2*(count+1), sizeof(double));
    ztop=(double*) calloc(2*(count+1), sizeof(double));
    plr=(double*) calloc(2*(count+1), sizeof(double));
    segm=(double*) calloc(2*(count+1), sizeof(double));
    aux=(double*) calloc(2*(count+1), sizeof(double));
    tgpsic=(double*) calloc(2*(count+1), sizeof(double));    
  };
};

/****** ÍÛÕÜÕÝâë ÒÝãâàÕÝÝÕÓÞ ßàÕÔáâÐÒÛÕÝØï ÜÞÔÕÛØ ******/

// ¿àÞØ×ÒÞÛìÝÐï ÓÛÐÔÚÐï ßÞÒÕàåÝÞáâì - áßÛÐÙÝclass any_spline{
protected:
  any_spline(){;};
};

typedef any_spline* p_any_spline;
typedef ObjectsAr<p_any_spline> any_splines;

// ¿àÞØ×ÒÞÛìÝÐï ÓàÐÝØæÐ àÞÒÝÞ ÜÕÖÔã ÔÒãÜï âÕÛÐÜØ
  
class any_border{
protected:
  border_flag bf;
  int iBodies[2];
  p_any_spline psSpline;
  any_border(any_border &bCopy);
  any_border(int ib[2],const p_any_spline ps,border_flag bf);
  any_border(const p_any_spline ps,border_flag bf);
public:
  int get_body(int ib) const;
  p_any_spline  spline()const{return psSpline;} ;
  void set_bodies(int i1,int i2){iBodies[0]=i1;iBodies[1]=i2;};
  border_flag kind() {return bf;};
  int another_body(int iB){
    if(iBodies[0]==iB)return iBodies[1];else
      if (iBodies[1]==iB)return iBodies[0];
    return -1;
  };
  ~any_border();
};

typedef any_border* p_any_border;
typedef ObjectsAr<p_any_border> any_borders;


// ÂÕÛÞ ÚÐÚ ÝÐÑÞà áÞÑáâÒÕÝÝëå ÓàÐÝØæ
class any_body{
protected:
  int iNumber;
  any_borders aBorders;
  IntAr aBordersNums;
  geophys_params gp;
  raytr_params rp;
  any_body(geophys_params _gp);
  any_body(any_body& bCopy);
  any_body(any_borders &ab,geophys_params gp);  
  void swap_borders(int,int);
  void gp2rp();
public:
  Logic numerated_borders()const{
    if(aBordersNums.count()>0)
      return TRUE;
    return FALSE;
  };
  void clear_borders();
  void set_number(int i){iNumber=i;};
  int n_of_borders() const {return aBorders.count();};
  void add_border(p_any_border pb){aBorders.addh(pb);};
  ~any_body();
  inline p_any_border get_border(int i)const ;
  geophys_params get_geophys() const {return gp;};
  void set_geophys(geophys_params GP) {gp=GP; gp2rp();};
  raytr_params get_raytr()const{return rp;};
  Logic contains_border(p_any_border pb)const;

  // ¾çÕÝÝÞ ÞßÐáÝÐï äãÝÚæØï, ÞÔÝÐÚÞ...
  void fix_numbers(const IntAr&aBorderNumbers);
  const IntAr& borders_numbers()const{return aBordersNums;};
  // ²ëÔÐÕâ áßØáÞÚ ÝÞÜÕàÞÒ âÕÛ-áÞáÕÔÕÙ (ÚàÞÜÕ -1, ÕááÝÞ)
  IntAr neighbors(int iNumber) const;

  // ·ÐÜÕÝØÛØ ÞÔÝã ÓàÐÝØæã ÝÐ ÔàãÓãî, ÕáÛØ çâÞ - ÚØÝãÛØ ØáÚÛîçÕÝØÕ
  void change_border(p_any_border,p_any_border); 

  // ¾âáÞàâØàÞÒÐâì ÓàÐÝØæë
  void sort_borders();
  int number()const{return iNumber;};

  void set_vp(double vp){gp.fVP=vp;gp2rp();};
  void set_kp(double kp){gp.fAn_KP=kp;gp2rp();};
  void set_vs(double vs){gp.fVS=vs;gp2rp();};
  void set_ks(double ks){gp.fAn_KS=ks;gp2rp();};  

  const any_borders& borders() const {return aBorders;};
};

// ¾çÕÝì ãÖ çÐáâÞ ÞÝÐ Òë×ëÒÐÕâáï...
inline p_any_border any_body::get_border (int i)const{
  // ¿àÞÒÕàÚØ ÝÐäØÓ, ßãáâì ObjectAr ßàÞÒÕàïÕâ...
#if 0
  if((i>=aBorders.count())||(i<0))
    return NULL;
#endif
  return aBorders(i);
};


typedef any_body* p_any_body;
typedef ObjectsAr<p_any_body> any_bodies;

// ÍÛÕÜÕÝâë ÒÝÕèÝÕÓÞ ßàÕÔáâÐÒÛÕÝØï ÜÞÔÕÛØ

// ÁßÛÐÙÝ
class any_ex_spline{
protected:
  any_ex_spline(){;};
};

// ³àÐÝØæÐ
class any_ex_border{
protected:
  border_flag bf;
  int iSpline;
  int iBodies[2];
  any_ex_border(int ib[2], int is, border_flag _bf);
  any_ex_border(){iSpline=-1; bf=NORMAL; iBodies[0]=iBodies[1]=-5;};
public:
  int get_body(int ib) ;
  int spline()const {return iSpline;};
  border_flag kind() {return bf;};
};

// ÂÕÛÞ
class any_ex_body{
protected:
  int iNumber;
  IntAr aBorders;
  geophys_params gp;
  raytr_params rp;
  any_ex_body(IntAr & ab, geophys_params const & gp, int in);
  any_ex_body(){;};
public:
  
  virtual int n_of_borders()const {return aBorders.count();};
  int make_body(const IntAr& ab, const geophys_params & gp,
        const raytr_params & rp,int in);
  geophys_params get_geophys()const{return gp;};
  raytr_params get_raytr()const{return rp;};
  int get_border(int i);
  int number()const{return iNumber;};
};

#endif /* any_elements.H */

/*  Melnikov George Y.  29-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 23-Nov-1999  */ 
/*          ...          */ 

class any_spline{
protected:
  any_spline(){;};
};

typedef any_spline* p_any_spline;
typedef ObjectsAr<p_any_spline> any_splines;

// ¿àÞØ×ÒÞÛìÝÐï ÓàÐÝØæÐ àÞÒÝÞ ÜÕÖÔã ÔÒãÜï âÕÛÐÜØ  
class any_border{
protected:
  border_flag bf;
  int iBodies[2];
  p_any_spline psSpline;
  any_border(any_border &bCopy);
  any_border(int ib[2],const p_any_spline ps,border_flag bf);
  any_border(const p_any_spline ps,border_flag bf);
public:
  int get_body(int ib) const;
  p_any_spline  spline()const{return psSpline;} ;
  void set_bodies(int i1,int i2){iBodies[0]=i1;iBodies[1]=i2;};
  border_flag kind() {return bf;};
  int another_body(int iB){
    if(iBodies[0]==iB)return iBodies[1];else
      if (iBodies[1]==iB)return iBodies[0];
    return -1;
  };
  ~any_border();
};

typedef any_border* p_any_border;
typedef ObjectsAr<p_any_border> any_borders;


// ÂÕÛÞ ÚÐÚ ÝÐÑÞà áÞÑáâÒÕÝÝëå ÓàÐÝØæ
class any_body{
protected:
  int iNumber;
  any_borders aBorders;
  IntAr aBordersNums;
  geophys_params gp;
  raytr_params rp;
  any_body(geophys_params _gp);
  any_body(any_body& bCopy);
  any_body(any_borders &ab,geophys_params gp);  
  void swap_borders(int,int);
  void gp2rp();
public:
  Logic numerated_borders()const{
    if(aBordersNums.count()>0)
      return TRUE;
    return FALSE;
  };
  void clear_borders();
  void set_number(int i){iNumber=i;};
  int n_of_borders() const {return aBorders.count();};
  void add_border(p_any_border pb){aBorders.addh(pb);};
  ~any_body();
  inline p_any_border get_border(int i)const ;
  geophys_params get_geophys() const {return gp;};
  void set_geophys(geophys_params GP) {gp=GP; gp2rp();};
  raytr_params get_raytr()const{return rp;};
  Logic contains_border(p_any_border pb)const;

  // ¾çÕÝÝÞ ÞßÐáÝÐï äãÝÚæØï, ÞÔÝÐÚÞ...
  void fix_numbers(const IntAr&aBorderNumbers);
  const IntAr& borders_numbers()const{return aBordersNums;};
  // ²ëÔÐÕâ áßØáÞÚ ÝÞÜÕàÞÒ âÕÛ-áÞáÕÔÕÙ (ÚàÞÜÕ -1, ÕááÝÞ)
  IntAr neighbors(int iNumber) const;

  // ·ÐÜÕÝØÛØ ÞÔÝã ÓàÐÝØæã ÝÐ ÔàãÓãî, ÕáÛØ çâÞ - ÚØÝãÛØ ØáÚÛîçÕÝØÕ
  void change_border(p_any_border,p_any_border); 

  // ¾âáÞàâØàÞÒÐâì ÓàÐÝØæë
  void sort_borders();
  int number()const{return iNumber;};

  void set_vp(double vp){gp.fVP=vp;gp2rp();};
  void set_kp(double kp){gp.fAn_KP=kp;gp2rp();};
  void set_vs(double vs){gp.fVS=vs;gp2rp();};
  void set_ks(double ks){gp.fAn_KS=ks;gp2rp();};  

  const any_borders& borders() const {return aBorders;};
};

// ¾çÕÝì ãÖ çÐáâÞ ÞÝÐ Òë×ëÒÐÕâáï...
inline p_any_border any_body::get_border (int i)const{
  // ¿àÞÒÕàÚØ ÝÐäØÓ, ßãáâì ObjectAr ßàÞÒÕàïÕâ...
#if 0
  if((i>=aBorders.count())||(i<0))
    return NULL;
#endif
  return aBorders(i);
};


typedef any_body* p_any_body;
typedef ObjectsAr<p_any_body> any_bodies;

// ÍÛÕÜÕÝâë ÒÝÕèÝÕÓÞ ßàÕÔáâÐÒÛÕÝØï ÜÞÔÕÛØ

// ÁßÛÐÙÝ
class any_ex_spline{
protected:
  any_ex_spline(){;};
};

// ³àÐÝØæÐ
class any_ex_border{
protected:
  border_flag bf;
  int iSpline;
  int iBodies[2];
  any_ex_border(int ib[2], int is, border_flag _bf);
  any_ex_border(){iSpline=-1; bf=NORMAL; iBodies[0]=iBodies[1]=-5;};
public:
  int get_body(int ib) ;
  int spline()const {return iSpline;};
  border_flag kind() {return bf;};
};

// ÂÕÛÞ
class any_ex_body{
protected:
  int iNumber;
  IntAr aBorders;
  geophys_params gp;
  raytr_params rp;
  any_ex_body(IntAr & ab, geophys_params const & gp, int in);
  any_ex_body(){;};
public:
  
  virtual int n_of_borders()const {return aBorders.count();};
  int make_body(const IntAr& ab, const geophys_params & gp,
        const raytr_params & rp,int in);
  geophys_params get_geophys()const{return gp;};
  raytr_params get_raytr()const{return rp;};
  int get_border(int i);
  int number()const{return iNumber;};
};

#endif /* any_elements.H */

/*  Melnikov George Y.  29-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 23-Nov-1999  */ 
/*          ...          */ 

  
class any_border{
protected:
  border_flag bf;
  int iBodies[2];
  p_any_spline psSpline;
  any_border(any_border &bCopy);
  any_border(int ib[2],const p_any_spline ps,border_flag bf);
  any_border(const p_any_spline ps,border_flag bf);
public:
  int get_body(int ib) const;
  p_any_spline  spline()const{return psSpline;} ;
  void set_bodies(int i1,int i2){iBodies[0]=i1;iBodies[1]=i2;};
  border_flag kind() {return bf;};
  int another_body(int iB){
    if(iBodies[0]==iB)return iBodies[1];else
      if (iBodies[1]==iB)return iBodies[0];
    return -1;
  };
  ~any_border();
};

typedef any_border* p_any_border;
typedef ObjectsAr<p_any_border> any_borders;


// ÂÕÛÞ ÚÐÚ ÝÐÑÞà áÞÑáâÒÕÝÝëå ÓàÐÝØæ
class any_body{
protected:
  int iNumber;
  any_borders aBorders;
  IntAr aBordersNums;
  geophys_params gp;
  raytr_params rp;
  any_body(geophys_params _gp);
  any_body(any_body& bCopy);
  any_body(any_borders &ab,geophys_params gp);  
  void swap_borders(int,int);
  void gp2rp();
public:
  Logic numerated_borders()const{
    if(aBordersNums.count()>0)
      return TRUE;
    return FALSE;
  };
  void clear_borders();
  void set_number(int i){iNumber=i;};
  int n_of_borders() const {return aBorders.count();};
  void add_border(p_any_border pb){aBorders.addh(pb);};
  ~any_body();
  inline p_any_border get_border(int i)const ;
  geophys_params get_geophys() const {return gp;};
  void set_geophys(geophys_params GP) {gp=GP; gp2rp();};
  raytr_params get_raytr()const{return rp;};
  Logic contains_border(p_any_border pb)const;

  // ¾çÕÝÝÞ ÞßÐáÝÐï äãÝÚæØï, ÞÔÝÐÚÞ...
  void fix_numbers(const IntAr&aBorderNumbers);
  const IntAr& borders_numbers()const{return aBordersNums;};
  // ²ëÔÐÕâ áßØáÞÚ ÝÞÜÕàÞÒ âÕÛ-áÞáÕÔÕÙ (ÚàÞÜÕ -1, ÕááÝÞ)
  IntAr neighbors(int iNumber) const;

  // ·ÐÜÕÝØÛØ ÞÔÝã ÓàÐÝØæã ÝÐ ÔàãÓãî, ÕáÛØ çâÞ - ÚØÝãÛØ ØáÚÛîçÕÝØÕ  void change_border(p_any_border,p_any_border); 

  // ¾âáÞàâØàÞÒÐâì ÓàÐÝØæë
  void sort_borders();
  int number()const{return iNumber;};

  void set_vp(double vp){gp.fVP=vp;gp2rp();};
  void set_kp(double kp){gp.fAn_KP=kp;gp2rp();};
  void set_vs(double vs){gp.fVS=vs;gp2rp();};
  void set_ks(double ks){gp.fAn_KS=ks;gp2rp();};  

  const any_borders& borders() const {return aBorders;};
};

// ¾çÕÝì ãÖ çÐáâÞ ÞÝÐ Òë×ëÒÐÕâáï...
inline p_any_border any_body::get_border (int i)const{
  // ¿àÞÒÕàÚØ ÝÐäØÓ, ßãáâì ObjectAr ßàÞÒÕàïÕâ...
#if 0
  if((i>=aBorders.count())||(i<0))
    return NULL;
#endif
  return aBorders(i);
};


typedef any_body* p_any_body;
typedef ObjectsAr<p_any_body> any_bodies;

// ÍÛÕÜÕÝâë ÒÝÕèÝÕÓÞ ßàÕÔáâÐÒÛÕÝØï ÜÞÔÕÛØ

// ÁßÛÐÙÝ
class any_ex_spline{
protected:
  any_ex_spline(){;};
};

// ³àÐÝØæÐ
class any_ex_border{
protected:
  border_flag bf;
  int iSpline;
  int iBodies[2];
  any_ex_border(int ib[2], int is, border_flag _bf);
  any_ex_border(){iSpline=-1; bf=NORMAL; iBodies[0]=iBodies[1]=-5;};
public:
  int get_body(int ib) ;
  int spline()const {return iSpline;};
  border_flag kind() {return bf;};
};

// ÂÕÛÞ
class any_ex_body{
protected:
  int iNumber;
  IntAr aBorders;
  geophys_params gp;
  raytr_params rp;
  any_ex_body(IntAr & ab, geophys_params const & gp, int in);
  any_ex_body(){;};
public:
  
  virtual int n_of_borders()const {return aBorders.count();};
  int make_body(const IntAr& ab, const geophys_params & gp,
        const raytr_params & rp,int in);
  geophys_params get_geophys()const{return gp;};
  raytr_params get_raytr()const{return rp;};
  int get_border(int i);
  int number()const{return iNumber;};
};

#endif /* any_elements.H */

/*  Melnikov George Y.  29-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 23-Nov-1999  */ 
/*          ...          */ 

  void change_border(p_any_border,p_any_border); 

  // ¾âáÞàâØàÞÒÐâì ÓàÐÝØæë  void sort_borders();
  int number()const{return iNumber;};

  void set_vp(double vp){gp.fVP=vp;gp2rp();};
  void set_kp(double kp){gp.fAn_KP=kp;gp2rp();};
  void set_vs(double vs){gp.fVS=vs;gp2rp();};
  void set_ks(double ks){gp.fAn_KS=ks;gp2rp();};  

  const any_borders& borders() const {return aBorders;};
};

// ¾çÕÝì ãÖ çÐáâÞ ÞÝÐ Òë×ëÒÐÕâáï...
inline p_any_border any_body::get_border (int i)const{
  // ¿àÞÒÕàÚØ ÝÐäØÓ, ßãáâì ObjectAr ßàÞÒÕàïÕâ...
#if 0
  if((i>=aBorders.count())||(i<0))
    return NULL;
#endif
  return aBorders(i);
};


typedef any_body* p_any_body;
typedef ObjectsAr<p_any_body> any_bodies;

// ÍÛÕÜÕÝâë ÒÝÕèÝÕÓÞ ßàÕÔáâÐÒÛÕÝØï ÜÞÔÕÛØ

// ÁßÛÐÙÝ
class any_ex_spline{
protected:
  any_ex_spline(){;};
};

// ³àÐÝØæÐ
class any_ex_border{
protected:
  border_flag bf;
  int iSpline;
  int iBodies[2];
  any_ex_border(int ib[2], int is, border_flag _bf);
  any_ex_border(){iSpline=-1; bf=NORMAL; iBodies[0]=iBodies[1]=-5;};
public:
  int get_body(int ib) ;
  int spline()const {return iSpline;};
  border_flag kind() {return bf;};
};

// ÂÕÛÞ
class any_ex_body{
protected:
  int iNumber;
  IntAr aBorders;
  geophys_params gp;
  raytr_params rp;
  any_ex_body(IntAr & ab, geophys_params const & gp, int in);
  any_ex_body(){;};
public:
  
  virtual int n_of_borders()const {return aBorders.count();};
  int make_body(const IntAr& ab, const geophys_params & gp,
        const raytr_params & rp,int in);
  geophys_params get_geophys()const{return gp;};
  raytr_params get_raytr()const{return rp;};
  int get_border(int i);
  int number()const{return iNumber;};
};

#endif /* any_elements.H */

/*  Melnikov George Y.  29-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 23-Nov-1999  */ 
/*          ...          */ 

  void sort_borders();
  int number()const{return iNumber;};

  void set_vp(double vp){gp.fVP=vp;gp2rp();};
  void set_kp(double kp){gp.fAn_KP=kp;gp2rp();};
  void set_vs(double vs){gp.fVS=vs;gp2rp();};
  void set_ks(double ks){gp.fAn_KS=ks;gp2rp();};  

  const any_borders& borders() const {return aBorders;};
};

// ¾çÕÝì ãÖ çÐáâÞ ÞÝÐ Òë×ëÒÐÕâáï...
inline p_any_border any_body::get_border (int i)const{
  // ¿àÞÒÕàÚØ ÝÐäØÓ, ßãáâì ObjectAr ßàÞÒÕàïÕâ...
#if 0
  if((i>=aBorders.count())||(i<0))
    return NULL;
#endif
  return aBorders(i);
};


typedef any_body* p_any_body;
typedef ObjectsAr<p_any_body> any_bodies;

// ÍÛÕÜÕÝâë ÒÝÕèÝÕÓÞ ßàÕÔáâÐÒÛÕÝØï ÜÞÔÕÛØ

// ÁßÛÐÙÝ
class any_ex_spline{
protected:
  any_ex_spline(){;};
};

// ³àÐÝØæÐclass any_ex_border{
protected:
  border_flag bf;
  int iSpline;
  int iBodies[2];
  any_ex_border(int ib[2], int is, border_flag _bf);
  any_ex_border(){iSpline=-1; bf=NORMAL; iBodies[0]=iBodies[1]=-5;};
public:
  int get_body(int ib) ;
  int spline()const {return iSpline;};
  border_flag kind() {return bf;};
};

// ÂÕÛÞ
class any_ex_body{
protected:
  int iNumber;
  IntAr aBorders;
  geophys_params gp;
  raytr_params rp;
  any_ex_body(IntAr & ab, geophys_params const & gp, int in);
  any_ex_body(){;};
public:
  
  virtual int n_of_borders()const {return aBorders.count();};
  int make_body(const IntAr& ab, const geophys_params & gp,
        const raytr_params & rp,int in);
  geophys_params get_geophys()const{return gp;};
  raytr_params get_raytr()const{return rp;};
  int get_border(int i);
  int number()const{return iNumber;};
};

#endif /* any_elements.H */

/*  Melnikov George Y.  29-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 23-Nov-1999  */ 
/*          ...          */ 

class any_ex_border{
protected:
  border_flag bf;
  int iSpline;
  int iBodies[2];
  any_ex_border(int ib[2], int is, border_flag _bf);
  any_ex_border(){iSpline=-1; bf=NORMAL; iBodies[0]=iBodies[1]=-5;};
public:
  int get_body(int ib) ;
  int spline()const {return iSpline;};
  border_flag kind() {return bf;};
};

// ÂÕÛÞ
class any_ex_body{
protected:
  int iNumber;
  IntAr aBorders;
  geophys_params gp;
  raytr_params rp;
  any_ex_body(IntAr & ab, geophys_params const & gp, int in);
  any_ex_body(){;};
public:
  
  virtual int n_of_borders()const {return aBorders.count();};
  int make_body(const IntAr& ab, const geophys_params & gp,
        const raytr_params & rp,int in);
  geophys_params get_geophys()const{return gp;};
  raytr_params get_raytr()const{return rp;};
  int get_border(int i);
  int number()const{return iNumber;};
};

#endif /* any_elements.H */

/*  Melnikov George Y.  29-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 23-Nov-1999  */ 
/*          ...          */