#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 r25674 = lambda1;
        float r25675 = theta;
        float r25676 = sin(r25675);
        float r25677 = delta;
        float r25678 = sin(r25677);
        float r25679 = r25676 * r25678;
        float r25680 = phi1;
        float r25681 = cos(r25680);
        float r25682 = r25679 * r25681;
        float r25683 = cos(r25677);
        float r25684 = sin(r25680);
        float r25685 = r25684 * r25683;
        float r25686 = r25681 * r25678;
        float r25687 = cos(r25675);
        float r25688 = r25686 * r25687;
        float r25689 = r25685 + r25688;
        float r25690 = asin(r25689);
        float r25691 = sin(r25690);
        float r25692 = r25684 * r25691;
        float r25693 = r25683 - r25692;
        float r25694 = atan2(r25682, r25693);
        float r25695 = r25674 + r25694;
        return r25695;
}

double f_id(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r25696 = lambda1;
        double r25697 = theta;
        double r25698 = sin(r25697);
        double r25699 = delta;
        double r25700 = sin(r25699);
        double r25701 = r25698 * r25700;
        double r25702 = phi1;
        double r25703 = cos(r25702);
        double r25704 = r25701 * r25703;
        double r25705 = cos(r25699);
        double r25706 = sin(r25702);
        double r25707 = r25706 * r25705;
        double r25708 = r25703 * r25700;
        double r25709 = cos(r25697);
        double r25710 = r25708 * r25709;
        double r25711 = r25707 + r25710;
        double r25712 = asin(r25711);
        double r25713 = sin(r25712);
        double r25714 = r25706 * r25713;
        double r25715 = r25705 - r25714;
        double r25716 = atan2(r25704, r25715);
        double r25717 = r25696 + r25716;
        return r25717;
}


double f_of(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) {
        float r25718 = delta;
        float r25719 = sin(r25718);
        float r25720 = phi1;
        float r25721 = cos(r25720);
        float r25722 = theta;
        float r25723 = sin(r25722);
        float r25724 = r25721 * r25723;
        float r25725 = r25719 * r25724;
        float r25726 = sin(r25720);
        float r25727 = -r25726;
        float r25728 = cos(r25722);
        float r25729 = r25721 * r25719;
        float r25730 = cos(r25718);
        float r25731 = r25726 * r25730;
        float r25732 = fma(r25728, r25729, r25731);
        float r25733 = asin(r25732);
        float r25734 = sin(r25733);
        float r25735 = fma(r25727, r25734, r25730);
        float r25736 = log1p(r25735);
        float r25737 = expm1(r25736);
        float r25738 = atan2(r25725, r25737);
        float r25739 = lambda1;
        float r25740 = r25738 + r25739;
        return r25740;
}

double f_od(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r25741 = delta;
        double r25742 = sin(r25741);
        double r25743 = phi1;
        double r25744 = cos(r25743);
        double r25745 = theta;
        double r25746 = sin(r25745);
        double r25747 = r25744 * r25746;
        double r25748 = r25742 * r25747;
        double r25749 = sin(r25743);
        double r25750 = -r25749;
        double r25751 = cos(r25745);
        double r25752 = r25744 * r25742;
        double r25753 = cos(r25741);
        double r25754 = r25749 * r25753;
        double r25755 = fma(r25751, r25752, r25754);
        double r25756 = asin(r25755);
        double r25757 = sin(r25756);
        double r25758 = fma(r25750, r25757, r25753);
        double r25759 = log1p(r25758);
        double r25760 = expm1(r25759);
        double r25761 = atan2(r25748, r25760);
        double r25762 = lambda1;
        double r25763 = r25761 + r25762;
        return r25763;
}

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 r25764, r25765, r25766, r25767, r25768, r25769, r25770, r25771, r25772, r25773, r25774, r25775, r25776, r25777, r25778, r25779, r25780, r25781, r25782, r25783, r25784, r25785;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r25764);
        mpfr_init(r25765);
        mpfr_init(r25766);
        mpfr_init(r25767);
        mpfr_init(r25768);
        mpfr_init(r25769);
        mpfr_init(r25770);
        mpfr_init(r25771);
        mpfr_init(r25772);
        mpfr_init(r25773);
        mpfr_init(r25774);
        mpfr_init(r25775);
        mpfr_init(r25776);
        mpfr_init(r25777);
        mpfr_init(r25778);
        mpfr_init(r25779);
        mpfr_init(r25780);
        mpfr_init(r25781);
        mpfr_init(r25782);
        mpfr_init(r25783);
        mpfr_init(r25784);
        mpfr_init(r25785);
}

double f_im(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r25764, lambda1, MPFR_RNDN);
        mpfr_set_d(r25765, theta, MPFR_RNDN);
        mpfr_sin(r25766, r25765, MPFR_RNDN);
        mpfr_set_d(r25767, delta, MPFR_RNDN);
        mpfr_sin(r25768, r25767, MPFR_RNDN);
        mpfr_mul(r25769, r25766, r25768, MPFR_RNDN);
        mpfr_set_d(r25770, phi1, MPFR_RNDN);
        mpfr_cos(r25771, r25770, MPFR_RNDN);
        mpfr_mul(r25772, r25769, r25771, MPFR_RNDN);
        mpfr_cos(r25773, r25767, MPFR_RNDN);
        mpfr_sin(r25774, r25770, MPFR_RNDN);
        mpfr_mul(r25775, r25774, r25773, MPFR_RNDN);
        mpfr_mul(r25776, r25771, r25768, MPFR_RNDN);
        mpfr_cos(r25777, r25765, MPFR_RNDN);
        mpfr_mul(r25778, r25776, r25777, MPFR_RNDN);
        mpfr_add(r25779, r25775, r25778, MPFR_RNDN);
        mpfr_asin(r25780, r25779, MPFR_RNDN);
        mpfr_sin(r25781, r25780, MPFR_RNDN);
        mpfr_mul(r25782, r25774, r25781, MPFR_RNDN);
        mpfr_sub(r25783, r25773, r25782, MPFR_RNDN);
        mpfr_atan2(r25784, r25772, r25783, MPFR_RNDN);
        mpfr_add(r25785, r25764, r25784, MPFR_RNDN);
        return mpfr_get_d(r25785, MPFR_RNDN);
}

static mpfr_t r25786, r25787, r25788, r25789, r25790, r25791, r25792, r25793, r25794, r25795, r25796, r25797, r25798, r25799, r25800, r25801, r25802, r25803, r25804, r25805, r25806, r25807, r25808;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25786);
        mpfr_init(r25787);
        mpfr_init(r25788);
        mpfr_init(r25789);
        mpfr_init(r25790);
        mpfr_init(r25791);
        mpfr_init(r25792);
        mpfr_init(r25793);
        mpfr_init(r25794);
        mpfr_init(r25795);
        mpfr_init(r25796);
        mpfr_init(r25797);
        mpfr_init(r25798);
        mpfr_init(r25799);
        mpfr_init(r25800);
        mpfr_init(r25801);
        mpfr_init(r25802);
        mpfr_init(r25803);
        mpfr_init(r25804);
        mpfr_init(r25805);
        mpfr_init(r25806);
        mpfr_init(r25807);
        mpfr_init(r25808);
}

double f_fm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r25786, delta, MPFR_RNDN);
        mpfr_sin(r25787, r25786, MPFR_RNDN);
        mpfr_set_d(r25788, phi1, MPFR_RNDN);
        mpfr_cos(r25789, r25788, MPFR_RNDN);
        mpfr_set_d(r25790, theta, MPFR_RNDN);
        mpfr_sin(r25791, r25790, MPFR_RNDN);
        mpfr_mul(r25792, r25789, r25791, MPFR_RNDN);
        mpfr_mul(r25793, r25787, r25792, MPFR_RNDN);
        mpfr_sin(r25794, r25788, MPFR_RNDN);
        mpfr_neg(r25795, r25794, MPFR_RNDN);
        mpfr_cos(r25796, r25790, MPFR_RNDN);
        mpfr_mul(r25797, r25789, r25787, MPFR_RNDN);
        mpfr_cos(r25798, r25786, MPFR_RNDN);
        mpfr_mul(r25799, r25794, r25798, MPFR_RNDN);
        mpfr_fma(r25800, r25796, r25797, r25799, MPFR_RNDN);
        mpfr_asin(r25801, r25800, MPFR_RNDN);
        mpfr_sin(r25802, r25801, MPFR_RNDN);
        mpfr_fma(r25803, r25795, r25802, r25798, MPFR_RNDN);
        mpfr_log1p(r25804, r25803, MPFR_RNDN);
        mpfr_expm1(r25805, r25804, MPFR_RNDN);
        mpfr_atan2(r25806, r25793, r25805, MPFR_RNDN);
        mpfr_set_d(r25807, lambda1, MPFR_RNDN);
        mpfr_add(r25808, r25806, r25807, MPFR_RNDN);
        return mpfr_get_d(r25808, MPFR_RNDN);
}

static mpfr_t r25809, r25810, r25811, r25812, r25813, r25814, r25815, r25816, r25817, r25818, r25819, r25820, r25821, r25822, r25823, r25824, r25825, r25826, r25827, r25828, r25829, r25830, r25831;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25809);
        mpfr_init(r25810);
        mpfr_init(r25811);
        mpfr_init(r25812);
        mpfr_init(r25813);
        mpfr_init(r25814);
        mpfr_init(r25815);
        mpfr_init(r25816);
        mpfr_init(r25817);
        mpfr_init(r25818);
        mpfr_init(r25819);
        mpfr_init(r25820);
        mpfr_init(r25821);
        mpfr_init(r25822);
        mpfr_init(r25823);
        mpfr_init(r25824);
        mpfr_init(r25825);
        mpfr_init(r25826);
        mpfr_init(r25827);
        mpfr_init(r25828);
        mpfr_init(r25829);
        mpfr_init(r25830);
        mpfr_init(r25831);
}

double f_dm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r25809, delta, MPFR_RNDN);
        mpfr_sin(r25810, r25809, MPFR_RNDN);
        mpfr_set_d(r25811, phi1, MPFR_RNDN);
        mpfr_cos(r25812, r25811, MPFR_RNDN);
        mpfr_set_d(r25813, theta, MPFR_RNDN);
        mpfr_sin(r25814, r25813, MPFR_RNDN);
        mpfr_mul(r25815, r25812, r25814, MPFR_RNDN);
        mpfr_mul(r25816, r25810, r25815, MPFR_RNDN);
        mpfr_sin(r25817, r25811, MPFR_RNDN);
        mpfr_neg(r25818, r25817, MPFR_RNDN);
        mpfr_cos(r25819, r25813, MPFR_RNDN);
        mpfr_mul(r25820, r25812, r25810, MPFR_RNDN);
        mpfr_cos(r25821, r25809, MPFR_RNDN);
        mpfr_mul(r25822, r25817, r25821, MPFR_RNDN);
        mpfr_fma(r25823, r25819, r25820, r25822, MPFR_RNDN);
        mpfr_asin(r25824, r25823, MPFR_RNDN);
        mpfr_sin(r25825, r25824, MPFR_RNDN);
        mpfr_fma(r25826, r25818, r25825, r25821, MPFR_RNDN);
        mpfr_log1p(r25827, r25826, MPFR_RNDN);
        mpfr_expm1(r25828, r25827, MPFR_RNDN);
        mpfr_atan2(r25829, r25816, r25828, MPFR_RNDN);
        mpfr_set_d(r25830, lambda1, MPFR_RNDN);
        mpfr_add(r25831, r25829, r25830, MPFR_RNDN);
        return mpfr_get_d(r25831, MPFR_RNDN);
}