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

char *name = "Equirectangular approximation to distance on a great circle";

double f_if(float R, float lambda1, float lambda2, float phi1, float phi2) {
        float r22058 = R;
        float r22059 = lambda1;
        float r22060 = lambda2;
        float r22061 = r22059 - r22060;
        float r22062 = phi1;
        float r22063 = phi2;
        float r22064 = r22062 + r22063;
        float r22065 = 2;
        float r22066 = r22064 / r22065;
        float r22067 = cos(r22066);
        float r22068 = r22061 * r22067;
        float r22069 = r22068 * r22068;
        float r22070 = r22062 - r22063;
        float r22071 = r22070 * r22070;
        float r22072 = r22069 + r22071;
        float r22073 = sqrt(r22072);
        float r22074 = r22058 * r22073;
        return r22074;
}

double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) {
        double r22075 = R;
        double r22076 = lambda1;
        double r22077 = lambda2;
        double r22078 = r22076 - r22077;
        double r22079 = phi1;
        double r22080 = phi2;
        double r22081 = r22079 + r22080;
        double r22082 = 2;
        double r22083 = r22081 / r22082;
        double r22084 = cos(r22083);
        double r22085 = r22078 * r22084;
        double r22086 = r22085 * r22085;
        double r22087 = r22079 - r22080;
        double r22088 = r22087 * r22087;
        double r22089 = r22086 + r22088;
        double r22090 = sqrt(r22089);
        double r22091 = r22075 * r22090;
        return r22091;
}


double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) {
        float r22092 = phi2;
        float r22093 = phi1;
        float r22094 = r22092 - r22093;
        float r22095 = 7.79840277025698e+121;
        bool r22096 = r22094 <= r22095;
        float r22097 = R;
        float r22098 = lambda1;
        float r22099 = lambda2;
        float r22100 = r22098 - r22099;
        float r22101 = r22093 + r22092;
        float r22102 = 2;
        float r22103 = r22101 / r22102;
        float r22104 = cos(r22103);
        float r22105 = r22100 * r22104;
        float r22106 = r22105 * r22105;
        float r22107 = r22093 - r22092;
        float r22108 = r22107 * r22107;
        float r22109 = r22106 + r22108;
        float r22110 = sqrt(r22109);
        float r22111 = r22097 * r22110;
        float r22112 = r22097 * r22094;
        float r22113 = r22096 ? r22111 : r22112;
        return r22113;
}

double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) {
        double r22114 = phi2;
        double r22115 = phi1;
        double r22116 = r22114 - r22115;
        double r22117 = 7.79840277025698e+121;
        bool r22118 = r22116 <= r22117;
        double r22119 = R;
        double r22120 = lambda1;
        double r22121 = lambda2;
        double r22122 = r22120 - r22121;
        double r22123 = r22115 + r22114;
        double r22124 = 2;
        double r22125 = r22123 / r22124;
        double r22126 = cos(r22125);
        double r22127 = r22122 * r22126;
        double r22128 = r22127 * r22127;
        double r22129 = r22115 - r22114;
        double r22130 = r22129 * r22129;
        double r22131 = r22128 + r22130;
        double r22132 = sqrt(r22131);
        double r22133 = r22119 * r22132;
        double r22134 = r22119 * r22116;
        double r22135 = r22118 ? r22133 : r22134;
        return r22135;
}

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 r22136, r22137, r22138, r22139, r22140, r22141, r22142, r22143, r22144, r22145, r22146, r22147, r22148, r22149, r22150, r22151, r22152;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r22136);
        mpfr_init(r22137);
        mpfr_init(r22138);
        mpfr_init(r22139);
        mpfr_init(r22140);
        mpfr_init(r22141);
        mpfr_init(r22142);
        mpfr_init_set_str(r22143, "2", 10, MPFR_RNDN);
        mpfr_init(r22144);
        mpfr_init(r22145);
        mpfr_init(r22146);
        mpfr_init(r22147);
        mpfr_init(r22148);
        mpfr_init(r22149);
        mpfr_init(r22150);
        mpfr_init(r22151);
        mpfr_init(r22152);
}

double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r22136, R, MPFR_RNDN);
        mpfr_set_d(r22137, lambda1, MPFR_RNDN);
        mpfr_set_d(r22138, lambda2, MPFR_RNDN);
        mpfr_sub(r22139, r22137, r22138, MPFR_RNDN);
        mpfr_set_d(r22140, phi1, MPFR_RNDN);
        mpfr_set_d(r22141, phi2, MPFR_RNDN);
        mpfr_add(r22142, r22140, r22141, MPFR_RNDN);
        ;
        mpfr_div(r22144, r22142, r22143, MPFR_RNDN);
        mpfr_cos(r22145, r22144, MPFR_RNDN);
        mpfr_mul(r22146, r22139, r22145, MPFR_RNDN);
        mpfr_mul(r22147, r22146, r22146, MPFR_RNDN);
        mpfr_sub(r22148, r22140, r22141, MPFR_RNDN);
        mpfr_mul(r22149, r22148, r22148, MPFR_RNDN);
        mpfr_add(r22150, r22147, r22149, MPFR_RNDN);
        mpfr_sqrt(r22151, r22150, MPFR_RNDN);
        mpfr_mul(r22152, r22136, r22151, MPFR_RNDN);
        return mpfr_get_d(r22152, MPFR_RNDN);
}

static mpfr_t r22153, r22154, r22155, r22156, r22157, r22158, r22159, r22160, r22161, r22162, r22163, r22164, r22165, r22166, r22167, r22168, r22169, r22170, r22171, r22172, r22173, r22174;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r22153);
        mpfr_init(r22154);
        mpfr_init(r22155);
        mpfr_init_set_str(r22156, "7.79840277025698e+121", 10, MPFR_RNDN);
        mpfr_init(r22157);
        mpfr_init(r22158);
        mpfr_init(r22159);
        mpfr_init(r22160);
        mpfr_init(r22161);
        mpfr_init(r22162);
        mpfr_init_set_str(r22163, "2", 10, MPFR_RNDN);
        mpfr_init(r22164);
        mpfr_init(r22165);
        mpfr_init(r22166);
        mpfr_init(r22167);
        mpfr_init(r22168);
        mpfr_init(r22169);
        mpfr_init(r22170);
        mpfr_init(r22171);
        mpfr_init(r22172);
        mpfr_init(r22173);
        mpfr_init(r22174);
}

double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r22153, phi2, MPFR_RNDN);
        mpfr_set_d(r22154, phi1, MPFR_RNDN);
        mpfr_sub(r22155, r22153, r22154, MPFR_RNDN);
        ;
        mpfr_set_si(r22157, mpfr_cmp(r22155, r22156) <= 0, MPFR_RNDN);
        mpfr_set_d(r22158, R, MPFR_RNDN);
        mpfr_set_d(r22159, lambda1, MPFR_RNDN);
        mpfr_set_d(r22160, lambda2, MPFR_RNDN);
        mpfr_sub(r22161, r22159, r22160, MPFR_RNDN);
        mpfr_add(r22162, r22154, r22153, MPFR_RNDN);
        ;
        mpfr_div(r22164, r22162, r22163, MPFR_RNDN);
        mpfr_cos(r22165, r22164, MPFR_RNDN);
        mpfr_mul(r22166, r22161, r22165, MPFR_RNDN);
        mpfr_mul(r22167, r22166, r22166, MPFR_RNDN);
        mpfr_sub(r22168, r22154, r22153, MPFR_RNDN);
        mpfr_mul(r22169, r22168, r22168, MPFR_RNDN);
        mpfr_add(r22170, r22167, r22169, MPFR_RNDN);
        mpfr_sqrt(r22171, r22170, MPFR_RNDN);
        mpfr_mul(r22172, r22158, r22171, MPFR_RNDN);
        mpfr_mul(r22173, r22158, r22155, MPFR_RNDN);
        if (mpfr_get_si(r22157, MPFR_RNDN)) { mpfr_set(r22174, r22172, MPFR_RNDN); } else { mpfr_set(r22174, r22173, MPFR_RNDN); };
        return mpfr_get_d(r22174, MPFR_RNDN);
}

static mpfr_t r22175, r22176, r22177, r22178, r22179, r22180, r22181, r22182, r22183, r22184, r22185, r22186, r22187, r22188, r22189, r22190, r22191, r22192, r22193, r22194, r22195, r22196;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r22175);
        mpfr_init(r22176);
        mpfr_init(r22177);
        mpfr_init_set_str(r22178, "7.79840277025698e+121", 10, MPFR_RNDN);
        mpfr_init(r22179);
        mpfr_init(r22180);
        mpfr_init(r22181);
        mpfr_init(r22182);
        mpfr_init(r22183);
        mpfr_init(r22184);
        mpfr_init_set_str(r22185, "2", 10, MPFR_RNDN);
        mpfr_init(r22186);
        mpfr_init(r22187);
        mpfr_init(r22188);
        mpfr_init(r22189);
        mpfr_init(r22190);
        mpfr_init(r22191);
        mpfr_init(r22192);
        mpfr_init(r22193);
        mpfr_init(r22194);
        mpfr_init(r22195);
        mpfr_init(r22196);
}

double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r22175, phi2, MPFR_RNDN);
        mpfr_set_d(r22176, phi1, MPFR_RNDN);
        mpfr_sub(r22177, r22175, r22176, MPFR_RNDN);
        ;
        mpfr_set_si(r22179, mpfr_cmp(r22177, r22178) <= 0, MPFR_RNDN);
        mpfr_set_d(r22180, R, MPFR_RNDN);
        mpfr_set_d(r22181, lambda1, MPFR_RNDN);
        mpfr_set_d(r22182, lambda2, MPFR_RNDN);
        mpfr_sub(r22183, r22181, r22182, MPFR_RNDN);
        mpfr_add(r22184, r22176, r22175, MPFR_RNDN);
        ;
        mpfr_div(r22186, r22184, r22185, MPFR_RNDN);
        mpfr_cos(r22187, r22186, MPFR_RNDN);
        mpfr_mul(r22188, r22183, r22187, MPFR_RNDN);
        mpfr_mul(r22189, r22188, r22188, MPFR_RNDN);
        mpfr_sub(r22190, r22176, r22175, MPFR_RNDN);
        mpfr_mul(r22191, r22190, r22190, MPFR_RNDN);
        mpfr_add(r22192, r22189, r22191, MPFR_RNDN);
        mpfr_sqrt(r22193, r22192, MPFR_RNDN);
        mpfr_mul(r22194, r22180, r22193, MPFR_RNDN);
        mpfr_mul(r22195, r22180, r22177, MPFR_RNDN);
        if (mpfr_get_si(r22179, MPFR_RNDN)) { mpfr_set(r22196, r22194, MPFR_RNDN); } else { mpfr_set(r22196, r22195, MPFR_RNDN); };
        return mpfr_get_d(r22196, MPFR_RNDN);
}