59#include <lal/XLALGSL.h>
60#include <lal/LALInspiral.h>
63#define UNUSED __attribute__ ((unused))
98 integrator->
input = input;
101 integrator->
type = gsl_odeiv_step_rkf45;
112 XLAL_CALLGSL( integrator->
control = gsl_odeiv_control_standard_new(1.0e-2, 1.0e-2, 1.0, 1.0) );
161 gsl_odeiv_system sys;
172 if ( !integrator->
input )
189 sys.params = &gslParams;
200 &t,
input->
h, &h, integrator->
y) );
203 if ( gslStatus != GSL_SUCCESS )
void XLALRungeKutta4Free(rk4GSLIntegrator *integrator)
static int derivativeGSLWrapper(REAL8 t, const REAL8 y[], REAL8 dydx[], void *params)
void LALRungeKutta4(LALStatus *status, REAL8Vector *yout, rk4GSLIntegrator *integrator, void *params)
rk4GSLIntegrator * XLALRungeKutta4Init(INT4 n, rk4In *input)
int XLALRungeKutta4(REAL8Vector *yout, rk4GSLIntegrator *integrator, void *params)
#define INITSTATUS(statusptr)
#define RETURN(statusptr)
#define XLAL_CALLGSL(statement)
#define XLAL_ERROR_VOID(...)
#define XLAL_ERROR_NULL(...)
int XLALPrintError(const char *fmt,...) _LAL_GCC_PRINTF_FORMAT_(1
#define XLAL_PRINT_DEPRECATION_WARNING(replacement)
Structure containing steps and controls for the GSL Runge-Kutta solver.
gsl_odeiv_control * control
const gsl_odeiv_step_type * type
gsl_odeiv_evolve * evolve
Structure used as an input to Runge-Kutta solver.