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

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


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

double f_od(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r26020 = lambda1;
        double r26021 = theta;
        double r26022 = sin(r26021);
        double r26023 = delta;
        double r26024 = sin(r26023);
        double r26025 = r26022 * r26024;
        double r26026 = phi1;
        double r26027 = cos(r26026);
        double r26028 = r26025 * r26027;
        double r26029 = cos(r26023);
        double r26030 = sin(r26026);
        double r26031 = r26029 * r26030;
        double r26032 = cos(r26021);
        double r26033 = r26032 * r26027;
        double r26034 = r26024 * r26033;
        double r26035 = r26031 + r26034;
        double r26036 = asin(r26035);
        double r26037 = sin(r26036);
        double r26038 = r26030 * r26037;
        double r26039 = r26029 - r26038;
        double r26040 = atan2(r26028, r26039);
        double r26041 = r26020 + r26040;
        return r26041;
}

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 r26042, r26043, r26044, r26045, r26046, r26047, r26048, r26049, r26050, r26051, r26052, r26053, r26054, r26055, r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        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);
        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);
}

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

static mpfr_t r26064, r26065, r26066, r26067, r26068, r26069, r26070, r26071, r26072, r26073, r26074, r26075, r26076, r26077, r26078, r26079, r26080, r26081, r26082, r26083, r26084, r26085;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26064);
        mpfr_init(r26065);
        mpfr_init(r26066);
        mpfr_init(r26067);
        mpfr_init(r26068);
        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(r26083);
        mpfr_init(r26084);
        mpfr_init(r26085);
}

double f_fm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r26064, lambda1, MPFR_RNDN);
        mpfr_set_d(r26065, theta, MPFR_RNDN);
        mpfr_sin(r26066, r26065, MPFR_RNDN);
        mpfr_set_d(r26067, delta, MPFR_RNDN);
        mpfr_sin(r26068, r26067, MPFR_RNDN);
        mpfr_mul(r26069, r26066, r26068, MPFR_RNDN);
        mpfr_set_d(r26070, phi1, MPFR_RNDN);
        mpfr_cos(r26071, r26070, MPFR_RNDN);
        mpfr_mul(r26072, r26069, r26071, MPFR_RNDN);
        mpfr_cos(r26073, r26067, MPFR_RNDN);
        mpfr_sin(r26074, r26070, MPFR_RNDN);
        mpfr_mul(r26075, r26073, r26074, MPFR_RNDN);
        mpfr_cos(r26076, r26065, MPFR_RNDN);
        mpfr_mul(r26077, r26076, r26071, MPFR_RNDN);
        mpfr_mul(r26078, r26068, r26077, MPFR_RNDN);
        mpfr_add(r26079, r26075, r26078, MPFR_RNDN);
        mpfr_asin(r26080, r26079, MPFR_RNDN);
        mpfr_sin(r26081, r26080, MPFR_RNDN);
        mpfr_mul(r26082, r26074, r26081, MPFR_RNDN);
        mpfr_sub(r26083, r26073, r26082, MPFR_RNDN);
        mpfr_atan2(r26084, r26072, r26083, MPFR_RNDN);
        mpfr_add(r26085, r26064, r26084, MPFR_RNDN);
        return mpfr_get_d(r26085, MPFR_RNDN);
}

static mpfr_t r26086, r26087, r26088, r26089, r26090, r26091, r26092, r26093, r26094, r26095, r26096, r26097, r26098, r26099, r26100, r26101, r26102, r26103, r26104, r26105, r26106, r26107;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26086);
        mpfr_init(r26087);
        mpfr_init(r26088);
        mpfr_init(r26089);
        mpfr_init(r26090);
        mpfr_init(r26091);
        mpfr_init(r26092);
        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(r26107);
}

double f_dm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r26086, lambda1, MPFR_RNDN);
        mpfr_set_d(r26087, theta, MPFR_RNDN);
        mpfr_sin(r26088, r26087, MPFR_RNDN);
        mpfr_set_d(r26089, delta, MPFR_RNDN);
        mpfr_sin(r26090, r26089, MPFR_RNDN);
        mpfr_mul(r26091, r26088, r26090, MPFR_RNDN);
        mpfr_set_d(r26092, phi1, MPFR_RNDN);
        mpfr_cos(r26093, r26092, MPFR_RNDN);
        mpfr_mul(r26094, r26091, r26093, MPFR_RNDN);
        mpfr_cos(r26095, r26089, MPFR_RNDN);
        mpfr_sin(r26096, r26092, MPFR_RNDN);
        mpfr_mul(r26097, r26095, r26096, MPFR_RNDN);
        mpfr_cos(r26098, r26087, MPFR_RNDN);
        mpfr_mul(r26099, r26098, r26093, MPFR_RNDN);
        mpfr_mul(r26100, r26090, r26099, MPFR_RNDN);
        mpfr_add(r26101, r26097, r26100, MPFR_RNDN);
        mpfr_asin(r26102, r26101, MPFR_RNDN);
        mpfr_sin(r26103, r26102, MPFR_RNDN);
        mpfr_mul(r26104, r26096, r26103, MPFR_RNDN);
        mpfr_sub(r26105, r26095, r26104, MPFR_RNDN);
        mpfr_atan2(r26106, r26094, r26105, MPFR_RNDN);
        mpfr_add(r26107, r26086, r26106, MPFR_RNDN);
        return mpfr_get_d(r26107, MPFR_RNDN);
}