#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 r9289 = lambda1;
        float r9290 = theta;
        float r9291 = sin(r9290);
        float r9292 = delta;
        float r9293 = sin(r9292);
        float r9294 = r9291 * r9293;
        float r9295 = phi1;
        float r9296 = cos(r9295);
        float r9297 = r9294 * r9296;
        float r9298 = cos(r9292);
        float r9299 = sin(r9295);
        float r9300 = r9299 * r9298;
        float r9301 = r9296 * r9293;
        float r9302 = cos(r9290);
        float r9303 = r9301 * r9302;
        float r9304 = r9300 + r9303;
        float r9305 = asin(r9304);
        float r9306 = sin(r9305);
        float r9307 = r9299 * r9306;
        float r9308 = r9298 - r9307;
        float r9309 = atan2(r9297, r9308);
        float r9310 = r9289 + r9309;
        return r9310;
}

double f_id(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r9311 = lambda1;
        double r9312 = theta;
        double r9313 = sin(r9312);
        double r9314 = delta;
        double r9315 = sin(r9314);
        double r9316 = r9313 * r9315;
        double r9317 = phi1;
        double r9318 = cos(r9317);
        double r9319 = r9316 * r9318;
        double r9320 = cos(r9314);
        double r9321 = sin(r9317);
        double r9322 = r9321 * r9320;
        double r9323 = r9318 * r9315;
        double r9324 = cos(r9312);
        double r9325 = r9323 * r9324;
        double r9326 = r9322 + r9325;
        double r9327 = asin(r9326);
        double r9328 = sin(r9327);
        double r9329 = r9321 * r9328;
        double r9330 = r9320 - r9329;
        double r9331 = atan2(r9319, r9330);
        double r9332 = r9311 + r9331;
        return r9332;
}


double f_of(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) {
        float r9333 = phi1;
        float r9334 = cos(r9333);
        float r9335 = theta;
        float r9336 = sin(r9335);
        float r9337 = r9334 * r9336;
        float r9338 = delta;
        float r9339 = sin(r9338);
        float r9340 = r9337 * r9339;
        float r9341 = cos(r9338);
        float r9342 = r9341 * r9341;
        float r9343 = cos(r9335);
        float r9344 = r9343 * r9334;
        float r9345 = r9344 * r9339;
        float r9346 = sin(r9333);
        float r9347 = r9345 * r9346;
        float r9348 = 2.0f;
        float r9349 = pow(r9346, r9348);
        float r9350 = r9349 * r9341;
        float r9351 = r9347 + r9350;
        float r9352 = r9351 * r9350;
        float r9353 = r9342 - r9352;
        float r9354 = r9351 * r9347;
        float r9355 = r9353 - r9354;
        float r9356 = r9341 + r9351;
        float r9357 = r9355 / r9356;
        float r9358 = atan2(r9340, r9357);
        float r9359 = lambda1;
        float r9360 = r9358 + r9359;
        return r9360;
}

double f_od(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r9361 = phi1;
        double r9362 = cos(r9361);
        double r9363 = theta;
        double r9364 = sin(r9363);
        double r9365 = r9362 * r9364;
        double r9366 = delta;
        double r9367 = sin(r9366);
        double r9368 = r9365 * r9367;
        double r9369 = cos(r9366);
        double r9370 = r9369 * r9369;
        double r9371 = cos(r9363);
        double r9372 = r9371 * r9362;
        double r9373 = r9372 * r9367;
        double r9374 = sin(r9361);
        double r9375 = r9373 * r9374;
        double r9376 = 2.0;
        double r9377 = pow(r9374, r9376);
        double r9378 = r9377 * r9369;
        double r9379 = r9375 + r9378;
        double r9380 = r9379 * r9378;
        double r9381 = r9370 - r9380;
        double r9382 = r9379 * r9375;
        double r9383 = r9381 - r9382;
        double r9384 = r9369 + r9379;
        double r9385 = r9383 / r9384;
        double r9386 = atan2(r9368, r9385);
        double r9387 = lambda1;
        double r9388 = r9386 + r9387;
        return r9388;
}

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 r9389, r9390, r9391, r9392, r9393, r9394, r9395, r9396, r9397, r9398, r9399, r9400, r9401, r9402, r9403, r9404, r9405, r9406, r9407, r9408, r9409, r9410;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9389);
        mpfr_init(r9390);
        mpfr_init(r9391);
        mpfr_init(r9392);
        mpfr_init(r9393);
        mpfr_init(r9394);
        mpfr_init(r9395);
        mpfr_init(r9396);
        mpfr_init(r9397);
        mpfr_init(r9398);
        mpfr_init(r9399);
        mpfr_init(r9400);
        mpfr_init(r9401);
        mpfr_init(r9402);
        mpfr_init(r9403);
        mpfr_init(r9404);
        mpfr_init(r9405);
        mpfr_init(r9406);
        mpfr_init(r9407);
        mpfr_init(r9408);
        mpfr_init(r9409);
        mpfr_init(r9410);
}

double f_im(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r9389, lambda1, MPFR_RNDN);
        mpfr_set_d(r9390, theta, MPFR_RNDN);
        mpfr_sin(r9391, r9390, MPFR_RNDN);
        mpfr_set_d(r9392, delta, MPFR_RNDN);
        mpfr_sin(r9393, r9392, MPFR_RNDN);
        mpfr_mul(r9394, r9391, r9393, MPFR_RNDN);
        mpfr_set_d(r9395, phi1, MPFR_RNDN);
        mpfr_cos(r9396, r9395, MPFR_RNDN);
        mpfr_mul(r9397, r9394, r9396, MPFR_RNDN);
        mpfr_cos(r9398, r9392, MPFR_RNDN);
        mpfr_sin(r9399, r9395, MPFR_RNDN);
        mpfr_mul(r9400, r9399, r9398, MPFR_RNDN);
        mpfr_mul(r9401, r9396, r9393, MPFR_RNDN);
        mpfr_cos(r9402, r9390, MPFR_RNDN);
        mpfr_mul(r9403, r9401, r9402, MPFR_RNDN);
        mpfr_add(r9404, r9400, r9403, MPFR_RNDN);
        mpfr_asin(r9405, r9404, MPFR_RNDN);
        mpfr_sin(r9406, r9405, MPFR_RNDN);
        mpfr_mul(r9407, r9399, r9406, MPFR_RNDN);
        mpfr_sub(r9408, r9398, r9407, MPFR_RNDN);
        mpfr_atan2(r9409, r9397, r9408, MPFR_RNDN);
        mpfr_add(r9410, r9389, r9409, MPFR_RNDN);
        return mpfr_get_d(r9410, MPFR_RNDN);
}

static mpfr_t r9411, r9412, r9413, r9414, r9415, r9416, r9417, r9418, r9419, r9420, r9421, r9422, r9423, r9424, r9425, r9426, r9427, r9428, r9429, r9430, r9431, r9432, r9433, r9434, r9435, r9436, r9437, r9438;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9411);
        mpfr_init(r9412);
        mpfr_init(r9413);
        mpfr_init(r9414);
        mpfr_init(r9415);
        mpfr_init(r9416);
        mpfr_init(r9417);
        mpfr_init(r9418);
        mpfr_init(r9419);
        mpfr_init(r9420);
        mpfr_init(r9421);
        mpfr_init(r9422);
        mpfr_init(r9423);
        mpfr_init(r9424);
        mpfr_init(r9425);
        mpfr_init_set_str(r9426, "2", 10, MPFR_RNDN);
        mpfr_init(r9427);
        mpfr_init(r9428);
        mpfr_init(r9429);
        mpfr_init(r9430);
        mpfr_init(r9431);
        mpfr_init(r9432);
        mpfr_init(r9433);
        mpfr_init(r9434);
        mpfr_init(r9435);
        mpfr_init(r9436);
        mpfr_init(r9437);
        mpfr_init(r9438);
}

double f_fm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r9411, phi1, MPFR_RNDN);
        mpfr_cos(r9412, r9411, MPFR_RNDN);
        mpfr_set_d(r9413, theta, MPFR_RNDN);
        mpfr_sin(r9414, r9413, MPFR_RNDN);
        mpfr_mul(r9415, r9412, r9414, MPFR_RNDN);
        mpfr_set_d(r9416, delta, MPFR_RNDN);
        mpfr_sin(r9417, r9416, MPFR_RNDN);
        mpfr_mul(r9418, r9415, r9417, MPFR_RNDN);
        mpfr_cos(r9419, r9416, MPFR_RNDN);
        mpfr_mul(r9420, r9419, r9419, MPFR_RNDN);
        mpfr_cos(r9421, r9413, MPFR_RNDN);
        mpfr_mul(r9422, r9421, r9412, MPFR_RNDN);
        mpfr_mul(r9423, r9422, r9417, MPFR_RNDN);
        mpfr_sin(r9424, r9411, MPFR_RNDN);
        mpfr_mul(r9425, r9423, r9424, MPFR_RNDN);
        ;
        mpfr_pow(r9427, r9424, r9426, MPFR_RNDN);
        mpfr_mul(r9428, r9427, r9419, MPFR_RNDN);
        mpfr_add(r9429, r9425, r9428, MPFR_RNDN);
        mpfr_mul(r9430, r9429, r9428, MPFR_RNDN);
        mpfr_sub(r9431, r9420, r9430, MPFR_RNDN);
        mpfr_mul(r9432, r9429, r9425, MPFR_RNDN);
        mpfr_sub(r9433, r9431, r9432, MPFR_RNDN);
        mpfr_add(r9434, r9419, r9429, MPFR_RNDN);
        mpfr_div(r9435, r9433, r9434, MPFR_RNDN);
        mpfr_atan2(r9436, r9418, r9435, MPFR_RNDN);
        mpfr_set_d(r9437, lambda1, MPFR_RNDN);
        mpfr_add(r9438, r9436, r9437, MPFR_RNDN);
        return mpfr_get_d(r9438, MPFR_RNDN);
}

static mpfr_t r9439, r9440, r9441, r9442, r9443, r9444, r9445, r9446, r9447, r9448, r9449, r9450, r9451, r9452, r9453, r9454, r9455, r9456, r9457, r9458, r9459, r9460, r9461, r9462, r9463, r9464, r9465, r9466;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9439);
        mpfr_init(r9440);
        mpfr_init(r9441);
        mpfr_init(r9442);
        mpfr_init(r9443);
        mpfr_init(r9444);
        mpfr_init(r9445);
        mpfr_init(r9446);
        mpfr_init(r9447);
        mpfr_init(r9448);
        mpfr_init(r9449);
        mpfr_init(r9450);
        mpfr_init(r9451);
        mpfr_init(r9452);
        mpfr_init(r9453);
        mpfr_init_set_str(r9454, "2", 10, MPFR_RNDN);
        mpfr_init(r9455);
        mpfr_init(r9456);
        mpfr_init(r9457);
        mpfr_init(r9458);
        mpfr_init(r9459);
        mpfr_init(r9460);
        mpfr_init(r9461);
        mpfr_init(r9462);
        mpfr_init(r9463);
        mpfr_init(r9464);
        mpfr_init(r9465);
        mpfr_init(r9466);
}

double f_dm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r9439, phi1, MPFR_RNDN);
        mpfr_cos(r9440, r9439, MPFR_RNDN);
        mpfr_set_d(r9441, theta, MPFR_RNDN);
        mpfr_sin(r9442, r9441, MPFR_RNDN);
        mpfr_mul(r9443, r9440, r9442, MPFR_RNDN);
        mpfr_set_d(r9444, delta, MPFR_RNDN);
        mpfr_sin(r9445, r9444, MPFR_RNDN);
        mpfr_mul(r9446, r9443, r9445, MPFR_RNDN);
        mpfr_cos(r9447, r9444, MPFR_RNDN);
        mpfr_mul(r9448, r9447, r9447, MPFR_RNDN);
        mpfr_cos(r9449, r9441, MPFR_RNDN);
        mpfr_mul(r9450, r9449, r9440, MPFR_RNDN);
        mpfr_mul(r9451, r9450, r9445, MPFR_RNDN);
        mpfr_sin(r9452, r9439, MPFR_RNDN);
        mpfr_mul(r9453, r9451, r9452, MPFR_RNDN);
        ;
        mpfr_pow(r9455, r9452, r9454, MPFR_RNDN);
        mpfr_mul(r9456, r9455, r9447, MPFR_RNDN);
        mpfr_add(r9457, r9453, r9456, MPFR_RNDN);
        mpfr_mul(r9458, r9457, r9456, MPFR_RNDN);
        mpfr_sub(r9459, r9448, r9458, MPFR_RNDN);
        mpfr_mul(r9460, r9457, r9453, MPFR_RNDN);
        mpfr_sub(r9461, r9459, r9460, MPFR_RNDN);
        mpfr_add(r9462, r9447, r9457, MPFR_RNDN);
        mpfr_div(r9463, r9461, r9462, MPFR_RNDN);
        mpfr_atan2(r9464, r9446, r9463, MPFR_RNDN);
        mpfr_set_d(r9465, lambda1, MPFR_RNDN);
        mpfr_add(r9466, r9464, r9465, MPFR_RNDN);
        return mpfr_get_d(r9466, MPFR_RNDN);
}