swh:1:snp:122bde0cb0e54f3d002c308e151c63f07e45e6be
Tip revision: b24da720877d30a8991fc430fabe1b321b62530b authored by Hanno Rein on 20 July 2023, 21:46:00 UTC
Updating version to 3.26.1
Updating version to 3.26.1
Tip revision: b24da72
tools.h
/**
* @file tools.h
* @brief Tools for creating distributions.
* @author Hanno Rein <hanno@hanno-rein.de>
*
* @section LICENSE
* Copyright (c) 2011 Hanno Rein, Shangfei Liu
*
* This file is part of rebound.
*
* rebound 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.
*
* rebound 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 rebound. If not, see <http://www.gnu.org/licenses/>.
*
*/
#ifndef TOOLS_H
#define TOOLS_H
#include <stdint.h>
struct reb_simulation;
struct reb_particles;
/**
* @brief Returns deltad/delta
* @details Note, there is a typo in Gozdziewski et al 2001.
* @param r REBOUND simulation to be considered.
*/
double reb_tools_megno_deltad_delta(struct reb_simulation* const r);
/**
* @brief Update MEGNO after a successful timestep by adding dY (=ddelta/delta*dt)
* @param r REBOUND simulation to be considered.
* @param dY Increment for MEGNO Y
*/
void reb_tools_megno_update(struct reb_simulation* r, double dY);
/**
* @brief Init random number generator based on time and process id.
*/
void reb_tools_init_srand(struct reb_simulation* r);
/**
* @brief Convert angles for orbit routines
*/
double reb_tools_M_to_E(double e, double M);
/**
* @brief Convert angles for orbit routines
*/
double reb_tools_M_to_f(double e, double M);
/**
* @brief Convert angles for orbit routines
*/
double reb_tools_E_to_f(double e, double M);
/**
* @brief Kepler solver in Pal coordinates
*/
void reb_tools_solve_kepler_pal(double h, double k, double lambda, double* p, double* q);
/**
* @brief Convert particle to Pal coordinates
*/
void reb_tools_particle_to_pal(double G, struct reb_particle p, struct reb_particle primary, double *a, double* lambda, double* k, double* h, double* ix, double* iy);
/**
* @brief internal function to handle outputs for the Fast Simulation Restarter.
*/
void reb_fsr_heartbeat(struct reb_simulation* const r);
#endif // TOOLS_H