GyotoXillverReflection.h

Go to the documentation of this file.

 
/*
  Copyright (c) 2017, 2018 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 __GyotoXillverReflection_H_
#define __GyotoXillverReflection_H_
 
#include <GyotoThinDisk.h>
 
 
namespace Gyoto{
  namespace Astrobj { class XillverReflection; }
}
 
class Gyoto::Astrobj::XillverReflection
: public Astrobj::ThinDisk,
  public Hook::Listener
{
  friend class Gyoto::SmartPointer<Gyoto::Astrobj::XillverReflection>;
 
 private:
  std::string filenameIllum_; 
  std::string filenameRefl_; 
  double * reflection_; 
 
  double * logxi_; 
  double * incl_; 
  double * freq_; 
  size_t nnu_; 
  size_t ni_; 
  size_t nxi_; 
 
  double * illumination_;
 
  double * radius_; 
  double * phi_; 
  size_t nr_; 
  size_t nphi_; 
 
  double aa_; 
  double lampradius_; 
  double timelampphizero_; 
 
  bool average_over_angle_; 
 
 protected:
 
 public:
  GYOTO_OBJECT;
  // fillProperty is overridden to remove leading "!" from FITS filename
  void fillProperty(Gyoto::FactoryMessenger *fmp, Property const &p) const;
  XillverReflection(); 
  XillverReflection(const XillverReflection& o); 
  virtual XillverReflection * clone() const ; 
  virtual ~XillverReflection() ; 
 
 public:
 
  void timelampphizero(double tt);
  double timelampphizero() const;
  void lampradius(double rr);
  double lampradius() const;
  void fileillumination(std::string const &f);
  std::string fileillumination() const ;
  void filereflection(std::string const &f);
  std::string filereflection() const ;
  
  void averageOverAngle(bool t);
  bool averageOverAngle()const;
 
  #ifdef GYOTO_USE_CFITSIO
  virtual void fitsReadIllum(std::string filename);
  virtual void fitsWriteIllum(std::string filename);
 
  virtual void fitsReadRefl(std::string filename);
  virtual void fitsWriteRefl(std::string filename);
#endif
 
  void setIllumination(double * pattern);
  void setReflection(double * pattern);
  
  
  virtual void copyIllumination(double const * const pattern = NULL,
                                size_t const naxes[2] = NULL);
  virtual double const * getIllumination() const;
  virtual void getIlluminationNaxes( size_t naxes[2] ) const ; 
 
  virtual void copyReflection(double const * const pattern = NULL,
                              size_t const naxes[3] = NULL);
  virtual double const * getReflection() const;
  virtual void getReflectionNaxes( size_t naxes[3] ) const ; 
 
  virtual void copyGridReflLogxi(double const * const pattern = NULL,
                                 size_t nxi = 0 );
  virtual double const * getGridReflLogxi() const; 
  virtual void copyGridReflIncl(double const * const pattern = NULL,
                                size_t ni = 0 );
  virtual double const * getGridReflIncl() const; 
  virtual void copyGridReflFreq(double const * const pattern = NULL,
                                size_t nnu = 0 );
  virtual double const * getGridReflFreq() const; 
  virtual void copyGridIllumRadius(double const * const pattern = NULL,
                                   size_t nr = 0 );
  virtual double const * getGridIllumRadius() const; 
  virtual void copyGridIllumPhi(double const * const pattern = NULL,
                                size_t nphi = 0 );
  virtual double const * getGridIllumPhi() const; 
  
 protected:
  
  void getIndicesRefl(size_t i[3], double const co[4], double logxi, double incl,
                      double nu=0.) const ;
  void getIndicesIllum(size_t i[3], double const co[4]) const ;
 
 public:
  
  virtual double emission(double nu_em, double dsem,
                          state_t const &_ph, double const _obj[8]=NULL) const;
 
  virtual void updateSpin() ;
  virtual void tell(Gyoto::Hook::Teller *msg);
  using Generic::metric;
  virtual void metric(SmartPointer<Metric::Generic>);
 
};
 
#endif