#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Destination given bearing on a great circle";

double f_if(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) {
        float r25819 = lambda1;
        float r25820 = theta;
        float r25821 = sin(r25820);
        float r25822 = delta;
        float r25823 = sin(r25822);
        float r25824 = r25821 * r25823;
        float r25825 = phi1;
        float r25826 = cos(r25825);
        float r25827 = r25824 * r25826;
        float r25828 = cos(r25822);
        float r25829 = sin(r25825);
        float r25830 = r25829 * r25828;
        float r25831 = r25826 * r25823;
        float r25832 = cos(r25820);
        float r25833 = r25831 * r25832;
        float r25834 = r25830 + r25833;
        float r25835 = asin(r25834);
        float r25836 = sin(r25835);
        float r25837 = r25829 * r25836;
        float r25838 = r25828 - r25837;
        float r25839 = atan2(r25827, r25838);
        float r25840 = r25819 + r25839;
        return r25840;
}

double f_id(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r25841 = lambda1;
        double r25842 = theta;
        double r25843 = sin(r25842);
        double r25844 = delta;
        double r25845 = sin(r25844);
        double r25846 = r25843 * r25845;
        double r25847 = phi1;
        double r25848 = cos(r25847);
        double r25849 = r25846 * r25848;
        double r25850 = cos(r25844);
        double r25851 = sin(r25847);
        double r25852 = r25851 * r25850;
        double r25853 = r25848 * r25845;
        double r25854 = cos(r25842);
        double r25855 = r25853 * r25854;
        double r25856 = r25852 + r25855;
        double r25857 = asin(r25856);
        double r25858 = sin(r25857);
        double r25859 = r25851 * r25858;
        double r25860 = r25850 - r25859;
        double r25861 = atan2(r25849, r25860);
        double r25862 = r25841 + r25861;
        return r25862;
}


double f_of(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) {
        float r25863 = lambda1;
        float r25864 = theta;
        float r25865 = sin(r25864);
        float r25866 = delta;
        float r25867 = sin(r25866);
        float r25868 = r25865 * r25867;
        float r25869 = phi1;
        float r25870 = cos(r25869);
        float r25871 = r25868 * r25870;
        float r25872 = cos(r25866);
        float r25873 = sin(r25869);
        float r25874 = r25873 * r25872;
        float r25875 = r25870 * r25867;
        float r25876 = cos(r25864);
        float r25877 = r25875 * r25876;
        float r25878 = r25874 + r25877;
        float r25879 = asin(r25878);
        float r25880 = sin(r25879);
        float r25881 = r25873 * r25880;
        float r25882 = r25872 - r25881;
        float r25883 = atan2(r25871, r25882);
        float r25884 = r25863 + r25883;
        return r25884;
}

double f_od(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r25885 = lambda1;
        double r25886 = theta;
        double r25887 = sin(r25886);
        double r25888 = delta;
        double r25889 = sin(r25888);
        double r25890 = r25887 * r25889;
        double r25891 = phi1;
        double r25892 = cos(r25891);
        double r25893 = r25890 * r25892;
        double r25894 = cos(r25888);
        double r25895 = sin(r25891);
        double r25896 = r25895 * r25894;
        double r25897 = r25892 * r25889;
        double r25898 = cos(r25886);
        double r25899 = r25897 * r25898;
        double r25900 = r25896 + r25899;
        double r25901 = asin(r25900);
        double r25902 = sin(r25901);
        double r25903 = r25895 * r25902;
        double r25904 = r25894 - r25903;
        double r25905 = atan2(r25893, r25904);
        double r25906 = r25885 + r25905;
        return r25906;
}

void mpfr_fmod2(mpfr_t r, mpfr_t n, mpfr_t d, mpfr_rnd_t rmd) {
        mpfr_fmod(r, n, d, rmd);
        if (mpfr_cmp_ui(r, 0) < 0) mpfr_add(r, r, d, rmd);
}


static mpfr_t r25907, r25908, r25909, r25910, r25911, r25912, r25913, r25914, r25915, r25916, r25917, r25918, r25919, r25920, r25921, r25922, r25923, r25924, r25925, r25926, r25927, r25928;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r25907);
        mpfr_init(r25908);
        mpfr_init(r25909);
        mpfr_init(r25910);
        mpfr_init(r25911);
        mpfr_init(r25912);
        mpfr_init(r25913);
        mpfr_init(r25914);
        mpfr_init(r25915);
        mpfr_init(r25916);
        mpfr_init(r25917);
        mpfr_init(r25918);
        mpfr_init(r25919);
        mpfr_init(r25920);
        mpfr_init(r25921);
        mpfr_init(r25922);
        mpfr_init(r25923);
        mpfr_init(r25924);
        mpfr_init(r25925);
        mpfr_init(r25926);
        mpfr_init(r25927);
        mpfr_init(r25928);
}

double f_im(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r25907, lambda1, MPFR_RNDN);
        mpfr_set_d(r25908, theta, MPFR_RNDN);
        mpfr_sin(r25909, r25908, MPFR_RNDN);
        mpfr_set_d(r25910, delta, MPFR_RNDN);
        mpfr_sin(r25911, r25910, MPFR_RNDN);
        mpfr_mul(r25912, r25909, r25911, MPFR_RNDN);
        mpfr_set_d(r25913, phi1, MPFR_RNDN);
        mpfr_cos(r25914, r25913, MPFR_RNDN);
        mpfr_mul(r25915, r25912, r25914, MPFR_RNDN);
        mpfr_cos(r25916, r25910, MPFR_RNDN);
        mpfr_sin(r25917, r25913, MPFR_RNDN);
        mpfr_mul(r25918, r25917, r25916, MPFR_RNDN);
        mpfr_mul(r25919, r25914, r25911, MPFR_RNDN);
        mpfr_cos(r25920, r25908, MPFR_RNDN);
        mpfr_mul(r25921, r25919, r25920, MPFR_RNDN);
        mpfr_add(r25922, r25918, r25921, MPFR_RNDN);
        mpfr_asin(r25923, r25922, MPFR_RNDN);
        mpfr_sin(r25924, r25923, MPFR_RNDN);
        mpfr_mul(r25925, r25917, r25924, MPFR_RNDN);
        mpfr_sub(r25926, r25916, r25925, MPFR_RNDN);
        mpfr_atan2(r25927, r25915, r25926, MPFR_RNDN);
        mpfr_add(r25928, r25907, r25927, MPFR_RNDN);
        return mpfr_get_d(r25928, MPFR_RNDN);
}

static mpfr_t r25929, r25930, r25931, r25932, r25933, r25934, r25935, r25936, r25937, r25938, r25939, r25940, r25941, r25942, r25943, r25944, r25945, r25946, r25947, r25948, r25949, r25950;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25929);
        mpfr_init(r25930);
        mpfr_init(r25931);
        mpfr_init(r25932);
        mpfr_init(r25933);
        mpfr_init(r25934);
        mpfr_init(r25935);
        mpfr_init(r25936);
        mpfr_init(r25937);
        mpfr_init(r25938);
        mpfr_init(r25939);
        mpfr_init(r25940);
        mpfr_init(r25941);
        mpfr_init(r25942);
        mpfr_init(r25943);
        mpfr_init(r25944);
        mpfr_init(r25945);
        mpfr_init(r25946);
        mpfr_init(r25947);
        mpfr_init(r25948);
        mpfr_init(r25949);
        mpfr_init(r25950);
}

double f_fm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r25929, lambda1, MPFR_RNDN);
        mpfr_set_d(r25930, theta, MPFR_RNDN);
        mpfr_sin(r25931, r25930, MPFR_RNDN);
        mpfr_set_d(r25932, delta, MPFR_RNDN);
        mpfr_sin(r25933, r25932, MPFR_RNDN);
        mpfr_mul(r25934, r25931, r25933, MPFR_RNDN);
        mpfr_set_d(r25935, phi1, MPFR_RNDN);
        mpfr_cos(r25936, r25935, MPFR_RNDN);
        mpfr_mul(r25937, r25934, r25936, MPFR_RNDN);
        mpfr_cos(r25938, r25932, MPFR_RNDN);
        mpfr_sin(r25939, r25935, MPFR_RNDN);
        mpfr_mul(r25940, r25939, r25938, MPFR_RNDN);
        mpfr_mul(r25941, r25936, r25933, MPFR_RNDN);
        mpfr_cos(r25942, r25930, MPFR_RNDN);
        mpfr_mul(r25943, r25941, r25942, MPFR_RNDN);
        mpfr_add(r25944, r25940, r25943, MPFR_RNDN);
        mpfr_asin(r25945, r25944, MPFR_RNDN);
        mpfr_sin(r25946, r25945, MPFR_RNDN);
        mpfr_mul(r25947, r25939, r25946, MPFR_RNDN);
        mpfr_sub(r25948, r25938, r25947, MPFR_RNDN);
        mpfr_atan2(r25949, r25937, r25948, MPFR_RNDN);
        mpfr_add(r25950, r25929, r25949, MPFR_RNDN);
        return mpfr_get_d(r25950, MPFR_RNDN);
}

static mpfr_t r25951, r25952, r25953, r25954, r25955, r25956, r25957, r25958, r25959, r25960, r25961, r25962, r25963, r25964, r25965, r25966, r25967, r25968, r25969, r25970, r25971, r25972;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25951);
        mpfr_init(r25952);
        mpfr_init(r25953);
        mpfr_init(r25954);
        mpfr_init(r25955);
        mpfr_init(r25956);
        mpfr_init(r25957);
        mpfr_init(r25958);
        mpfr_init(r25959);
        mpfr_init(r25960);
        mpfr_init(r25961);
        mpfr_init(r25962);
        mpfr_init(r25963);
        mpfr_init(r25964);
        mpfr_init(r25965);
        mpfr_init(r25966);
        mpfr_init(r25967);
        mpfr_init(r25968);
        mpfr_init(r25969);
        mpfr_init(r25970);
        mpfr_init(r25971);
        mpfr_init(r25972);
}

double f_dm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r25951, lambda1, MPFR_RNDN);
        mpfr_set_d(r25952, theta, MPFR_RNDN);
        mpfr_sin(r25953, r25952, MPFR_RNDN);
        mpfr_set_d(r25954, delta, MPFR_RNDN);
        mpfr_sin(r25955, r25954, MPFR_RNDN);
        mpfr_mul(r25956, r25953, r25955, MPFR_RNDN);
        mpfr_set_d(r25957, phi1, MPFR_RNDN);
        mpfr_cos(r25958, r25957, MPFR_RNDN);
        mpfr_mul(r25959, r25956, r25958, MPFR_RNDN);
        mpfr_cos(r25960, r25954, MPFR_RNDN);
        mpfr_sin(r25961, r25957, MPFR_RNDN);
        mpfr_mul(r25962, r25961, r25960, MPFR_RNDN);
        mpfr_mul(r25963, r25958, r25955, MPFR_RNDN);
        mpfr_cos(r25964, r25952, MPFR_RNDN);
        mpfr_mul(r25965, r25963, r25964, MPFR_RNDN);
        mpfr_add(r25966, r25962, r25965, MPFR_RNDN);
        mpfr_asin(r25967, r25966, MPFR_RNDN);
        mpfr_sin(r25968, r25967, MPFR_RNDN);
        mpfr_mul(r25969, r25961, r25968, MPFR_RNDN);
        mpfr_sub(r25970, r25960, r25969, MPFR_RNDN);
        mpfr_atan2(r25971, r25959, r25970, MPFR_RNDN);
        mpfr_add(r25972, r25951, r25971, MPFR_RNDN);
        return mpfr_get_d(r25972, MPFR_RNDN);
}