#include "orsa_orbit_gsl.h"#include <gsl/gsl_roots.h>#include <gsl/gsl_errno.h>#include <gsl/gsl_siman.h>#include <gsl/gsl_multimin.h>#include <gsl/gsl_rng.h>#include <gsl/gsl_multifit_nlin.h>#include <gsl/gsl_blas.h>#include <gsl/gsl_diff.h>#include <gsl/gsl_deriv.h>#include <gsl/gsl_linalg.h>#include <gsl/gsl_cblas.h>#include <gsl/gsl_randist.h>#include <gsl/gsl_eigen.h>#include "orsa_file.h"#include "orsa_error.h"#include <sys/time.h>#include <time.h>#include <string.h>#include <algorithm>#include <limits>Include dependency graph for orsa_orbit_gsl.cc:
Go to the source code of this file.
Namespaces | |
| namespace | orsa |
Defines | |
| #define | N_TRIES 10 |
| #define | ITERS_FIXED_T 3 |
| #define | STEP_SIZE 1.0 |
| #define | K 1.0e5 |
| #define | T_INITIAL 0.01 |
| #define | MU_T 1.005 |
| #define | T_MIN 1.0e-6 |
| #define | FIT(i) gsl_vector_get(s->x, i) |
| #define | ERR(i) sqrt(gsl_matrix_get(covar,i,i)) |
Functions | |
| double | orsa::Weight (const Observation &obs1, const Observation &obs2, const double optimal_dt=FromUnits(1, DAY)) |
| double | orsa::Weight (vector< Observation > &obs, const double optimal_dt=FromUnits(1, DAY)) |
| void | orsa::BuildObservationTriplets (const vector< Observation > &obs, vector< ThreeObservations > &triplets, const double optimal_dt=FromUnits(1, DAY)) |
| double | orsa::MOID_f (const gsl_vector *v, void *params) |
| double | orsa::MOID (const Orbit &, const Orbit &, Vector &, Vector &) |
| Minimal Orbit Intersection Distance between two orbits. | |
| double | orsa::MOID2RB_f (const gsl_vector *v, void *params) |
| double | orsa::MOID2RB (const Vector &, const Vector &, const Orbit &, const Orbit &, Vector &, Vector &) |
| Minimal Orbit Intersection Distance between two orbits with different reference bodies. | |
| double | orsa::E1 (void *p) |
| double | orsa::M1 (void *p1, void *p2) |
| void | orsa::S1 (const gsl_rng *r, void *p, double step_size) |
| int | orsa::least_sq_f (const gsl_vector *v, void *params, gsl_vector *f) |
| double | orsa::least_sq_diff_f (double x, void *params) |
| int | orsa::least_sq_df (const gsl_vector *v, void *params, gsl_matrix *J) |
| int | orsa::least_sq_fdf (const gsl_vector *v, void *params, gsl_vector *f, gsl_matrix *J) |
| void | orsa::OrbitDifferentialCorrectionsLeastSquares (OrbitWithCovarianceMatrixGSL &orbit, const vector< Observation > &obs) |
| int | orsa::gauss_v_f (const gsl_vector *v, void *params, gsl_vector *f) |
| double | orsa::gauss_v_diff_f (double x, void *params) |
| int | orsa::gauss_v_df (const gsl_vector *v, void *params, gsl_matrix *J) |
| int | orsa::gauss_v_fdf (const gsl_vector *v, void *params, gsl_vector *f, gsl_matrix *J) |
| void | orsa::gauss_v (const Vector &r, Vector &v, const vector< Observation > &obs) |
| OrbitWithCovarianceMatrixGSL | orsa::Compute (const vector< Observation > &obs) |
| double | orsa::poly_8 (double x, void *params) |
| double | orsa::poly_8_df (double x, void *params) |
| void | orsa::poly_8_fdf (double x, void *params, double *y, double *dy) |
| void | orsa::poly_8_gsl_solve (poly_8_params ¶ms, vector< poly_8_solution > &solutions) |
| void | orsa::Compute_Gauss (const vector< Observation > &obs_in, vector< PreliminaryOrbit > &preliminary_orbits) |
| void | orsa::Compute_TestMethod (const vector< Observation > &obs_in, vector< PreliminaryOrbit > &preliminary_orbits) |
| double | orsa::CloseApproaches_gsl_f (const gsl_vector *v, void *params) |
| void | orsa::SearchCloseApproaches (const Evolution *evol, const unsigned int obj_index, const unsigned int index, std::vector< CloseApproach > &clapp, const double distance_threshold, const double time_accuracy) |
Variables | |
| const gsl_siman_params_t | orsa::params = {N_TRIES, ITERS_FIXED_T, STEP_SIZE, K, T_INITIAL, MU_T, T_MIN} |
| const double | orsa::d_a = 0.1 |
| const double | orsa::d_e = 0.01 |
| const double | orsa::d_i = 0.1*(pi/180) |
| const double | orsa::d_omega_node = 1.0*(pi/180) |
| const double | orsa::d_omega_pericenter = 1.0*(pi/180) |
| const double | orsa::d_M = 5.0*(pi/180) |
| const double | orsa::c_a = secure_pow(1.0/d_a, 2) |
| const double | orsa::c_e = secure_pow(1.0/d_e, 2) |
| const double | orsa::c_i = secure_pow(1.0/d_i, 2) |
| const double | orsa::c_omega_node = secure_pow(1.0/d_omega_node, 2) |
| const double | orsa::c_omega_pericenter = secure_pow(1.0/d_omega_pericenter, 2) |
| const double | orsa::c_M = secure_pow(1.0/d_M, 2) |
| #define ERR | ( | i | ) | sqrt(gsl_matrix_get(covar,i,i)) |
Referenced by orsa::OrbitDifferentialCorrectionsLeastSquares().
| #define FIT | ( | i | ) | gsl_vector_get(s->x, i) |
Referenced by orsa::OrbitDifferentialCorrectionsLeastSquares().
| #define ITERS_FIXED_T 3 |
Definition at line 551 of file orsa_orbit_gsl.cc.
| #define K 1.0e5 |
Definition at line 557 of file orsa_orbit_gsl.cc.
| #define MU_T 1.005 |
Definition at line 563 of file orsa_orbit_gsl.cc.
| #define N_TRIES 10 |
Definition at line 548 of file orsa_orbit_gsl.cc.
| #define STEP_SIZE 1.0 |
Definition at line 554 of file orsa_orbit_gsl.cc.
| #define T_INITIAL 0.01 |
Definition at line 560 of file orsa_orbit_gsl.cc.
| #define T_MIN 1.0e-6 |
Definition at line 564 of file orsa_orbit_gsl.cc.
1.4.7