GyotoBalasinGrumiller.h

/*
    Copyright 2025 Filipe Costa, Frédéric 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 __BalasinGrumillerheader_H_
#define __BalasinGrumillerheader_H_

#include <iostream> // for basic I/O
#include <cmath>    // for mathematical functions
#include <fstream>
#include <string>
#include <vector>

#include <GyotoSmartPointer.h>
#include <GyotoObject.h>
#include <GyotoAstrobj.h>
#include <GyotoHooks.h>
#include <GyotoDefs.h>
#include "GyotoWorldline.h"
#include <GyotoMetric.h> // Include Gyoto Metric base class

#ifdef GYOTO_USE_XERCES
#include <GyotoRegister.h>
#endif


namespace Gyoto {
  namespace Metric { class BalasinGrumiller; }
}

class Gyoto::Metric::BalasinGrumiller : public Gyoto::Metric::Generic 
{
friend class Gyoto::SmartPointer<Gyoto::Metric::BalasinGrumiller>;
protected:
 double V0value_ ; 
 double Rvalue_ ;
 double r0value_ ;
public:
GYOTO_OBJECT;
BalasinGrumiller(); 
virtual BalasinGrumiller* clone() const ;

 // Mutators / assignment
  // ---------------------
 public:
  // default operator= is fine
  void V0value(const double V0value); 
  void Rvalue(const double Rvalue); 
  void r0value(const double r0value); 

  // Accessors
  // ---------
 public:
  double V0value() const ; 
  double Rvalue() const ; 
  double r0value() const ; 

  void gmunu(double g[4][4], const double x[4]) const ;
  int christoffel(double dst[4][4][4], const double x[4]) const ;

  // Those two are implemented as examples.
 // double gmunu(const double x[4], int mu, int nu) const ;
 // double christoffel(const double coord[4],
//  const int alpha, const int mu, const int nu) const ;
 };

#endif