GyotoThinDiskGridIntensity.h

Go to the documentation of this file.

/*
    Copyright 2019-2021 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 __GyotoThinDiskGridIntensity_H_ 
#define __GyotoThinDiskGridIntensity_H_ 

#include <iostream>
#include <fstream>
#include <iomanip>

namespace Gyoto{
  namespace Astrobj { class ThinDiskGridIntensity; }
  class GridData2D;
}

#include <GyotoThinDisk.h>
#include <GyotoGridData2D.h>

class Gyoto::Astrobj::ThinDiskGridIntensity
: public Astrobj::ThinDisk,
  public GridData2D,
  public Hook::Listener
{
  friend class Gyoto::SmartPointer<Gyoto::Astrobj::ThinDiskGridIntensity>;
 private:
  std::string filename_; 

  double * intensity_; 
  double * time_array_; 
  double deltat_;

 public:
  GYOTO_OBJECT;
  GYOTO_OBJECT_THREAD_SAFETY;
  
  // Constructors - Destructor
  // -------------------------
  ThinDiskGridIntensity(); 
  
  ThinDiskGridIntensity(const ThinDiskGridIntensity& ) ;
  virtual ThinDiskGridIntensity* clone () const; 
  
  virtual ~ThinDiskGridIntensity() ;                        
  
  // Accessors
  // ---------
 public:
  void file(std::string const &f) ;
  std::string file() const;
  /*
    timeTranslation shifts the value of GridData2D::tmin_ and tmax_,
    allowing to scan the full simulation without having to change
    the value of the Screen observation time (which is typically
    not provided in M unit in the XML). 
    Choosing a negative timeTranslation, i.e. performing tmin_,tmax_-=dt, 
    amounts to increasing the Screen observation time by the same value, 
    tobs+=dt.

   */
  void timeTranslation_inMunit(double const dt) ;
  double timeTranslation_inMunit() const ;
  //void dt(double dd);
  void copyIntensity(double const *const intensity,
                     size_t const naxes[3]);
  double const * getIntensity() const;
  void copyTimeArray(double const *const time_array, size_t const ntimes);
  double const * getTimeArray() const;
 public:
  using Generic::metric;
  std::vector<size_t> fitsRead(std::string filename) ;
  virtual double emission(double nu_em, double dsem,
                          state_t const &c_ph,double const c_obj[8]=NULL) const;
  virtual void getVelocity(double const pos[4], double vel[4]) ;



};

#endif