u3d_body.hpp

/* u3d_body.hpp */
/* $Id: u3d_body.hpp,v 1.1 2008/05/13 12:19:18 sad Exp $ */
#ifndef __u3d_body_hpp
#define __u3d_body_hpp

#include <mth/u3d_point.hpp>
//#include <u3d_defines.hpp>
#include "u3d_model_defines.hpp"


//============================================================================
class U3dBody
{
public:
  
  U3dBody();
  
  U3dBody(U3dSignAr signatures, double vp, double vs, 
      double kzp, double kzs, double density, double an_kp=1.0, double an_ks=1.0, double R=1.0);

  U3dBody(U3dSignAr signatures, U3dBodyVels vels, 
      double kzp, double kzs, double density, double an_kp=1.0, double an_ks=1.0, double R=1.0);

  U3dBody(U3dBody*);
  
  ~U3dBody();
  
  double getCPV() const;
  double getCSV() const;
  
  double getPV(U3dPoint) const;
  double getSV(U3dPoint) const;
  
  double getPKZ() const;
  double getSKZ() const;
  
  double getDen() const ;

  double getAnKP() const {  return m_an_kp; }

  double getAnKS() const {  return m_an_ks; }

  double getRParam() const { return m_R; }
  
  U3dBodyVels getVels() const;
  
  U3dSignAr getSignature() const { return m_signatures; };
  
  bool compareSign(U3dSignAr pointSign) const;
  
  bool constParams() const;
  
  
  void setAnParams(double an_kp, double an_ks, double an_R);
  
  void setAnKP(double an_kp) { m_an_kp = an_kp; }
  
  void setAnKS(double an_ks) { m_an_ks = an_ks; }
  
  void setRParam(double R) { m_R = R; }
  
  
  void getAnParams(double& an_kp, double& an_ks, double& an_R);
    
  double getAnKP() { return m_an_kp; }
  
  double getAnKS() { return m_an_ks; }
  
  double getRParam() { return m_R; }
  
  
protected: 
  
  bool m_ready;            
  bool m_constParams;      
  double m_pV;             
  double m_sV;             
  double m_pKz;            
  double m_sKz;            
  double m_density;        
  double m_an_kp;          
  double m_an_ks;          
  double m_R;              
  U3dSignAr m_signatures;  
  U3dBodyVels m_vels;      
};

typedef std::vector<U3dBody*> U3dBodiesAr; 

#endif /* u3d_body.hpp */