p_elements.H

/* p_elements.H */
/* $Id: p_elements.H,v 2.15 2001/10/09 14:54:38 hoh Exp $ */
#ifndef __p_elements_H
#define __p_elements_H


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

#include "any_elements.H"
/*
   ÍÛÕÜÕÝâë ÜÞÔÕÛØ á ßÛÞáÚØÜØ ÓàÐÝØæÐÜØ: ÓàÐÝØæë, âÕÛÐ, áßÛÐÙÝë.

   ¿àÕÔÞáâÐÒÛïîâ ÝÕÞÑåÞÔØÜãî äãÝÚæØÞÝÐÛìÝÞáâì Ø ßàØÞâÚàëÒÐîâáï ÔÛï ßÞÛì×ÞÒÐâÕÛï. */

class pl_model;
class any_model;
class pl_spline;

// ÍÛÕÜÕÝâë ÒÝãâàÕÝÝÕÓÞ ßàÕÔáâÐÒÛÕÝØï
class pl_intern_spline:public any_spline{

// ³ÛÐÔÚÐï ßÞÒÕàåÝÞáâì - ßÛÞáÚÞáâì
protected:
  double a,b,c,d;  // ºÞíääØæØÕÝâë ÚÐÝÞÝØçÕáÚÞÓÞ ãàÐÒÝÕÝØï ßÛÞáÚÞáâØ
  pl_intern_spline(); // ¸áßÞÛì×ãÕÜ ßàØ ßÞàÞÖÔÕÝØØ ÔÕâÕÙ

public:
  ~pl_intern_spline(){;};
  pl_intern_spline(double A,double B,double C, double D);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  // (ÑÞÛÕÕ ÝÞÒëÙ, ÝÞ áâÐàëÙ ßÞÚÐ âàÞÓÐâì ÝÕ ÑãÔÕÜ)
  point2D local(const point& pMain,
        const point& p0) const;

  // ¿ÕàÕÒÞÔ Ò ÓÛÞÑÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point global(const point2D& pMain,
           const point& p0) const;

  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};

  double distance(const point3d &pWhere);

  // ´àã×ìï 
friend class pl_model;
};

typedef pl_intern_spline* p_pl_intern_spline;
typedef ObjectsAr<p_pl_intern_spline> pl_intern_splines;


// ³àÐÝØæÐ á áÞÞâÒÕâáâÒãîéØÜØ ÓÕÞäØ×ØçÕáÚØÜØ ßÐàÐÜÕâàÐÜØ
class pl_overlying_spline:public pl_intern_spline{
protected:
  geophys_params gp;
  point3d pInternal; // ÂÞçÚÐ, ÛÕÖÐéÐï ÝÐ ÓàÐÝØæÕ, ÞÑàÐ×ÞÒÐÝÝÞÙ íâØÜ áßÛÐÙÝÞÜ
public:
  pl_overlying_spline(const geophys_params &gp,const point3d& p,
              double dFi,double dAz);
  geophys_params get_geophys(){return gp;};
  point3d get_internal_point(){return pInternal;};
};


typedef pl_overlying_spline* p_pl_overlying_spline;
typedef ObjectsAr<p_pl_overlying_spline> pl_overlying_splines;

// ³àÐÝØæÐ - ÜÝÞÓÞãÓÞÛìÝØÚ
// ÂàÕÑãÕâáï ÝÕÚÞâÞàëÙ áßÞáÞÑ ÞßàÕÔÕÛÕÝØï æÕÛÞáâÝÞáâØ ÓàÐÝØæë...
class pl_intern_border:public any_border{
protected:
  aPoints aTops; // ²ÕàèØÝë ÜÝÞÓÞãÓÞÛìÝØÚÐ, âàÕÑÞÒÐÝØÕ: aTops[0]=aTops[_last_].
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
  pl_intern_border(pl_intern_border &bCopy);

  // ÂÞçÚØ, ÞßàÕÔÕÛïîéØÕ ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ ÝÐ ßÐßÕ-áßÛÐÙÝÕ
  Logic points_created;
  point p0;

  void init_local_coords(); // ²ëçØáÛïÕâ âÞçÚã p0
public:
  ~pl_intern_border();
  pl_intern_border(const aPoints & p, int ib[2],p_pl_intern_spline ps,
           border_flag bf);
  pl_intern_border(int ib[2],p_pl_intern_spline psSpline,border_flag bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag _bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag,aSegments3D&);

  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;
  point get_a_top(int i)const{return aTops(i);};
  // ¾ßÐáÝÐï äãÝÚæØï, ÞÔÝÐÚÞ...
  int add_a_top(const point& pt);
  int n_of_tops()const{return aTops.count()-1;};
  virtual int contains(const point &pWhat)const;
  // ¿Þ ÞâàÕ×ÚÐÜ - ßÞáÛÕÔÞÒÐâÕÛìÝÞáâì âÞçÕÚ.
  int fix_points(aSegments3D&);
  const aPoints& points()const{return aTops;};
  point closest_border_point(const point&)const;

  // ¿ÕàÕÒÞÔ Ø× ÞÔÝØå ÚÞÞàÔØÝÐâ Ò ÔàãÓØÕ
  point2D local(const point3d&);
  point3d global(const point2D&);
  point P0(){return p0;};
  point middle()const;
  friend class pl_model;
};

typedef pl_intern_border* p_pl_intern_border;
#define AnyType p_pl_intern_border
#define AnyAr pl_intern_borders
#include <mix/AnyAr.H>


class pl_fict_border:public pl_intern_border,public pl_intern_spline{
public:
  pl_fict_border(const aPoints & p, int ib[2],
         double A,double B, double C, double D):
   pl_intern_border(p,ib,NULL,TEMP),pl_intern_spline(A,B,C,D){};
};

typedef pl_fict_border* p_pl_fict_border;
typedef ObjectsAr<p_pl_fict_border> pl_fict_borders;

// ÂÕÛÞ - ÝÐÑÞà ÜÝÞÓÞãÓÞÛìÝØÚÞÒ
class pl_intern_body:public any_body{
protected:
  pl_intern_body(pl_intern_body &bCopy);
  pl_intern_splines* get_bodies_splines()const;
public:
  pl_intern_body(geophys_params gp);
  pl_intern_body(any_borders &ab,geophys_params gp);
  ~pl_intern_body(){;};

  // ´ÞÑÐÒÛïÕÜ Ò íâÞâ ÚÛÐáá, âÐÚ ÚÐÚ âàÕÑãÕâáï ßÕàÕáÕçÕÝØÕ á ßÛÞáÚÞáâïÜØ
  virtual int contains(const point &pWhat)const;

  // ÁÞ×ÔÐÝØÕ äØÚâØÒÝÞÙ ÓàÐÝØæë ÒÝãâàØ âÕÛÐ ßÞ ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâØ
  int fict_border (const pl_spline& sSlice,p_pl_fict_border& bNew)const;

friend class pl_model;  
};

typedef pl_intern_body* p_pl_intern_body;
typedef ObjectsAr<p_pl_intern_body> pl_intern_bodies;

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

class pl_spline:public any_ex_spline{
protected:
  double a,b,c,d;

public:
  pl_spline();
  pl_spline(double A,double B,double C, double D);

  // ßÛÞáÚÞáâì ßÞ âàÕÜ âÞçÚÐÜ
  pl_spline(const point& p1,const point& p2,const point& p3);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  void set_coef(double A,double B,double C, double D);
  void get_coef(double &A,double &B,double &C, double &D);
  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};
  double distance(const point3d &pWhere);
};

typedef pl_spline* p_pl_spline;
typedef ObjectsAr<p_pl_spline> p_pl_splines;

class pl_border:public any_ex_border{
protected:
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
public:
  pl_border();
  pl_border(const aPoints & at, int ib[2], int is, border_flag bf);
  pl_border(int ib[2], int is, border_flag bf);
  point get_a_top(int i)const {return point(aTops_X(i),aTops_Y(i),aTops_Z(i));};
  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;

  // °ÝÐÛÞÓ ÚÞÝáâàãÚâÞàÐ ÔÛï ßÕàÕ×ÐßØáØ áÞ×ÔÐÝÝÞÓÞ ÞÑêÕÚâÐ
  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞ
class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

 */

class pl_model;
class any_model;
class pl_spline;

// ÍÛÕÜÕÝâë ÒÝãâàÕÝÝÕÓÞ ßàÕÔáâÐÒÛÕÝØï
class pl_intern_spline:public any_spline{

// ³ÛÐÔÚÐï ßÞÒÕàåÝÞáâì - ßÛÞáÚÞáâì
protected:
  double a,b,c,d;  // ºÞíääØæØÕÝâë ÚÐÝÞÝØçÕáÚÞÓÞ ãàÐÒÝÕÝØï ßÛÞáÚÞáâØ
  pl_intern_spline(); // ¸áßÞÛì×ãÕÜ ßàØ ßÞàÞÖÔÕÝØØ ÔÕâÕÙ

public:
  ~pl_intern_spline(){;};
  pl_intern_spline(double A,double B,double C, double D);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  // (ÑÞÛÕÕ ÝÞÒëÙ, ÝÞ áâÐàëÙ ßÞÚÐ âàÞÓÐâì ÝÕ ÑãÔÕÜ)
  point2D local(const point& pMain,
        const point& p0) const;

  // ¿ÕàÕÒÞÔ Ò ÓÛÞÑÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point global(const point2D& pMain,
           const point& p0) const;

  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};

  double distance(const point3d &pWhere);

  // ´àã×ìï 
friend class pl_model;
};

typedef pl_intern_spline* p_pl_intern_spline;
typedef ObjectsAr<p_pl_intern_spline> pl_intern_splines;


// ³àÐÝØæÐ á áÞÞâÒÕâáâÒãîéØÜØ ÓÕÞäØ×ØçÕáÚØÜØ ßÐàÐÜÕâàÐÜØ
class pl_overlying_spline:public pl_intern_spline{
protected:
  geophys_params gp;
  point3d pInternal; // ÂÞçÚÐ, ÛÕÖÐéÐï ÝÐ ÓàÐÝØæÕ, ÞÑàÐ×ÞÒÐÝÝÞÙ íâØÜ áßÛÐÙÝÞÜ
public:
  pl_overlying_spline(const geophys_params &gp,const point3d& p,
              double dFi,double dAz);
  geophys_params get_geophys(){return gp;};
  point3d get_internal_point(){return pInternal;};
};


typedef pl_overlying_spline* p_pl_overlying_spline;
typedef ObjectsAr<p_pl_overlying_spline> pl_overlying_splines;

// ³àÐÝØæÐ - ÜÝÞÓÞãÓÞÛìÝØÚ
// ÂàÕÑãÕâáï ÝÕÚÞâÞàëÙ áßÞáÞÑ ÞßàÕÔÕÛÕÝØï æÕÛÞáâÝÞáâØ ÓàÐÝØæë...
class pl_intern_border:public any_border{
protected:
  aPoints aTops; // ²ÕàèØÝë ÜÝÞÓÞãÓÞÛìÝØÚÐ, âàÕÑÞÒÐÝØÕ: aTops[0]=aTops[_last_].
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
  pl_intern_border(pl_intern_border &bCopy);

  // ÂÞçÚØ, ÞßàÕÔÕÛïîéØÕ ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ ÝÐ ßÐßÕ-áßÛÐÙÝÕ
  Logic points_created;
  point p0;

  void init_local_coords(); // ²ëçØáÛïÕâ âÞçÚã p0
public:
  ~pl_intern_border();
  pl_intern_border(const aPoints & p, int ib[2],p_pl_intern_spline ps,
           border_flag bf);
  pl_intern_border(int ib[2],p_pl_intern_spline psSpline,border_flag bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag _bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag,aSegments3D&);

  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;
  point get_a_top(int i)const{return aTops(i);};
  // ¾ßÐáÝÐï äãÝÚæØï, ÞÔÝÐÚÞ...
  int add_a_top(const point& pt);
  int n_of_tops()const{return aTops.count()-1;};
  virtual int contains(const point &pWhat)const;
  // ¿Þ ÞâàÕ×ÚÐÜ - ßÞáÛÕÔÞÒÐâÕÛìÝÞáâì âÞçÕÚ.
  int fix_points(aSegments3D&);
  const aPoints& points()const{return aTops;};
  point closest_border_point(const point&)const;

  // ¿ÕàÕÒÞÔ Ø× ÞÔÝØå ÚÞÞàÔØÝÐâ Ò ÔàãÓØÕ
  point2D local(const point3d&);
  point3d global(const point2D&);
  point P0(){return p0;};
  point middle()const;
  friend class pl_model;
};

typedef pl_intern_border* p_pl_intern_border;
#define AnyType p_pl_intern_border
#define AnyAr pl_intern_borders
#include <mix/AnyAr.H>


class pl_fict_border:public pl_intern_border,public pl_intern_spline{
public:
  pl_fict_border(const aPoints & p, int ib[2],
         double A,double B, double C, double D):
   pl_intern_border(p,ib,NULL,TEMP),pl_intern_spline(A,B,C,D){};
};

typedef pl_fict_border* p_pl_fict_border;
typedef ObjectsAr<p_pl_fict_border> pl_fict_borders;

// ÂÕÛÞ - ÝÐÑÞà ÜÝÞÓÞãÓÞÛìÝØÚÞÒ
class pl_intern_body:public any_body{
protected:
  pl_intern_body(pl_intern_body &bCopy);
  pl_intern_splines* get_bodies_splines()const;
public:
  pl_intern_body(geophys_params gp);
  pl_intern_body(any_borders &ab,geophys_params gp);
  ~pl_intern_body(){;};

  // ´ÞÑÐÒÛïÕÜ Ò íâÞâ ÚÛÐáá, âÐÚ ÚÐÚ âàÕÑãÕâáï ßÕàÕáÕçÕÝØÕ á ßÛÞáÚÞáâïÜØ
  virtual int contains(const point &pWhat)const;

  // ÁÞ×ÔÐÝØÕ äØÚâØÒÝÞÙ ÓàÐÝØæë ÒÝãâàØ âÕÛÐ ßÞ ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâØ
  int fict_border (const pl_spline& sSlice,p_pl_fict_border& bNew)const;

friend class pl_model;  
};

typedef pl_intern_body* p_pl_intern_body;
typedef ObjectsAr<p_pl_intern_body> pl_intern_bodies;

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

class pl_spline:public any_ex_spline{
protected:
  double a,b,c,d;

public:
  pl_spline();
  pl_spline(double A,double B,double C, double D);

  // ßÛÞáÚÞáâì ßÞ âàÕÜ âÞçÚÐÜ
  pl_spline(const point& p1,const point& p2,const point& p3);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  void set_coef(double A,double B,double C, double D);
  void get_coef(double &A,double &B,double &C, double &D);
  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};
  double distance(const point3d &pWhere);
};

typedef pl_spline* p_pl_spline;
typedef ObjectsAr<p_pl_spline> p_pl_splines;

class pl_border:public any_ex_border{
protected:
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
public:
  pl_border();
  pl_border(const aPoints & at, int ib[2], int is, border_flag bf);
  pl_border(int ib[2], int is, border_flag bf);
  point get_a_top(int i)const {return point(aTops_X(i),aTops_Y(i),aTops_Z(i));};
  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;

  // °ÝÐÛÞÓ ÚÞÝáâàãÚâÞàÐ ÔÛï ßÕàÕ×ÐßØáØ áÞ×ÔÐÝÝÞÓÞ ÞÑêÕÚâÐ
  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞ
class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ  // (ÑÞÛÕÕ ÝÞÒëÙ, ÝÞ áâÐàëÙ ßÞÚÐ âàÞÓÐâì ÝÕ ÑãÔÕÜ)
  point2D local(const point& pMain,
        const point& p0) const;

  // ¿ÕàÕÒÞÔ Ò ÓÛÞÑÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point global(const point2D& pMain,
           const point& p0) const;

  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};

  double distance(const point3d &pWhere);

  // ´àã×ìï 
friend class pl_model;
};

typedef pl_intern_spline* p_pl_intern_spline;
typedef ObjectsAr<p_pl_intern_spline> pl_intern_splines;


// ³àÐÝØæÐ á áÞÞâÒÕâáâÒãîéØÜØ ÓÕÞäØ×ØçÕáÚØÜØ ßÐàÐÜÕâàÐÜØ
class pl_overlying_spline:public pl_intern_spline{
protected:
  geophys_params gp;
  point3d pInternal; // ÂÞçÚÐ, ÛÕÖÐéÐï ÝÐ ÓàÐÝØæÕ, ÞÑàÐ×ÞÒÐÝÝÞÙ íâØÜ áßÛÐÙÝÞÜ
public:
  pl_overlying_spline(const geophys_params &gp,const point3d& p,
              double dFi,double dAz);
  geophys_params get_geophys(){return gp;};
  point3d get_internal_point(){return pInternal;};
};


typedef pl_overlying_spline* p_pl_overlying_spline;
typedef ObjectsAr<p_pl_overlying_spline> pl_overlying_splines;

// ³àÐÝØæÐ - ÜÝÞÓÞãÓÞÛìÝØÚ
// ÂàÕÑãÕâáï ÝÕÚÞâÞàëÙ áßÞáÞÑ ÞßàÕÔÕÛÕÝØï æÕÛÞáâÝÞáâØ ÓàÐÝØæë...
class pl_intern_border:public any_border{
protected:
  aPoints aTops; // ²ÕàèØÝë ÜÝÞÓÞãÓÞÛìÝØÚÐ, âàÕÑÞÒÐÝØÕ: aTops[0]=aTops[_last_].
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
  pl_intern_border(pl_intern_border &bCopy);

  // ÂÞçÚØ, ÞßàÕÔÕÛïîéØÕ ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ ÝÐ ßÐßÕ-áßÛÐÙÝÕ
  Logic points_created;
  point p0;

  void init_local_coords(); // ²ëçØáÛïÕâ âÞçÚã p0
public:
  ~pl_intern_border();
  pl_intern_border(const aPoints & p, int ib[2],p_pl_intern_spline ps,
           border_flag bf);
  pl_intern_border(int ib[2],p_pl_intern_spline psSpline,border_flag bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag _bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag,aSegments3D&);

  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;
  point get_a_top(int i)const{return aTops(i);};
  // ¾ßÐáÝÐï äãÝÚæØï, ÞÔÝÐÚÞ...
  int add_a_top(const point& pt);
  int n_of_tops()const{return aTops.count()-1;};
  virtual int contains(const point &pWhat)const;
  // ¿Þ ÞâàÕ×ÚÐÜ - ßÞáÛÕÔÞÒÐâÕÛìÝÞáâì âÞçÕÚ.
  int fix_points(aSegments3D&);
  const aPoints& points()const{return aTops;};
  point closest_border_point(const point&)const;

  // ¿ÕàÕÒÞÔ Ø× ÞÔÝØå ÚÞÞàÔØÝÐâ Ò ÔàãÓØÕ
  point2D local(const point3d&);
  point3d global(const point2D&);
  point P0(){return p0;};
  point middle()const;
  friend class pl_model;
};

typedef pl_intern_border* p_pl_intern_border;
#define AnyType p_pl_intern_border
#define AnyAr pl_intern_borders
#include <mix/AnyAr.H>


class pl_fict_border:public pl_intern_border,public pl_intern_spline{
public:
  pl_fict_border(const aPoints & p, int ib[2],
         double A,double B, double C, double D):
   pl_intern_border(p,ib,NULL,TEMP),pl_intern_spline(A,B,C,D){};
};

typedef pl_fict_border* p_pl_fict_border;
typedef ObjectsAr<p_pl_fict_border> pl_fict_borders;

// ÂÕÛÞ - ÝÐÑÞà ÜÝÞÓÞãÓÞÛìÝØÚÞÒ
class pl_intern_body:public any_body{
protected:
  pl_intern_body(pl_intern_body &bCopy);
  pl_intern_splines* get_bodies_splines()const;
public:
  pl_intern_body(geophys_params gp);
  pl_intern_body(any_borders &ab,geophys_params gp);
  ~pl_intern_body(){;};

  // ´ÞÑÐÒÛïÕÜ Ò íâÞâ ÚÛÐáá, âÐÚ ÚÐÚ âàÕÑãÕâáï ßÕàÕáÕçÕÝØÕ á ßÛÞáÚÞáâïÜØ
  virtual int contains(const point &pWhat)const;

  // ÁÞ×ÔÐÝØÕ äØÚâØÒÝÞÙ ÓàÐÝØæë ÒÝãâàØ âÕÛÐ ßÞ ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâØ
  int fict_border (const pl_spline& sSlice,p_pl_fict_border& bNew)const;

friend class pl_model;  
};

typedef pl_intern_body* p_pl_intern_body;
typedef ObjectsAr<p_pl_intern_body> pl_intern_bodies;

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

class pl_spline:public any_ex_spline{
protected:
  double a,b,c,d;

public:
  pl_spline();
  pl_spline(double A,double B,double C, double D);

  // ßÛÞáÚÞáâì ßÞ âàÕÜ âÞçÚÐÜ
  pl_spline(const point& p1,const point& p2,const point& p3);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  void set_coef(double A,double B,double C, double D);
  void get_coef(double &A,double &B,double &C, double &D);
  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};
  double distance(const point3d &pWhere);
};

typedef pl_spline* p_pl_spline;
typedef ObjectsAr<p_pl_spline> p_pl_splines;

class pl_border:public any_ex_border{
protected:
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
public:
  pl_border();
  pl_border(const aPoints & at, int ib[2], int is, border_flag bf);
  pl_border(int ib[2], int is, border_flag bf);
  point get_a_top(int i)const {return point(aTops_X(i),aTops_Y(i),aTops_Z(i));};
  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;

  // °ÝÐÛÞÓ ÚÞÝáâàãÚâÞàÐ ÔÛï ßÕàÕ×ÐßØáØ áÞ×ÔÐÝÝÞÓÞ ÞÑêÕÚâÐ
  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞ
class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

  // (ÑÞÛÕÕ ÝÞÒëÙ, ÝÞ áâÐàëÙ ßÞÚÐ âàÞÓÐâì ÝÕ ÑãÔÕÜ)
  point2D local(const point& pMain,
        const point& p0) const;

  // ¿ÕàÕÒÞÔ Ò ÓÛÞÑÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ  point global(const point2D& pMain,
           const point& p0) const;

  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};

  double distance(const point3d &pWhere);

  // ´àã×ìï 
friend class pl_model;
};

typedef pl_intern_spline* p_pl_intern_spline;
typedef ObjectsAr<p_pl_intern_spline> pl_intern_splines;


// ³àÐÝØæÐ á áÞÞâÒÕâáâÒãîéØÜØ ÓÕÞäØ×ØçÕáÚØÜØ ßÐàÐÜÕâàÐÜØ
class pl_overlying_spline:public pl_intern_spline{
protected:
  geophys_params gp;
  point3d pInternal; // ÂÞçÚÐ, ÛÕÖÐéÐï ÝÐ ÓàÐÝØæÕ, ÞÑàÐ×ÞÒÐÝÝÞÙ íâØÜ áßÛÐÙÝÞÜ
public:
  pl_overlying_spline(const geophys_params &gp,const point3d& p,
              double dFi,double dAz);
  geophys_params get_geophys(){return gp;};
  point3d get_internal_point(){return pInternal;};
};


typedef pl_overlying_spline* p_pl_overlying_spline;
typedef ObjectsAr<p_pl_overlying_spline> pl_overlying_splines;

// ³àÐÝØæÐ - ÜÝÞÓÞãÓÞÛìÝØÚ
// ÂàÕÑãÕâáï ÝÕÚÞâÞàëÙ áßÞáÞÑ ÞßàÕÔÕÛÕÝØï æÕÛÞáâÝÞáâØ ÓàÐÝØæë...
class pl_intern_border:public any_border{
protected:
  aPoints aTops; // ²ÕàèØÝë ÜÝÞÓÞãÓÞÛìÝØÚÐ, âàÕÑÞÒÐÝØÕ: aTops[0]=aTops[_last_].
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
  pl_intern_border(pl_intern_border &bCopy);

  // ÂÞçÚØ, ÞßàÕÔÕÛïîéØÕ ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ ÝÐ ßÐßÕ-áßÛÐÙÝÕ
  Logic points_created;
  point p0;

  void init_local_coords(); // ²ëçØáÛïÕâ âÞçÚã p0
public:
  ~pl_intern_border();
  pl_intern_border(const aPoints & p, int ib[2],p_pl_intern_spline ps,
           border_flag bf);
  pl_intern_border(int ib[2],p_pl_intern_spline psSpline,border_flag bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag _bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag,aSegments3D&);

  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;
  point get_a_top(int i)const{return aTops(i);};
  // ¾ßÐáÝÐï äãÝÚæØï, ÞÔÝÐÚÞ...
  int add_a_top(const point& pt);
  int n_of_tops()const{return aTops.count()-1;};
  virtual int contains(const point &pWhat)const;
  // ¿Þ ÞâàÕ×ÚÐÜ - ßÞáÛÕÔÞÒÐâÕÛìÝÞáâì âÞçÕÚ.
  int fix_points(aSegments3D&);
  const aPoints& points()const{return aTops;};
  point closest_border_point(const point&)const;

  // ¿ÕàÕÒÞÔ Ø× ÞÔÝØå ÚÞÞàÔØÝÐâ Ò ÔàãÓØÕ
  point2D local(const point3d&);
  point3d global(const point2D&);
  point P0(){return p0;};
  point middle()const;
  friend class pl_model;
};

typedef pl_intern_border* p_pl_intern_border;
#define AnyType p_pl_intern_border
#define AnyAr pl_intern_borders
#include <mix/AnyAr.H>


class pl_fict_border:public pl_intern_border,public pl_intern_spline{
public:
  pl_fict_border(const aPoints & p, int ib[2],
         double A,double B, double C, double D):
   pl_intern_border(p,ib,NULL,TEMP),pl_intern_spline(A,B,C,D){};
};

typedef pl_fict_border* p_pl_fict_border;
typedef ObjectsAr<p_pl_fict_border> pl_fict_borders;

// ÂÕÛÞ - ÝÐÑÞà ÜÝÞÓÞãÓÞÛìÝØÚÞÒ
class pl_intern_body:public any_body{
protected:
  pl_intern_body(pl_intern_body &bCopy);
  pl_intern_splines* get_bodies_splines()const;
public:
  pl_intern_body(geophys_params gp);
  pl_intern_body(any_borders &ab,geophys_params gp);
  ~pl_intern_body(){;};

  // ´ÞÑÐÒÛïÕÜ Ò íâÞâ ÚÛÐáá, âÐÚ ÚÐÚ âàÕÑãÕâáï ßÕàÕáÕçÕÝØÕ á ßÛÞáÚÞáâïÜØ
  virtual int contains(const point &pWhat)const;

  // ÁÞ×ÔÐÝØÕ äØÚâØÒÝÞÙ ÓàÐÝØæë ÒÝãâàØ âÕÛÐ ßÞ ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâØ
  int fict_border (const pl_spline& sSlice,p_pl_fict_border& bNew)const;

friend class pl_model;  
};

typedef pl_intern_body* p_pl_intern_body;
typedef ObjectsAr<p_pl_intern_body> pl_intern_bodies;

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

class pl_spline:public any_ex_spline{
protected:
  double a,b,c,d;

public:
  pl_spline();
  pl_spline(double A,double B,double C, double D);

  // ßÛÞáÚÞáâì ßÞ âàÕÜ âÞçÚÐÜ
  pl_spline(const point& p1,const point& p2,const point& p3);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  void set_coef(double A,double B,double C, double D);
  void get_coef(double &A,double &B,double &C, double &D);
  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};
  double distance(const point3d &pWhere);
};

typedef pl_spline* p_pl_spline;
typedef ObjectsAr<p_pl_spline> p_pl_splines;

class pl_border:public any_ex_border{
protected:
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
public:
  pl_border();
  pl_border(const aPoints & at, int ib[2], int is, border_flag bf);
  pl_border(int ib[2], int is, border_flag bf);
  point get_a_top(int i)const {return point(aTops_X(i),aTops_Y(i),aTops_Z(i));};
  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;

  // °ÝÐÛÞÓ ÚÞÝáâàãÚâÞàÐ ÔÛï ßÕàÕ×ÐßØáØ áÞ×ÔÐÝÝÞÓÞ ÞÑêÕÚâÐ
  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞ
class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

  point global(const point2D& pMain,
           const point& p0) const;

  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};

  double distance(const point3d &pWhere);

  // ´àã×ìï friend class pl_model;
};

typedef pl_intern_spline* p_pl_intern_spline;
typedef ObjectsAr<p_pl_intern_spline> pl_intern_splines;


// ³àÐÝØæÐ á áÞÞâÒÕâáâÒãîéØÜØ ÓÕÞäØ×ØçÕáÚØÜØ ßÐàÐÜÕâàÐÜØ
class pl_overlying_spline:public pl_intern_spline{
protected:
  geophys_params gp;
  point3d pInternal; // ÂÞçÚÐ, ÛÕÖÐéÐï ÝÐ ÓàÐÝØæÕ, ÞÑàÐ×ÞÒÐÝÝÞÙ íâØÜ áßÛÐÙÝÞÜ
public:
  pl_overlying_spline(const geophys_params &gp,const point3d& p,
              double dFi,double dAz);
  geophys_params get_geophys(){return gp;};
  point3d get_internal_point(){return pInternal;};
};


typedef pl_overlying_spline* p_pl_overlying_spline;
typedef ObjectsAr<p_pl_overlying_spline> pl_overlying_splines;

// ³àÐÝØæÐ - ÜÝÞÓÞãÓÞÛìÝØÚ
// ÂàÕÑãÕâáï ÝÕÚÞâÞàëÙ áßÞáÞÑ ÞßàÕÔÕÛÕÝØï æÕÛÞáâÝÞáâØ ÓàÐÝØæë...
class pl_intern_border:public any_border{
protected:
  aPoints aTops; // ²ÕàèØÝë ÜÝÞÓÞãÓÞÛìÝØÚÐ, âàÕÑÞÒÐÝØÕ: aTops[0]=aTops[_last_].
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
  pl_intern_border(pl_intern_border &bCopy);

  // ÂÞçÚØ, ÞßàÕÔÕÛïîéØÕ ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ ÝÐ ßÐßÕ-áßÛÐÙÝÕ
  Logic points_created;
  point p0;

  void init_local_coords(); // ²ëçØáÛïÕâ âÞçÚã p0
public:
  ~pl_intern_border();
  pl_intern_border(const aPoints & p, int ib[2],p_pl_intern_spline ps,
           border_flag bf);
  pl_intern_border(int ib[2],p_pl_intern_spline psSpline,border_flag bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag _bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag,aSegments3D&);

  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;
  point get_a_top(int i)const{return aTops(i);};
  // ¾ßÐáÝÐï äãÝÚæØï, ÞÔÝÐÚÞ...
  int add_a_top(const point& pt);
  int n_of_tops()const{return aTops.count()-1;};
  virtual int contains(const point &pWhat)const;
  // ¿Þ ÞâàÕ×ÚÐÜ - ßÞáÛÕÔÞÒÐâÕÛìÝÞáâì âÞçÕÚ.
  int fix_points(aSegments3D&);
  const aPoints& points()const{return aTops;};
  point closest_border_point(const point&)const;

  // ¿ÕàÕÒÞÔ Ø× ÞÔÝØå ÚÞÞàÔØÝÐâ Ò ÔàãÓØÕ
  point2D local(const point3d&);
  point3d global(const point2D&);
  point P0(){return p0;};
  point middle()const;
  friend class pl_model;
};

typedef pl_intern_border* p_pl_intern_border;
#define AnyType p_pl_intern_border
#define AnyAr pl_intern_borders
#include <mix/AnyAr.H>


class pl_fict_border:public pl_intern_border,public pl_intern_spline{
public:
  pl_fict_border(const aPoints & p, int ib[2],
         double A,double B, double C, double D):
   pl_intern_border(p,ib,NULL,TEMP),pl_intern_spline(A,B,C,D){};
};

typedef pl_fict_border* p_pl_fict_border;
typedef ObjectsAr<p_pl_fict_border> pl_fict_borders;

// ÂÕÛÞ - ÝÐÑÞà ÜÝÞÓÞãÓÞÛìÝØÚÞÒ
class pl_intern_body:public any_body{
protected:
  pl_intern_body(pl_intern_body &bCopy);
  pl_intern_splines* get_bodies_splines()const;
public:
  pl_intern_body(geophys_params gp);
  pl_intern_body(any_borders &ab,geophys_params gp);
  ~pl_intern_body(){;};

  // ´ÞÑÐÒÛïÕÜ Ò íâÞâ ÚÛÐáá, âÐÚ ÚÐÚ âàÕÑãÕâáï ßÕàÕáÕçÕÝØÕ á ßÛÞáÚÞáâïÜØ
  virtual int contains(const point &pWhat)const;

  // ÁÞ×ÔÐÝØÕ äØÚâØÒÝÞÙ ÓàÐÝØæë ÒÝãâàØ âÕÛÐ ßÞ ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâØ
  int fict_border (const pl_spline& sSlice,p_pl_fict_border& bNew)const;

friend class pl_model;  
};

typedef pl_intern_body* p_pl_intern_body;
typedef ObjectsAr<p_pl_intern_body> pl_intern_bodies;

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

class pl_spline:public any_ex_spline{
protected:
  double a,b,c,d;

public:
  pl_spline();
  pl_spline(double A,double B,double C, double D);

  // ßÛÞáÚÞáâì ßÞ âàÕÜ âÞçÚÐÜ
  pl_spline(const point& p1,const point& p2,const point& p3);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  void set_coef(double A,double B,double C, double D);
  void get_coef(double &A,double &B,double &C, double &D);
  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};
  double distance(const point3d &pWhere);
};

typedef pl_spline* p_pl_spline;
typedef ObjectsAr<p_pl_spline> p_pl_splines;

class pl_border:public any_ex_border{
protected:
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
public:
  pl_border();
  pl_border(const aPoints & at, int ib[2], int is, border_flag bf);
  pl_border(int ib[2], int is, border_flag bf);
  point get_a_top(int i)const {return point(aTops_X(i),aTops_Y(i),aTops_Z(i));};
  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;

  // °ÝÐÛÞÓ ÚÞÝáâàãÚâÞàÐ ÔÛï ßÕàÕ×ÐßØáØ áÞ×ÔÐÝÝÞÓÞ ÞÑêÕÚâÐ
  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞ
class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

friend class pl_model;
};

typedef pl_intern_spline* p_pl_intern_spline;
typedef ObjectsAr<p_pl_intern_spline> pl_intern_splines;


// ³àÐÝØæÐ á áÞÞâÒÕâáâÒãîéØÜØ ÓÕÞäØ×ØçÕáÚØÜØ ßÐàÐÜÕâàÐÜØclass pl_overlying_spline:public pl_intern_spline{
protected:
  geophys_params gp;
  point3d pInternal; // ÂÞçÚÐ, ÛÕÖÐéÐï ÝÐ ÓàÐÝØæÕ, ÞÑàÐ×ÞÒÐÝÝÞÙ íâØÜ áßÛÐÙÝÞÜ
public:
  pl_overlying_spline(const geophys_params &gp,const point3d& p,
              double dFi,double dAz);
  geophys_params get_geophys(){return gp;};
  point3d get_internal_point(){return pInternal;};
};


typedef pl_overlying_spline* p_pl_overlying_spline;
typedef ObjectsAr<p_pl_overlying_spline> pl_overlying_splines;

// ³àÐÝØæÐ - ÜÝÞÓÞãÓÞÛìÝØÚ
// ÂàÕÑãÕâáï ÝÕÚÞâÞàëÙ áßÞáÞÑ ÞßàÕÔÕÛÕÝØï æÕÛÞáâÝÞáâØ ÓàÐÝØæë...
class pl_intern_border:public any_border{
protected:
  aPoints aTops; // ²ÕàèØÝë ÜÝÞÓÞãÓÞÛìÝØÚÐ, âàÕÑÞÒÐÝØÕ: aTops[0]=aTops[_last_].
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
  pl_intern_border(pl_intern_border &bCopy);

  // ÂÞçÚØ, ÞßàÕÔÕÛïîéØÕ ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ ÝÐ ßÐßÕ-áßÛÐÙÝÕ
  Logic points_created;
  point p0;

  void init_local_coords(); // ²ëçØáÛïÕâ âÞçÚã p0
public:
  ~pl_intern_border();
  pl_intern_border(const aPoints & p, int ib[2],p_pl_intern_spline ps,
           border_flag bf);
  pl_intern_border(int ib[2],p_pl_intern_spline psSpline,border_flag bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag _bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag,aSegments3D&);

  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;
  point get_a_top(int i)const{return aTops(i);};
  // ¾ßÐáÝÐï äãÝÚæØï, ÞÔÝÐÚÞ...
  int add_a_top(const point& pt);
  int n_of_tops()const{return aTops.count()-1;};
  virtual int contains(const point &pWhat)const;
  // ¿Þ ÞâàÕ×ÚÐÜ - ßÞáÛÕÔÞÒÐâÕÛìÝÞáâì âÞçÕÚ.
  int fix_points(aSegments3D&);
  const aPoints& points()const{return aTops;};
  point closest_border_point(const point&)const;

  // ¿ÕàÕÒÞÔ Ø× ÞÔÝØå ÚÞÞàÔØÝÐâ Ò ÔàãÓØÕ
  point2D local(const point3d&);
  point3d global(const point2D&);
  point P0(){return p0;};
  point middle()const;
  friend class pl_model;
};

typedef pl_intern_border* p_pl_intern_border;
#define AnyType p_pl_intern_border
#define AnyAr pl_intern_borders
#include <mix/AnyAr.H>


class pl_fict_border:public pl_intern_border,public pl_intern_spline{
public:
  pl_fict_border(const aPoints & p, int ib[2],
         double A,double B, double C, double D):
   pl_intern_border(p,ib,NULL,TEMP),pl_intern_spline(A,B,C,D){};
};

typedef pl_fict_border* p_pl_fict_border;
typedef ObjectsAr<p_pl_fict_border> pl_fict_borders;

// ÂÕÛÞ - ÝÐÑÞà ÜÝÞÓÞãÓÞÛìÝØÚÞÒ
class pl_intern_body:public any_body{
protected:
  pl_intern_body(pl_intern_body &bCopy);
  pl_intern_splines* get_bodies_splines()const;
public:
  pl_intern_body(geophys_params gp);
  pl_intern_body(any_borders &ab,geophys_params gp);
  ~pl_intern_body(){;};

  // ´ÞÑÐÒÛïÕÜ Ò íâÞâ ÚÛÐáá, âÐÚ ÚÐÚ âàÕÑãÕâáï ßÕàÕáÕçÕÝØÕ á ßÛÞáÚÞáâïÜØ
  virtual int contains(const point &pWhat)const;

  // ÁÞ×ÔÐÝØÕ äØÚâØÒÝÞÙ ÓàÐÝØæë ÒÝãâàØ âÕÛÐ ßÞ ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâØ
  int fict_border (const pl_spline& sSlice,p_pl_fict_border& bNew)const;

friend class pl_model;  
};

typedef pl_intern_body* p_pl_intern_body;
typedef ObjectsAr<p_pl_intern_body> pl_intern_bodies;

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

class pl_spline:public any_ex_spline{
protected:
  double a,b,c,d;

public:
  pl_spline();
  pl_spline(double A,double B,double C, double D);

  // ßÛÞáÚÞáâì ßÞ âàÕÜ âÞçÚÐÜ
  pl_spline(const point& p1,const point& p2,const point& p3);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  void set_coef(double A,double B,double C, double D);
  void get_coef(double &A,double &B,double &C, double &D);
  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};
  double distance(const point3d &pWhere);
};

typedef pl_spline* p_pl_spline;
typedef ObjectsAr<p_pl_spline> p_pl_splines;

class pl_border:public any_ex_border{
protected:
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
public:
  pl_border();
  pl_border(const aPoints & at, int ib[2], int is, border_flag bf);
  pl_border(int ib[2], int is, border_flag bf);
  point get_a_top(int i)const {return point(aTops_X(i),aTops_Y(i),aTops_Z(i));};
  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;

  // °ÝÐÛÞÓ ÚÞÝáâàãÚâÞàÐ ÔÛï ßÕàÕ×ÐßØáØ áÞ×ÔÐÝÝÞÓÞ ÞÑêÕÚâÐ
  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞ
class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

class pl_overlying_spline:public pl_intern_spline{
protected:
  geophys_params gp;
  point3d pInternal; // ÂÞçÚÐ, ÛÕÖÐéÐï ÝÐ ÓàÐÝØæÕ, ÞÑàÐ×ÞÒÐÝÝÞÙ íâØÜ áßÛÐÙÝÞÜpublic:
  pl_overlying_spline(const geophys_params &gp,const point3d& p,
              double dFi,double dAz);
  geophys_params get_geophys(){return gp;};
  point3d get_internal_point(){return pInternal;};
};


typedef pl_overlying_spline* p_pl_overlying_spline;
typedef ObjectsAr<p_pl_overlying_spline> pl_overlying_splines;

// ³àÐÝØæÐ - ÜÝÞÓÞãÓÞÛìÝØÚ
// ÂàÕÑãÕâáï ÝÕÚÞâÞàëÙ áßÞáÞÑ ÞßàÕÔÕÛÕÝØï æÕÛÞáâÝÞáâØ ÓàÐÝØæë...
class pl_intern_border:public any_border{
protected:
  aPoints aTops; // ²ÕàèØÝë ÜÝÞÓÞãÓÞÛìÝØÚÐ, âàÕÑÞÒÐÝØÕ: aTops[0]=aTops[_last_].
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
  pl_intern_border(pl_intern_border &bCopy);

  // ÂÞçÚØ, ÞßàÕÔÕÛïîéØÕ ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ ÝÐ ßÐßÕ-áßÛÐÙÝÕ
  Logic points_created;
  point p0;

  void init_local_coords(); // ²ëçØáÛïÕâ âÞçÚã p0
public:
  ~pl_intern_border();
  pl_intern_border(const aPoints & p, int ib[2],p_pl_intern_spline ps,
           border_flag bf);
  pl_intern_border(int ib[2],p_pl_intern_spline psSpline,border_flag bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag _bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag,aSegments3D&);

  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;
  point get_a_top(int i)const{return aTops(i);};
  // ¾ßÐáÝÐï äãÝÚæØï, ÞÔÝÐÚÞ...
  int add_a_top(const point& pt);
  int n_of_tops()const{return aTops.count()-1;};
  virtual int contains(const point &pWhat)const;
  // ¿Þ ÞâàÕ×ÚÐÜ - ßÞáÛÕÔÞÒÐâÕÛìÝÞáâì âÞçÕÚ.
  int fix_points(aSegments3D&);
  const aPoints& points()const{return aTops;};
  point closest_border_point(const point&)const;

  // ¿ÕàÕÒÞÔ Ø× ÞÔÝØå ÚÞÞàÔØÝÐâ Ò ÔàãÓØÕ
  point2D local(const point3d&);
  point3d global(const point2D&);
  point P0(){return p0;};
  point middle()const;
  friend class pl_model;
};

typedef pl_intern_border* p_pl_intern_border;
#define AnyType p_pl_intern_border
#define AnyAr pl_intern_borders
#include <mix/AnyAr.H>


class pl_fict_border:public pl_intern_border,public pl_intern_spline{
public:
  pl_fict_border(const aPoints & p, int ib[2],
         double A,double B, double C, double D):
   pl_intern_border(p,ib,NULL,TEMP),pl_intern_spline(A,B,C,D){};
};

typedef pl_fict_border* p_pl_fict_border;
typedef ObjectsAr<p_pl_fict_border> pl_fict_borders;

// ÂÕÛÞ - ÝÐÑÞà ÜÝÞÓÞãÓÞÛìÝØÚÞÒ
class pl_intern_body:public any_body{
protected:
  pl_intern_body(pl_intern_body &bCopy);
  pl_intern_splines* get_bodies_splines()const;
public:
  pl_intern_body(geophys_params gp);
  pl_intern_body(any_borders &ab,geophys_params gp);
  ~pl_intern_body(){;};

  // ´ÞÑÐÒÛïÕÜ Ò íâÞâ ÚÛÐáá, âÐÚ ÚÐÚ âàÕÑãÕâáï ßÕàÕáÕçÕÝØÕ á ßÛÞáÚÞáâïÜØ
  virtual int contains(const point &pWhat)const;

  // ÁÞ×ÔÐÝØÕ äØÚâØÒÝÞÙ ÓàÐÝØæë ÒÝãâàØ âÕÛÐ ßÞ ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâØ
  int fict_border (const pl_spline& sSlice,p_pl_fict_border& bNew)const;

friend class pl_model;  
};

typedef pl_intern_body* p_pl_intern_body;
typedef ObjectsAr<p_pl_intern_body> pl_intern_bodies;

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

class pl_spline:public any_ex_spline{
protected:
  double a,b,c,d;

public:
  pl_spline();
  pl_spline(double A,double B,double C, double D);

  // ßÛÞáÚÞáâì ßÞ âàÕÜ âÞçÚÐÜ
  pl_spline(const point& p1,const point& p2,const point& p3);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  void set_coef(double A,double B,double C, double D);
  void get_coef(double &A,double &B,double &C, double &D);
  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};
  double distance(const point3d &pWhere);
};

typedef pl_spline* p_pl_spline;
typedef ObjectsAr<p_pl_spline> p_pl_splines;

class pl_border:public any_ex_border{
protected:
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
public:
  pl_border();
  pl_border(const aPoints & at, int ib[2], int is, border_flag bf);
  pl_border(int ib[2], int is, border_flag bf);
  point get_a_top(int i)const {return point(aTops_X(i),aTops_Y(i),aTops_Z(i));};
  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;

  // °ÝÐÛÞÓ ÚÞÝáâàãÚâÞàÐ ÔÛï ßÕàÕ×ÐßØáØ áÞ×ÔÐÝÝÞÓÞ ÞÑêÕÚâÐ
  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞ
class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

public:
  pl_overlying_spline(const geophys_params &gp,const point3d& p,
              double dFi,double dAz);
  geophys_params get_geophys(){return gp;};
  point3d get_internal_point(){return pInternal;};
};


typedef pl_overlying_spline* p_pl_overlying_spline;
typedef ObjectsAr<p_pl_overlying_spline> pl_overlying_splines;

// ³àÐÝØæÐ - ÜÝÞÓÞãÓÞÛìÝØÚ// ÂàÕÑãÕâáï ÝÕÚÞâÞàëÙ áßÞáÞÑ ÞßàÕÔÕÛÕÝØï æÕÛÞáâÝÞáâØ ÓàÐÝØæë...
class pl_intern_border:public any_border{
protected:
  aPoints aTops; // ²ÕàèØÝë ÜÝÞÓÞãÓÞÛìÝØÚÐ, âàÕÑÞÒÐÝØÕ: aTops[0]=aTops[_last_].
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
  pl_intern_border(pl_intern_border &bCopy);

  // ÂÞçÚØ, ÞßàÕÔÕÛïîéØÕ ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ ÝÐ ßÐßÕ-áßÛÐÙÝÕ
  Logic points_created;
  point p0;

  void init_local_coords(); // ²ëçØáÛïÕâ âÞçÚã p0
public:
  ~pl_intern_border();
  pl_intern_border(const aPoints & p, int ib[2],p_pl_intern_spline ps,
           border_flag bf);
  pl_intern_border(int ib[2],p_pl_intern_spline psSpline,border_flag bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag _bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag,aSegments3D&);

  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;
  point get_a_top(int i)const{return aTops(i);};
  // ¾ßÐáÝÐï äãÝÚæØï, ÞÔÝÐÚÞ...
  int add_a_top(const point& pt);
  int n_of_tops()const{return aTops.count()-1;};
  virtual int contains(const point &pWhat)const;
  // ¿Þ ÞâàÕ×ÚÐÜ - ßÞáÛÕÔÞÒÐâÕÛìÝÞáâì âÞçÕÚ.
  int fix_points(aSegments3D&);
  const aPoints& points()const{return aTops;};
  point closest_border_point(const point&)const;

  // ¿ÕàÕÒÞÔ Ø× ÞÔÝØå ÚÞÞàÔØÝÐâ Ò ÔàãÓØÕ
  point2D local(const point3d&);
  point3d global(const point2D&);
  point P0(){return p0;};
  point middle()const;
  friend class pl_model;
};

typedef pl_intern_border* p_pl_intern_border;
#define AnyType p_pl_intern_border
#define AnyAr pl_intern_borders
#include <mix/AnyAr.H>


class pl_fict_border:public pl_intern_border,public pl_intern_spline{
public:
  pl_fict_border(const aPoints & p, int ib[2],
         double A,double B, double C, double D):
   pl_intern_border(p,ib,NULL,TEMP),pl_intern_spline(A,B,C,D){};
};

typedef pl_fict_border* p_pl_fict_border;
typedef ObjectsAr<p_pl_fict_border> pl_fict_borders;

// ÂÕÛÞ - ÝÐÑÞà ÜÝÞÓÞãÓÞÛìÝØÚÞÒ
class pl_intern_body:public any_body{
protected:
  pl_intern_body(pl_intern_body &bCopy);
  pl_intern_splines* get_bodies_splines()const;
public:
  pl_intern_body(geophys_params gp);
  pl_intern_body(any_borders &ab,geophys_params gp);
  ~pl_intern_body(){;};

  // ´ÞÑÐÒÛïÕÜ Ò íâÞâ ÚÛÐáá, âÐÚ ÚÐÚ âàÕÑãÕâáï ßÕàÕáÕçÕÝØÕ á ßÛÞáÚÞáâïÜØ
  virtual int contains(const point &pWhat)const;

  // ÁÞ×ÔÐÝØÕ äØÚâØÒÝÞÙ ÓàÐÝØæë ÒÝãâàØ âÕÛÐ ßÞ ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâØ
  int fict_border (const pl_spline& sSlice,p_pl_fict_border& bNew)const;

friend class pl_model;  
};

typedef pl_intern_body* p_pl_intern_body;
typedef ObjectsAr<p_pl_intern_body> pl_intern_bodies;

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

class pl_spline:public any_ex_spline{
protected:
  double a,b,c,d;

public:
  pl_spline();
  pl_spline(double A,double B,double C, double D);

  // ßÛÞáÚÞáâì ßÞ âàÕÜ âÞçÚÐÜ
  pl_spline(const point& p1,const point& p2,const point& p3);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  void set_coef(double A,double B,double C, double D);
  void get_coef(double &A,double &B,double &C, double &D);
  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};
  double distance(const point3d &pWhere);
};

typedef pl_spline* p_pl_spline;
typedef ObjectsAr<p_pl_spline> p_pl_splines;

class pl_border:public any_ex_border{
protected:
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
public:
  pl_border();
  pl_border(const aPoints & at, int ib[2], int is, border_flag bf);
  pl_border(int ib[2], int is, border_flag bf);
  point get_a_top(int i)const {return point(aTops_X(i),aTops_Y(i),aTops_Z(i));};
  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;

  // °ÝÐÛÞÓ ÚÞÝáâàãÚâÞàÐ ÔÛï ßÕàÕ×ÐßØáØ áÞ×ÔÐÝÝÞÓÞ ÞÑêÕÚâÐ
  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞ
class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

// ÂàÕÑãÕâáï ÝÕÚÞâÞàëÙ áßÞáÞÑ ÞßàÕÔÕÛÕÝØï æÕÛÞáâÝÞáâØ ÓàÐÝØæë...
class pl_intern_border:public any_border{
protected:
  aPoints aTops; // ²ÕàèØÝë ÜÝÞÓÞãÓÞÛìÝØÚÐ, âàÕÑÞÒÐÝØÕ: aTops[0]=aTops[_last_].
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
  pl_intern_border(pl_intern_border &bCopy);

  // ÂÞçÚØ, ÞßàÕÔÕÛïîéØÕ ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ ÝÐ ßÐßÕ-áßÛÐÙÝÕ
  Logic points_created;
  point p0;

  void init_local_coords(); // ²ëçØáÛïÕâ âÞçÚã p0
public:
  ~pl_intern_border();
  pl_intern_border(const aPoints & p, int ib[2],p_pl_intern_spline ps,
           border_flag bf);
  pl_intern_border(int ib[2],p_pl_intern_spline psSpline,border_flag bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag _bf);
  pl_intern_border(p_pl_intern_spline ps,border_flag,aSegments3D&);

  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;
  point get_a_top(int i)const{return aTops(i);};
  // ¾ßÐáÝÐï äãÝÚæØï, ÞÔÝÐÚÞ...
  int add_a_top(const point& pt);
  int n_of_tops()const{return aTops.count()-1;};
  virtual int contains(const point &pWhat)const;
  // ¿Þ ÞâàÕ×ÚÐÜ - ßÞáÛÕÔÞÒÐâÕÛìÝÞáâì âÞçÕÚ.
  int fix_points(aSegments3D&);
  const aPoints& points()const{return aTops;};
  point closest_border_point(const point&)const;

  // ¿ÕàÕÒÞÔ Ø× ÞÔÝØå ÚÞÞàÔØÝÐâ Ò ÔàãÓØÕ  point2D local(const point3d&);
  point3d global(const point2D&);
  point P0(){return p0;};
  point middle()const;
  friend class pl_model;
};

typedef pl_intern_border* p_pl_intern_border;
#define AnyType p_pl_intern_border
#define AnyAr pl_intern_borders
#include <mix/AnyAr.H>


class pl_fict_border:public pl_intern_border,public pl_intern_spline{
public:
  pl_fict_border(const aPoints & p, int ib[2],
         double A,double B, double C, double D):
   pl_intern_border(p,ib,NULL,TEMP),pl_intern_spline(A,B,C,D){};
};

typedef pl_fict_border* p_pl_fict_border;
typedef ObjectsAr<p_pl_fict_border> pl_fict_borders;

// ÂÕÛÞ - ÝÐÑÞà ÜÝÞÓÞãÓÞÛìÝØÚÞÒ
class pl_intern_body:public any_body{
protected:
  pl_intern_body(pl_intern_body &bCopy);
  pl_intern_splines* get_bodies_splines()const;
public:
  pl_intern_body(geophys_params gp);
  pl_intern_body(any_borders &ab,geophys_params gp);
  ~pl_intern_body(){;};

  // ´ÞÑÐÒÛïÕÜ Ò íâÞâ ÚÛÐáá, âÐÚ ÚÐÚ âàÕÑãÕâáï ßÕàÕáÕçÕÝØÕ á ßÛÞáÚÞáâïÜØ
  virtual int contains(const point &pWhat)const;

  // ÁÞ×ÔÐÝØÕ äØÚâØÒÝÞÙ ÓàÐÝØæë ÒÝãâàØ âÕÛÐ ßÞ ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâØ
  int fict_border (const pl_spline& sSlice,p_pl_fict_border& bNew)const;

friend class pl_model;  
};

typedef pl_intern_body* p_pl_intern_body;
typedef ObjectsAr<p_pl_intern_body> pl_intern_bodies;

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

class pl_spline:public any_ex_spline{
protected:
  double a,b,c,d;

public:
  pl_spline();
  pl_spline(double A,double B,double C, double D);

  // ßÛÞáÚÞáâì ßÞ âàÕÜ âÞçÚÐÜ
  pl_spline(const point& p1,const point& p2,const point& p3);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  void set_coef(double A,double B,double C, double D);
  void get_coef(double &A,double &B,double &C, double &D);
  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};
  double distance(const point3d &pWhere);
};

typedef pl_spline* p_pl_spline;
typedef ObjectsAr<p_pl_spline> p_pl_splines;

class pl_border:public any_ex_border{
protected:
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
public:
  pl_border();
  pl_border(const aPoints & at, int ib[2], int is, border_flag bf);
  pl_border(int ib[2], int is, border_flag bf);
  point get_a_top(int i)const {return point(aTops_X(i),aTops_Y(i),aTops_Z(i));};
  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;

  // °ÝÐÛÞÓ ÚÞÝáâàãÚâÞàÐ ÔÛï ßÕàÕ×ÐßØáØ áÞ×ÔÐÝÝÞÓÞ ÞÑêÕÚâÐ
  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞ
class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

  point2D local(const point3d&);
  point3d global(const point2D&);
  point P0(){return p0;};
  point middle()const;
  friend class pl_model;
};

typedef pl_intern_border* p_pl_intern_border;
#define AnyType p_pl_intern_border
#define AnyAr pl_intern_borders
#include <mix/AnyAr.H>


class pl_fict_border:public pl_intern_border,public pl_intern_spline{
public:
  pl_fict_border(const aPoints & p, int ib[2],
         double A,double B, double C, double D):
   pl_intern_border(p,ib,NULL,TEMP),pl_intern_spline(A,B,C,D){};
};

typedef pl_fict_border* p_pl_fict_border;
typedef ObjectsAr<p_pl_fict_border> pl_fict_borders;

// ÂÕÛÞ - ÝÐÑÞà ÜÝÞÓÞãÓÞÛìÝØÚÞÒclass pl_intern_body:public any_body{
protected:
  pl_intern_body(pl_intern_body &bCopy);
  pl_intern_splines* get_bodies_splines()const;
public:
  pl_intern_body(geophys_params gp);
  pl_intern_body(any_borders &ab,geophys_params gp);
  ~pl_intern_body(){;};

  // ´ÞÑÐÒÛïÕÜ Ò íâÞâ ÚÛÐáá, âÐÚ ÚÐÚ âàÕÑãÕâáï ßÕàÕáÕçÕÝØÕ á ßÛÞáÚÞáâïÜØ
  virtual int contains(const point &pWhat)const;

  // ÁÞ×ÔÐÝØÕ äØÚâØÒÝÞÙ ÓàÐÝØæë ÒÝãâàØ âÕÛÐ ßÞ ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâØ
  int fict_border (const pl_spline& sSlice,p_pl_fict_border& bNew)const;

friend class pl_model;  
};

typedef pl_intern_body* p_pl_intern_body;
typedef ObjectsAr<p_pl_intern_body> pl_intern_bodies;

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

class pl_spline:public any_ex_spline{
protected:
  double a,b,c,d;

public:
  pl_spline();
  pl_spline(double A,double B,double C, double D);

  // ßÛÞáÚÞáâì ßÞ âàÕÜ âÞçÚÐÜ
  pl_spline(const point& p1,const point& p2,const point& p3);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  void set_coef(double A,double B,double C, double D);
  void get_coef(double &A,double &B,double &C, double &D);
  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};
  double distance(const point3d &pWhere);
};

typedef pl_spline* p_pl_spline;
typedef ObjectsAr<p_pl_spline> p_pl_splines;

class pl_border:public any_ex_border{
protected:
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
public:
  pl_border();
  pl_border(const aPoints & at, int ib[2], int is, border_flag bf);
  pl_border(int ib[2], int is, border_flag bf);
  point get_a_top(int i)const {return point(aTops_X(i),aTops_Y(i),aTops_Z(i));};
  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;

  // °ÝÐÛÞÓ ÚÞÝáâàãÚâÞàÐ ÔÛï ßÕàÕ×ÐßØáØ áÞ×ÔÐÝÝÞÓÞ ÞÑêÕÚâÐ
  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞ
class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

class pl_intern_body:public any_body{
protected:
  pl_intern_body(pl_intern_body &bCopy);
  pl_intern_splines* get_bodies_splines()const;
public:
  pl_intern_body(geophys_params gp);
  pl_intern_body(any_borders &ab,geophys_params gp);
  ~pl_intern_body(){;};

  // ´ÞÑÐÒÛïÕÜ Ò íâÞâ ÚÛÐáá, âÐÚ ÚÐÚ âàÕÑãÕâáï ßÕàÕáÕçÕÝØÕ á ßÛÞáÚÞáâïÜØ
  virtual int contains(const point &pWhat)const;

  // ÁÞ×ÔÐÝØÕ äØÚâØÒÝÞÙ ÓàÐÝØæë ÒÝãâàØ âÕÛÐ ßÞ ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâØ
  int fict_border (const pl_spline& sSlice,p_pl_fict_border& bNew)const;

friend class pl_model;  
};

typedef pl_intern_body* p_pl_intern_body;
typedef ObjectsAr<p_pl_intern_body> pl_intern_bodies;

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

class pl_spline:public any_ex_spline{
protected:
  double a,b,c,d;

public:
  pl_spline();
  pl_spline(double A,double B,double C, double D);

  // ßÛÞáÚÞáâì ßÞ âàÕÜ âÞçÚÐÜ  pl_spline(const point& p1,const point& p2,const point& p3);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ
  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  void set_coef(double A,double B,double C, double D);
  void get_coef(double &A,double &B,double &C, double &D);
  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};
  double distance(const point3d &pWhere);
};

typedef pl_spline* p_pl_spline;
typedef ObjectsAr<p_pl_spline> p_pl_splines;

class pl_border:public any_ex_border{
protected:
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
public:
  pl_border();
  pl_border(const aPoints & at, int ib[2], int is, border_flag bf);
  pl_border(int ib[2], int is, border_flag bf);
  point get_a_top(int i)const {return point(aTops_X(i),aTops_Y(i),aTops_Z(i));};
  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;

  // °ÝÐÛÞÓ ÚÞÝáâàãÚâÞàÐ ÔÛï ßÕàÕ×ÐßØáØ áÞ×ÔÐÝÝÞÓÞ ÞÑêÕÚâÐ
  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞ
class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

  pl_spline(const point& p1,const point& p2,const point& p3);

  // ¿ÕàÕáÕçÕÝØÕ á ÞâàÕ×ÚÞÜ ÒåÞÔÝëå âÞçÕÚ
  point cross(const point& p1,const point& p2) const;

  // ¿ÕàÕÒÞÔ Ò ÛÞÚÐÛìÝãî áØáâÕÜã ÚÞÞàÔØÝÐâ  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  void set_coef(double A,double B,double C, double D);
  void get_coef(double &A,double &B,double &C, double &D);
  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};
  double distance(const point3d &pWhere);
};

typedef pl_spline* p_pl_spline;
typedef ObjectsAr<p_pl_spline> p_pl_splines;

class pl_border:public any_ex_border{
protected:
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
public:
  pl_border();
  pl_border(const aPoints & at, int ib[2], int is, border_flag bf);
  pl_border(int ib[2], int is, border_flag bf);
  point get_a_top(int i)const {return point(aTops_X(i),aTops_Y(i),aTops_Z(i));};
  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;

  // °ÝÐÛÞÓ ÚÞÝáâàãÚâÞàÐ ÔÛï ßÕàÕ×ÐßØáØ áÞ×ÔÐÝÝÞÓÞ ÞÑêÕÚâÐ
  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞ
class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

  point2D local(const point& pMain,
        const point& p0,const point& p1,const point& p2) const;

  // Á ÚÐÚÞÙ áâÞàÞÝë Þâ ßÛÞáÚÞáâØ ÛÕÖØâ âÞçÚÐ
  int sign(const point & p) const;

  void set_coef(double A,double B,double C, double D);
  void get_coef(double &A,double &B,double &C, double &D);
  double A() const {return a;};
  double B() const {return b;};
  double C() const {return c;};
  double D() const {return d;};
  double distance(const point3d &pWhere);
};

typedef pl_spline* p_pl_spline;
typedef ObjectsAr<p_pl_spline> p_pl_splines;

class pl_border:public any_ex_border{
protected:
  DoubleAr aTops_X; 
  DoubleAr aTops_Y; 
  DoubleAr aTops_Z; 
public:
  pl_border();
  pl_border(const aPoints & at, int ib[2], int is, border_flag bf);
  pl_border(int ib[2], int is, border_flag bf);
  point get_a_top(int i)const {return point(aTops_X(i),aTops_Y(i),aTops_Z(i));};
  const DoubleAr& get_x()const;
  const DoubleAr& get_y()const;
  const DoubleAr& get_z()const;

  // °ÝÐÛÞÓ ÚÞÝáâàãÚâÞàÐ ÔÛï ßÕàÕ×ÐßØáØ áÞ×ÔÐÝÝÞÓÞ ÞÑêÕÚâÐ  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞ
class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

  int make_border(DoubleAr&ax,DoubleAr&ay,DoubleAr&az,
          int ib[2], int is, border_flag bf);

  int n_of_tops()const{return aTops_Z.count();};
  int get_body(int ib);
  int get_spline()const {return iSpline;};
};

// ¿ÞçâØ âÞ ÖÕ, çâÞ Ø ßàÞØ×ÒÞÛìÝÞÕ âÕÛÞclass pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */ 

class pl_body:public any_ex_body{
public:
  pl_body():any_ex_body(){;};
  pl_body(IntAr & ab, const geophys_params & gp, int in);
};


/* ÄãÝÚæØï-ÛîÜßÕÝ...*/
void get_splines_weights(const point& pLast,
             const pl_intern_splines & paSplines,
             const any_borders & paBorders,
             DoubleAr&daWeights);
#endif /* p_elements.H */
/*  Melnikov George Y.  16-Sep-1998  */
/*          ...          */ 
/*  Melnikov George Yu. 31-Jan-2000  */