#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 r26002 = R;
        float r26003 = 2;
        float r26004 = phi1;
        float r26005 = phi2;
        float r26006 = r26004 - r26005;
        float r26007 = r26006 / r26003;
        float r26008 = sin(r26007);
        float r26009 = pow(r26008, r26003);
        float r26010 = cos(r26004);
        float r26011 = cos(r26005);
        float r26012 = r26010 * r26011;
        float r26013 = lambda1;
        float r26014 = lambda2;
        float r26015 = r26013 - r26014;
        float r26016 = r26015 / r26003;
        float r26017 = sin(r26016);
        float r26018 = r26012 * r26017;
        float r26019 = r26018 * r26017;
        float r26020 = r26009 + r26019;
        float r26021 = sqrt(r26020);
        float r26022 = 1;
        float r26023 = r26022 - r26020;
        float r26024 = sqrt(r26023);
        float r26025 = atan2(r26021, r26024);
        float r26026 = r26003 * r26025;
        float r26027 = r26002 * r26026;
        return r26027;
}

double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) {
        double r26028 = R;
        double r26029 = 2;
        double r26030 = phi1;
        double r26031 = phi2;
        double r26032 = r26030 - r26031;
        double r26033 = r26032 / r26029;
        double r26034 = sin(r26033);
        double r26035 = pow(r26034, r26029);
        double r26036 = cos(r26030);
        double r26037 = cos(r26031);
        double r26038 = r26036 * r26037;
        double r26039 = lambda1;
        double r26040 = lambda2;
        double r26041 = r26039 - r26040;
        double r26042 = r26041 / r26029;
        double r26043 = sin(r26042);
        double r26044 = r26038 * r26043;
        double r26045 = r26044 * r26043;
        double r26046 = r26035 + r26045;
        double r26047 = sqrt(r26046);
        double r26048 = 1;
        double r26049 = r26048 - r26046;
        double r26050 = sqrt(r26049);
        double r26051 = atan2(r26047, r26050);
        double r26052 = r26029 * r26051;
        double r26053 = r26028 * r26052;
        return r26053;
}


double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) {
        float r26054 = R;
        float r26055 = 2;
        float r26056 = phi1;
        float r26057 = phi2;
        float r26058 = r26056 - r26057;
        float r26059 = r26058 / r26055;
        float r26060 = sin(r26059);
        float r26061 = pow(r26060, r26055);
        float r26062 = cos(r26056);
        float r26063 = cos(r26057);
        float r26064 = r26062 * r26063;
        float r26065 = lambda1;
        float r26066 = lambda2;
        float r26067 = r26065 - r26066;
        float r26068 = r26067 / r26055;
        float r26069 = sin(r26068);
        float r26070 = r26064 * r26069;
        float r26071 = log1p(r26069);
        float r26072 = expm1(r26071);
        float r26073 = r26070 * r26072;
        float r26074 = r26061 + r26073;
        float r26075 = sqrt(r26074);
        float r26076 = 1;
        float r26077 = r26070 * r26069;
        float r26078 = r26061 + r26077;
        float r26079 = r26076 - r26078;
        float r26080 = sqrt(r26079);
        float r26081 = atan2(r26075, r26080);
        float r26082 = r26055 * r26081;
        float r26083 = r26054 * r26082;
        return r26083;
}

double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) {
        double r26084 = R;
        double r26085 = 2;
        double r26086 = phi1;
        double r26087 = phi2;
        double r26088 = r26086 - r26087;
        double r26089 = r26088 / r26085;
        double r26090 = sin(r26089);
        double r26091 = pow(r26090, r26085);
        double r26092 = cos(r26086);
        double r26093 = cos(r26087);
        double r26094 = r26092 * r26093;
        double r26095 = lambda1;
        double r26096 = lambda2;
        double r26097 = r26095 - r26096;
        double r26098 = r26097 / r26085;
        double r26099 = sin(r26098);
        double r26100 = r26094 * r26099;
        double r26101 = log1p(r26099);
        double r26102 = expm1(r26101);
        double r26103 = r26100 * r26102;
        double r26104 = r26091 + r26103;
        double r26105 = sqrt(r26104);
        double r26106 = 1;
        double r26107 = r26100 * r26099;
        double r26108 = r26091 + r26107;
        double r26109 = r26106 - r26108;
        double r26110 = sqrt(r26109);
        double r26111 = atan2(r26105, r26110);
        double r26112 = r26085 * r26111;
        double r26113 = r26084 * r26112;
        return r26113;
}

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 r26114, r26115, r26116, r26117, r26118, r26119, r26120, r26121, r26122, r26123, r26124, r26125, r26126, r26127, r26128, r26129, r26130, r26131, r26132, r26133, r26134, r26135, r26136, r26137, r26138, r26139;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26114);
        mpfr_init_set_str(r26115, "2", 10, MPFR_RNDN);
        mpfr_init(r26116);
        mpfr_init(r26117);
        mpfr_init(r26118);
        mpfr_init(r26119);
        mpfr_init(r26120);
        mpfr_init(r26121);
        mpfr_init(r26122);
        mpfr_init(r26123);
        mpfr_init(r26124);
        mpfr_init(r26125);
        mpfr_init(r26126);
        mpfr_init(r26127);
        mpfr_init(r26128);
        mpfr_init(r26129);
        mpfr_init(r26130);
        mpfr_init(r26131);
        mpfr_init(r26132);
        mpfr_init(r26133);
        mpfr_init_set_str(r26134, "1", 10, MPFR_RNDN);
        mpfr_init(r26135);
        mpfr_init(r26136);
        mpfr_init(r26137);
        mpfr_init(r26138);
        mpfr_init(r26139);
}

double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r26114, R, MPFR_RNDN);
        ;
        mpfr_set_d(r26116, phi1, MPFR_RNDN);
        mpfr_set_d(r26117, phi2, MPFR_RNDN);
        mpfr_sub(r26118, r26116, r26117, MPFR_RNDN);
        mpfr_div(r26119, r26118, r26115, MPFR_RNDN);
        mpfr_sin(r26120, r26119, MPFR_RNDN);
        mpfr_pow(r26121, r26120, r26115, MPFR_RNDN);
        mpfr_cos(r26122, r26116, MPFR_RNDN);
        mpfr_cos(r26123, r26117, MPFR_RNDN);
        mpfr_mul(r26124, r26122, r26123, MPFR_RNDN);
        mpfr_set_d(r26125, lambda1, MPFR_RNDN);
        mpfr_set_d(r26126, lambda2, MPFR_RNDN);
        mpfr_sub(r26127, r26125, r26126, MPFR_RNDN);
        mpfr_div(r26128, r26127, r26115, MPFR_RNDN);
        mpfr_sin(r26129, r26128, MPFR_RNDN);
        mpfr_mul(r26130, r26124, r26129, MPFR_RNDN);
        mpfr_mul(r26131, r26130, r26129, MPFR_RNDN);
        mpfr_add(r26132, r26121, r26131, MPFR_RNDN);
        mpfr_sqrt(r26133, r26132, MPFR_RNDN);
        ;
        mpfr_sub(r26135, r26134, r26132, MPFR_RNDN);
        mpfr_sqrt(r26136, r26135, MPFR_RNDN);
        mpfr_atan2(r26137, r26133, r26136, MPFR_RNDN);
        mpfr_mul(r26138, r26115, r26137, MPFR_RNDN);
        mpfr_mul(r26139, r26114, r26138, MPFR_RNDN);
        return mpfr_get_d(r26139, MPFR_RNDN);
}

static mpfr_t r26140, r26141, r26142, r26143, r26144, r26145, r26146, r26147, r26148, r26149, r26150, r26151, r26152, r26153, r26154, r26155, r26156, r26157, r26158, r26159, r26160, r26161, r26162, r26163, r26164, r26165, r26166, r26167, r26168, r26169;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26140);
        mpfr_init_set_str(r26141, "2", 10, MPFR_RNDN);
        mpfr_init(r26142);
        mpfr_init(r26143);
        mpfr_init(r26144);
        mpfr_init(r26145);
        mpfr_init(r26146);
        mpfr_init(r26147);
        mpfr_init(r26148);
        mpfr_init(r26149);
        mpfr_init(r26150);
        mpfr_init(r26151);
        mpfr_init(r26152);
        mpfr_init(r26153);
        mpfr_init(r26154);
        mpfr_init(r26155);
        mpfr_init(r26156);
        mpfr_init(r26157);
        mpfr_init(r26158);
        mpfr_init(r26159);
        mpfr_init(r26160);
        mpfr_init(r26161);
        mpfr_init_set_str(r26162, "1", 10, MPFR_RNDN);
        mpfr_init(r26163);
        mpfr_init(r26164);
        mpfr_init(r26165);
        mpfr_init(r26166);
        mpfr_init(r26167);
        mpfr_init(r26168);
        mpfr_init(r26169);
}

double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r26140, R, MPFR_RNDN);
        ;
        mpfr_set_d(r26142, phi1, MPFR_RNDN);
        mpfr_set_d(r26143, phi2, MPFR_RNDN);
        mpfr_sub(r26144, r26142, r26143, MPFR_RNDN);
        mpfr_div(r26145, r26144, r26141, MPFR_RNDN);
        mpfr_sin(r26146, r26145, MPFR_RNDN);
        mpfr_pow(r26147, r26146, r26141, MPFR_RNDN);
        mpfr_cos(r26148, r26142, MPFR_RNDN);
        mpfr_cos(r26149, r26143, MPFR_RNDN);
        mpfr_mul(r26150, r26148, r26149, MPFR_RNDN);
        mpfr_set_d(r26151, lambda1, MPFR_RNDN);
        mpfr_set_d(r26152, lambda2, MPFR_RNDN);
        mpfr_sub(r26153, r26151, r26152, MPFR_RNDN);
        mpfr_div(r26154, r26153, r26141, MPFR_RNDN);
        mpfr_sin(r26155, r26154, MPFR_RNDN);
        mpfr_mul(r26156, r26150, r26155, MPFR_RNDN);
        mpfr_log1p(r26157, r26155, MPFR_RNDN);
        mpfr_expm1(r26158, r26157, MPFR_RNDN);
        mpfr_mul(r26159, r26156, r26158, MPFR_RNDN);
        mpfr_add(r26160, r26147, r26159, MPFR_RNDN);
        mpfr_sqrt(r26161, r26160, MPFR_RNDN);
        ;
        mpfr_mul(r26163, r26156, r26155, MPFR_RNDN);
        mpfr_add(r26164, r26147, r26163, MPFR_RNDN);
        mpfr_sub(r26165, r26162, r26164, MPFR_RNDN);
        mpfr_sqrt(r26166, r26165, MPFR_RNDN);
        mpfr_atan2(r26167, r26161, r26166, MPFR_RNDN);
        mpfr_mul(r26168, r26141, r26167, MPFR_RNDN);
        mpfr_mul(r26169, r26140, r26168, MPFR_RNDN);
        return mpfr_get_d(r26169, MPFR_RNDN);
}

static mpfr_t r26170, r26171, r26172, r26173, r26174, r26175, r26176, r26177, r26178, r26179, r26180, r26181, r26182, r26183, r26184, r26185, r26186, r26187, r26188, r26189, r26190, r26191, r26192, r26193, r26194, r26195, r26196, r26197, r26198, r26199;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26170);
        mpfr_init_set_str(r26171, "2", 10, MPFR_RNDN);
        mpfr_init(r26172);
        mpfr_init(r26173);
        mpfr_init(r26174);
        mpfr_init(r26175);
        mpfr_init(r26176);
        mpfr_init(r26177);
        mpfr_init(r26178);
        mpfr_init(r26179);
        mpfr_init(r26180);
        mpfr_init(r26181);
        mpfr_init(r26182);
        mpfr_init(r26183);
        mpfr_init(r26184);
        mpfr_init(r26185);
        mpfr_init(r26186);
        mpfr_init(r26187);
        mpfr_init(r26188);
        mpfr_init(r26189);
        mpfr_init(r26190);
        mpfr_init(r26191);
        mpfr_init_set_str(r26192, "1", 10, MPFR_RNDN);
        mpfr_init(r26193);
        mpfr_init(r26194);
        mpfr_init(r26195);
        mpfr_init(r26196);
        mpfr_init(r26197);
        mpfr_init(r26198);
        mpfr_init(r26199);
}

double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r26170, R, MPFR_RNDN);
        ;
        mpfr_set_d(r26172, phi1, MPFR_RNDN);
        mpfr_set_d(r26173, phi2, MPFR_RNDN);
        mpfr_sub(r26174, r26172, r26173, MPFR_RNDN);
        mpfr_div(r26175, r26174, r26171, MPFR_RNDN);
        mpfr_sin(r26176, r26175, MPFR_RNDN);
        mpfr_pow(r26177, r26176, r26171, MPFR_RNDN);
        mpfr_cos(r26178, r26172, MPFR_RNDN);
        mpfr_cos(r26179, r26173, MPFR_RNDN);
        mpfr_mul(r26180, r26178, r26179, MPFR_RNDN);
        mpfr_set_d(r26181, lambda1, MPFR_RNDN);
        mpfr_set_d(r26182, lambda2, MPFR_RNDN);
        mpfr_sub(r26183, r26181, r26182, MPFR_RNDN);
        mpfr_div(r26184, r26183, r26171, MPFR_RNDN);
        mpfr_sin(r26185, r26184, MPFR_RNDN);
        mpfr_mul(r26186, r26180, r26185, MPFR_RNDN);
        mpfr_log1p(r26187, r26185, MPFR_RNDN);
        mpfr_expm1(r26188, r26187, MPFR_RNDN);
        mpfr_mul(r26189, r26186, r26188, MPFR_RNDN);
        mpfr_add(r26190, r26177, r26189, MPFR_RNDN);
        mpfr_sqrt(r26191, r26190, MPFR_RNDN);
        ;
        mpfr_mul(r26193, r26186, r26185, MPFR_RNDN);
        mpfr_add(r26194, r26177, r26193, MPFR_RNDN);
        mpfr_sub(r26195, r26192, r26194, MPFR_RNDN);
        mpfr_sqrt(r26196, r26195, MPFR_RNDN);
        mpfr_atan2(r26197, r26191, r26196, MPFR_RNDN);
        mpfr_mul(r26198, r26171, r26197, MPFR_RNDN);
        mpfr_mul(r26199, r26170, r26198, MPFR_RNDN);
        return mpfr_get_d(r26199, MPFR_RNDN);
}