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

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

double f_if(float lambda1, float lambda2, float phi1, float phi2) {
        float r28274 = lambda1;
        float r28275 = phi2;
        float r28276 = cos(r28275);
        float r28277 = lambda2;
        float r28278 = r28274 - r28277;
        float r28279 = sin(r28278);
        float r28280 = r28276 * r28279;
        float r28281 = phi1;
        float r28282 = cos(r28281);
        float r28283 = cos(r28278);
        float r28284 = r28276 * r28283;
        float r28285 = r28282 + r28284;
        float r28286 = atan2(r28280, r28285);
        float r28287 = r28274 + r28286;
        return r28287;
}

double f_id(double lambda1, double lambda2, double phi1, double phi2) {
        double r28288 = lambda1;
        double r28289 = phi2;
        double r28290 = cos(r28289);
        double r28291 = lambda2;
        double r28292 = r28288 - r28291;
        double r28293 = sin(r28292);
        double r28294 = r28290 * r28293;
        double r28295 = phi1;
        double r28296 = cos(r28295);
        double r28297 = cos(r28292);
        double r28298 = r28290 * r28297;
        double r28299 = r28296 + r28298;
        double r28300 = atan2(r28294, r28299);
        double r28301 = r28288 + r28300;
        return r28301;
}


double f_of(float lambda1, float lambda2, float phi1, float phi2) {
        float r28302 = lambda1;
        float r28303 = phi2;
        float r28304 = cos(r28303);
        float r28305 = lambda2;
        float r28306 = r28302 - r28305;
        float r28307 = sin(r28306);
        float r28308 = r28304 * r28307;
        float r28309 = phi1;
        float r28310 = cos(r28309);
        float r28311 = exp(r28310);
        float r28312 = cos(r28306);
        float r28313 = r28304 * r28312;
        float r28314 = exp(r28313);
        float r28315 = r28311 * r28314;
        float r28316 = log(r28315);
        float r28317 = atan2(r28308, r28316);
        float r28318 = r28302 + r28317;
        return r28318;
}

double f_od(double lambda1, double lambda2, double phi1, double phi2) {
        double r28319 = lambda1;
        double r28320 = phi2;
        double r28321 = cos(r28320);
        double r28322 = lambda2;
        double r28323 = r28319 - r28322;
        double r28324 = sin(r28323);
        double r28325 = r28321 * r28324;
        double r28326 = phi1;
        double r28327 = cos(r28326);
        double r28328 = exp(r28327);
        double r28329 = cos(r28323);
        double r28330 = r28321 * r28329;
        double r28331 = exp(r28330);
        double r28332 = r28328 * r28331;
        double r28333 = log(r28332);
        double r28334 = atan2(r28325, r28333);
        double r28335 = r28319 + r28334;
        return r28335;
}

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 r28336, r28337, r28338, r28339, r28340, r28341, r28342, r28343, r28344, r28345, r28346, r28347, r28348, r28349;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r28336);
        mpfr_init(r28337);
        mpfr_init(r28338);
        mpfr_init(r28339);
        mpfr_init(r28340);
        mpfr_init(r28341);
        mpfr_init(r28342);
        mpfr_init(r28343);
        mpfr_init(r28344);
        mpfr_init(r28345);
        mpfr_init(r28346);
        mpfr_init(r28347);
        mpfr_init(r28348);
        mpfr_init(r28349);
}

double f_im(double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r28336, lambda1, MPFR_RNDN);
        mpfr_set_d(r28337, phi2, MPFR_RNDN);
        mpfr_cos(r28338, r28337, MPFR_RNDN);
        mpfr_set_d(r28339, lambda2, MPFR_RNDN);
        mpfr_sub(r28340, r28336, r28339, MPFR_RNDN);
        mpfr_sin(r28341, r28340, MPFR_RNDN);
        mpfr_mul(r28342, r28338, r28341, MPFR_RNDN);
        mpfr_set_d(r28343, phi1, MPFR_RNDN);
        mpfr_cos(r28344, r28343, MPFR_RNDN);
        mpfr_cos(r28345, r28340, MPFR_RNDN);
        mpfr_mul(r28346, r28338, r28345, MPFR_RNDN);
        mpfr_add(r28347, r28344, r28346, MPFR_RNDN);
        mpfr_atan2(r28348, r28342, r28347, MPFR_RNDN);
        mpfr_add(r28349, r28336, r28348, MPFR_RNDN);
        return mpfr_get_d(r28349, MPFR_RNDN);
}

static mpfr_t r28350, r28351, r28352, r28353, r28354, r28355, r28356, r28357, r28358, r28359, r28360, r28361, r28362, r28363, r28364, r28365, r28366;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r28350);
        mpfr_init(r28351);
        mpfr_init(r28352);
        mpfr_init(r28353);
        mpfr_init(r28354);
        mpfr_init(r28355);
        mpfr_init(r28356);
        mpfr_init(r28357);
        mpfr_init(r28358);
        mpfr_init(r28359);
        mpfr_init(r28360);
        mpfr_init(r28361);
        mpfr_init(r28362);
        mpfr_init(r28363);
        mpfr_init(r28364);
        mpfr_init(r28365);
        mpfr_init(r28366);
}

double f_fm(double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r28350, lambda1, MPFR_RNDN);
        mpfr_set_d(r28351, phi2, MPFR_RNDN);
        mpfr_cos(r28352, r28351, MPFR_RNDN);
        mpfr_set_d(r28353, lambda2, MPFR_RNDN);
        mpfr_sub(r28354, r28350, r28353, MPFR_RNDN);
        mpfr_sin(r28355, r28354, MPFR_RNDN);
        mpfr_mul(r28356, r28352, r28355, MPFR_RNDN);
        mpfr_set_d(r28357, phi1, MPFR_RNDN);
        mpfr_cos(r28358, r28357, MPFR_RNDN);
        mpfr_exp(r28359, r28358, MPFR_RNDN);
        mpfr_cos(r28360, r28354, MPFR_RNDN);
        mpfr_mul(r28361, r28352, r28360, MPFR_RNDN);
        mpfr_exp(r28362, r28361, MPFR_RNDN);
        mpfr_mul(r28363, r28359, r28362, MPFR_RNDN);
        mpfr_log(r28364, r28363, MPFR_RNDN);
        mpfr_atan2(r28365, r28356, r28364, MPFR_RNDN);
        mpfr_add(r28366, r28350, r28365, MPFR_RNDN);
        return mpfr_get_d(r28366, MPFR_RNDN);
}

static mpfr_t r28367, r28368, r28369, r28370, r28371, r28372, r28373, r28374, r28375, r28376, r28377, r28378, r28379, r28380, r28381, r28382, r28383;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r28367);
        mpfr_init(r28368);
        mpfr_init(r28369);
        mpfr_init(r28370);
        mpfr_init(r28371);
        mpfr_init(r28372);
        mpfr_init(r28373);
        mpfr_init(r28374);
        mpfr_init(r28375);
        mpfr_init(r28376);
        mpfr_init(r28377);
        mpfr_init(r28378);
        mpfr_init(r28379);
        mpfr_init(r28380);
        mpfr_init(r28381);
        mpfr_init(r28382);
        mpfr_init(r28383);
}

double f_dm(double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r28367, lambda1, MPFR_RNDN);
        mpfr_set_d(r28368, phi2, MPFR_RNDN);
        mpfr_cos(r28369, r28368, MPFR_RNDN);
        mpfr_set_d(r28370, lambda2, MPFR_RNDN);
        mpfr_sub(r28371, r28367, r28370, MPFR_RNDN);
        mpfr_sin(r28372, r28371, MPFR_RNDN);
        mpfr_mul(r28373, r28369, r28372, MPFR_RNDN);
        mpfr_set_d(r28374, phi1, MPFR_RNDN);
        mpfr_cos(r28375, r28374, MPFR_RNDN);
        mpfr_exp(r28376, r28375, MPFR_RNDN);
        mpfr_cos(r28377, r28371, MPFR_RNDN);
        mpfr_mul(r28378, r28369, r28377, MPFR_RNDN);
        mpfr_exp(r28379, r28378, MPFR_RNDN);
        mpfr_mul(r28380, r28376, r28379, MPFR_RNDN);
        mpfr_log(r28381, r28380, MPFR_RNDN);
        mpfr_atan2(r28382, r28373, r28381, MPFR_RNDN);
        mpfr_add(r28383, r28367, r28382, MPFR_RNDN);
        return mpfr_get_d(r28383, MPFR_RNDN);
}