integrators.c File Reference

#include "libgravix2/integrators.h"
#include <assert.h>
#include <math.h>
#include "libgravix2/config.h"
#include "libgravix2/helpers.h"
#include "libgravix2/pot.h"

Include dependency graph for integrators.c:

Functions
void	grvx_integration_step (struct GrvxQP *qp, struct GrvxQP *e, double h, const struct GrvxPlanets *planets)
	Single integration step. More...
unsigned	grvx_integration_loop (struct GrvxQP *qp, double h, unsigned n, const struct GrvxPlanets *planets)
	Multiple integration steps. More...

Function Documentation

grvx_integration_loop()

unsigned grvx_integration_loop	(	struct GrvxQP *	qp,
		double	h,
		unsigned	n,
		const struct GrvxPlanets *	planets
	)

Multiple integration steps.

qp is advanced by the integrator for a maximum of n steps. The propagation is stopped prematurely if the distance to a planet becomes too small. (The threshold is controlled via GrvxConfig.min_dist.)

Parameters

qp	Phase space.
h	Step size passed to grvx_integration_step().
n	Maximum number of integration steps.
planets	Planets handle.

Returns

Number of processed integration steps. (Smaller than n if integration was stopped prematurely.)

Here is the call graph for this function:

Here is the caller graph for this function:

grvx_integration_step()

void grvx_integration_step	(	struct GrvxQP *	qp,
		struct GrvxQP *	e,
		double	h,
		const struct GrvxPlanets *	planets
	)

Single integration step.

qp is advanced by the integrator for a single step. Accumulation errors are compensated by tracking them in e. Initialize e with zeros for the first iteration.

Parameters

qp	Phase space.
e	Accumulation error.
h	Step size.
planets	Planets handle.

Here is the caller graph for this function: