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

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

double f_if(float lambda1, float lambda2, float phi1, float phi2) {
        float r26285 = lambda1;
        float r26286 = lambda2;
        float r26287 = r26285 - r26286;
        float r26288 = sin(r26287);
        float r26289 = phi2;
        float r26290 = cos(r26289);
        float r26291 = r26288 * r26290;
        float r26292 = phi1;
        float r26293 = cos(r26292);
        float r26294 = sin(r26289);
        float r26295 = r26293 * r26294;
        float r26296 = sin(r26292);
        float r26297 = r26296 * r26290;
        float r26298 = cos(r26287);
        float r26299 = r26297 * r26298;
        float r26300 = r26295 - r26299;
        float r26301 = atan2(r26291, r26300);
        return r26301;
}

double f_id(double lambda1, double lambda2, double phi1, double phi2) {
        double r26302 = lambda1;
        double r26303 = lambda2;
        double r26304 = r26302 - r26303;
        double r26305 = sin(r26304);
        double r26306 = phi2;
        double r26307 = cos(r26306);
        double r26308 = r26305 * r26307;
        double r26309 = phi1;
        double r26310 = cos(r26309);
        double r26311 = sin(r26306);
        double r26312 = r26310 * r26311;
        double r26313 = sin(r26309);
        double r26314 = r26313 * r26307;
        double r26315 = cos(r26304);
        double r26316 = r26314 * r26315;
        double r26317 = r26312 - r26316;
        double r26318 = atan2(r26308, r26317);
        return r26318;
}


double f_of(float lambda1, float lambda2, float phi1, float phi2) {
        float r26319 = lambda1;
        float r26320 = sin(r26319);
        float r26321 = lambda2;
        float r26322 = cos(r26321);
        float r26323 = r26320 * r26322;
        float r26324 = cos(r26319);
        float r26325 = -r26321;
        float r26326 = sin(r26325);
        float r26327 = r26324 * r26326;
        float r26328 = r26323 + r26327;
        float r26329 = phi2;
        float r26330 = cos(r26329);
        float r26331 = r26328 * r26330;
        float r26332 = phi1;
        float r26333 = cos(r26332);
        float r26334 = sin(r26329);
        float r26335 = r26333 * r26334;
        float r26336 = sin(r26332);
        float r26337 = r26336 * r26330;
        float r26338 = r26324 * r26322;
        float r26339 = r26337 * r26338;
        float r26340 = sin(r26321);
        float r26341 = r26320 * r26340;
        float r26342 = r26337 * r26341;
        float r26343 = exp(r26342);
        float r26344 = log(r26343);
        float r26345 = r26339 + r26344;
        float r26346 = r26335 - r26345;
        float r26347 = atan2(r26331, r26346);
        return r26347;
}

double f_od(double lambda1, double lambda2, double phi1, double phi2) {
        double r26348 = lambda1;
        double r26349 = sin(r26348);
        double r26350 = lambda2;
        double r26351 = cos(r26350);
        double r26352 = r26349 * r26351;
        double r26353 = cos(r26348);
        double r26354 = -r26350;
        double r26355 = sin(r26354);
        double r26356 = r26353 * r26355;
        double r26357 = r26352 + r26356;
        double r26358 = phi2;
        double r26359 = cos(r26358);
        double r26360 = r26357 * r26359;
        double r26361 = phi1;
        double r26362 = cos(r26361);
        double r26363 = sin(r26358);
        double r26364 = r26362 * r26363;
        double r26365 = sin(r26361);
        double r26366 = r26365 * r26359;
        double r26367 = r26353 * r26351;
        double r26368 = r26366 * r26367;
        double r26369 = sin(r26350);
        double r26370 = r26349 * r26369;
        double r26371 = r26366 * r26370;
        double r26372 = exp(r26371);
        double r26373 = log(r26372);
        double r26374 = r26368 + r26373;
        double r26375 = r26364 - r26374;
        double r26376 = atan2(r26360, r26375);
        return r26376;
}

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 r26377, r26378, r26379, r26380, r26381, r26382, r26383, r26384, r26385, r26386, r26387, r26388, r26389, r26390, r26391, r26392, r26393;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26377);
        mpfr_init(r26378);
        mpfr_init(r26379);
        mpfr_init(r26380);
        mpfr_init(r26381);
        mpfr_init(r26382);
        mpfr_init(r26383);
        mpfr_init(r26384);
        mpfr_init(r26385);
        mpfr_init(r26386);
        mpfr_init(r26387);
        mpfr_init(r26388);
        mpfr_init(r26389);
        mpfr_init(r26390);
        mpfr_init(r26391);
        mpfr_init(r26392);
        mpfr_init(r26393);
}

double f_im(double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r26377, lambda1, MPFR_RNDN);
        mpfr_set_d(r26378, lambda2, MPFR_RNDN);
        mpfr_sub(r26379, r26377, r26378, MPFR_RNDN);
        mpfr_sin(r26380, r26379, MPFR_RNDN);
        mpfr_set_d(r26381, phi2, MPFR_RNDN);
        mpfr_cos(r26382, r26381, MPFR_RNDN);
        mpfr_mul(r26383, r26380, r26382, MPFR_RNDN);
        mpfr_set_d(r26384, phi1, MPFR_RNDN);
        mpfr_cos(r26385, r26384, MPFR_RNDN);
        mpfr_sin(r26386, r26381, MPFR_RNDN);
        mpfr_mul(r26387, r26385, r26386, MPFR_RNDN);
        mpfr_sin(r26388, r26384, MPFR_RNDN);
        mpfr_mul(r26389, r26388, r26382, MPFR_RNDN);
        mpfr_cos(r26390, r26379, MPFR_RNDN);
        mpfr_mul(r26391, r26389, r26390, MPFR_RNDN);
        mpfr_sub(r26392, r26387, r26391, MPFR_RNDN);
        mpfr_atan2(r26393, r26383, r26392, MPFR_RNDN);
        return mpfr_get_d(r26393, MPFR_RNDN);
}

static mpfr_t r26394, r26395, r26396, r26397, r26398, r26399, r26400, r26401, r26402, r26403, r26404, r26405, r26406, r26407, r26408, r26409, r26410, r26411, r26412, r26413, r26414, r26415, r26416, r26417, r26418, r26419, r26420, r26421, r26422;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26394);
        mpfr_init(r26395);
        mpfr_init(r26396);
        mpfr_init(r26397);
        mpfr_init(r26398);
        mpfr_init(r26399);
        mpfr_init(r26400);
        mpfr_init(r26401);
        mpfr_init(r26402);
        mpfr_init(r26403);
        mpfr_init(r26404);
        mpfr_init(r26405);
        mpfr_init(r26406);
        mpfr_init(r26407);
        mpfr_init(r26408);
        mpfr_init(r26409);
        mpfr_init(r26410);
        mpfr_init(r26411);
        mpfr_init(r26412);
        mpfr_init(r26413);
        mpfr_init(r26414);
        mpfr_init(r26415);
        mpfr_init(r26416);
        mpfr_init(r26417);
        mpfr_init(r26418);
        mpfr_init(r26419);
        mpfr_init(r26420);
        mpfr_init(r26421);
        mpfr_init(r26422);
}

double f_fm(double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r26394, lambda1, MPFR_RNDN);
        mpfr_sin(r26395, r26394, MPFR_RNDN);
        mpfr_set_d(r26396, lambda2, MPFR_RNDN);
        mpfr_cos(r26397, r26396, MPFR_RNDN);
        mpfr_mul(r26398, r26395, r26397, MPFR_RNDN);
        mpfr_cos(r26399, r26394, MPFR_RNDN);
        mpfr_neg(r26400, r26396, MPFR_RNDN);
        mpfr_sin(r26401, r26400, MPFR_RNDN);
        mpfr_mul(r26402, r26399, r26401, MPFR_RNDN);
        mpfr_add(r26403, r26398, r26402, MPFR_RNDN);
        mpfr_set_d(r26404, phi2, MPFR_RNDN);
        mpfr_cos(r26405, r26404, MPFR_RNDN);
        mpfr_mul(r26406, r26403, r26405, MPFR_RNDN);
        mpfr_set_d(r26407, phi1, MPFR_RNDN);
        mpfr_cos(r26408, r26407, MPFR_RNDN);
        mpfr_sin(r26409, r26404, MPFR_RNDN);
        mpfr_mul(r26410, r26408, r26409, MPFR_RNDN);
        mpfr_sin(r26411, r26407, MPFR_RNDN);
        mpfr_mul(r26412, r26411, r26405, MPFR_RNDN);
        mpfr_mul(r26413, r26399, r26397, MPFR_RNDN);
        mpfr_mul(r26414, r26412, r26413, MPFR_RNDN);
        mpfr_sin(r26415, r26396, MPFR_RNDN);
        mpfr_mul(r26416, r26395, r26415, MPFR_RNDN);
        mpfr_mul(r26417, r26412, r26416, MPFR_RNDN);
        mpfr_exp(r26418, r26417, MPFR_RNDN);
        mpfr_log(r26419, r26418, MPFR_RNDN);
        mpfr_add(r26420, r26414, r26419, MPFR_RNDN);
        mpfr_sub(r26421, r26410, r26420, MPFR_RNDN);
        mpfr_atan2(r26422, r26406, r26421, MPFR_RNDN);
        return mpfr_get_d(r26422, MPFR_RNDN);
}

static mpfr_t r26423, r26424, r26425, r26426, r26427, r26428, r26429, r26430, r26431, r26432, r26433, r26434, r26435, r26436, r26437, r26438, r26439, r26440, r26441, r26442, r26443, r26444, r26445, r26446, r26447, r26448, r26449, r26450, r26451;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26423);
        mpfr_init(r26424);
        mpfr_init(r26425);
        mpfr_init(r26426);
        mpfr_init(r26427);
        mpfr_init(r26428);
        mpfr_init(r26429);
        mpfr_init(r26430);
        mpfr_init(r26431);
        mpfr_init(r26432);
        mpfr_init(r26433);
        mpfr_init(r26434);
        mpfr_init(r26435);
        mpfr_init(r26436);
        mpfr_init(r26437);
        mpfr_init(r26438);
        mpfr_init(r26439);
        mpfr_init(r26440);
        mpfr_init(r26441);
        mpfr_init(r26442);
        mpfr_init(r26443);
        mpfr_init(r26444);
        mpfr_init(r26445);
        mpfr_init(r26446);
        mpfr_init(r26447);
        mpfr_init(r26448);
        mpfr_init(r26449);
        mpfr_init(r26450);
        mpfr_init(r26451);
}

double f_dm(double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r26423, lambda1, MPFR_RNDN);
        mpfr_sin(r26424, r26423, MPFR_RNDN);
        mpfr_set_d(r26425, lambda2, MPFR_RNDN);
        mpfr_cos(r26426, r26425, MPFR_RNDN);
        mpfr_mul(r26427, r26424, r26426, MPFR_RNDN);
        mpfr_cos(r26428, r26423, MPFR_RNDN);
        mpfr_neg(r26429, r26425, MPFR_RNDN);
        mpfr_sin(r26430, r26429, MPFR_RNDN);
        mpfr_mul(r26431, r26428, r26430, MPFR_RNDN);
        mpfr_add(r26432, r26427, r26431, MPFR_RNDN);
        mpfr_set_d(r26433, phi2, MPFR_RNDN);
        mpfr_cos(r26434, r26433, MPFR_RNDN);
        mpfr_mul(r26435, r26432, r26434, MPFR_RNDN);
        mpfr_set_d(r26436, phi1, MPFR_RNDN);
        mpfr_cos(r26437, r26436, MPFR_RNDN);
        mpfr_sin(r26438, r26433, MPFR_RNDN);
        mpfr_mul(r26439, r26437, r26438, MPFR_RNDN);
        mpfr_sin(r26440, r26436, MPFR_RNDN);
        mpfr_mul(r26441, r26440, r26434, MPFR_RNDN);
        mpfr_mul(r26442, r26428, r26426, MPFR_RNDN);
        mpfr_mul(r26443, r26441, r26442, MPFR_RNDN);
        mpfr_sin(r26444, r26425, MPFR_RNDN);
        mpfr_mul(r26445, r26424, r26444, MPFR_RNDN);
        mpfr_mul(r26446, r26441, r26445, MPFR_RNDN);
        mpfr_exp(r26447, r26446, MPFR_RNDN);
        mpfr_log(r26448, r26447, MPFR_RNDN);
        mpfr_add(r26449, r26443, r26448, MPFR_RNDN);
        mpfr_sub(r26450, r26439, r26449, MPFR_RNDN);
        mpfr_atan2(r26451, r26435, r26450, MPFR_RNDN);
        return mpfr_get_d(r26451, MPFR_RNDN);
}