vel2dmod.hpp

// vel2dmod.hpp
// $Id: vel2dmod.hpp,v 1.25 2008/04/09 08:52:33 sasha Exp $
//-----------------------------------------------------------------------------
#ifndef __2dModelH
#define __2dModelH

#include "vel2dmod_global.hpp"
#include "grid2d.hpp"
#include "os.hpp"
#include "vel2dmod_service.hpp"


//==============================================================================
class Grid2d;



//==============================================================================
class Vel2dModel
{
 public:
  
  Vel2dModel();
  
  Vel2dModel(double *x, double *z, double *ang, double *vp, double *vs, 
         double *pkz, double *skz,double *den, 
         double *an_kp, double *an_ks, double *R, int n, 
         double xMin, double xMax, double zMin, double zMax, 
         double mod_x, double mod_y, double mod_az, int approx_type = 0, double grid_cell_size = RT_GRIDSIZE);
  
  Vel2dModel(BodyInfoAr& bodyInfo, double xMin, double xMax, double zMin, double zMax, 
         double mod_x, double mod_y, double mod_az, int approx_type = 0);
  
  Vel2dModel(Vel2dModel*);

  virtual ~Vel2dModel();
  
  int getApproxType() {
    return m_approx_type;
  }
  
  void     AddGridPtr(Grid2d *gridPtr);  

  void     CheckModelStatus();                  

  int      DefineBody(double x, double z);      

  int      ExactDefineBody(double x, double z); 
  
  int      ExactDefineBound(double x, double z, double eps=DEFINE_BOUND_EPS, bool cry=true);

  int      GetBoundNum(double x, double z);             //??????????

  double   GetVel(double x, double z, WAVE_TYPE vPS, int bodyNum=-1);
                                           
  double   GetVel(double x, double z, WAVE_TYPE vPS, double angle, int bodyNum=-1);
            
  void     GetVelGr(double x, double z, double &vk, double &vka, WAVE_TYPE wt);

  double   GetVelMax(WAVE_TYPE wt);

  double   AnyBound(int bNum, double x, bool &flag, double z=-1);
  
  bool     GetParams(double x, double z, double &vP, double &vS, double &den, double &anKP, double &anKS, double &R);
  bool     GetParams(double x, double z, double &vP, double &vS, double &den) {
    double  anKP, anKS, R;
    return GetParams(x, z, vP, vS, den, anKP, anKS, R);
  }
  bool     GetParams2(double x, double z, double &vPx, double &vPz, double &vSx, double &vSz, double &den, double &R);

  bool     BoundInt(int bNum, double x1, double z1, double x2, double z2, double &ix, double &iz,
            double &iAngle, double &curvity, bool ANY=false);
    
  bool     IsOuterBound(int bNum); 

  /*
  void     ReplaceCurve(int cNum, int num, double *xx, double *zz, double rc, double sp);
  void     ReplaceCurve(int cNum, double rc, double sp, double *xx, double *zz, 
                double *dx, double *dz, double *curv, bool *flag, int nn);
  */

  int      GetSegNum(int bn, int n);
  int      GetSegNum(int bn, double x, double z);
  
  bool     GetLegitimValues(int bn, bool st_end, double &zmin, double &zmax);

  void     RebuildModel(int bn, double *x, double *z, bool *flag, int num, double rc=-1, double sp=-1, bool *im_flag=NULL);
  
  void     GenerateGrid(double step=-1);  
  void     RegenerateGrid(double xmn, double zmn, double xmx, double zmx);
  double   GetGridCellSize();

  void     AddWellInfo(ObservSys *os);
  bool     GetMod1dFlag();

  bool     IsInsideOfModelArea(double x, double z);
  void     CleanOS();
  void     SetModelXYAZ(double m_x, double m_y, double m_az);
  void     GetModelXYAZ(double &m_x, double &m_y, double &m_az);
  
  virtual bool     RemoveBound(int boundN);
  
  virtual bool     AddBound(double *x, double *z, bool *flag, int p_num, double rc=-1, double sp=-1);
  
  bool     CombineBounds(Bound *top_b, Bound *bot_b);
  
  Bound*   SplitBounds(Bound *bound, double z);
  
  void     TestModelEntity();
  
  bool     checkInside(double x, double z);
  
  bool     belongPerimeter(double x, double z);
  
  BodiesAr Bodies;      
  BoundsAr Bounds;      
  char modelName[16];
  bool ready;            
  double xMin, xMax, zMin, zMax; 
  Grid2d *pGrid;  
  bool    grid_defined;
  
  virtual void    SetPBLimits(int n, double lmin, double lmax) 
  {n=0; lmin=0; lmax=0;}; 
  virtual void    SetPBComment(int n, const char* str) 
  {n=0; if (0x0 == str) return;};         
  virtual void    SetPBValue(int n, double val) 
  {n=0; val=0;};                
  virtual void    AddPBValue(int n, double val) 
  {n=0; val=0;};                
  double GetModelAz() { return modelAZ; };                         
  void   GetModelXY(double &x, double &y) { x=modelX; y=modelY; }; 
  bool LineInt(int bNum, double p, double coord1,double coord2,bool finiteFlag,
           bool xzFlag, double &result);

  bool ExportToTextFiles(const char *files_path);
  
protected:

  bool     PrepareModel();                       
  void     DefineBoundPos();                     
  bool TestLine(int bNum, double p, double coord1, double coord2, bool finiteFlag, bool xzFlag);

  void RebuildOuterBound(int bn, double x0, double z0, double x1, double z1, bool st_end);
  void BoundsSegmentation();

  bool    Model1dFlag; 
  
  int m_approx_type;
  
  double  modelAZ, modelX, modelY;
};


//-----------------------------------------------------------------------------
#endif