#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 r26508 = lambda1;
        float r26509 = lambda2;
        float r26510 = r26508 - r26509;
        float r26511 = sin(r26510);
        float r26512 = phi2;
        float r26513 = cos(r26512);
        float r26514 = r26511 * r26513;
        float r26515 = phi1;
        float r26516 = cos(r26515);
        float r26517 = sin(r26512);
        float r26518 = r26516 * r26517;
        float r26519 = sin(r26515);
        float r26520 = r26519 * r26513;
        float r26521 = cos(r26510);
        float r26522 = r26520 * r26521;
        float r26523 = r26518 - r26522;
        float r26524 = atan2(r26514, r26523);
        return r26524;
}

double f_id(double lambda1, double lambda2, double phi1, double phi2) {
        double r26525 = lambda1;
        double r26526 = lambda2;
        double r26527 = r26525 - r26526;
        double r26528 = sin(r26527);
        double r26529 = phi2;
        double r26530 = cos(r26529);
        double r26531 = r26528 * r26530;
        double r26532 = phi1;
        double r26533 = cos(r26532);
        double r26534 = sin(r26529);
        double r26535 = r26533 * r26534;
        double r26536 = sin(r26532);
        double r26537 = r26536 * r26530;
        double r26538 = cos(r26527);
        double r26539 = r26537 * r26538;
        double r26540 = r26535 - r26539;
        double r26541 = atan2(r26531, r26540);
        return r26541;
}


double f_of(float lambda1, float lambda2, float phi1, float phi2) {
        float r26542 = lambda1;
        float r26543 = sin(r26542);
        float r26544 = lambda2;
        float r26545 = cos(r26544);
        float r26546 = r26543 * r26545;
        float r26547 = cos(r26542);
        float r26548 = -r26544;
        float r26549 = sin(r26548);
        float r26550 = r26547 * r26549;
        float r26551 = log1p(r26550);
        float r26552 = expm1(r26551);
        float r26553 = r26546 + r26552;
        float r26554 = phi2;
        float r26555 = cos(r26554);
        float r26556 = r26553 * r26555;
        float r26557 = phi1;
        float r26558 = cos(r26557);
        float r26559 = sin(r26554);
        float r26560 = r26558 * r26559;
        float r26561 = sin(r26557);
        float r26562 = r26561 * r26555;
        float r26563 = r26547 * r26545;
        float r26564 = r26562 * r26563;
        float r26565 = sin(r26544);
        float r26566 = r26543 * r26565;
        float r26567 = r26562 * r26566;
        float r26568 = r26564 + r26567;
        float r26569 = r26560 - r26568;
        float r26570 = atan2(r26556, r26569);
        return r26570;
}

double f_od(double lambda1, double lambda2, double phi1, double phi2) {
        double r26571 = lambda1;
        double r26572 = sin(r26571);
        double r26573 = lambda2;
        double r26574 = cos(r26573);
        double r26575 = r26572 * r26574;
        double r26576 = cos(r26571);
        double r26577 = -r26573;
        double r26578 = sin(r26577);
        double r26579 = r26576 * r26578;
        double r26580 = log1p(r26579);
        double r26581 = expm1(r26580);
        double r26582 = r26575 + r26581;
        double r26583 = phi2;
        double r26584 = cos(r26583);
        double r26585 = r26582 * r26584;
        double r26586 = phi1;
        double r26587 = cos(r26586);
        double r26588 = sin(r26583);
        double r26589 = r26587 * r26588;
        double r26590 = sin(r26586);
        double r26591 = r26590 * r26584;
        double r26592 = r26576 * r26574;
        double r26593 = r26591 * r26592;
        double r26594 = sin(r26573);
        double r26595 = r26572 * r26594;
        double r26596 = r26591 * r26595;
        double r26597 = r26593 + r26596;
        double r26598 = r26589 - r26597;
        double r26599 = atan2(r26585, r26598);
        return r26599;
}

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 r26600, r26601, r26602, r26603, r26604, r26605, r26606, r26607, r26608, r26609, r26610, r26611, r26612, r26613, r26614, r26615, r26616;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26600);
        mpfr_init(r26601);
        mpfr_init(r26602);
        mpfr_init(r26603);
        mpfr_init(r26604);
        mpfr_init(r26605);
        mpfr_init(r26606);
        mpfr_init(r26607);
        mpfr_init(r26608);
        mpfr_init(r26609);
        mpfr_init(r26610);
        mpfr_init(r26611);
        mpfr_init(r26612);
        mpfr_init(r26613);
        mpfr_init(r26614);
        mpfr_init(r26615);
        mpfr_init(r26616);
}

double f_im(double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r26600, lambda1, MPFR_RNDN);
        mpfr_set_d(r26601, lambda2, MPFR_RNDN);
        mpfr_sub(r26602, r26600, r26601, MPFR_RNDN);
        mpfr_sin(r26603, r26602, MPFR_RNDN);
        mpfr_set_d(r26604, phi2, MPFR_RNDN);
        mpfr_cos(r26605, r26604, MPFR_RNDN);
        mpfr_mul(r26606, r26603, r26605, MPFR_RNDN);
        mpfr_set_d(r26607, phi1, MPFR_RNDN);
        mpfr_cos(r26608, r26607, MPFR_RNDN);
        mpfr_sin(r26609, r26604, MPFR_RNDN);
        mpfr_mul(r26610, r26608, r26609, MPFR_RNDN);
        mpfr_sin(r26611, r26607, MPFR_RNDN);
        mpfr_mul(r26612, r26611, r26605, MPFR_RNDN);
        mpfr_cos(r26613, r26602, MPFR_RNDN);
        mpfr_mul(r26614, r26612, r26613, MPFR_RNDN);
        mpfr_sub(r26615, r26610, r26614, MPFR_RNDN);
        mpfr_atan2(r26616, r26606, r26615, MPFR_RNDN);
        return mpfr_get_d(r26616, MPFR_RNDN);
}

static mpfr_t r26617, r26618, r26619, r26620, r26621, r26622, r26623, r26624, r26625, r26626, r26627, r26628, r26629, r26630, r26631, r26632, r26633, r26634, r26635, r26636, r26637, r26638, r26639, r26640, r26641, r26642, r26643, r26644, r26645;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26617);
        mpfr_init(r26618);
        mpfr_init(r26619);
        mpfr_init(r26620);
        mpfr_init(r26621);
        mpfr_init(r26622);
        mpfr_init(r26623);
        mpfr_init(r26624);
        mpfr_init(r26625);
        mpfr_init(r26626);
        mpfr_init(r26627);
        mpfr_init(r26628);
        mpfr_init(r26629);
        mpfr_init(r26630);
        mpfr_init(r26631);
        mpfr_init(r26632);
        mpfr_init(r26633);
        mpfr_init(r26634);
        mpfr_init(r26635);
        mpfr_init(r26636);
        mpfr_init(r26637);
        mpfr_init(r26638);
        mpfr_init(r26639);
        mpfr_init(r26640);
        mpfr_init(r26641);
        mpfr_init(r26642);
        mpfr_init(r26643);
        mpfr_init(r26644);
        mpfr_init(r26645);
}

double f_fm(double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r26617, lambda1, MPFR_RNDN);
        mpfr_sin(r26618, r26617, MPFR_RNDN);
        mpfr_set_d(r26619, lambda2, MPFR_RNDN);
        mpfr_cos(r26620, r26619, MPFR_RNDN);
        mpfr_mul(r26621, r26618, r26620, MPFR_RNDN);
        mpfr_cos(r26622, r26617, MPFR_RNDN);
        mpfr_neg(r26623, r26619, MPFR_RNDN);
        mpfr_sin(r26624, r26623, MPFR_RNDN);
        mpfr_mul(r26625, r26622, r26624, MPFR_RNDN);
        mpfr_log1p(r26626, r26625, MPFR_RNDN);
        mpfr_expm1(r26627, r26626, MPFR_RNDN);
        mpfr_add(r26628, r26621, r26627, MPFR_RNDN);
        mpfr_set_d(r26629, phi2, MPFR_RNDN);
        mpfr_cos(r26630, r26629, MPFR_RNDN);
        mpfr_mul(r26631, r26628, r26630, MPFR_RNDN);
        mpfr_set_d(r26632, phi1, MPFR_RNDN);
        mpfr_cos(r26633, r26632, MPFR_RNDN);
        mpfr_sin(r26634, r26629, MPFR_RNDN);
        mpfr_mul(r26635, r26633, r26634, MPFR_RNDN);
        mpfr_sin(r26636, r26632, MPFR_RNDN);
        mpfr_mul(r26637, r26636, r26630, MPFR_RNDN);
        mpfr_mul(r26638, r26622, r26620, MPFR_RNDN);
        mpfr_mul(r26639, r26637, r26638, MPFR_RNDN);
        mpfr_sin(r26640, r26619, MPFR_RNDN);
        mpfr_mul(r26641, r26618, r26640, MPFR_RNDN);
        mpfr_mul(r26642, r26637, r26641, MPFR_RNDN);
        mpfr_add(r26643, r26639, r26642, MPFR_RNDN);
        mpfr_sub(r26644, r26635, r26643, MPFR_RNDN);
        mpfr_atan2(r26645, r26631, r26644, MPFR_RNDN);
        return mpfr_get_d(r26645, MPFR_RNDN);
}

static mpfr_t r26646, r26647, r26648, r26649, r26650, r26651, r26652, r26653, r26654, r26655, r26656, r26657, r26658, r26659, r26660, r26661, r26662, r26663, r26664, r26665, r26666, r26667, r26668, r26669, r26670, r26671, r26672, r26673, r26674;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26646);
        mpfr_init(r26647);
        mpfr_init(r26648);
        mpfr_init(r26649);
        mpfr_init(r26650);
        mpfr_init(r26651);
        mpfr_init(r26652);
        mpfr_init(r26653);
        mpfr_init(r26654);
        mpfr_init(r26655);
        mpfr_init(r26656);
        mpfr_init(r26657);
        mpfr_init(r26658);
        mpfr_init(r26659);
        mpfr_init(r26660);
        mpfr_init(r26661);
        mpfr_init(r26662);
        mpfr_init(r26663);
        mpfr_init(r26664);
        mpfr_init(r26665);
        mpfr_init(r26666);
        mpfr_init(r26667);
        mpfr_init(r26668);
        mpfr_init(r26669);
        mpfr_init(r26670);
        mpfr_init(r26671);
        mpfr_init(r26672);
        mpfr_init(r26673);
        mpfr_init(r26674);
}

double f_dm(double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r26646, lambda1, MPFR_RNDN);
        mpfr_sin(r26647, r26646, MPFR_RNDN);
        mpfr_set_d(r26648, lambda2, MPFR_RNDN);
        mpfr_cos(r26649, r26648, MPFR_RNDN);
        mpfr_mul(r26650, r26647, r26649, MPFR_RNDN);
        mpfr_cos(r26651, r26646, MPFR_RNDN);
        mpfr_neg(r26652, r26648, MPFR_RNDN);
        mpfr_sin(r26653, r26652, MPFR_RNDN);
        mpfr_mul(r26654, r26651, r26653, MPFR_RNDN);
        mpfr_log1p(r26655, r26654, MPFR_RNDN);
        mpfr_expm1(r26656, r26655, MPFR_RNDN);
        mpfr_add(r26657, r26650, r26656, MPFR_RNDN);
        mpfr_set_d(r26658, phi2, MPFR_RNDN);
        mpfr_cos(r26659, r26658, MPFR_RNDN);
        mpfr_mul(r26660, r26657, r26659, MPFR_RNDN);
        mpfr_set_d(r26661, phi1, MPFR_RNDN);
        mpfr_cos(r26662, r26661, MPFR_RNDN);
        mpfr_sin(r26663, r26658, MPFR_RNDN);
        mpfr_mul(r26664, r26662, r26663, MPFR_RNDN);
        mpfr_sin(r26665, r26661, MPFR_RNDN);
        mpfr_mul(r26666, r26665, r26659, MPFR_RNDN);
        mpfr_mul(r26667, r26651, r26649, MPFR_RNDN);
        mpfr_mul(r26668, r26666, r26667, MPFR_RNDN);
        mpfr_sin(r26669, r26648, MPFR_RNDN);
        mpfr_mul(r26670, r26647, r26669, MPFR_RNDN);
        mpfr_mul(r26671, r26666, r26670, MPFR_RNDN);
        mpfr_add(r26672, r26668, r26671, MPFR_RNDN);
        mpfr_sub(r26673, r26664, r26672, MPFR_RNDN);
        mpfr_atan2(r26674, r26660, r26673, MPFR_RNDN);
        return mpfr_get_d(r26674, MPFR_RNDN);
}