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