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,
                         const std::string & prefix = "");
#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