#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 r25955 = lambda1;
        float r25956 = theta;
        float r25957 = sin(r25956);
        float r25958 = delta;
        float r25959 = sin(r25958);
        float r25960 = r25957 * r25959;
        float r25961 = phi1;
        float r25962 = cos(r25961);
        float r25963 = r25960 * r25962;
        float r25964 = cos(r25958);
        float r25965 = sin(r25961);
        float r25966 = r25965 * r25964;
        float r25967 = r25962 * r25959;
        float r25968 = cos(r25956);
        float r25969 = r25967 * r25968;
        float r25970 = r25966 + r25969;
        float r25971 = asin(r25970);
        float r25972 = sin(r25971);
        float r25973 = r25965 * r25972;
        float r25974 = r25964 - r25973;
        float r25975 = atan2(r25963, r25974);
        float r25976 = r25955 + r25975;
        return r25976;
}

double f_id(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r25977 = lambda1;
        double r25978 = theta;
        double r25979 = sin(r25978);
        double r25980 = delta;
        double r25981 = sin(r25980);
        double r25982 = r25979 * r25981;
        double r25983 = phi1;
        double r25984 = cos(r25983);
        double r25985 = r25982 * r25984;
        double r25986 = cos(r25980);
        double r25987 = sin(r25983);
        double r25988 = r25987 * r25986;
        double r25989 = r25984 * r25981;
        double r25990 = cos(r25978);
        double r25991 = r25989 * r25990;
        double r25992 = r25988 + r25991;
        double r25993 = asin(r25992);
        double r25994 = sin(r25993);
        double r25995 = r25987 * r25994;
        double r25996 = r25986 - r25995;
        double r25997 = atan2(r25985, r25996);
        double r25998 = r25977 + r25997;
        return r25998;
}


double f_of(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) {
        float r25999 = lambda1;
        float r26000 = delta;
        float r26001 = sin(r26000);
        float r26002 = phi1;
        float r26003 = cos(r26002);
        float r26004 = theta;
        float r26005 = sin(r26004);
        float r26006 = r26003 * r26005;
        float r26007 = r26001 * r26006;
        float r26008 = r26003 * r26003;
        float r26009 = cos(r26000);
        float r26010 = r26008 * r26009;
        float r26011 = r26002 + r26002;
        float r26012 = sin(r26011);
        float r26013 = 0;
        float r26014 = r26012 + r26013;
        float r26015 = 2;
        float r26016 = r26014 / r26015;
        float r26017 = cos(r26004);
        float r26018 = r26017 * r26001;
        float r26019 = r26016 * r26018;
        float r26020 = r26010 - r26019;
        float r26021 = atan2(r26007, r26020);
        float r26022 = r25999 + r26021;
        return r26022;
}

double f_od(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r26023 = lambda1;
        double r26024 = delta;
        double r26025 = sin(r26024);
        double r26026 = phi1;
        double r26027 = cos(r26026);
        double r26028 = theta;
        double r26029 = sin(r26028);
        double r26030 = r26027 * r26029;
        double r26031 = r26025 * r26030;
        double r26032 = r26027 * r26027;
        double r26033 = cos(r26024);
        double r26034 = r26032 * r26033;
        double r26035 = r26026 + r26026;
        double r26036 = sin(r26035);
        double r26037 = 0;
        double r26038 = r26036 + r26037;
        double r26039 = 2;
        double r26040 = r26038 / r26039;
        double r26041 = cos(r26028);
        double r26042 = r26041 * r26025;
        double r26043 = r26040 * r26042;
        double r26044 = r26034 - r26043;
        double r26045 = atan2(r26031, r26044);
        double r26046 = r26023 + r26045;
        return r26046;
}

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 r26047, r26048, r26049, r26050, r26051, r26052, r26053, r26054, r26055, r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063, r26064, r26065, r26066, r26067, r26068;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26047);
        mpfr_init(r26048);
        mpfr_init(r26049);
        mpfr_init(r26050);
        mpfr_init(r26051);
        mpfr_init(r26052);
        mpfr_init(r26053);
        mpfr_init(r26054);
        mpfr_init(r26055);
        mpfr_init(r26056);
        mpfr_init(r26057);
        mpfr_init(r26058);
        mpfr_init(r26059);
        mpfr_init(r26060);
        mpfr_init(r26061);
        mpfr_init(r26062);
        mpfr_init(r26063);
        mpfr_init(r26064);
        mpfr_init(r26065);
        mpfr_init(r26066);
        mpfr_init(r26067);
        mpfr_init(r26068);
}

double f_im(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r26047, lambda1, MPFR_RNDN);
        mpfr_set_d(r26048, theta, MPFR_RNDN);
        mpfr_sin(r26049, r26048, MPFR_RNDN);
        mpfr_set_d(r26050, delta, MPFR_RNDN);
        mpfr_sin(r26051, r26050, MPFR_RNDN);
        mpfr_mul(r26052, r26049, r26051, MPFR_RNDN);
        mpfr_set_d(r26053, phi1, MPFR_RNDN);
        mpfr_cos(r26054, r26053, MPFR_RNDN);
        mpfr_mul(r26055, r26052, r26054, MPFR_RNDN);
        mpfr_cos(r26056, r26050, MPFR_RNDN);
        mpfr_sin(r26057, r26053, MPFR_RNDN);
        mpfr_mul(r26058, r26057, r26056, MPFR_RNDN);
        mpfr_mul(r26059, r26054, r26051, MPFR_RNDN);
        mpfr_cos(r26060, r26048, MPFR_RNDN);
        mpfr_mul(r26061, r26059, r26060, MPFR_RNDN);
        mpfr_add(r26062, r26058, r26061, MPFR_RNDN);
        mpfr_asin(r26063, r26062, MPFR_RNDN);
        mpfr_sin(r26064, r26063, MPFR_RNDN);
        mpfr_mul(r26065, r26057, r26064, MPFR_RNDN);
        mpfr_sub(r26066, r26056, r26065, MPFR_RNDN);
        mpfr_atan2(r26067, r26055, r26066, MPFR_RNDN);
        mpfr_add(r26068, r26047, r26067, MPFR_RNDN);
        return mpfr_get_d(r26068, MPFR_RNDN);
}

static mpfr_t r26069, r26070, r26071, r26072, r26073, r26074, r26075, r26076, r26077, r26078, r26079, r26080, r26081, r26082, r26083, r26084, r26085, r26086, r26087, r26088, r26089, r26090, r26091, r26092;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26069);
        mpfr_init(r26070);
        mpfr_init(r26071);
        mpfr_init(r26072);
        mpfr_init(r26073);
        mpfr_init(r26074);
        mpfr_init(r26075);
        mpfr_init(r26076);
        mpfr_init(r26077);
        mpfr_init(r26078);
        mpfr_init(r26079);
        mpfr_init(r26080);
        mpfr_init(r26081);
        mpfr_init(r26082);
        mpfr_init_set_str(r26083, "0", 10, MPFR_RNDN);
        mpfr_init(r26084);
        mpfr_init_set_str(r26085, "2", 10, MPFR_RNDN);
        mpfr_init(r26086);
        mpfr_init(r26087);
        mpfr_init(r26088);
        mpfr_init(r26089);
        mpfr_init(r26090);
        mpfr_init(r26091);
        mpfr_init(r26092);
}

double f_fm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r26069, lambda1, MPFR_RNDN);
        mpfr_set_d(r26070, delta, MPFR_RNDN);
        mpfr_sin(r26071, r26070, MPFR_RNDN);
        mpfr_set_d(r26072, phi1, MPFR_RNDN);
        mpfr_cos(r26073, r26072, MPFR_RNDN);
        mpfr_set_d(r26074, theta, MPFR_RNDN);
        mpfr_sin(r26075, r26074, MPFR_RNDN);
        mpfr_mul(r26076, r26073, r26075, MPFR_RNDN);
        mpfr_mul(r26077, r26071, r26076, MPFR_RNDN);
        mpfr_mul(r26078, r26073, r26073, MPFR_RNDN);
        mpfr_cos(r26079, r26070, MPFR_RNDN);
        mpfr_mul(r26080, r26078, r26079, MPFR_RNDN);
        mpfr_add(r26081, r26072, r26072, MPFR_RNDN);
        mpfr_sin(r26082, r26081, MPFR_RNDN);
        ;
        mpfr_add(r26084, r26082, r26083, MPFR_RNDN);
        ;
        mpfr_div(r26086, r26084, r26085, MPFR_RNDN);
        mpfr_cos(r26087, r26074, MPFR_RNDN);
        mpfr_mul(r26088, r26087, r26071, MPFR_RNDN);
        mpfr_mul(r26089, r26086, r26088, MPFR_RNDN);
        mpfr_sub(r26090, r26080, r26089, MPFR_RNDN);
        mpfr_atan2(r26091, r26077, r26090, MPFR_RNDN);
        mpfr_add(r26092, r26069, r26091, MPFR_RNDN);
        return mpfr_get_d(r26092, MPFR_RNDN);
}

static mpfr_t r26093, r26094, r26095, r26096, r26097, r26098, r26099, r26100, r26101, r26102, r26103, r26104, r26105, r26106, r26107, r26108, r26109, r26110, r26111, r26112, r26113, r26114, r26115, r26116;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26093);
        mpfr_init(r26094);
        mpfr_init(r26095);
        mpfr_init(r26096);
        mpfr_init(r26097);
        mpfr_init(r26098);
        mpfr_init(r26099);
        mpfr_init(r26100);
        mpfr_init(r26101);
        mpfr_init(r26102);
        mpfr_init(r26103);
        mpfr_init(r26104);
        mpfr_init(r26105);
        mpfr_init(r26106);
        mpfr_init_set_str(r26107, "0", 10, MPFR_RNDN);
        mpfr_init(r26108);
        mpfr_init_set_str(r26109, "2", 10, MPFR_RNDN);
        mpfr_init(r26110);
        mpfr_init(r26111);
        mpfr_init(r26112);
        mpfr_init(r26113);
        mpfr_init(r26114);
        mpfr_init(r26115);
        mpfr_init(r26116);
}

double f_dm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r26093, lambda1, MPFR_RNDN);
        mpfr_set_d(r26094, delta, MPFR_RNDN);
        mpfr_sin(r26095, r26094, MPFR_RNDN);
        mpfr_set_d(r26096, phi1, MPFR_RNDN);
        mpfr_cos(r26097, r26096, MPFR_RNDN);
        mpfr_set_d(r26098, theta, MPFR_RNDN);
        mpfr_sin(r26099, r26098, MPFR_RNDN);
        mpfr_mul(r26100, r26097, r26099, MPFR_RNDN);
        mpfr_mul(r26101, r26095, r26100, MPFR_RNDN);
        mpfr_mul(r26102, r26097, r26097, MPFR_RNDN);
        mpfr_cos(r26103, r26094, MPFR_RNDN);
        mpfr_mul(r26104, r26102, r26103, MPFR_RNDN);
        mpfr_add(r26105, r26096, r26096, MPFR_RNDN);
        mpfr_sin(r26106, r26105, MPFR_RNDN);
        ;
        mpfr_add(r26108, r26106, r26107, MPFR_RNDN);
        ;
        mpfr_div(r26110, r26108, r26109, MPFR_RNDN);
        mpfr_cos(r26111, r26098, MPFR_RNDN);
        mpfr_mul(r26112, r26111, r26095, MPFR_RNDN);
        mpfr_mul(r26113, r26110, r26112, MPFR_RNDN);
        mpfr_sub(r26114, r26104, r26113, MPFR_RNDN);
        mpfr_atan2(r26115, r26101, r26114, MPFR_RNDN);
        mpfr_add(r26116, r26093, r26115, MPFR_RNDN);
        return mpfr_get_d(r26116, MPFR_RNDN);
}