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

char *name = "Destination given bearing on a great circle";

double f_if(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) {
        float r20539 = lambda1;
        float r20540 = theta;
        float r20541 = sin(r20540);
        float r20542 = delta;
        float r20543 = sin(r20542);
        float r20544 = r20541 * r20543;
        float r20545 = phi1;
        float r20546 = cos(r20545);
        float r20547 = r20544 * r20546;
        float r20548 = cos(r20542);
        float r20549 = sin(r20545);
        float r20550 = r20549 * r20548;
        float r20551 = r20546 * r20543;
        float r20552 = cos(r20540);
        float r20553 = r20551 * r20552;
        float r20554 = r20550 + r20553;
        float r20555 = asin(r20554);
        float r20556 = sin(r20555);
        float r20557 = r20549 * r20556;
        float r20558 = r20548 - r20557;
        float r20559 = atan2(r20547, r20558);
        float r20560 = r20539 + r20559;
        return r20560;
}

double f_id(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r20561 = lambda1;
        double r20562 = theta;
        double r20563 = sin(r20562);
        double r20564 = delta;
        double r20565 = sin(r20564);
        double r20566 = r20563 * r20565;
        double r20567 = phi1;
        double r20568 = cos(r20567);
        double r20569 = r20566 * r20568;
        double r20570 = cos(r20564);
        double r20571 = sin(r20567);
        double r20572 = r20571 * r20570;
        double r20573 = r20568 * r20565;
        double r20574 = cos(r20562);
        double r20575 = r20573 * r20574;
        double r20576 = r20572 + r20575;
        double r20577 = asin(r20576);
        double r20578 = sin(r20577);
        double r20579 = r20571 * r20578;
        double r20580 = r20570 - r20579;
        double r20581 = atan2(r20569, r20580);
        double r20582 = r20561 + r20581;
        return r20582;
}


double f_of(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) {
        float r20583 = lambda1;
        float r20584 = theta;
        float r20585 = sin(r20584);
        float r20586 = delta;
        float r20587 = sin(r20586);
        float r20588 = r20585 * r20587;
        float r20589 = phi1;
        float r20590 = cos(r20589);
        float r20591 = r20588 * r20590;
        float r20592 = cos(r20586);
        float r20593 = sin(r20589);
        float r20594 = r20592 * r20593;
        float r20595 = cos(r20584);
        float r20596 = r20587 * r20590;
        float r20597 = r20595 * r20596;
        float r20598 = r20594 + r20597;
        float r20599 = asin(r20598);
        float r20600 = sin(r20599);
        float r20601 = 3;
        float r20602 = pow(r20593, r20601);
        float r20603 = r20600 * r20602;
        float r20604 = r20600 * r20600;
        float r20605 = r20603 * r20604;
        float r20606 = cbrt(r20605);
        float r20607 = r20592 - r20606;
        float r20608 = atan2(r20591, r20607);
        float r20609 = r20583 + r20608;
        return r20609;
}

double f_od(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r20610 = lambda1;
        double r20611 = theta;
        double r20612 = sin(r20611);
        double r20613 = delta;
        double r20614 = sin(r20613);
        double r20615 = r20612 * r20614;
        double r20616 = phi1;
        double r20617 = cos(r20616);
        double r20618 = r20615 * r20617;
        double r20619 = cos(r20613);
        double r20620 = sin(r20616);
        double r20621 = r20619 * r20620;
        double r20622 = cos(r20611);
        double r20623 = r20614 * r20617;
        double r20624 = r20622 * r20623;
        double r20625 = r20621 + r20624;
        double r20626 = asin(r20625);
        double r20627 = sin(r20626);
        double r20628 = 3;
        double r20629 = pow(r20620, r20628);
        double r20630 = r20627 * r20629;
        double r20631 = r20627 * r20627;
        double r20632 = r20630 * r20631;
        double r20633 = cbrt(r20632);
        double r20634 = r20619 - r20633;
        double r20635 = atan2(r20618, r20634);
        double r20636 = r20610 + r20635;
        return r20636;
}

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 r20637, r20638, r20639, r20640, r20641, r20642, r20643, r20644, r20645, r20646, r20647, r20648, r20649, r20650, r20651, r20652, r20653, r20654, r20655, r20656, r20657, r20658;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r20637);
        mpfr_init(r20638);
        mpfr_init(r20639);
        mpfr_init(r20640);
        mpfr_init(r20641);
        mpfr_init(r20642);
        mpfr_init(r20643);
        mpfr_init(r20644);
        mpfr_init(r20645);
        mpfr_init(r20646);
        mpfr_init(r20647);
        mpfr_init(r20648);
        mpfr_init(r20649);
        mpfr_init(r20650);
        mpfr_init(r20651);
        mpfr_init(r20652);
        mpfr_init(r20653);
        mpfr_init(r20654);
        mpfr_init(r20655);
        mpfr_init(r20656);
        mpfr_init(r20657);
        mpfr_init(r20658);
}

double f_im(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r20637, lambda1, MPFR_RNDN);
        mpfr_set_d(r20638, theta, MPFR_RNDN);
        mpfr_sin(r20639, r20638, MPFR_RNDN);
        mpfr_set_d(r20640, delta, MPFR_RNDN);
        mpfr_sin(r20641, r20640, MPFR_RNDN);
        mpfr_mul(r20642, r20639, r20641, MPFR_RNDN);
        mpfr_set_d(r20643, phi1, MPFR_RNDN);
        mpfr_cos(r20644, r20643, MPFR_RNDN);
        mpfr_mul(r20645, r20642, r20644, MPFR_RNDN);
        mpfr_cos(r20646, r20640, MPFR_RNDN);
        mpfr_sin(r20647, r20643, MPFR_RNDN);
        mpfr_mul(r20648, r20647, r20646, MPFR_RNDN);
        mpfr_mul(r20649, r20644, r20641, MPFR_RNDN);
        mpfr_cos(r20650, r20638, MPFR_RNDN);
        mpfr_mul(r20651, r20649, r20650, MPFR_RNDN);
        mpfr_add(r20652, r20648, r20651, MPFR_RNDN);
        mpfr_asin(r20653, r20652, MPFR_RNDN);
        mpfr_sin(r20654, r20653, MPFR_RNDN);
        mpfr_mul(r20655, r20647, r20654, MPFR_RNDN);
        mpfr_sub(r20656, r20646, r20655, MPFR_RNDN);
        mpfr_atan2(r20657, r20645, r20656, MPFR_RNDN);
        mpfr_add(r20658, r20637, r20657, MPFR_RNDN);
        return mpfr_get_d(r20658, MPFR_RNDN);
}

static mpfr_t r20659, r20660, r20661, r20662, r20663, r20664, r20665, r20666, r20667, r20668, r20669, r20670, r20671, r20672, r20673, r20674, r20675, r20676, r20677, r20678, r20679, r20680, r20681, r20682, r20683, r20684, r20685;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20659);
        mpfr_init(r20660);
        mpfr_init(r20661);
        mpfr_init(r20662);
        mpfr_init(r20663);
        mpfr_init(r20664);
        mpfr_init(r20665);
        mpfr_init(r20666);
        mpfr_init(r20667);
        mpfr_init(r20668);
        mpfr_init(r20669);
        mpfr_init(r20670);
        mpfr_init(r20671);
        mpfr_init(r20672);
        mpfr_init(r20673);
        mpfr_init(r20674);
        mpfr_init(r20675);
        mpfr_init(r20676);
        mpfr_init_set_str(r20677, "3", 10, MPFR_RNDN);
        mpfr_init(r20678);
        mpfr_init(r20679);
        mpfr_init(r20680);
        mpfr_init(r20681);
        mpfr_init(r20682);
        mpfr_init(r20683);
        mpfr_init(r20684);
        mpfr_init(r20685);
}

double f_fm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r20659, lambda1, MPFR_RNDN);
        mpfr_set_d(r20660, theta, MPFR_RNDN);
        mpfr_sin(r20661, r20660, MPFR_RNDN);
        mpfr_set_d(r20662, delta, MPFR_RNDN);
        mpfr_sin(r20663, r20662, MPFR_RNDN);
        mpfr_mul(r20664, r20661, r20663, MPFR_RNDN);
        mpfr_set_d(r20665, phi1, MPFR_RNDN);
        mpfr_cos(r20666, r20665, MPFR_RNDN);
        mpfr_mul(r20667, r20664, r20666, MPFR_RNDN);
        mpfr_cos(r20668, r20662, MPFR_RNDN);
        mpfr_sin(r20669, r20665, MPFR_RNDN);
        mpfr_mul(r20670, r20668, r20669, MPFR_RNDN);
        mpfr_cos(r20671, r20660, MPFR_RNDN);
        mpfr_mul(r20672, r20663, r20666, MPFR_RNDN);
        mpfr_mul(r20673, r20671, r20672, MPFR_RNDN);
        mpfr_add(r20674, r20670, r20673, MPFR_RNDN);
        mpfr_asin(r20675, r20674, MPFR_RNDN);
        mpfr_sin(r20676, r20675, MPFR_RNDN);
        ;
        mpfr_pow(r20678, r20669, r20677, MPFR_RNDN);
        mpfr_mul(r20679, r20676, r20678, MPFR_RNDN);
        mpfr_mul(r20680, r20676, r20676, MPFR_RNDN);
        mpfr_mul(r20681, r20679, r20680, MPFR_RNDN);
        mpfr_cbrt(r20682, r20681, MPFR_RNDN);
        mpfr_sub(r20683, r20668, r20682, MPFR_RNDN);
        mpfr_atan2(r20684, r20667, r20683, MPFR_RNDN);
        mpfr_add(r20685, r20659, r20684, MPFR_RNDN);
        return mpfr_get_d(r20685, MPFR_RNDN);
}

static mpfr_t r20686, r20687, r20688, r20689, r20690, r20691, r20692, r20693, r20694, r20695, r20696, r20697, r20698, r20699, r20700, r20701, r20702, r20703, r20704, r20705, r20706, r20707, r20708, r20709, r20710, r20711, r20712;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20686);
        mpfr_init(r20687);
        mpfr_init(r20688);
        mpfr_init(r20689);
        mpfr_init(r20690);
        mpfr_init(r20691);
        mpfr_init(r20692);
        mpfr_init(r20693);
        mpfr_init(r20694);
        mpfr_init(r20695);
        mpfr_init(r20696);
        mpfr_init(r20697);
        mpfr_init(r20698);
        mpfr_init(r20699);
        mpfr_init(r20700);
        mpfr_init(r20701);
        mpfr_init(r20702);
        mpfr_init(r20703);
        mpfr_init_set_str(r20704, "3", 10, MPFR_RNDN);
        mpfr_init(r20705);
        mpfr_init(r20706);
        mpfr_init(r20707);
        mpfr_init(r20708);
        mpfr_init(r20709);
        mpfr_init(r20710);
        mpfr_init(r20711);
        mpfr_init(r20712);
}

double f_dm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r20686, lambda1, MPFR_RNDN);
        mpfr_set_d(r20687, theta, MPFR_RNDN);
        mpfr_sin(r20688, r20687, MPFR_RNDN);
        mpfr_set_d(r20689, delta, MPFR_RNDN);
        mpfr_sin(r20690, r20689, MPFR_RNDN);
        mpfr_mul(r20691, r20688, r20690, MPFR_RNDN);
        mpfr_set_d(r20692, phi1, MPFR_RNDN);
        mpfr_cos(r20693, r20692, MPFR_RNDN);
        mpfr_mul(r20694, r20691, r20693, MPFR_RNDN);
        mpfr_cos(r20695, r20689, MPFR_RNDN);
        mpfr_sin(r20696, r20692, MPFR_RNDN);
        mpfr_mul(r20697, r20695, r20696, MPFR_RNDN);
        mpfr_cos(r20698, r20687, MPFR_RNDN);
        mpfr_mul(r20699, r20690, r20693, MPFR_RNDN);
        mpfr_mul(r20700, r20698, r20699, MPFR_RNDN);
        mpfr_add(r20701, r20697, r20700, MPFR_RNDN);
        mpfr_asin(r20702, r20701, MPFR_RNDN);
        mpfr_sin(r20703, r20702, MPFR_RNDN);
        ;
        mpfr_pow(r20705, r20696, r20704, MPFR_RNDN);
        mpfr_mul(r20706, r20703, r20705, MPFR_RNDN);
        mpfr_mul(r20707, r20703, r20703, MPFR_RNDN);
        mpfr_mul(r20708, r20706, r20707, MPFR_RNDN);
        mpfr_cbrt(r20709, r20708, MPFR_RNDN);
        mpfr_sub(r20710, r20695, r20709, MPFR_RNDN);
        mpfr_atan2(r20711, r20694, r20710, MPFR_RNDN);
        mpfr_add(r20712, r20686, r20711, MPFR_RNDN);
        return mpfr_get_d(r20712, MPFR_RNDN);
}