GyotoDeformedTorus.h

Go to the documentation of this file.

/*
 *   Copyright (c) 2013, 2018 Frederic Vincent & Thibaut Paumard
 *
 *
 */

#ifndef __DeformedTorus_h
#define __DeformedTorus_h

#include <GyotoAstrobj.h>
#include <GyotoStandardAstrobj.h>
#include <GyotoKerrBL.h>


namespace Gyoto {
  namespace Astrobj {
    class DeformedTorus;
  };
  class FactoryMessenger;
  namespace Spectrum {
    class Generic;
  }
};

class Gyoto::Astrobj::DeformedTorus
: public Gyoto::Astrobj::Standard {
  friend class Gyoto::SmartPointer<Gyoto::Astrobj::DeformedTorus>;

 private:
  SmartPointer<Gyoto::Metric::KerrBL> gg_;
  SmartPointer<Spectrum::Generic> spectrum_;
  double c_;
  unsigned long mode_;
  double param_beta_;
  double param_beta_st_;
  double param_eta_;
  enum perturb_t {RadialTranslation=1,
                  VerticalTranslation=2,
                  Rotation=3,
                  Expansion=4,
                  RadialShear=5,
                  VerticalShear=6,
                  PureShear=7};
  perturb_t perturb_kind_;
 public:
  GYOTO_OBJECT;
  DeformedTorus();
  DeformedTorus(const DeformedTorus &o);
  virtual ~DeformedTorus();
  virtual DeformedTorus * clone() const ;

  // Standard accessors
  GYOTO_OBJECT_ACCESSORS(SmartPointer<Spectrum::Generic>, spectrum);
  GYOTO_OBJECT_ACCESSORS(double, largeRadius);
  GYOTO_OBJECT_ACCESSORS(double, beta);
  GYOTO_OBJECT_ACCESSORS(double, betaSt);
  GYOTO_OBJECT_ACCESSORS(double, eta);
  GYOTO_OBJECT_ACCESSORS(unsigned long, mode);
  GYOTO_OBJECT_ACCESSORS_STRING(perturbKind);

  using Generic::metric;
  virtual void metric(Gyoto::SmartPointer<Gyoto::Metric::Generic>);
  virtual double operator()(double const coord[4]) ;
  virtual void getVelocity(double const pos[4], double vel[4]) ;
  /*virtual int Impact(Gyoto::Photon* ph, size_t index,
    Astrobj::Properties *data=NULL);*/
  double emission(double nuem,double,state_t const &,double const *) const;
#endif
};