GyotoKerrKS.h

Go to the documentation of this file.

 
/*
    Copyright 2011-2015 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 __GyotoKerrKS_H_
#define __GyotoKerrKS_H_ 
 
namespace Gyoto {
  namespace Metric { class KerrKS; }
}
 
#include <GyotoMetric.h>
#include <GyotoWorldline.h>
#ifdef GYOTO_USE_XERCES
#include <GyotoRegister.h>
#endif
 
class Gyoto::Metric::KerrKS
: public Metric::Generic
{
  friend class Gyoto::SmartPointer<Gyoto::Metric::KerrKS>;
  
  // Data : 
  // -----
 
 protected:
  double spin_ ;  
  double a2_;     
  double rsink_;  
  double drhor_;  
 
  // Constructors - Destructor
  // -------------------------
 public: 
  GYOTO_OBJECT;
  KerrKS(); 
  virtual KerrKS* clone () const;         
  
  // Mutators / assignment
  // ---------------------
 public:
  // default operator= is fine
  void spin(const double spin); 
 
  // Accessors
  // ---------
 public:
  double spin() const ; 
  void horizonSecurity(double drhor);
  double horizonSecurity() const;
  
  double gmunu(const double x[4],
                       int alpha, int beta) const ;
 
  void gmunu(double g[4][4], const double pos[4]) const;
 
  void gmunu_up(double gup[4][4], const double pos[4]) const;
 
  void jacobian(double dst[4][4][4], const double x[4]) const ;
 
  using Generic::christoffel;
  int christoffel(double dst[4][4][4], const double x[4]) const ;
  int christoffel(double dst[4][4][4], const double pos[4], double gup[4][4], double jac[4][4][4]) const ;
 
  virtual void circularVelocity(double const pos[4], double vel [4],
                                double dir=1.) const ;
 
  virtual int isStopCondition(double const * const coord) const;
 
  virtual int setParameter(std::string name,
                           std::string content,
                           std::string unit);
};
 
#endif