#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 r25897 = lambda1;
        float r25898 = theta;
        float r25899 = sin(r25898);
        float r25900 = delta;
        float r25901 = sin(r25900);
        float r25902 = r25899 * r25901;
        float r25903 = phi1;
        float r25904 = cos(r25903);
        float r25905 = r25902 * r25904;
        float r25906 = cos(r25900);
        float r25907 = sin(r25903);
        float r25908 = r25907 * r25906;
        float r25909 = r25904 * r25901;
        float r25910 = cos(r25898);
        float r25911 = r25909 * r25910;
        float r25912 = r25908 + r25911;
        float r25913 = asin(r25912);
        float r25914 = sin(r25913);
        float r25915 = r25907 * r25914;
        float r25916 = r25906 - r25915;
        float r25917 = atan2(r25905, r25916);
        float r25918 = r25897 + r25917;
        return r25918;
}

double f_id(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r25919 = lambda1;
        double r25920 = theta;
        double r25921 = sin(r25920);
        double r25922 = delta;
        double r25923 = sin(r25922);
        double r25924 = r25921 * r25923;
        double r25925 = phi1;
        double r25926 = cos(r25925);
        double r25927 = r25924 * r25926;
        double r25928 = cos(r25922);
        double r25929 = sin(r25925);
        double r25930 = r25929 * r25928;
        double r25931 = r25926 * r25923;
        double r25932 = cos(r25920);
        double r25933 = r25931 * r25932;
        double r25934 = r25930 + r25933;
        double r25935 = asin(r25934);
        double r25936 = sin(r25935);
        double r25937 = r25929 * r25936;
        double r25938 = r25928 - r25937;
        double r25939 = atan2(r25927, r25938);
        double r25940 = r25919 + r25939;
        return r25940;
}


double f_of(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) {
        float r25941 = delta;
        float r25942 = sin(r25941);
        float r25943 = phi1;
        float r25944 = cos(r25943);
        float r25945 = theta;
        float r25946 = sin(r25945);
        float r25947 = r25944 * r25946;
        float r25948 = r25942 * r25947;
        float r25949 = sin(r25943);
        float r25950 = -r25949;
        float r25951 = cos(r25945);
        float r25952 = r25944 * r25942;
        float r25953 = cos(r25941);
        float r25954 = r25949 * r25953;
        float r25955 = fma(r25951, r25952, r25954);
        float r25956 = asin(r25955);
        float r25957 = sin(r25956);
        float r25958 = fma(r25950, r25957, r25953);
        float r25959 = exp(r25958);
        float r25960 = log(r25959);
        float r25961 = atan2(r25948, r25960);
        float r25962 = lambda1;
        float r25963 = r25961 + r25962;
        return r25963;
}

double f_od(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r25964 = delta;
        double r25965 = sin(r25964);
        double r25966 = phi1;
        double r25967 = cos(r25966);
        double r25968 = theta;
        double r25969 = sin(r25968);
        double r25970 = r25967 * r25969;
        double r25971 = r25965 * r25970;
        double r25972 = sin(r25966);
        double r25973 = -r25972;
        double r25974 = cos(r25968);
        double r25975 = r25967 * r25965;
        double r25976 = cos(r25964);
        double r25977 = r25972 * r25976;
        double r25978 = fma(r25974, r25975, r25977);
        double r25979 = asin(r25978);
        double r25980 = sin(r25979);
        double r25981 = fma(r25973, r25980, r25976);
        double r25982 = exp(r25981);
        double r25983 = log(r25982);
        double r25984 = atan2(r25971, r25983);
        double r25985 = lambda1;
        double r25986 = r25984 + r25985;
        return r25986;
}

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 r25987, r25988, r25989, r25990, r25991, r25992, r25993, r25994, r25995, r25996, r25997, r25998, r25999, r26000, r26001, r26002, r26003, r26004, r26005, r26006, r26007, r26008;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r25987);
        mpfr_init(r25988);
        mpfr_init(r25989);
        mpfr_init(r25990);
        mpfr_init(r25991);
        mpfr_init(r25992);
        mpfr_init(r25993);
        mpfr_init(r25994);
        mpfr_init(r25995);
        mpfr_init(r25996);
        mpfr_init(r25997);
        mpfr_init(r25998);
        mpfr_init(r25999);
        mpfr_init(r26000);
        mpfr_init(r26001);
        mpfr_init(r26002);
        mpfr_init(r26003);
        mpfr_init(r26004);
        mpfr_init(r26005);
        mpfr_init(r26006);
        mpfr_init(r26007);
        mpfr_init(r26008);
}

double f_im(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r25987, lambda1, MPFR_RNDN);
        mpfr_set_d(r25988, theta, MPFR_RNDN);
        mpfr_sin(r25989, r25988, MPFR_RNDN);
        mpfr_set_d(r25990, delta, MPFR_RNDN);
        mpfr_sin(r25991, r25990, MPFR_RNDN);
        mpfr_mul(r25992, r25989, r25991, MPFR_RNDN);
        mpfr_set_d(r25993, phi1, MPFR_RNDN);
        mpfr_cos(r25994, r25993, MPFR_RNDN);
        mpfr_mul(r25995, r25992, r25994, MPFR_RNDN);
        mpfr_cos(r25996, r25990, MPFR_RNDN);
        mpfr_sin(r25997, r25993, MPFR_RNDN);
        mpfr_mul(r25998, r25997, r25996, MPFR_RNDN);
        mpfr_mul(r25999, r25994, r25991, MPFR_RNDN);
        mpfr_cos(r26000, r25988, MPFR_RNDN);
        mpfr_mul(r26001, r25999, r26000, MPFR_RNDN);
        mpfr_add(r26002, r25998, r26001, MPFR_RNDN);
        mpfr_asin(r26003, r26002, MPFR_RNDN);
        mpfr_sin(r26004, r26003, MPFR_RNDN);
        mpfr_mul(r26005, r25997, r26004, MPFR_RNDN);
        mpfr_sub(r26006, r25996, r26005, MPFR_RNDN);
        mpfr_atan2(r26007, r25995, r26006, MPFR_RNDN);
        mpfr_add(r26008, r25987, r26007, MPFR_RNDN);
        return mpfr_get_d(r26008, MPFR_RNDN);
}

static mpfr_t r26009, r26010, r26011, r26012, r26013, r26014, r26015, r26016, r26017, r26018, r26019, r26020, r26021, r26022, r26023, r26024, r26025, r26026, r26027, r26028, r26029, r26030, r26031;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26009);
        mpfr_init(r26010);
        mpfr_init(r26011);
        mpfr_init(r26012);
        mpfr_init(r26013);
        mpfr_init(r26014);
        mpfr_init(r26015);
        mpfr_init(r26016);
        mpfr_init(r26017);
        mpfr_init(r26018);
        mpfr_init(r26019);
        mpfr_init(r26020);
        mpfr_init(r26021);
        mpfr_init(r26022);
        mpfr_init(r26023);
        mpfr_init(r26024);
        mpfr_init(r26025);
        mpfr_init(r26026);
        mpfr_init(r26027);
        mpfr_init(r26028);
        mpfr_init(r26029);
        mpfr_init(r26030);
        mpfr_init(r26031);
}

double f_fm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r26009, delta, MPFR_RNDN);
        mpfr_sin(r26010, r26009, MPFR_RNDN);
        mpfr_set_d(r26011, phi1, MPFR_RNDN);
        mpfr_cos(r26012, r26011, MPFR_RNDN);
        mpfr_set_d(r26013, theta, MPFR_RNDN);
        mpfr_sin(r26014, r26013, MPFR_RNDN);
        mpfr_mul(r26015, r26012, r26014, MPFR_RNDN);
        mpfr_mul(r26016, r26010, r26015, MPFR_RNDN);
        mpfr_sin(r26017, r26011, MPFR_RNDN);
        mpfr_neg(r26018, r26017, MPFR_RNDN);
        mpfr_cos(r26019, r26013, MPFR_RNDN);
        mpfr_mul(r26020, r26012, r26010, MPFR_RNDN);
        mpfr_cos(r26021, r26009, MPFR_RNDN);
        mpfr_mul(r26022, r26017, r26021, MPFR_RNDN);
        mpfr_fma(r26023, r26019, r26020, r26022, MPFR_RNDN);
        mpfr_asin(r26024, r26023, MPFR_RNDN);
        mpfr_sin(r26025, r26024, MPFR_RNDN);
        mpfr_fma(r26026, r26018, r26025, r26021, MPFR_RNDN);
        mpfr_exp(r26027, r26026, MPFR_RNDN);
        mpfr_log(r26028, r26027, MPFR_RNDN);
        mpfr_atan2(r26029, r26016, r26028, MPFR_RNDN);
        mpfr_set_d(r26030, lambda1, MPFR_RNDN);
        mpfr_add(r26031, r26029, r26030, MPFR_RNDN);
        return mpfr_get_d(r26031, MPFR_RNDN);
}

static mpfr_t r26032, r26033, r26034, r26035, r26036, r26037, r26038, r26039, r26040, r26041, r26042, r26043, r26044, r26045, r26046, r26047, r26048, r26049, r26050, r26051, r26052, r26053, r26054;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26032);
        mpfr_init(r26033);
        mpfr_init(r26034);
        mpfr_init(r26035);
        mpfr_init(r26036);
        mpfr_init(r26037);
        mpfr_init(r26038);
        mpfr_init(r26039);
        mpfr_init(r26040);
        mpfr_init(r26041);
        mpfr_init(r26042);
        mpfr_init(r26043);
        mpfr_init(r26044);
        mpfr_init(r26045);
        mpfr_init(r26046);
        mpfr_init(r26047);
        mpfr_init(r26048);
        mpfr_init(r26049);
        mpfr_init(r26050);
        mpfr_init(r26051);
        mpfr_init(r26052);
        mpfr_init(r26053);
        mpfr_init(r26054);
}

double f_dm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r26032, delta, MPFR_RNDN);
        mpfr_sin(r26033, r26032, MPFR_RNDN);
        mpfr_set_d(r26034, phi1, MPFR_RNDN);
        mpfr_cos(r26035, r26034, MPFR_RNDN);
        mpfr_set_d(r26036, theta, MPFR_RNDN);
        mpfr_sin(r26037, r26036, MPFR_RNDN);
        mpfr_mul(r26038, r26035, r26037, MPFR_RNDN);
        mpfr_mul(r26039, r26033, r26038, MPFR_RNDN);
        mpfr_sin(r26040, r26034, MPFR_RNDN);
        mpfr_neg(r26041, r26040, MPFR_RNDN);
        mpfr_cos(r26042, r26036, MPFR_RNDN);
        mpfr_mul(r26043, r26035, r26033, MPFR_RNDN);
        mpfr_cos(r26044, r26032, MPFR_RNDN);
        mpfr_mul(r26045, r26040, r26044, MPFR_RNDN);
        mpfr_fma(r26046, r26042, r26043, r26045, MPFR_RNDN);
        mpfr_asin(r26047, r26046, MPFR_RNDN);
        mpfr_sin(r26048, r26047, MPFR_RNDN);
        mpfr_fma(r26049, r26041, r26048, r26044, MPFR_RNDN);
        mpfr_exp(r26050, r26049, MPFR_RNDN);
        mpfr_log(r26051, r26050, MPFR_RNDN);
        mpfr_atan2(r26052, r26039, r26051, MPFR_RNDN);
        mpfr_set_d(r26053, lambda1, MPFR_RNDN);
        mpfr_add(r26054, r26052, r26053, MPFR_RNDN);
        return mpfr_get_d(r26054, MPFR_RNDN);
}