GyotoDisk3D.h

Go to the documentation of this file.

 
/*
    Copyright 2011 Frederic Vincent, Thibaut Paumard
 
    This file is part of Gyoto.
 
    Gyoto is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
 
    Gyoto is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.
 
    You should have received a copy of the GNU General Public License
    along with Gyoto.  If not, see <http://www.gnu.org/licenses/>.
 */
 
#ifndef __GyotoDisk3D_H_ 
#define __GyotoDisk3D_H_ 
 
#include <iostream>
#include <fstream>
#include <iomanip>
 
namespace Gyoto{
  namespace Astrobj { class Disk3D; }
}
 
class Gyoto::Astrobj::Disk3D : public Gyoto::Astrobj::Generic {
  friend class Gyoto::SmartPointer<Gyoto::Astrobj::Disk3D>;
 private:
  std::string filename_; 
  double * emissquant_; 
 
  double * opacity_; 
 
  double * velocity_; 
 
  double dnu_; 
  double nu0_; 
  size_t nnu_; 
 
  double dphi_; 
  double phimin_;
  size_t nphi_; 
  double phimax_; 
 
  size_t repeat_phi_; 
  //double phi0_==0, phi max is always 2*M_PI
 
  double dz_; 
  double zmin_; 
  size_t nz_; 
  double zmax_; 
 
  double dr_; 
  double rin_; 
  size_t nr_; 
  double rout_; 
 
  int zsym_; 
 
  double tPattern_; 
  double omegaPattern_; 
 
  // Constructors - Destructor
  // -------------------------
 public:
  GYOTO_OBJECT;
  // fillProperty is overridden to remove leading "!" from FITS filename
  void fillProperty(Gyoto::FactoryMessenger *fmp, Property const &p) const;
 
  Disk3D(); 
  
  Disk3D(const Disk3D& ) ;
  virtual Disk3D* clone () const; 
  
  virtual ~Disk3D() ;                        
  
  // Accessors
  // ---------
 public:
 
#ifdef GYOTO_USE_CFITSIO
  virtual void fitsRead(std::string filename_);
 
  virtual void fitsWrite(std::string filename_);
#endif
 
  void file(std::string const &f);
  std::string file() const;
  void zsym(bool t);
  bool zsym() const;
  void tPattern(double t);
  double tPattern() const;
  void omegaPattern(double t);
  double omegaPattern() const;
 
 
 
  void setEmissquant(double * pattern);
 
  void opacity(double * pattern);
 
 
  void setVelocity(double * pattern);
 
 
  virtual void copyEmissquant(double const * const pattern = NULL,
                              size_t const naxes[4] = NULL);
 
  virtual double const * getEmissquant() const;
 
  virtual void getEmissquantNaxes( size_t naxes[4] ) const ;
 
  virtual void copyOpacity(double const * const pattern = NULL,
                              size_t const naxes[4] = NULL);
 
  virtual double const * opacity() const;
 
 
  virtual void copyVelocity(double const * const pattern = NULL,
                              size_t const naxes[3] = NULL);
  virtual double const * getVelocity() const;
 
  virtual void repeatPhi(size_t n);
  virtual size_t repeatPhi() const;
 
  virtual void nu0(double freq);
  virtual double nu0() const;
 
  virtual void dnu(double dfreq);
  virtual double dnu() const;
 
  void rin(double rrin);
  double rin() const;
 
  void rout(double rout);
  double rout() const;
 
  void zmin(double zmin);
  double zmin() const;
 
  void zmax(double zmax);
  double zmax() const;
 
  void phimin(double phimin);
  double phimin() const;
 
  void phimax(double phimax);
  double phimax() const;
 
 protected:
  void getIndices(size_t i[4], double const co[4], double nu=0.) const ;
 
 public:
  int Impact(Photon *ph, size_t index, Astrobj::Properties *data);
 
 
  virtual void getVelocity(double const pos[4], double vel[4])  ;
 
};
 
#endif