p_model.H

/* p_model.H */
/* $Id: p_model.H,v 2.26 2005/07/18 07:15:51 nikita Exp $ */
#ifndef __p_model_H
#define __p_model_H

#include <mix/Progress.hpp>
#include "any_elements.H"
#include "any_model.H"
#include "p_elements.H"
#include "db_managers.H"
/*
   ºÞÝÚàÕâÝÐï ÜÞÔÕÛì: áßÛÐÙÝë - ßÛÞáÚÞáâØ, Ð ÓàÐÝØæë âÕÛ - ÜÝÞÓÞãÓÞÛìÝØÚØ.
   ´ÞßÞÛÝØâÕÛìÝëÕ ÒÞ×ÜÞÖÝÞáâØ ßÞ áàÐÒÝÕÝØî á ÑÐ×ÞÒëÜ ÚÛÐááÞÜ -
   - ÞâÞÑàÐÖÕÝØÕ íÛÕÜÕÝâÞÒ ÜÞÔÕÛØ ßÞÛì×ÞÒÐâÕÛî (Ð ÝÕ âÞÛìÚÞ Øå ÚÞÛØçÕáâÒÐ)
   - ÒÜÕáâÞ any_border(s) Ø ßÞÔÞÑÝëå - Øå ÔÕâØ p_...
// ÚÛÐááë Ø× p_elements.H
class point;
class point2D;
class aSegments2D;
class pl_intern_spline;
class pl_intern_border;
class pl_intern_body;
class pl_spline;
class pl_border;
class pl_body;
class pl_overlying_splines;
 */

class pl_model:public any_model{
protected:
  db_volume_manager *pdbvTables;
  db_results_manager *pdbmTables;
  IntAr* piaInclinometry;
  aPoints* papInclinometry;
  pl_intern_borders modelBorders; //nikita 22.04.05
  Logic parallel_structure_flag;
  //Parameters for parallel model
  parallel_model_params *pPMP;

  int read_from_DB(const dbItem &dbiWhere);  // ÁçØâÐÛØ ÞâÚãÔÐ ÝÐÔÞ
  // ½ÐÙâØ ßÞ âÞçÚÕ âÕÛÞ  int get_body(const point &pWhere, pl_intern_body * &pb) const;
  // ÁÒï×ëÒÐÝØÕ ÝÐÑÞàÐ ßÛÞáÚÞáâÕÙ Ò áØáâÕÜã âÕÛ Ø ÓàÐÝØæ.
  virtual int correlate(const pl_overlying_splines& apOSplines,
            pl_intern_body* bInitial,const IntAr& iaSequence);

  // ´Ûï ØáßÞÛì×ÞÒÐÝØï Ò ÚÞÝáâàãÚâÞàÕ ßÞ ÑÐ×Õ ÔÐÝÝëå
  virtual pl_intern_border* new_border(pl_intern_spline*,border_flag);
  virtual pl_intern_spline* new_spline(double,double,double,double);
  virtual pl_intern_body* new_body(geophys_params&);
  void add_breakpoints (const aPoints&apInclIn,int iI,int iCBody,int iTB);
  Logic cross_border (const point &p1,const point & p2,point&pRes,int iB);
  void incl_pass_body (const point &p1,const point & p2,int iC,int&iN);

  /**Read segments defining borders of model volume */
  ErrCode readModelVolumeSegs(dbHandle);//nikita 25.04.05
  dbNameStruct dbsModGlass;
  int write_bodies()const;
  // ¿ãáâÞÙ, çâÞÑë Ø× ÒÕàåÝØå ãàÞÒÝÕÙ ÝÕ Òë×ëÒÐâì ÝØçÕÓÞ
  pl_model(){
    pdbmTables=NULL;
    piaInclinometry=NULL;
    papInclinometry=NULL;
    pdbvTables=NULL;
    pPMP=NULL;
  };

public:
  pl_model(const dbItem &dbiWhere);
  pl_model(const dbNameStruct &dbh1D, const dbItem &dbiBorders,
       char* sC,IntAr *iaSequence=NULL);
  pl_model(pl_model& am);   // ºÞÝáâàãÚâÞà ÚÞßØàÞÒÐÝØï... ½Õ ÑëÒÐÕâ
  virtual ~pl_model();

  unsigned get_model_bounds_num();
  Logic get_model_bound(unsigned n, aPoints* points);

  // ¾âÞÑàÐÖÕÝØÕ íÛÕÜÕÝâÞÒ ÜÞÔÕÛØ
  int get_border(int i, pl_border &bord) const;
  // Ø× âÕÛÐ ßÞÔ ÝÞÜÕàÞÜ n1 ÓàÐÝØæã á ÝÞÜÕàÞÜ n2
  int get_border(int n1,int n2, pl_border &pb) const;
  int get_body(int i, pl_body &bod) const;
  // ÂÕÛÞ ßÞ áÞÔÕàÖÐéÕÙáï Ò ÝÕÜ âÞçÚÕ
  int get_body(const point &pWhere, pl_body &pb) const;
  int get_spline(int i, pl_spline &spl) const;

  // ½ÞÜÕà ÓàÐÝØæë Ò ÜÞÔÕÛØ ßÞ ÝÞÜÕàã ÓàÐÝØæë Ò âÕÛÕ
  int get_border_num(int ibod,int iborloc,int& iborglob) const;

  // ÄãÝÚæØï ÞâÞÑàÐÖÕÝØï ßÕàÕáÕçÕÝØï ßÛÞáÚÞáâÕÙ ÜÞÔÕÛØ á ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâìî
  int slice(const pl_spline& sSlice,         
        // »ÞÚÐÛìÝÐï áØáâÕÜÐ ÚÞÞàÔØÝÐâ
	    const point& p0,const point& p1,const point& p2,
        // ÁîÔÐ ×ÐáãÝÕÜ àÕ×ãÛìâÐâ
        aSegments2D & aSegments,
        // ÁîÔÐ ×ÐáãÝÕÜ ÓàÐÝØæë ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
        aSegments2D & aBorder,
        // Here push numbers of bounds to get possibility 
        // select bounds in 3DToolKit program interactively (inserted by nikita) 
        // this parameter is NULL by default to avoid possible uncompatibility with other UNIVERS programs  
        IntAr* borderNumbers = NULL)const;

  // ¾ßàÕÔÕÛÕÝØÕ âÞÓÞ, çâÞ âÞçÚÐ ÝÐåÞÔØâáï ÒÝãâàØ ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
  Logic in_model_area(point p)const;
  // ·ÐßØáì
  //  virtual int write_to_DB(const dbItem &dbiWhere)const;
  virtual int write_to_DB(Progress* zp=NULL)const;
  void get_shots_data(aPoints& apSPs,IntAr& iaDT,IntAr& iaUT){
    pdbmTables->get_shots_data(apSPs,iaDT,iaUT);
  };
  void delete_modeling_data(const point& pSP);

  // ¿Þ ØÝÚÛØÝÞÜÕâàØØ - ÝÐÑÞà âÞçÕÚ ßÕàÕáÕçÕÝØï ÓàÐÝØæ ÜÞÔÕÛØ ßÞ íâÞÙ ØÝÚÛ.
  const aPoints& install_inclinometry(const aPoints& apInclIn);
  // ¸ÝÚÛØÝÞÜÕâàØï ×ÐÔÐÕâáï ßÐàÞÙ ÚÞÞàÔØÝÐâ, ÒÕàâØÚÐÛìÝÐï ßÐÛÚÐ
  const aPoints& install_vertical_inclinometry(double X, double Y);

  // ¿ÞÛãçØâì ÔÐÝÝëÕ Þ ÜÞÔÕÛØ ßÞ ^ ÚàØÒÞÙ
  void get_vp(DoubleAr &)const;
  void get_kp(DoubleAr &)const;
  void get_vs(DoubleAr &)const;
  void get_ks(DoubleAr &)const;

  void set_vp(const DoubleAr &);
  void set_kp(const DoubleAr &);
  void set_vs(const DoubleAr &);
  void set_ks(const DoubleAr &);

  Logic get_parallel_flag()const {return parallel_structure_flag;};
  void set_parallel_flag(Logic Flag){parallel_structure_flag=Flag;};
  parallel_model_params* get_parallel_model_params() {return pPMP;};
  
  void fix_dbStruct(dbNameStruct& dbnsIn)const{pdbvTables->fix_dbStruct(dbnsIn);};
  void prepare_tables(point p ,int i1,int i2){
    pdbmTables->prepare_tables(p,i1,i2,this);
  };
  p_pl_spline get_point_direction(const point3d &pWhere)const;
  
  /* New improvements by Garry */
  Logic check_structure(); /* Checking of model structure for possibility to
                  raytracing as in parallel structured model */
  Logic get_layers_inclination(double& a, double& b, double& c);
  Logic set_params_for_parallel();
  Logic analyse_gp(int LeadingWaveType);

  virtual const char* class_id(){return "pl_model";};
  virtual Logic isDerivedFrom(const char* class_name);

};

#endif /* p_model.H */

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

  int get_body(const point &pWhere, pl_intern_body * &pb) const;
  // ÁÒï×ëÒÐÝØÕ ÝÐÑÞàÐ ßÛÞáÚÞáâÕÙ Ò áØáâÕÜã âÕÛ Ø ÓàÐÝØæ.
  virtual int correlate(const pl_overlying_splines& apOSplines,
            pl_intern_body* bInitial,const IntAr& iaSequence);

  // ´Ûï ØáßÞÛì×ÞÒÐÝØï Ò ÚÞÝáâàãÚâÞàÕ ßÞ ÑÐ×Õ ÔÐÝÝëå  virtual pl_intern_border* new_border(pl_intern_spline*,border_flag);
  virtual pl_intern_spline* new_spline(double,double,double,double);
  virtual pl_intern_body* new_body(geophys_params&);
  void add_breakpoints (const aPoints&apInclIn,int iI,int iCBody,int iTB);
  Logic cross_border (const point &p1,const point & p2,point&pRes,int iB);
  void incl_pass_body (const point &p1,const point & p2,int iC,int&iN);

  /**Read segments defining borders of model volume */
  ErrCode readModelVolumeSegs(dbHandle);//nikita 25.04.05
  dbNameStruct dbsModGlass;
  int write_bodies()const;
  // ¿ãáâÞÙ, çâÞÑë Ø× ÒÕàåÝØå ãàÞÒÝÕÙ ÝÕ Òë×ëÒÐâì ÝØçÕÓÞ
  pl_model(){
    pdbmTables=NULL;
    piaInclinometry=NULL;
    papInclinometry=NULL;
    pdbvTables=NULL;
    pPMP=NULL;
  };

public:
  pl_model(const dbItem &dbiWhere);
  pl_model(const dbNameStruct &dbh1D, const dbItem &dbiBorders,
       char* sC,IntAr *iaSequence=NULL);
  pl_model(pl_model& am);   // ºÞÝáâàãÚâÞà ÚÞßØàÞÒÐÝØï... ½Õ ÑëÒÐÕâ
  virtual ~pl_model();

  unsigned get_model_bounds_num();
  Logic get_model_bound(unsigned n, aPoints* points);

  // ¾âÞÑàÐÖÕÝØÕ íÛÕÜÕÝâÞÒ ÜÞÔÕÛØ
  int get_border(int i, pl_border &bord) const;
  // Ø× âÕÛÐ ßÞÔ ÝÞÜÕàÞÜ n1 ÓàÐÝØæã á ÝÞÜÕàÞÜ n2
  int get_border(int n1,int n2, pl_border &pb) const;
  int get_body(int i, pl_body &bod) const;
  // ÂÕÛÞ ßÞ áÞÔÕàÖÐéÕÙáï Ò ÝÕÜ âÞçÚÕ
  int get_body(const point &pWhere, pl_body &pb) const;
  int get_spline(int i, pl_spline &spl) const;

  // ½ÞÜÕà ÓàÐÝØæë Ò ÜÞÔÕÛØ ßÞ ÝÞÜÕàã ÓàÐÝØæë Ò âÕÛÕ
  int get_border_num(int ibod,int iborloc,int& iborglob) const;

  // ÄãÝÚæØï ÞâÞÑàÐÖÕÝØï ßÕàÕáÕçÕÝØï ßÛÞáÚÞáâÕÙ ÜÞÔÕÛØ á ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâìî
  int slice(const pl_spline& sSlice,         
        // »ÞÚÐÛìÝÐï áØáâÕÜÐ ÚÞÞàÔØÝÐâ
	    const point& p0,const point& p1,const point& p2,
        // ÁîÔÐ ×ÐáãÝÕÜ àÕ×ãÛìâÐâ
        aSegments2D & aSegments,
        // ÁîÔÐ ×ÐáãÝÕÜ ÓàÐÝØæë ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
        aSegments2D & aBorder,
        // Here push numbers of bounds to get possibility 
        // select bounds in 3DToolKit program interactively (inserted by nikita) 
        // this parameter is NULL by default to avoid possible uncompatibility with other UNIVERS programs  
        IntAr* borderNumbers = NULL)const;

  // ¾ßàÕÔÕÛÕÝØÕ âÞÓÞ, çâÞ âÞçÚÐ ÝÐåÞÔØâáï ÒÝãâàØ ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
  Logic in_model_area(point p)const;
  // ·ÐßØáì
  //  virtual int write_to_DB(const dbItem &dbiWhere)const;
  virtual int write_to_DB(Progress* zp=NULL)const;
  void get_shots_data(aPoints& apSPs,IntAr& iaDT,IntAr& iaUT){
    pdbmTables->get_shots_data(apSPs,iaDT,iaUT);
  };
  void delete_modeling_data(const point& pSP);

  // ¿Þ ØÝÚÛØÝÞÜÕâàØØ - ÝÐÑÞà âÞçÕÚ ßÕàÕáÕçÕÝØï ÓàÐÝØæ ÜÞÔÕÛØ ßÞ íâÞÙ ØÝÚÛ.
  const aPoints& install_inclinometry(const aPoints& apInclIn);
  // ¸ÝÚÛØÝÞÜÕâàØï ×ÐÔÐÕâáï ßÐàÞÙ ÚÞÞàÔØÝÐâ, ÒÕàâØÚÐÛìÝÐï ßÐÛÚÐ
  const aPoints& install_vertical_inclinometry(double X, double Y);

  // ¿ÞÛãçØâì ÔÐÝÝëÕ Þ ÜÞÔÕÛØ ßÞ ^ ÚàØÒÞÙ
  void get_vp(DoubleAr &)const;
  void get_kp(DoubleAr &)const;
  void get_vs(DoubleAr &)const;
  void get_ks(DoubleAr &)const;

  void set_vp(const DoubleAr &);
  void set_kp(const DoubleAr &);
  void set_vs(const DoubleAr &);
  void set_ks(const DoubleAr &);

  Logic get_parallel_flag()const {return parallel_structure_flag;};
  void set_parallel_flag(Logic Flag){parallel_structure_flag=Flag;};
  parallel_model_params* get_parallel_model_params() {return pPMP;};
  
  void fix_dbStruct(dbNameStruct& dbnsIn)const{pdbvTables->fix_dbStruct(dbnsIn);};
  void prepare_tables(point p ,int i1,int i2){
    pdbmTables->prepare_tables(p,i1,i2,this);
  };
  p_pl_spline get_point_direction(const point3d &pWhere)const;
  
  /* New improvements by Garry */
  Logic check_structure(); /* Checking of model structure for possibility to
                  raytracing as in parallel structured model */
  Logic get_layers_inclination(double& a, double& b, double& c);
  Logic set_params_for_parallel();
  Logic analyse_gp(int LeadingWaveType);

  virtual const char* class_id(){return "pl_model";};
  virtual Logic isDerivedFrom(const char* class_name);

};

#endif /* p_model.H */

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

  virtual pl_intern_border* new_border(pl_intern_spline*,border_flag);
  virtual pl_intern_spline* new_spline(double,double,double,double);
  virtual pl_intern_body* new_body(geophys_params&);
  void add_breakpoints (const aPoints&apInclIn,int iI,int iCBody,int iTB);
  Logic cross_border (const point &p1,const point & p2,point&pRes,int iB);
  void incl_pass_body (const point &p1,const point & p2,int iC,int&iN);

  ErrCode readModelVolumeSegs(dbHandle);//nikita 25.04.05
  dbNameStruct dbsModGlass;
  int write_bodies()const;
  // ¿ãáâÞÙ, çâÞÑë Ø× ÒÕàåÝØå ãàÞÒÝÕÙ ÝÕ Òë×ëÒÐâì ÝØçÕÓÞ  pl_model(){
    pdbmTables=NULL;
    piaInclinometry=NULL;
    papInclinometry=NULL;
    pdbvTables=NULL;
    pPMP=NULL;
  };

public:
  pl_model(const dbItem &dbiWhere);
  pl_model(const dbNameStruct &dbh1D, const dbItem &dbiBorders,
       char* sC,IntAr *iaSequence=NULL);
  pl_model(pl_model& am);   // ºÞÝáâàãÚâÞà ÚÞßØàÞÒÐÝØï... ½Õ ÑëÒÐÕâ
  virtual ~pl_model();

  unsigned get_model_bounds_num();
  Logic get_model_bound(unsigned n, aPoints* points);

  // ¾âÞÑàÐÖÕÝØÕ íÛÕÜÕÝâÞÒ ÜÞÔÕÛØ
  int get_border(int i, pl_border &bord) const;
  // Ø× âÕÛÐ ßÞÔ ÝÞÜÕàÞÜ n1 ÓàÐÝØæã á ÝÞÜÕàÞÜ n2
  int get_border(int n1,int n2, pl_border &pb) const;
  int get_body(int i, pl_body &bod) const;
  // ÂÕÛÞ ßÞ áÞÔÕàÖÐéÕÙáï Ò ÝÕÜ âÞçÚÕ
  int get_body(const point &pWhere, pl_body &pb) const;
  int get_spline(int i, pl_spline &spl) const;

  // ½ÞÜÕà ÓàÐÝØæë Ò ÜÞÔÕÛØ ßÞ ÝÞÜÕàã ÓàÐÝØæë Ò âÕÛÕ
  int get_border_num(int ibod,int iborloc,int& iborglob) const;

  // ÄãÝÚæØï ÞâÞÑàÐÖÕÝØï ßÕàÕáÕçÕÝØï ßÛÞáÚÞáâÕÙ ÜÞÔÕÛØ á ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâìî
  int slice(const pl_spline& sSlice,         
        // »ÞÚÐÛìÝÐï áØáâÕÜÐ ÚÞÞàÔØÝÐâ
	    const point& p0,const point& p1,const point& p2,
        // ÁîÔÐ ×ÐáãÝÕÜ àÕ×ãÛìâÐâ
        aSegments2D & aSegments,
        // ÁîÔÐ ×ÐáãÝÕÜ ÓàÐÝØæë ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
        aSegments2D & aBorder,
        // Here push numbers of bounds to get possibility 
        // select bounds in 3DToolKit program interactively (inserted by nikita) 
        // this parameter is NULL by default to avoid possible uncompatibility with other UNIVERS programs  
        IntAr* borderNumbers = NULL)const;

  // ¾ßàÕÔÕÛÕÝØÕ âÞÓÞ, çâÞ âÞçÚÐ ÝÐåÞÔØâáï ÒÝãâàØ ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
  Logic in_model_area(point p)const;
  // ·ÐßØáì
  //  virtual int write_to_DB(const dbItem &dbiWhere)const;
  virtual int write_to_DB(Progress* zp=NULL)const;
  void get_shots_data(aPoints& apSPs,IntAr& iaDT,IntAr& iaUT){
    pdbmTables->get_shots_data(apSPs,iaDT,iaUT);
  };
  void delete_modeling_data(const point& pSP);

  // ¿Þ ØÝÚÛØÝÞÜÕâàØØ - ÝÐÑÞà âÞçÕÚ ßÕàÕáÕçÕÝØï ÓàÐÝØæ ÜÞÔÕÛØ ßÞ íâÞÙ ØÝÚÛ.
  const aPoints& install_inclinometry(const aPoints& apInclIn);
  // ¸ÝÚÛØÝÞÜÕâàØï ×ÐÔÐÕâáï ßÐàÞÙ ÚÞÞàÔØÝÐâ, ÒÕàâØÚÐÛìÝÐï ßÐÛÚÐ
  const aPoints& install_vertical_inclinometry(double X, double Y);

  // ¿ÞÛãçØâì ÔÐÝÝëÕ Þ ÜÞÔÕÛØ ßÞ ^ ÚàØÒÞÙ
  void get_vp(DoubleAr &)const;
  void get_kp(DoubleAr &)const;
  void get_vs(DoubleAr &)const;
  void get_ks(DoubleAr &)const;

  void set_vp(const DoubleAr &);
  void set_kp(const DoubleAr &);
  void set_vs(const DoubleAr &);
  void set_ks(const DoubleAr &);

  Logic get_parallel_flag()const {return parallel_structure_flag;};
  void set_parallel_flag(Logic Flag){parallel_structure_flag=Flag;};
  parallel_model_params* get_parallel_model_params() {return pPMP;};
  
  void fix_dbStruct(dbNameStruct& dbnsIn)const{pdbvTables->fix_dbStruct(dbnsIn);};
  void prepare_tables(point p ,int i1,int i2){
    pdbmTables->prepare_tables(p,i1,i2,this);
  };
  p_pl_spline get_point_direction(const point3d &pWhere)const;
  
  /* New improvements by Garry */
  Logic check_structure(); /* Checking of model structure for possibility to
                  raytracing as in parallel structured model */
  Logic get_layers_inclination(double& a, double& b, double& c);
  Logic set_params_for_parallel();
  Logic analyse_gp(int LeadingWaveType);

  virtual const char* class_id(){return "pl_model";};
  virtual Logic isDerivedFrom(const char* class_name);

};

#endif /* p_model.H */

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

  pl_model(){
    pdbmTables=NULL;
    piaInclinometry=NULL;
    papInclinometry=NULL;
    pdbvTables=NULL;
    pPMP=NULL;
  };

public:
  pl_model(const dbItem &dbiWhere);
  pl_model(const dbNameStruct &dbh1D, const dbItem &dbiBorders,
       char* sC,IntAr *iaSequence=NULL);
  pl_model(pl_model& am);   // ºÞÝáâàãÚâÞà ÚÞßØàÞÒÐÝØï... ½Õ ÑëÒÐÕâ
  virtual ~pl_model();

  unsigned get_model_bounds_num();
  Logic get_model_bound(unsigned n, aPoints* points);

  // ¾âÞÑàÐÖÕÝØÕ íÛÕÜÕÝâÞÒ ÜÞÔÕÛØ
  int get_border(int i, pl_border &bord) const;
  // Ø× âÕÛÐ ßÞÔ ÝÞÜÕàÞÜ n1 ÓàÐÝØæã á ÝÞÜÕàÞÜ n2
  int get_border(int n1,int n2, pl_border &pb) const;
  int get_body(int i, pl_body &bod) const;
  // ÂÕÛÞ ßÞ áÞÔÕàÖÐéÕÙáï Ò ÝÕÜ âÞçÚÕ
  int get_body(const point &pWhere, pl_body &pb) const;
  int get_spline(int i, pl_spline &spl) const;

  // ½ÞÜÕà ÓàÐÝØæë Ò ÜÞÔÕÛØ ßÞ ÝÞÜÕàã ÓàÐÝØæë Ò âÕÛÕ  int get_border_num(int ibod,int iborloc,int& iborglob) const;

  // ÄãÝÚæØï ÞâÞÑàÐÖÕÝØï ßÕàÕáÕçÕÝØï ßÛÞáÚÞáâÕÙ ÜÞÔÕÛØ á ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâìî
  int slice(const pl_spline& sSlice,         
        // »ÞÚÐÛìÝÐï áØáâÕÜÐ ÚÞÞàÔØÝÐâ
	    const point& p0,const point& p1,const point& p2,
        // ÁîÔÐ ×ÐáãÝÕÜ àÕ×ãÛìâÐâ
        aSegments2D & aSegments,
        // ÁîÔÐ ×ÐáãÝÕÜ ÓàÐÝØæë ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
        aSegments2D & aBorder,
        // Here push numbers of bounds to get possibility 
        // select bounds in 3DToolKit program interactively (inserted by nikita) 
        // this parameter is NULL by default to avoid possible uncompatibility with other UNIVERS programs  
        IntAr* borderNumbers = NULL)const;

  // ¾ßàÕÔÕÛÕÝØÕ âÞÓÞ, çâÞ âÞçÚÐ ÝÐåÞÔØâáï ÒÝãâàØ ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
  Logic in_model_area(point p)const;
  // ·ÐßØáì
  //  virtual int write_to_DB(const dbItem &dbiWhere)const;
  virtual int write_to_DB(Progress* zp=NULL)const;
  void get_shots_data(aPoints& apSPs,IntAr& iaDT,IntAr& iaUT){
    pdbmTables->get_shots_data(apSPs,iaDT,iaUT);
  };
  void delete_modeling_data(const point& pSP);

  // ¿Þ ØÝÚÛØÝÞÜÕâàØØ - ÝÐÑÞà âÞçÕÚ ßÕàÕáÕçÕÝØï ÓàÐÝØæ ÜÞÔÕÛØ ßÞ íâÞÙ ØÝÚÛ.
  const aPoints& install_inclinometry(const aPoints& apInclIn);
  // ¸ÝÚÛØÝÞÜÕâàØï ×ÐÔÐÕâáï ßÐàÞÙ ÚÞÞàÔØÝÐâ, ÒÕàâØÚÐÛìÝÐï ßÐÛÚÐ
  const aPoints& install_vertical_inclinometry(double X, double Y);

  // ¿ÞÛãçØâì ÔÐÝÝëÕ Þ ÜÞÔÕÛØ ßÞ ^ ÚàØÒÞÙ
  void get_vp(DoubleAr &)const;
  void get_kp(DoubleAr &)const;
  void get_vs(DoubleAr &)const;
  void get_ks(DoubleAr &)const;

  void set_vp(const DoubleAr &);
  void set_kp(const DoubleAr &);
  void set_vs(const DoubleAr &);
  void set_ks(const DoubleAr &);

  Logic get_parallel_flag()const {return parallel_structure_flag;};
  void set_parallel_flag(Logic Flag){parallel_structure_flag=Flag;};
  parallel_model_params* get_parallel_model_params() {return pPMP;};
  
  void fix_dbStruct(dbNameStruct& dbnsIn)const{pdbvTables->fix_dbStruct(dbnsIn);};
  void prepare_tables(point p ,int i1,int i2){
    pdbmTables->prepare_tables(p,i1,i2,this);
  };
  p_pl_spline get_point_direction(const point3d &pWhere)const;
  
  /* New improvements by Garry */
  Logic check_structure(); /* Checking of model structure for possibility to
                  raytracing as in parallel structured model */
  Logic get_layers_inclination(double& a, double& b, double& c);
  Logic set_params_for_parallel();
  Logic analyse_gp(int LeadingWaveType);

  virtual const char* class_id(){return "pl_model";};
  virtual Logic isDerivedFrom(const char* class_name);

};

#endif /* p_model.H */

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

  int get_border_num(int ibod,int iborloc,int& iborglob) const;

  // ÄãÝÚæØï ÞâÞÑàÐÖÕÝØï ßÕàÕáÕçÕÝØï ßÛÞáÚÞáâÕÙ ÜÞÔÕÛØ á ×ÐÔÐÝÝÞÙ ßÛÞáÚÞáâìî  int slice(const pl_spline& sSlice,         
        // »ÞÚÐÛìÝÐï áØáâÕÜÐ ÚÞÞàÔØÝÐâ
	    const point& p0,const point& p1,const point& p2,
        // ÁîÔÐ ×ÐáãÝÕÜ àÕ×ãÛìâÐâ
        aSegments2D & aSegments,
        // ÁîÔÐ ×ÐáãÝÕÜ ÓàÐÝØæë ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
        aSegments2D & aBorder,
        // Here push numbers of bounds to get possibility 
        // select bounds in 3DToolKit program interactively (inserted by nikita) 
        // this parameter is NULL by default to avoid possible uncompatibility with other UNIVERS programs  
        IntAr* borderNumbers = NULL)const;

  // ¾ßàÕÔÕÛÕÝØÕ âÞÓÞ, çâÞ âÞçÚÐ ÝÐåÞÔØâáï ÒÝãâàØ ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
  Logic in_model_area(point p)const;
  // ·ÐßØáì
  //  virtual int write_to_DB(const dbItem &dbiWhere)const;
  virtual int write_to_DB(Progress* zp=NULL)const;
  void get_shots_data(aPoints& apSPs,IntAr& iaDT,IntAr& iaUT){
    pdbmTables->get_shots_data(apSPs,iaDT,iaUT);
  };
  void delete_modeling_data(const point& pSP);

  // ¿Þ ØÝÚÛØÝÞÜÕâàØØ - ÝÐÑÞà âÞçÕÚ ßÕàÕáÕçÕÝØï ÓàÐÝØæ ÜÞÔÕÛØ ßÞ íâÞÙ ØÝÚÛ.
  const aPoints& install_inclinometry(const aPoints& apInclIn);
  // ¸ÝÚÛØÝÞÜÕâàØï ×ÐÔÐÕâáï ßÐàÞÙ ÚÞÞàÔØÝÐâ, ÒÕàâØÚÐÛìÝÐï ßÐÛÚÐ
  const aPoints& install_vertical_inclinometry(double X, double Y);

  // ¿ÞÛãçØâì ÔÐÝÝëÕ Þ ÜÞÔÕÛØ ßÞ ^ ÚàØÒÞÙ
  void get_vp(DoubleAr &)const;
  void get_kp(DoubleAr &)const;
  void get_vs(DoubleAr &)const;
  void get_ks(DoubleAr &)const;

  void set_vp(const DoubleAr &);
  void set_kp(const DoubleAr &);
  void set_vs(const DoubleAr &);
  void set_ks(const DoubleAr &);

  Logic get_parallel_flag()const {return parallel_structure_flag;};
  void set_parallel_flag(Logic Flag){parallel_structure_flag=Flag;};
  parallel_model_params* get_parallel_model_params() {return pPMP;};
  
  void fix_dbStruct(dbNameStruct& dbnsIn)const{pdbvTables->fix_dbStruct(dbnsIn);};
  void prepare_tables(point p ,int i1,int i2){
    pdbmTables->prepare_tables(p,i1,i2,this);
  };
  p_pl_spline get_point_direction(const point3d &pWhere)const;
  
  /* New improvements by Garry */
  Logic check_structure(); /* Checking of model structure for possibility to
                  raytracing as in parallel structured model */
  Logic get_layers_inclination(double& a, double& b, double& c);
  Logic set_params_for_parallel();
  Logic analyse_gp(int LeadingWaveType);

  virtual const char* class_id(){return "pl_model";};
  virtual Logic isDerivedFrom(const char* class_name);

};

#endif /* p_model.H */

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

  int slice(const pl_spline& sSlice,            
        // »ÞÚÐÛìÝÐï áØáâÕÜÐ ÚÞÞàÔØÝÐâ	    const point& p0,const point& p1,const point& p2,
        // ÁîÔÐ ×ÐáãÝÕÜ àÕ×ãÛìâÐâ
        aSegments2D & aSegments,
        // ÁîÔÐ ×ÐáãÝÕÜ ÓàÐÝØæë ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
        aSegments2D & aBorder,
        // Here push numbers of bounds to get possibility 
        // select bounds in 3DToolKit program interactively (inserted by nikita) 
        // this parameter is NULL by default to avoid possible uncompatibility with other UNIVERS programs  
        IntAr* borderNumbers = NULL)const;

  // ¾ßàÕÔÕÛÕÝØÕ âÞÓÞ, çâÞ âÞçÚÐ ÝÐåÞÔØâáï ÒÝãâàØ ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
  Logic in_model_area(point p)const;
  // ·ÐßØáì
  //  virtual int write_to_DB(const dbItem &dbiWhere)const;
  virtual int write_to_DB(Progress* zp=NULL)const;
  void get_shots_data(aPoints& apSPs,IntAr& iaDT,IntAr& iaUT){
    pdbmTables->get_shots_data(apSPs,iaDT,iaUT);
  };
  void delete_modeling_data(const point& pSP);

  // ¿Þ ØÝÚÛØÝÞÜÕâàØØ - ÝÐÑÞà âÞçÕÚ ßÕàÕáÕçÕÝØï ÓàÐÝØæ ÜÞÔÕÛØ ßÞ íâÞÙ ØÝÚÛ.
  const aPoints& install_inclinometry(const aPoints& apInclIn);
  // ¸ÝÚÛØÝÞÜÕâàØï ×ÐÔÐÕâáï ßÐàÞÙ ÚÞÞàÔØÝÐâ, ÒÕàâØÚÐÛìÝÐï ßÐÛÚÐ
  const aPoints& install_vertical_inclinometry(double X, double Y);

  // ¿ÞÛãçØâì ÔÐÝÝëÕ Þ ÜÞÔÕÛØ ßÞ ^ ÚàØÒÞÙ
  void get_vp(DoubleAr &)const;
  void get_kp(DoubleAr &)const;
  void get_vs(DoubleAr &)const;
  void get_ks(DoubleAr &)const;

  void set_vp(const DoubleAr &);
  void set_kp(const DoubleAr &);
  void set_vs(const DoubleAr &);
  void set_ks(const DoubleAr &);

  Logic get_parallel_flag()const {return parallel_structure_flag;};
  void set_parallel_flag(Logic Flag){parallel_structure_flag=Flag;};
  parallel_model_params* get_parallel_model_params() {return pPMP;};
  
  void fix_dbStruct(dbNameStruct& dbnsIn)const{pdbvTables->fix_dbStruct(dbnsIn);};
  void prepare_tables(point p ,int i1,int i2){
    pdbmTables->prepare_tables(p,i1,i2,this);
  };
  p_pl_spline get_point_direction(const point3d &pWhere)const;
  
  /* New improvements by Garry */
  Logic check_structure(); /* Checking of model structure for possibility to
                  raytracing as in parallel structured model */
  Logic get_layers_inclination(double& a, double& b, double& c);
  Logic set_params_for_parallel();
  Logic analyse_gp(int LeadingWaveType);

  virtual const char* class_id(){return "pl_model";};
  virtual Logic isDerivedFrom(const char* class_name);

};

#endif /* p_model.H */

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

        const point& p0,const point& p1,const point& p2,
        // ÁîÔÐ ×ÐáãÝÕÜ àÕ×ãÛìâÐâ
        aSegments2D & aSegments,
        // ÁîÔÐ ×ÐáãÝÕÜ ÓàÐÝØæë ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
        aSegments2D & aBorder,
        // Here push numbers of bounds to get possibility 
        // select bounds in 3DToolKit program interactively (inserted by nikita) 
        // this parameter is NULL by default to avoid possible uncompatibility with other UNIVERS programs  
        IntAr* borderNumbers = NULL)const;

  // ¾ßàÕÔÕÛÕÝØÕ âÞÓÞ, çâÞ âÞçÚÐ ÝÐåÞÔØâáï ÒÝãâàØ ÞÑÛÐáâØ ÜÞÔÕÛØàÞÒÐÝØï
  Logic in_model_area(point p)const;
  // ·ÐßØáì  //  virtual int write_to_DB(const dbItem &dbiWhere)const;
  virtual int write_to_DB(Progress* zp=NULL)const;
  void get_shots_data(aPoints& apSPs,IntAr& iaDT,IntAr& iaUT){
    pdbmTables->get_shots_data(apSPs,iaDT,iaUT);
  };
  void delete_modeling_data(const point& pSP);

  // ¿Þ ØÝÚÛØÝÞÜÕâàØØ - ÝÐÑÞà âÞçÕÚ ßÕàÕáÕçÕÝØï ÓàÐÝØæ ÜÞÔÕÛØ ßÞ íâÞÙ ØÝÚÛ.
  const aPoints& install_inclinometry(const aPoints& apInclIn);
  // ¸ÝÚÛØÝÞÜÕâàØï ×ÐÔÐÕâáï ßÐàÞÙ ÚÞÞàÔØÝÐâ, ÒÕàâØÚÐÛìÝÐï ßÐÛÚÐ
  const aPoints& install_vertical_inclinometry(double X, double Y);

  // ¿ÞÛãçØâì ÔÐÝÝëÕ Þ ÜÞÔÕÛØ ßÞ ^ ÚàØÒÞÙ
  void get_vp(DoubleAr &)const;
  void get_kp(DoubleAr &)const;
  void get_vs(DoubleAr &)const;
  void get_ks(DoubleAr &)const;

  void set_vp(const DoubleAr &);
  void set_kp(const DoubleAr &);
  void set_vs(const DoubleAr &);
  void set_ks(const DoubleAr &);

  Logic get_parallel_flag()const {return parallel_structure_flag;};
  void set_parallel_flag(Logic Flag){parallel_structure_flag=Flag;};
  parallel_model_params* get_parallel_model_params() {return pPMP;};
  
  void fix_dbStruct(dbNameStruct& dbnsIn)const{pdbvTables->fix_dbStruct(dbnsIn);};
  void prepare_tables(point p ,int i1,int i2){
    pdbmTables->prepare_tables(p,i1,i2,this);
  };
  p_pl_spline get_point_direction(const point3d &pWhere)const;
  
  /* New improvements by Garry */
  Logic check_structure(); /* Checking of model structure for possibility to
                  raytracing as in parallel structured model */
  Logic get_layers_inclination(double& a, double& b, double& c);
  Logic set_params_for_parallel();
  Logic analyse_gp(int LeadingWaveType);

  virtual const char* class_id(){return "pl_model";};
  virtual Logic isDerivedFrom(const char* class_name);

};

#endif /* p_model.H */

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

  //  virtual int write_to_DB(const dbItem &dbiWhere)const;
  virtual int write_to_DB(Progress* zp=NULL)const;
  void get_shots_data(aPoints& apSPs,IntAr& iaDT,IntAr& iaUT){
    pdbmTables->get_shots_data(apSPs,iaDT,iaUT);
  };
  void delete_modeling_data(const point& pSP);

  // ¿Þ ØÝÚÛØÝÞÜÕâàØØ - ÝÐÑÞà âÞçÕÚ ßÕàÕáÕçÕÝØï ÓàÐÝØæ ÜÞÔÕÛØ ßÞ íâÞÙ ØÝÚÛ.
  const aPoints& install_inclinometry(const aPoints& apInclIn);
  // ¸ÝÚÛØÝÞÜÕâàØï ×ÐÔÐÕâáï ßÐàÞÙ ÚÞÞàÔØÝÐâ, ÒÕàâØÚÐÛìÝÐï ßÐÛÚÐ
  const aPoints& install_vertical_inclinometry(double X, double Y);

  // ¿ÞÛãçØâì ÔÐÝÝëÕ Þ ÜÞÔÕÛØ ßÞ ^ ÚàØÒÞÙ
  void get_vp(DoubleAr &)const;
  void get_kp(DoubleAr &)const;
  void get_vs(DoubleAr &)const;
  void get_ks(DoubleAr &)const;

  void set_vp(const DoubleAr &);
  void set_kp(const DoubleAr &);
  void set_vs(const DoubleAr &);
  void set_ks(const DoubleAr &);

  Logic get_parallel_flag()const {return parallel_structure_flag;};
  void set_parallel_flag(Logic Flag){parallel_structure_flag=Flag;};
  parallel_model_params* get_parallel_model_params() {return pPMP;};
  
  void fix_dbStruct(dbNameStruct& dbnsIn)const{pdbvTables->fix_dbStruct(dbnsIn);};
  void prepare_tables(point p ,int i1,int i2){
    pdbmTables->prepare_tables(p,i1,i2,this);
  };
  p_pl_spline get_point_direction(const point3d &pWhere)const;
  
  /* New improvements by Garry */
  Logic check_structure(); /* Checking of model structure for possibility to
                  raytracing as in parallel structured model */
  Logic get_layers_inclination(double& a, double& b, double& c);
  Logic set_params_for_parallel();
  Logic analyse_gp(int LeadingWaveType);

  virtual const char* class_id(){return "pl_model";};
  virtual Logic isDerivedFrom(const char* class_name);

};

#endif /* p_model.H */

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