CInterface.h

Go to the documentation of this file.

#include "../Core/SharedLibraryExportMacros.h"
#include "DiskBinarySystem.h"
#include "OrbitSolver.h"
#include "../Star/CInterface.h"
#include "../Planet/CInterface.h"

extern "C" {
    LIB_PUBLIC extern const int LOCKED_SURFACE_SPIN_EVOL_MODE;

    LIB_PUBLIC extern const int BINARY_EVOL_MODE;

    LIB_PUBLIC extern const int SINGLE_EVOL_MODE;

    LIB_PUBLIC extern const int TABULATION_EVOL_MODE;

    LIB_PUBLIC extern const double NaN;

    struct LIB_PUBLIC DiskBinarySystem;

    struct LIB_PUBLIC OrbitSolver;

    struct LIB_PUBLIC BrokenPowerlawPhaseLagZone;

    struct LIB_PUBLIC EccentricityExpansionCoefficients;

    LIB_PUBLIC void prepare_eccentricity_expansion(
        const char *filename,
        double precision,
        bool pre_load,
        bool disable_precision_fail
    );

    LIB_PUBLIC void set_zone_dissipation(
        BrokenPowerlawPhaseLagZone *zone,

        unsigned num_tidal_frequency_breaks,

        unsigned num_spin_frequency_breaks,

        double *tidal_frequency_breaks,

        double *spin_frequency_breaks,

        double *tidal_frequency_powers,

        double *spin_frequency_powers,

        double reference_phase_lag,

        double inertial_mode_enhancement,

        double inertial_mode_sharpness
    );

    LIB_PUBLIC DiskBinarySystem *create_star_planet_system(
        EvolvingStar *star,

        CPlanet *planet,

        double initial_semimajor,

        double initial_eccentricity,

        double initial_inclination,

        double disk_lock_frequency,

        double disk_dissipation_age,

        double secondary_formation_age
    );

    LIB_PUBLIC DiskBinarySystem *create_star_star_system(
        EvolvingStar *primary,

        EvolvingStar *secondary,

        double initial_semimajor,

        double initial_eccentricity,

        double initial_inclination,

        double disk_lock_frequency,

        double disk_dissipation_age,

        double secondary_formation_age
    );

    LIB_PUBLIC DiskBinarySystem *create_planet_planet_system(
        CPlanet *primary,

        CPlanet *secondary,

        double initial_semimajor,

        double initial_eccentricity,

        double initial_inclination,

        double disk_lock_frequency,

        double disk_dissipation_age
    );

    LIB_PUBLIC void destroy_binary(
        DiskBinarySystem *system
    );

    LIB_PUBLIC void configure_planet(
        CPlanet *planet,

        double age,

        double companion_mass,

        double semimajor,

        double eccentricity,

        const double *spin_angmom,

        const double *inclination,

        const double *periapsis,

        bool locked_surface,

        bool zero_outer_inclination,

        bool zero_outer_periapsis
    );

    LIB_PUBLIC void configure_star(
        EvolvingStar *star,

        double age,

        double companion_mass,

        double semimajor,

        double eccentricity,

        const double *spin_angmom,

        const double *inclination,

        const double *periapsis,

        bool locked_surface,

        bool zero_outer_inclination,

        bool zero_outer_periapsis
    );

    LIB_PUBLIC void configure_system(
        DiskBinarySystem *system,

        double age,

        double semimajor,

        double eccentricity,

        const double *spin_angmom,

        const double *inclination,

        const double *periapsis,

        int evolution_mode
    );

    LIB_PUBLIC OrbitSolver *evolve_system(
        DiskBinarySystem *system,

        double final_age,

        double max_time_step,

        double precision,

        double *required_ages,

        unsigned num_required_ages,

        bool print_progress,

        //values results in infinite timeout.
        double max_runtime,

        unsigned max_time_steps
    );

    LIB_PUBLIC void destroy_solver(
        OrbitSolver *solver
    );

    LIB_PUBLIC unsigned num_evolution_steps(
        OrbitSolver *solver
    );

    LIB_PUBLIC void get_star_planet_evolution(
        const OrbitSolver *solver,

        const DiskBinarySystem *system,

        const EvolvingStar *star,

        const CPlanet *planet,

        double *age,

        double *semimajor,

        double *eccentricity,

        double *envelope_inclination,

        double *core_inclination,

        double *envelope_periapsis,

        double *core_periapsis,

        double *envelope_angmom,

        double *core_angmom,

        double *planet_inclination,

        double *planet_periapsis,

        double *planet_angmom,

        int *evolution_mode,

        bool *wind_saturation,

        double *semimajor_rate,

        double *eccentricity_rate,

        double *envelope_inclination_rate,

        double *core_inclination_rate,

        double *envelope_periapsis_rate,

        double *core_periapsis_rate,

        double *envelope_angmom_rate,

        double *core_angmom_rate,

        double *planet_inclination_rate,

        double *planet_periapsis_rate,

        double *planet_angmom_rate
    );

    LIB_PUBLIC void get_star_star_evolution(
        const OrbitSolver *solver,

        const DiskBinarySystem *system,

        const EvolvingStar *primary,

        const EvolvingStar *secondary,

        double *age,

        double *semimajor,

        double *eccentricity,

        double *primary_envelope_inclination,

        double *primary_core_inclination,

        double *primary_envelope_periapsis,

        double *primary_core_periapsis,

        double *primary_envelope_angmom,

        double *primary_core_angmom,

        double *secondary_envelope_inclination,

        double *secondary_core_inclination,

        double *secondary_envelope_periapsis,

        double *secondary_core_periapsis,

        double *secondary_envelope_angmom,

        double *secondary_core_angmom,

        int *evolution_mode,

        bool *primary_wind_saturation,

        bool *secondary_wind_saturation,

        double *semimajor_rate,

        double *eccentricity_rate,

        double *primary_envelope_inclination_rate,

        double *primary_core_inclination_rate,

        double *primary_envelope_periapsis_rate,

        double *primary_core_periapsis_rate,

        double *primary_envelope_angmom_rate,

        double *primary_core_angmom_rate,

        double *secondary_envelope_inclination_rate,

        double *secondary_core_inclination_rate,

        double *secondary_envelope_periapsis_rate,

        double *secondary_core_periapsis_rate,

        double *secondary_envelope_angmom_rate,

        double *secondary_core_angmom_rate
    );

    LIB_PUBLIC void get_planet_planet_evolution(
        const OrbitSolver *solver,

        const DiskBinarySystem *system,

        const CPlanet *primary,

        const CPlanet *secondary,

        double *age,

        double *semimajor,

        double *eccentricity,

        double *primary_inclination,

        double *primary_periapsis,

        double *primary_angmom,

        double *secondary_inclination,

        double *secondary_periapsis,

        double *secondary_angmom,

        int *evolution_mode,

        double *semimajor_rate,

        double *eccentricity_rate,

        double *primary_inclination_rate,

        double *primary_periapsis_rate,

        double *primary_angmom_rate,

        double *secondary_inclination_rate,

        double *secondary_periapsis_rate,

        double *secondary_angmom_rate
    );

    LIB_PUBLIC void get_star_planet_final_state(
        const OrbitSolver *solver,

        const DiskBinarySystem *system,

        const EvolvingStar *star,

        const CPlanet *planet,

        double *age,

        double *semimajor,

        double *eccentricity,

        double *envelope_inclination,

        double *core_inclination,

        double *envelope_periapsis,

        double *core_periapsis,

        double *envelope_angmom,

        double *core_angmom,

        double *planet_inclination,

        double *planet_periapsis,

        double *planet_angmom,

        int *evolution_mode,

        bool *wind_saturation
    );

    LIB_PUBLIC void get_star_star_final_state(
        const OrbitSolver *solver,

        const DiskBinarySystem *system,

        const EvolvingStar *primary,

        const EvolvingStar *secondary,

        double *age,

        double *semimajor,

        double *eccentricity,

        double *primary_envelope_inclination,

        double *primary_core_inclination,

        double *primary_envelope_periapsis,

        double *primary_core_periapsis,

        double *primary_envelope_angmom,

        double *primary_core_angmom,

        double *secondary_envelope_inclination,

        double *secondary_core_inclination,

        double *secondary_envelope_periapsis,

        double *secondary_core_periapsis,

        double *secondary_envelope_angmom,

        double *secondary_core_angmom,

        int *evolution_mode,

        bool *primary_wind_saturation,

        bool *secondary_wind_saturation
    );

    LIB_PUBLIC double get_expansion_coeff_precision(
        int m,

        int s
    );

    LIB_PUBLIC double evaluate_expansion_coeff(
        int m,

        int s,

        double e,

        bool deriv
    );

}//End Extern "C"