#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 r10375 = lambda1;
        float r10376 = lambda2;
        float r10377 = r10375 - r10376;
        float r10378 = sin(r10377);
        float r10379 = phi2;
        float r10380 = cos(r10379);
        float r10381 = r10378 * r10380;
        float r10382 = phi1;
        float r10383 = cos(r10382);
        float r10384 = sin(r10379);
        float r10385 = r10383 * r10384;
        float r10386 = sin(r10382);
        float r10387 = r10386 * r10380;
        float r10388 = cos(r10377);
        float r10389 = r10387 * r10388;
        float r10390 = r10385 - r10389;
        float r10391 = atan2(r10381, r10390);
        return r10391;
}

double f_id(double lambda1, double lambda2, double phi1, double phi2) {
        double r10392 = lambda1;
        double r10393 = lambda2;
        double r10394 = r10392 - r10393;
        double r10395 = sin(r10394);
        double r10396 = phi2;
        double r10397 = cos(r10396);
        double r10398 = r10395 * r10397;
        double r10399 = phi1;
        double r10400 = cos(r10399);
        double r10401 = sin(r10396);
        double r10402 = r10400 * r10401;
        double r10403 = sin(r10399);
        double r10404 = r10403 * r10397;
        double r10405 = cos(r10394);
        double r10406 = r10404 * r10405;
        double r10407 = r10402 - r10406;
        double r10408 = atan2(r10398, r10407);
        return r10408;
}


double f_of(float lambda1, float lambda2, float phi1, float phi2) {
        float r10409 = phi2;
        float r10410 = cos(r10409);
        float r10411 = lambda2;
        float r10412 = cos(r10411);
        float r10413 = lambda1;
        float r10414 = sin(r10413);
        float r10415 = r10412 * r10414;
        float r10416 = cos(r10413);
        float r10417 = sin(r10411);
        float r10418 = r10416 * r10417;
        float r10419 = r10415 - r10418;
        float r10420 = r10410 * r10419;
        float r10421 = phi1;
        float r10422 = sin(r10421);
        float r10423 = -r10422;
        float r10424 = r10423 * r10410;
        float r10425 = r10412 * r10416;
        float r10426 = fma(r10417, r10414, r10425);
        float r10427 = sin(r10409);
        float r10428 = cos(r10421);
        float r10429 = r10427 * r10428;
        float r10430 = fma(r10424, r10426, r10429);
        float r10431 = log1p(r10430);
        float r10432 = expm1(r10431);
        float r10433 = atan2(r10420, r10432);
        return r10433;
}

double f_od(double lambda1, double lambda2, double phi1, double phi2) {
        double r10434 = phi2;
        double r10435 = cos(r10434);
        double r10436 = lambda2;
        double r10437 = cos(r10436);
        double r10438 = lambda1;
        double r10439 = sin(r10438);
        double r10440 = r10437 * r10439;
        double r10441 = cos(r10438);
        double r10442 = sin(r10436);
        double r10443 = r10441 * r10442;
        double r10444 = r10440 - r10443;
        double r10445 = r10435 * r10444;
        double r10446 = phi1;
        double r10447 = sin(r10446);
        double r10448 = -r10447;
        double r10449 = r10448 * r10435;
        double r10450 = r10437 * r10441;
        double r10451 = fma(r10442, r10439, r10450);
        double r10452 = sin(r10434);
        double r10453 = cos(r10446);
        double r10454 = r10452 * r10453;
        double r10455 = fma(r10449, r10451, r10454);
        double r10456 = log1p(r10455);
        double r10457 = expm1(r10456);
        double r10458 = atan2(r10445, r10457);
        return r10458;
}

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 r10459, r10460, r10461, r10462, r10463, r10464, r10465, r10466, r10467, r10468, r10469, r10470, r10471, r10472, r10473, r10474, r10475;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r10459);
        mpfr_init(r10460);
        mpfr_init(r10461);
        mpfr_init(r10462);
        mpfr_init(r10463);
        mpfr_init(r10464);
        mpfr_init(r10465);
        mpfr_init(r10466);
        mpfr_init(r10467);
        mpfr_init(r10468);
        mpfr_init(r10469);
        mpfr_init(r10470);
        mpfr_init(r10471);
        mpfr_init(r10472);
        mpfr_init(r10473);
        mpfr_init(r10474);
        mpfr_init(r10475);
}

double f_im(double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r10459, lambda1, MPFR_RNDN);
        mpfr_set_d(r10460, lambda2, MPFR_RNDN);
        mpfr_sub(r10461, r10459, r10460, MPFR_RNDN);
        mpfr_sin(r10462, r10461, MPFR_RNDN);
        mpfr_set_d(r10463, phi2, MPFR_RNDN);
        mpfr_cos(r10464, r10463, MPFR_RNDN);
        mpfr_mul(r10465, r10462, r10464, MPFR_RNDN);
        mpfr_set_d(r10466, phi1, MPFR_RNDN);
        mpfr_cos(r10467, r10466, MPFR_RNDN);
        mpfr_sin(r10468, r10463, MPFR_RNDN);
        mpfr_mul(r10469, r10467, r10468, MPFR_RNDN);
        mpfr_sin(r10470, r10466, MPFR_RNDN);
        mpfr_mul(r10471, r10470, r10464, MPFR_RNDN);
        mpfr_cos(r10472, r10461, MPFR_RNDN);
        mpfr_mul(r10473, r10471, r10472, MPFR_RNDN);
        mpfr_sub(r10474, r10469, r10473, MPFR_RNDN);
        mpfr_atan2(r10475, r10465, r10474, MPFR_RNDN);
        return mpfr_get_d(r10475, MPFR_RNDN);
}

static mpfr_t r10476, r10477, r10478, r10479, r10480, r10481, r10482, r10483, r10484, r10485, r10486, r10487, r10488, r10489, r10490, r10491, r10492, r10493, r10494, r10495, r10496, r10497, r10498, r10499, r10500;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r10476);
        mpfr_init(r10477);
        mpfr_init(r10478);
        mpfr_init(r10479);
        mpfr_init(r10480);
        mpfr_init(r10481);
        mpfr_init(r10482);
        mpfr_init(r10483);
        mpfr_init(r10484);
        mpfr_init(r10485);
        mpfr_init(r10486);
        mpfr_init(r10487);
        mpfr_init(r10488);
        mpfr_init(r10489);
        mpfr_init(r10490);
        mpfr_init(r10491);
        mpfr_init(r10492);
        mpfr_init(r10493);
        mpfr_init(r10494);
        mpfr_init(r10495);
        mpfr_init(r10496);
        mpfr_init(r10497);
        mpfr_init(r10498);
        mpfr_init(r10499);
        mpfr_init(r10500);
}

double f_fm(double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r10476, phi2, MPFR_RNDN);
        mpfr_cos(r10477, r10476, MPFR_RNDN);
        mpfr_set_d(r10478, lambda2, MPFR_RNDN);
        mpfr_cos(r10479, r10478, MPFR_RNDN);
        mpfr_set_d(r10480, lambda1, MPFR_RNDN);
        mpfr_sin(r10481, r10480, MPFR_RNDN);
        mpfr_mul(r10482, r10479, r10481, MPFR_RNDN);
        mpfr_cos(r10483, r10480, MPFR_RNDN);
        mpfr_sin(r10484, r10478, MPFR_RNDN);
        mpfr_mul(r10485, r10483, r10484, MPFR_RNDN);
        mpfr_sub(r10486, r10482, r10485, MPFR_RNDN);
        mpfr_mul(r10487, r10477, r10486, MPFR_RNDN);
        mpfr_set_d(r10488, phi1, MPFR_RNDN);
        mpfr_sin(r10489, r10488, MPFR_RNDN);
        mpfr_neg(r10490, r10489, MPFR_RNDN);
        mpfr_mul(r10491, r10490, r10477, MPFR_RNDN);
        mpfr_mul(r10492, r10479, r10483, MPFR_RNDN);
        mpfr_fma(r10493, r10484, r10481, r10492, MPFR_RNDN);
        mpfr_sin(r10494, r10476, MPFR_RNDN);
        mpfr_cos(r10495, r10488, MPFR_RNDN);
        mpfr_mul(r10496, r10494, r10495, MPFR_RNDN);
        mpfr_fma(r10497, r10491, r10493, r10496, MPFR_RNDN);
        mpfr_log1p(r10498, r10497, MPFR_RNDN);
        mpfr_expm1(r10499, r10498, MPFR_RNDN);
        mpfr_atan2(r10500, r10487, r10499, MPFR_RNDN);
        return mpfr_get_d(r10500, MPFR_RNDN);
}

static mpfr_t r10501, r10502, r10503, r10504, r10505, r10506, r10507, r10508, r10509, r10510, r10511, r10512, r10513, r10514, r10515, r10516, r10517, r10518, r10519, r10520, r10521, r10522, r10523, r10524, r10525;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r10501);
        mpfr_init(r10502);
        mpfr_init(r10503);
        mpfr_init(r10504);
        mpfr_init(r10505);
        mpfr_init(r10506);
        mpfr_init(r10507);
        mpfr_init(r10508);
        mpfr_init(r10509);
        mpfr_init(r10510);
        mpfr_init(r10511);
        mpfr_init(r10512);
        mpfr_init(r10513);
        mpfr_init(r10514);
        mpfr_init(r10515);
        mpfr_init(r10516);
        mpfr_init(r10517);
        mpfr_init(r10518);
        mpfr_init(r10519);
        mpfr_init(r10520);
        mpfr_init(r10521);
        mpfr_init(r10522);
        mpfr_init(r10523);
        mpfr_init(r10524);
        mpfr_init(r10525);
}

double f_dm(double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r10501, phi2, MPFR_RNDN);
        mpfr_cos(r10502, r10501, MPFR_RNDN);
        mpfr_set_d(r10503, lambda2, MPFR_RNDN);
        mpfr_cos(r10504, r10503, MPFR_RNDN);
        mpfr_set_d(r10505, lambda1, MPFR_RNDN);
        mpfr_sin(r10506, r10505, MPFR_RNDN);
        mpfr_mul(r10507, r10504, r10506, MPFR_RNDN);
        mpfr_cos(r10508, r10505, MPFR_RNDN);
        mpfr_sin(r10509, r10503, MPFR_RNDN);
        mpfr_mul(r10510, r10508, r10509, MPFR_RNDN);
        mpfr_sub(r10511, r10507, r10510, MPFR_RNDN);
        mpfr_mul(r10512, r10502, r10511, MPFR_RNDN);
        mpfr_set_d(r10513, phi1, MPFR_RNDN);
        mpfr_sin(r10514, r10513, MPFR_RNDN);
        mpfr_neg(r10515, r10514, MPFR_RNDN);
        mpfr_mul(r10516, r10515, r10502, MPFR_RNDN);
        mpfr_mul(r10517, r10504, r10508, MPFR_RNDN);
        mpfr_fma(r10518, r10509, r10506, r10517, MPFR_RNDN);
        mpfr_sin(r10519, r10501, MPFR_RNDN);
        mpfr_cos(r10520, r10513, MPFR_RNDN);
        mpfr_mul(r10521, r10519, r10520, MPFR_RNDN);
        mpfr_fma(r10522, r10516, r10518, r10521, MPFR_RNDN);
        mpfr_log1p(r10523, r10522, MPFR_RNDN);
        mpfr_expm1(r10524, r10523, MPFR_RNDN);
        mpfr_atan2(r10525, r10512, r10524, MPFR_RNDN);
        return mpfr_get_d(r10525, MPFR_RNDN);
}