geodesics.hh

Go to the documentation of this file.

#ifndef GEODESICS_HH
#define GEODESICS_HH

#include "cvode/cvode.h"
#include "nvector/nvector_serial.h"  /* serial N_Vector types, fcts., macros */

void printfutil(double *y);
void f_util(double* y, double Jspin, double mu);

class geodesic{
public:
   geodesic(double r0, double phi0, double E0, double L0, 
      int dir, double Jspin, double mu=0, double dt=1);
   ~geodesic();
   
   bool integrate(double D_t, int& N, double rmin=0, bool stop_alr=false);
   //returns "success" and the variables below:
   double* r, *pr, *phi, *tau, *t, *Et, *Lt, *hplus, *hcross;
      
private:
   N_Vector y;
   realtype udata[2], tret;      //udata[0] = Jspin, udata[1] = mu
   realtype* ydata;
   void* cvode_mem;

};

#endif