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

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

double f_if(float R, float lambda1, float lambda2, float phi1, float phi2) {
        float r9536 = R;
        float r9537 = 2;
        float r9538 = phi1;
        float r9539 = phi2;
        float r9540 = r9538 - r9539;
        float r9541 = r9540 / r9537;
        float r9542 = sin(r9541);
        float r9543 = pow(r9542, r9537);
        float r9544 = cos(r9538);
        float r9545 = cos(r9539);
        float r9546 = r9544 * r9545;
        float r9547 = lambda1;
        float r9548 = lambda2;
        float r9549 = r9547 - r9548;
        float r9550 = r9549 / r9537;
        float r9551 = sin(r9550);
        float r9552 = r9546 * r9551;
        float r9553 = r9552 * r9551;
        float r9554 = r9543 + r9553;
        float r9555 = sqrt(r9554);
        float r9556 = 1;
        float r9557 = r9556 - r9554;
        float r9558 = sqrt(r9557);
        float r9559 = atan2(r9555, r9558);
        float r9560 = r9537 * r9559;
        float r9561 = r9536 * r9560;
        return r9561;
}

double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) {
        double r9562 = R;
        double r9563 = 2;
        double r9564 = phi1;
        double r9565 = phi2;
        double r9566 = r9564 - r9565;
        double r9567 = r9566 / r9563;
        double r9568 = sin(r9567);
        double r9569 = pow(r9568, r9563);
        double r9570 = cos(r9564);
        double r9571 = cos(r9565);
        double r9572 = r9570 * r9571;
        double r9573 = lambda1;
        double r9574 = lambda2;
        double r9575 = r9573 - r9574;
        double r9576 = r9575 / r9563;
        double r9577 = sin(r9576);
        double r9578 = r9572 * r9577;
        double r9579 = r9578 * r9577;
        double r9580 = r9569 + r9579;
        double r9581 = sqrt(r9580);
        double r9582 = 1;
        double r9583 = r9582 - r9580;
        double r9584 = sqrt(r9583);
        double r9585 = atan2(r9581, r9584);
        double r9586 = r9563 * r9585;
        double r9587 = r9562 * r9586;
        return r9587;
}


double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) {
        float r9588 = R;
        float r9589 = 2;
        float r9590 = phi1;
        float r9591 = phi2;
        float r9592 = r9590 - r9591;
        float r9593 = r9592 / r9589;
        float r9594 = sin(r9593);
        float r9595 = pow(r9594, r9589);
        float r9596 = cos(r9590);
        float r9597 = cos(r9591);
        float r9598 = r9596 * r9597;
        float r9599 = lambda1;
        float r9600 = lambda2;
        float r9601 = r9599 - r9600;
        float r9602 = r9601 / r9589;
        float r9603 = sin(r9602);
        float r9604 = r9598 * r9603;
        float r9605 = expm1(r9603);
        float r9606 = log1p(r9605);
        float r9607 = r9604 * r9606;
        float r9608 = r9595 + r9607;
        float r9609 = sqrt(r9608);
        float r9610 = 1;
        float r9611 = r9610 - r9608;
        float r9612 = sqrt(r9611);
        float r9613 = atan2(r9609, r9612);
        float r9614 = r9589 * r9613;
        float r9615 = r9588 * r9614;
        return r9615;
}

double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) {
        double r9616 = R;
        double r9617 = 2;
        double r9618 = phi1;
        double r9619 = phi2;
        double r9620 = r9618 - r9619;
        double r9621 = r9620 / r9617;
        double r9622 = sin(r9621);
        double r9623 = pow(r9622, r9617);
        double r9624 = cos(r9618);
        double r9625 = cos(r9619);
        double r9626 = r9624 * r9625;
        double r9627 = lambda1;
        double r9628 = lambda2;
        double r9629 = r9627 - r9628;
        double r9630 = r9629 / r9617;
        double r9631 = sin(r9630);
        double r9632 = r9626 * r9631;
        double r9633 = expm1(r9631);
        double r9634 = log1p(r9633);
        double r9635 = r9632 * r9634;
        double r9636 = r9623 + r9635;
        double r9637 = sqrt(r9636);
        double r9638 = 1;
        double r9639 = r9638 - r9636;
        double r9640 = sqrt(r9639);
        double r9641 = atan2(r9637, r9640);
        double r9642 = r9617 * r9641;
        double r9643 = r9616 * r9642;
        return r9643;
}

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 r9644, r9645, r9646, r9647, r9648, r9649, r9650, r9651, r9652, r9653, r9654, r9655, r9656, r9657, r9658, r9659, r9660, r9661, r9662, r9663, r9664, r9665, r9666, r9667, r9668, r9669;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9644);
        mpfr_init_set_str(r9645, "2", 10, MPFR_RNDN);
        mpfr_init(r9646);
        mpfr_init(r9647);
        mpfr_init(r9648);
        mpfr_init(r9649);
        mpfr_init(r9650);
        mpfr_init(r9651);
        mpfr_init(r9652);
        mpfr_init(r9653);
        mpfr_init(r9654);
        mpfr_init(r9655);
        mpfr_init(r9656);
        mpfr_init(r9657);
        mpfr_init(r9658);
        mpfr_init(r9659);
        mpfr_init(r9660);
        mpfr_init(r9661);
        mpfr_init(r9662);
        mpfr_init(r9663);
        mpfr_init_set_str(r9664, "1", 10, MPFR_RNDN);
        mpfr_init(r9665);
        mpfr_init(r9666);
        mpfr_init(r9667);
        mpfr_init(r9668);
        mpfr_init(r9669);
}

double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r9644, R, MPFR_RNDN);
        ;
        mpfr_set_d(r9646, phi1, MPFR_RNDN);
        mpfr_set_d(r9647, phi2, MPFR_RNDN);
        mpfr_sub(r9648, r9646, r9647, MPFR_RNDN);
        mpfr_div(r9649, r9648, r9645, MPFR_RNDN);
        mpfr_sin(r9650, r9649, MPFR_RNDN);
        mpfr_pow(r9651, r9650, r9645, MPFR_RNDN);
        mpfr_cos(r9652, r9646, MPFR_RNDN);
        mpfr_cos(r9653, r9647, MPFR_RNDN);
        mpfr_mul(r9654, r9652, r9653, MPFR_RNDN);
        mpfr_set_d(r9655, lambda1, MPFR_RNDN);
        mpfr_set_d(r9656, lambda2, MPFR_RNDN);
        mpfr_sub(r9657, r9655, r9656, MPFR_RNDN);
        mpfr_div(r9658, r9657, r9645, MPFR_RNDN);
        mpfr_sin(r9659, r9658, MPFR_RNDN);
        mpfr_mul(r9660, r9654, r9659, MPFR_RNDN);
        mpfr_mul(r9661, r9660, r9659, MPFR_RNDN);
        mpfr_add(r9662, r9651, r9661, MPFR_RNDN);
        mpfr_sqrt(r9663, r9662, MPFR_RNDN);
        ;
        mpfr_sub(r9665, r9664, r9662, MPFR_RNDN);
        mpfr_sqrt(r9666, r9665, MPFR_RNDN);
        mpfr_atan2(r9667, r9663, r9666, MPFR_RNDN);
        mpfr_mul(r9668, r9645, r9667, MPFR_RNDN);
        mpfr_mul(r9669, r9644, r9668, MPFR_RNDN);
        return mpfr_get_d(r9669, MPFR_RNDN);
}

static mpfr_t r9670, r9671, r9672, r9673, r9674, r9675, r9676, r9677, r9678, r9679, r9680, r9681, r9682, r9683, r9684, r9685, r9686, r9687, r9688, r9689, r9690, r9691, r9692, r9693, r9694, r9695, r9696, r9697;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9670);
        mpfr_init_set_str(r9671, "2", 10, MPFR_RNDN);
        mpfr_init(r9672);
        mpfr_init(r9673);
        mpfr_init(r9674);
        mpfr_init(r9675);
        mpfr_init(r9676);
        mpfr_init(r9677);
        mpfr_init(r9678);
        mpfr_init(r9679);
        mpfr_init(r9680);
        mpfr_init(r9681);
        mpfr_init(r9682);
        mpfr_init(r9683);
        mpfr_init(r9684);
        mpfr_init(r9685);
        mpfr_init(r9686);
        mpfr_init(r9687);
        mpfr_init(r9688);
        mpfr_init(r9689);
        mpfr_init(r9690);
        mpfr_init(r9691);
        mpfr_init_set_str(r9692, "1", 10, MPFR_RNDN);
        mpfr_init(r9693);
        mpfr_init(r9694);
        mpfr_init(r9695);
        mpfr_init(r9696);
        mpfr_init(r9697);
}

double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r9670, R, MPFR_RNDN);
        ;
        mpfr_set_d(r9672, phi1, MPFR_RNDN);
        mpfr_set_d(r9673, phi2, MPFR_RNDN);
        mpfr_sub(r9674, r9672, r9673, MPFR_RNDN);
        mpfr_div(r9675, r9674, r9671, MPFR_RNDN);
        mpfr_sin(r9676, r9675, MPFR_RNDN);
        mpfr_pow(r9677, r9676, r9671, MPFR_RNDN);
        mpfr_cos(r9678, r9672, MPFR_RNDN);
        mpfr_cos(r9679, r9673, MPFR_RNDN);
        mpfr_mul(r9680, r9678, r9679, MPFR_RNDN);
        mpfr_set_d(r9681, lambda1, MPFR_RNDN);
        mpfr_set_d(r9682, lambda2, MPFR_RNDN);
        mpfr_sub(r9683, r9681, r9682, MPFR_RNDN);
        mpfr_div(r9684, r9683, r9671, MPFR_RNDN);
        mpfr_sin(r9685, r9684, MPFR_RNDN);
        mpfr_mul(r9686, r9680, r9685, MPFR_RNDN);
        mpfr_expm1(r9687, r9685, MPFR_RNDN);
        mpfr_log1p(r9688, r9687, MPFR_RNDN);
        mpfr_mul(r9689, r9686, r9688, MPFR_RNDN);
        mpfr_add(r9690, r9677, r9689, MPFR_RNDN);
        mpfr_sqrt(r9691, r9690, MPFR_RNDN);
        ;
        mpfr_sub(r9693, r9692, r9690, MPFR_RNDN);
        mpfr_sqrt(r9694, r9693, MPFR_RNDN);
        mpfr_atan2(r9695, r9691, r9694, MPFR_RNDN);
        mpfr_mul(r9696, r9671, r9695, MPFR_RNDN);
        mpfr_mul(r9697, r9670, r9696, MPFR_RNDN);
        return mpfr_get_d(r9697, MPFR_RNDN);
}

static mpfr_t r9698, r9699, r9700, r9701, r9702, r9703, r9704, r9705, r9706, r9707, r9708, r9709, r9710, r9711, r9712, r9713, r9714, r9715, r9716, r9717, r9718, r9719, r9720, r9721, r9722, r9723, r9724, r9725;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9698);
        mpfr_init_set_str(r9699, "2", 10, MPFR_RNDN);
        mpfr_init(r9700);
        mpfr_init(r9701);
        mpfr_init(r9702);
        mpfr_init(r9703);
        mpfr_init(r9704);
        mpfr_init(r9705);
        mpfr_init(r9706);
        mpfr_init(r9707);
        mpfr_init(r9708);
        mpfr_init(r9709);
        mpfr_init(r9710);
        mpfr_init(r9711);
        mpfr_init(r9712);
        mpfr_init(r9713);
        mpfr_init(r9714);
        mpfr_init(r9715);
        mpfr_init(r9716);
        mpfr_init(r9717);
        mpfr_init(r9718);
        mpfr_init(r9719);
        mpfr_init_set_str(r9720, "1", 10, MPFR_RNDN);
        mpfr_init(r9721);
        mpfr_init(r9722);
        mpfr_init(r9723);
        mpfr_init(r9724);
        mpfr_init(r9725);
}

double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r9698, R, MPFR_RNDN);
        ;
        mpfr_set_d(r9700, phi1, MPFR_RNDN);
        mpfr_set_d(r9701, phi2, MPFR_RNDN);
        mpfr_sub(r9702, r9700, r9701, MPFR_RNDN);
        mpfr_div(r9703, r9702, r9699, MPFR_RNDN);
        mpfr_sin(r9704, r9703, MPFR_RNDN);
        mpfr_pow(r9705, r9704, r9699, MPFR_RNDN);
        mpfr_cos(r9706, r9700, MPFR_RNDN);
        mpfr_cos(r9707, r9701, MPFR_RNDN);
        mpfr_mul(r9708, r9706, r9707, MPFR_RNDN);
        mpfr_set_d(r9709, lambda1, MPFR_RNDN);
        mpfr_set_d(r9710, lambda2, MPFR_RNDN);
        mpfr_sub(r9711, r9709, r9710, MPFR_RNDN);
        mpfr_div(r9712, r9711, r9699, MPFR_RNDN);
        mpfr_sin(r9713, r9712, MPFR_RNDN);
        mpfr_mul(r9714, r9708, r9713, MPFR_RNDN);
        mpfr_expm1(r9715, r9713, MPFR_RNDN);
        mpfr_log1p(r9716, r9715, MPFR_RNDN);
        mpfr_mul(r9717, r9714, r9716, MPFR_RNDN);
        mpfr_add(r9718, r9705, r9717, MPFR_RNDN);
        mpfr_sqrt(r9719, r9718, MPFR_RNDN);
        ;
        mpfr_sub(r9721, r9720, r9718, MPFR_RNDN);
        mpfr_sqrt(r9722, r9721, MPFR_RNDN);
        mpfr_atan2(r9723, r9719, r9722, MPFR_RNDN);
        mpfr_mul(r9724, r9699, r9723, MPFR_RNDN);
        mpfr_mul(r9725, r9698, r9724, MPFR_RNDN);
        return mpfr_get_d(r9725, MPFR_RNDN);
}