#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Distance on a great circle";

double f_if(float R, float lambda1, float lambda2, float phi1, float phi2) {
        float r18968 = R;
        float r18969 = 2.0f;
        float r18970 = phi1;
        float r18971 = phi2;
        float r18972 = r18970 - r18971;
        float r18973 = r18972 / r18969;
        float r18974 = sin(r18973);
        float r18975 = r18974 * r18974;
        float r18976 = cos(r18970);
        float r18977 = cos(r18971);
        float r18978 = r18976 * r18977;
        float r18979 = lambda1;
        float r18980 = lambda2;
        float r18981 = r18979 - r18980;
        float r18982 = r18981 / r18969;
        float r18983 = sin(r18982);
        float r18984 = r18978 * r18983;
        float r18985 = r18984 * r18983;
        float r18986 = r18975 + r18985;
        float r18987 = sqrt(r18986);
        float r18988 = 1.0f;
        float r18989 = r18988 - r18986;
        float r18990 = sqrt(r18989);
        float r18991 = atan2(r18987, r18990);
        float r18992 = r18969 * r18991;
        float r18993 = r18968 * r18992;
        return r18993;
}

double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) {
        double r18994 = R;
        double r18995 = 2.0;
        double r18996 = phi1;
        double r18997 = phi2;
        double r18998 = r18996 - r18997;
        double r18999 = r18998 / r18995;
        double r19000 = sin(r18999);
        double r19001 = r19000 * r19000;
        double r19002 = cos(r18996);
        double r19003 = cos(r18997);
        double r19004 = r19002 * r19003;
        double r19005 = lambda1;
        double r19006 = lambda2;
        double r19007 = r19005 - r19006;
        double r19008 = r19007 / r18995;
        double r19009 = sin(r19008);
        double r19010 = r19004 * r19009;
        double r19011 = r19010 * r19009;
        double r19012 = r19001 + r19011;
        double r19013 = sqrt(r19012);
        double r19014 = 1.0;
        double r19015 = r19014 - r19012;
        double r19016 = sqrt(r19015);
        double r19017 = atan2(r19013, r19016);
        double r19018 = r18995 * r19017;
        double r19019 = r18994 * r19018;
        return r19019;
}


double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) {
        float r19020 = R;
        float r19021 = 2.0f;
        float r19022 = phi1;
        float r19023 = r19022 / r19021;
        float r19024 = sin(r19023);
        float r19025 = phi2;
        float r19026 = r19025 / r19021;
        float r19027 = cos(r19026);
        float r19028 = r19024 * r19027;
        float r19029 = cos(r19023);
        float r19030 = sin(r19026);
        float r19031 = r19029 * r19030;
        float r19032 = r19028 - r19031;
        float r19033 = r19032 * r19032;
        float r19034 = cos(r19022);
        float r19035 = cos(r19025);
        float r19036 = r19034 * r19035;
        float r19037 = lambda1;
        float r19038 = lambda2;
        float r19039 = r19037 - r19038;
        float r19040 = r19039 / r19021;
        float r19041 = sin(r19040);
        float r19042 = r19036 * r19041;
        float r19043 = r19042 * r19041;
        float r19044 = r19033 + r19043;
        float r19045 = sqrt(r19044);
        float r19046 = 1.0f;
        float r19047 = exp(r19041);
        float r19048 = log(r19047);
        float r19049 = r19042 * r19048;
        float r19050 = r19033 + r19049;
        float r19051 = r19046 - r19050;
        float r19052 = sqrt(r19051);
        float r19053 = atan2(r19045, r19052);
        float r19054 = r19021 * r19053;
        float r19055 = r19020 * r19054;
        return r19055;
}

double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) {
        double r19056 = R;
        double r19057 = 2.0;
        double r19058 = phi1;
        double r19059 = r19058 / r19057;
        double r19060 = sin(r19059);
        double r19061 = phi2;
        double r19062 = r19061 / r19057;
        double r19063 = cos(r19062);
        double r19064 = r19060 * r19063;
        double r19065 = cos(r19059);
        double r19066 = sin(r19062);
        double r19067 = r19065 * r19066;
        double r19068 = r19064 - r19067;
        double r19069 = r19068 * r19068;
        double r19070 = cos(r19058);
        double r19071 = cos(r19061);
        double r19072 = r19070 * r19071;
        double r19073 = lambda1;
        double r19074 = lambda2;
        double r19075 = r19073 - r19074;
        double r19076 = r19075 / r19057;
        double r19077 = sin(r19076);
        double r19078 = r19072 * r19077;
        double r19079 = r19078 * r19077;
        double r19080 = r19069 + r19079;
        double r19081 = sqrt(r19080);
        double r19082 = 1.0;
        double r19083 = exp(r19077);
        double r19084 = log(r19083);
        double r19085 = r19078 * r19084;
        double r19086 = r19069 + r19085;
        double r19087 = r19082 - r19086;
        double r19088 = sqrt(r19087);
        double r19089 = atan2(r19081, r19088);
        double r19090 = r19057 * r19089;
        double r19091 = r19056 * r19090;
        return r19091;
}

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 r19092, r19093, r19094, r19095, r19096, r19097, r19098, r19099, r19100, r19101, r19102, r19103, r19104, r19105, r19106, r19107, r19108, r19109, r19110, r19111, r19112, r19113, r19114, r19115, r19116, r19117;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r19092);
        mpfr_init_set_str(r19093, "2", 10, MPFR_RNDN);
        mpfr_init(r19094);
        mpfr_init(r19095);
        mpfr_init(r19096);
        mpfr_init(r19097);
        mpfr_init(r19098);
        mpfr_init(r19099);
        mpfr_init(r19100);
        mpfr_init(r19101);
        mpfr_init(r19102);
        mpfr_init(r19103);
        mpfr_init(r19104);
        mpfr_init(r19105);
        mpfr_init(r19106);
        mpfr_init(r19107);
        mpfr_init(r19108);
        mpfr_init(r19109);
        mpfr_init(r19110);
        mpfr_init(r19111);
        mpfr_init_set_str(r19112, "1", 10, MPFR_RNDN);
        mpfr_init(r19113);
        mpfr_init(r19114);
        mpfr_init(r19115);
        mpfr_init(r19116);
        mpfr_init(r19117);
}

double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r19092, R, MPFR_RNDN);
        ;
        mpfr_set_d(r19094, phi1, MPFR_RNDN);
        mpfr_set_d(r19095, phi2, MPFR_RNDN);
        mpfr_sub(r19096, r19094, r19095, MPFR_RNDN);
        mpfr_div(r19097, r19096, r19093, MPFR_RNDN);
        mpfr_sin(r19098, r19097, MPFR_RNDN);
        mpfr_sqr(r19099, r19098, MPFR_RNDN);
        mpfr_cos(r19100, r19094, MPFR_RNDN);
        mpfr_cos(r19101, r19095, MPFR_RNDN);
        mpfr_mul(r19102, r19100, r19101, MPFR_RNDN);
        mpfr_set_d(r19103, lambda1, MPFR_RNDN);
        mpfr_set_d(r19104, lambda2, MPFR_RNDN);
        mpfr_sub(r19105, r19103, r19104, MPFR_RNDN);
        mpfr_div(r19106, r19105, r19093, MPFR_RNDN);
        mpfr_sin(r19107, r19106, MPFR_RNDN);
        mpfr_mul(r19108, r19102, r19107, MPFR_RNDN);
        mpfr_mul(r19109, r19108, r19107, MPFR_RNDN);
        mpfr_add(r19110, r19099, r19109, MPFR_RNDN);
        mpfr_sqrt(r19111, r19110, MPFR_RNDN);
        ;
        mpfr_sub(r19113, r19112, r19110, MPFR_RNDN);
        mpfr_sqrt(r19114, r19113, MPFR_RNDN);
        mpfr_atan2(r19115, r19111, r19114, MPFR_RNDN);
        mpfr_mul(r19116, r19093, r19115, MPFR_RNDN);
        mpfr_mul(r19117, r19092, r19116, MPFR_RNDN);
        return mpfr_get_d(r19117, MPFR_RNDN);
}

static mpfr_t r19118, r19119, r19120, r19121, r19122, r19123, r19124, r19125, r19126, r19127, r19128, r19129, r19130, r19131, r19132, r19133, r19134, r19135, r19136, r19137, r19138, r19139, r19140, r19141, r19142, r19143, r19144, r19145, r19146, r19147, r19148, r19149, r19150, r19151, r19152, r19153;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r19118);
        mpfr_init_set_str(r19119, "2", 10, MPFR_RNDN);
        mpfr_init(r19120);
        mpfr_init(r19121);
        mpfr_init(r19122);
        mpfr_init(r19123);
        mpfr_init(r19124);
        mpfr_init(r19125);
        mpfr_init(r19126);
        mpfr_init(r19127);
        mpfr_init(r19128);
        mpfr_init(r19129);
        mpfr_init(r19130);
        mpfr_init(r19131);
        mpfr_init(r19132);
        mpfr_init(r19133);
        mpfr_init(r19134);
        mpfr_init(r19135);
        mpfr_init(r19136);
        mpfr_init(r19137);
        mpfr_init(r19138);
        mpfr_init(r19139);
        mpfr_init(r19140);
        mpfr_init(r19141);
        mpfr_init(r19142);
        mpfr_init(r19143);
        mpfr_init_set_str(r19144, "1", 10, MPFR_RNDN);
        mpfr_init(r19145);
        mpfr_init(r19146);
        mpfr_init(r19147);
        mpfr_init(r19148);
        mpfr_init(r19149);
        mpfr_init(r19150);
        mpfr_init(r19151);
        mpfr_init(r19152);
        mpfr_init(r19153);
}

double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r19118, R, MPFR_RNDN);
        ;
        mpfr_set_d(r19120, phi1, MPFR_RNDN);
        mpfr_div(r19121, r19120, r19119, MPFR_RNDN);
        mpfr_sin(r19122, r19121, MPFR_RNDN);
        mpfr_set_d(r19123, phi2, MPFR_RNDN);
        mpfr_div(r19124, r19123, r19119, MPFR_RNDN);
        mpfr_cos(r19125, r19124, MPFR_RNDN);
        mpfr_mul(r19126, r19122, r19125, MPFR_RNDN);
        mpfr_cos(r19127, r19121, MPFR_RNDN);
        mpfr_sin(r19128, r19124, MPFR_RNDN);
        mpfr_mul(r19129, r19127, r19128, MPFR_RNDN);
        mpfr_sub(r19130, r19126, r19129, MPFR_RNDN);
        mpfr_sqr(r19131, r19130, MPFR_RNDN);
        mpfr_cos(r19132, r19120, MPFR_RNDN);
        mpfr_cos(r19133, r19123, MPFR_RNDN);
        mpfr_mul(r19134, r19132, r19133, MPFR_RNDN);
        mpfr_set_d(r19135, lambda1, MPFR_RNDN);
        mpfr_set_d(r19136, lambda2, MPFR_RNDN);
        mpfr_sub(r19137, r19135, r19136, MPFR_RNDN);
        mpfr_div(r19138, r19137, r19119, MPFR_RNDN);
        mpfr_sin(r19139, r19138, MPFR_RNDN);
        mpfr_mul(r19140, r19134, r19139, MPFR_RNDN);
        mpfr_mul(r19141, r19140, r19139, MPFR_RNDN);
        mpfr_add(r19142, r19131, r19141, MPFR_RNDN);
        mpfr_sqrt(r19143, r19142, MPFR_RNDN);
        ;
        mpfr_exp(r19145, r19139, MPFR_RNDN);
        mpfr_log(r19146, r19145, MPFR_RNDN);
        mpfr_mul(r19147, r19140, r19146, MPFR_RNDN);
        mpfr_add(r19148, r19131, r19147, MPFR_RNDN);
        mpfr_sub(r19149, r19144, r19148, MPFR_RNDN);
        mpfr_sqrt(r19150, r19149, MPFR_RNDN);
        mpfr_atan2(r19151, r19143, r19150, MPFR_RNDN);
        mpfr_mul(r19152, r19119, r19151, MPFR_RNDN);
        mpfr_mul(r19153, r19118, r19152, MPFR_RNDN);
        return mpfr_get_d(r19153, MPFR_RNDN);
}

static mpfr_t r19154, r19155, r19156, r19157, r19158, r19159, r19160, r19161, r19162, r19163, r19164, r19165, r19166, r19167, r19168, r19169, r19170, r19171, r19172, r19173, r19174, r19175, r19176, r19177, r19178, r19179, r19180, r19181, r19182, r19183, r19184, r19185, r19186, r19187, r19188, r19189;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r19154);
        mpfr_init_set_str(r19155, "2", 10, MPFR_RNDN);
        mpfr_init(r19156);
        mpfr_init(r19157);
        mpfr_init(r19158);
        mpfr_init(r19159);
        mpfr_init(r19160);
        mpfr_init(r19161);
        mpfr_init(r19162);
        mpfr_init(r19163);
        mpfr_init(r19164);
        mpfr_init(r19165);
        mpfr_init(r19166);
        mpfr_init(r19167);
        mpfr_init(r19168);
        mpfr_init(r19169);
        mpfr_init(r19170);
        mpfr_init(r19171);
        mpfr_init(r19172);
        mpfr_init(r19173);
        mpfr_init(r19174);
        mpfr_init(r19175);
        mpfr_init(r19176);
        mpfr_init(r19177);
        mpfr_init(r19178);
        mpfr_init(r19179);
        mpfr_init_set_str(r19180, "1", 10, MPFR_RNDN);
        mpfr_init(r19181);
        mpfr_init(r19182);
        mpfr_init(r19183);
        mpfr_init(r19184);
        mpfr_init(r19185);
        mpfr_init(r19186);
        mpfr_init(r19187);
        mpfr_init(r19188);
        mpfr_init(r19189);
}

double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r19154, R, MPFR_RNDN);
        ;
        mpfr_set_d(r19156, phi1, MPFR_RNDN);
        mpfr_div(r19157, r19156, r19155, MPFR_RNDN);
        mpfr_sin(r19158, r19157, MPFR_RNDN);
        mpfr_set_d(r19159, phi2, MPFR_RNDN);
        mpfr_div(r19160, r19159, r19155, MPFR_RNDN);
        mpfr_cos(r19161, r19160, MPFR_RNDN);
        mpfr_mul(r19162, r19158, r19161, MPFR_RNDN);
        mpfr_cos(r19163, r19157, MPFR_RNDN);
        mpfr_sin(r19164, r19160, MPFR_RNDN);
        mpfr_mul(r19165, r19163, r19164, MPFR_RNDN);
        mpfr_sub(r19166, r19162, r19165, MPFR_RNDN);
        mpfr_sqr(r19167, r19166, MPFR_RNDN);
        mpfr_cos(r19168, r19156, MPFR_RNDN);
        mpfr_cos(r19169, r19159, MPFR_RNDN);
        mpfr_mul(r19170, r19168, r19169, MPFR_RNDN);
        mpfr_set_d(r19171, lambda1, MPFR_RNDN);
        mpfr_set_d(r19172, lambda2, MPFR_RNDN);
        mpfr_sub(r19173, r19171, r19172, MPFR_RNDN);
        mpfr_div(r19174, r19173, r19155, MPFR_RNDN);
        mpfr_sin(r19175, r19174, MPFR_RNDN);
        mpfr_mul(r19176, r19170, r19175, MPFR_RNDN);
        mpfr_mul(r19177, r19176, r19175, MPFR_RNDN);
        mpfr_add(r19178, r19167, r19177, MPFR_RNDN);
        mpfr_sqrt(r19179, r19178, MPFR_RNDN);
        ;
        mpfr_exp(r19181, r19175, MPFR_RNDN);
        mpfr_log(r19182, r19181, MPFR_RNDN);
        mpfr_mul(r19183, r19176, r19182, MPFR_RNDN);
        mpfr_add(r19184, r19167, r19183, MPFR_RNDN);
        mpfr_sub(r19185, r19180, r19184, MPFR_RNDN);
        mpfr_sqrt(r19186, r19185, MPFR_RNDN);
        mpfr_atan2(r19187, r19179, r19186, MPFR_RNDN);
        mpfr_mul(r19188, r19155, r19187, MPFR_RNDN);
        mpfr_mul(r19189, r19154, r19188, MPFR_RNDN);
        return mpfr_get_d(r19189, MPFR_RNDN);
}