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

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


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

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

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

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

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

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

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

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

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

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

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