#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 r26200 = lambda1;
        float r26201 = theta;
        float r26202 = sin(r26201);
        float r26203 = delta;
        float r26204 = sin(r26203);
        float r26205 = r26202 * r26204;
        float r26206 = phi1;
        float r26207 = cos(r26206);
        float r26208 = r26205 * r26207;
        float r26209 = cos(r26203);
        float r26210 = sin(r26206);
        float r26211 = r26210 * r26209;
        float r26212 = r26207 * r26204;
        float r26213 = cos(r26201);
        float r26214 = r26212 * r26213;
        float r26215 = r26211 + r26214;
        float r26216 = asin(r26215);
        float r26217 = sin(r26216);
        float r26218 = r26210 * r26217;
        float r26219 = r26209 - r26218;
        float r26220 = atan2(r26208, r26219);
        float r26221 = r26200 + r26220;
        return r26221;
}

double f_id(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r26222 = lambda1;
        double r26223 = theta;
        double r26224 = sin(r26223);
        double r26225 = delta;
        double r26226 = sin(r26225);
        double r26227 = r26224 * r26226;
        double r26228 = phi1;
        double r26229 = cos(r26228);
        double r26230 = r26227 * r26229;
        double r26231 = cos(r26225);
        double r26232 = sin(r26228);
        double r26233 = r26232 * r26231;
        double r26234 = r26229 * r26226;
        double r26235 = cos(r26223);
        double r26236 = r26234 * r26235;
        double r26237 = r26233 + r26236;
        double r26238 = asin(r26237);
        double r26239 = sin(r26238);
        double r26240 = r26232 * r26239;
        double r26241 = r26231 - r26240;
        double r26242 = atan2(r26230, r26241);
        double r26243 = r26222 + r26242;
        return r26243;
}


double f_of(float lambda1, float phi1, float __attribute__((unused)) phi2, float delta, float theta) {
        float r26244 = lambda1;
        float r26245 = phi1;
        float r26246 = cos(r26245);
        float r26247 = delta;
        float r26248 = sin(r26247);
        float r26249 = r26246 * r26248;
        float r26250 = expm1(r26249);
        float r26251 = log1p(r26250);
        float r26252 = theta;
        float r26253 = sin(r26252);
        float r26254 = r26251 * r26253;
        float r26255 = sin(r26245);
        float r26256 = -r26255;
        float r26257 = cos(r26252);
        float r26258 = r26257 * r26248;
        float r26259 = cos(r26247);
        float r26260 = r26259 * r26255;
        float r26261 = fma(r26246, r26258, r26260);
        float r26262 = fma(r26256, r26261, r26259);
        float r26263 = atan2(r26254, r26262);
        float r26264 = r26244 + r26263;
        return r26264;
}

double f_od(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        double r26265 = lambda1;
        double r26266 = phi1;
        double r26267 = cos(r26266);
        double r26268 = delta;
        double r26269 = sin(r26268);
        double r26270 = r26267 * r26269;
        double r26271 = expm1(r26270);
        double r26272 = log1p(r26271);
        double r26273 = theta;
        double r26274 = sin(r26273);
        double r26275 = r26272 * r26274;
        double r26276 = sin(r26266);
        double r26277 = -r26276;
        double r26278 = cos(r26273);
        double r26279 = r26278 * r26269;
        double r26280 = cos(r26268);
        double r26281 = r26280 * r26276;
        double r26282 = fma(r26267, r26279, r26281);
        double r26283 = fma(r26277, r26282, r26280);
        double r26284 = atan2(r26275, r26283);
        double r26285 = r26265 + r26284;
        return r26285;
}

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 r26286, r26287, r26288, r26289, r26290, r26291, r26292, r26293, r26294, r26295, r26296, r26297, r26298, r26299, r26300, r26301, r26302, r26303, r26304, r26305, r26306, r26307;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26286);
        mpfr_init(r26287);
        mpfr_init(r26288);
        mpfr_init(r26289);
        mpfr_init(r26290);
        mpfr_init(r26291);
        mpfr_init(r26292);
        mpfr_init(r26293);
        mpfr_init(r26294);
        mpfr_init(r26295);
        mpfr_init(r26296);
        mpfr_init(r26297);
        mpfr_init(r26298);
        mpfr_init(r26299);
        mpfr_init(r26300);
        mpfr_init(r26301);
        mpfr_init(r26302);
        mpfr_init(r26303);
        mpfr_init(r26304);
        mpfr_init(r26305);
        mpfr_init(r26306);
        mpfr_init(r26307);
}

double f_im(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r26286, lambda1, MPFR_RNDN);
        mpfr_set_d(r26287, theta, MPFR_RNDN);
        mpfr_sin(r26288, r26287, MPFR_RNDN);
        mpfr_set_d(r26289, delta, MPFR_RNDN);
        mpfr_sin(r26290, r26289, MPFR_RNDN);
        mpfr_mul(r26291, r26288, r26290, MPFR_RNDN);
        mpfr_set_d(r26292, phi1, MPFR_RNDN);
        mpfr_cos(r26293, r26292, MPFR_RNDN);
        mpfr_mul(r26294, r26291, r26293, MPFR_RNDN);
        mpfr_cos(r26295, r26289, MPFR_RNDN);
        mpfr_sin(r26296, r26292, MPFR_RNDN);
        mpfr_mul(r26297, r26296, r26295, MPFR_RNDN);
        mpfr_mul(r26298, r26293, r26290, MPFR_RNDN);
        mpfr_cos(r26299, r26287, MPFR_RNDN);
        mpfr_mul(r26300, r26298, r26299, MPFR_RNDN);
        mpfr_add(r26301, r26297, r26300, MPFR_RNDN);
        mpfr_asin(r26302, r26301, MPFR_RNDN);
        mpfr_sin(r26303, r26302, MPFR_RNDN);
        mpfr_mul(r26304, r26296, r26303, MPFR_RNDN);
        mpfr_sub(r26305, r26295, r26304, MPFR_RNDN);
        mpfr_atan2(r26306, r26294, r26305, MPFR_RNDN);
        mpfr_add(r26307, r26286, r26306, MPFR_RNDN);
        return mpfr_get_d(r26307, MPFR_RNDN);
}

static mpfr_t r26308, r26309, r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318, r26319, r26320, r26321, r26322, r26323, r26324, r26325, r26326, r26327, r26328;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26308);
        mpfr_init(r26309);
        mpfr_init(r26310);
        mpfr_init(r26311);
        mpfr_init(r26312);
        mpfr_init(r26313);
        mpfr_init(r26314);
        mpfr_init(r26315);
        mpfr_init(r26316);
        mpfr_init(r26317);
        mpfr_init(r26318);
        mpfr_init(r26319);
        mpfr_init(r26320);
        mpfr_init(r26321);
        mpfr_init(r26322);
        mpfr_init(r26323);
        mpfr_init(r26324);
        mpfr_init(r26325);
        mpfr_init(r26326);
        mpfr_init(r26327);
        mpfr_init(r26328);
}

double f_fm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r26308, lambda1, MPFR_RNDN);
        mpfr_set_d(r26309, phi1, MPFR_RNDN);
        mpfr_cos(r26310, r26309, MPFR_RNDN);
        mpfr_set_d(r26311, delta, MPFR_RNDN);
        mpfr_sin(r26312, r26311, MPFR_RNDN);
        mpfr_mul(r26313, r26310, r26312, MPFR_RNDN);
        mpfr_expm1(r26314, r26313, MPFR_RNDN);
        mpfr_log1p(r26315, r26314, MPFR_RNDN);
        mpfr_set_d(r26316, theta, MPFR_RNDN);
        mpfr_sin(r26317, r26316, MPFR_RNDN);
        mpfr_mul(r26318, r26315, r26317, MPFR_RNDN);
        mpfr_sin(r26319, r26309, MPFR_RNDN);
        mpfr_neg(r26320, r26319, MPFR_RNDN);
        mpfr_cos(r26321, r26316, MPFR_RNDN);
        mpfr_mul(r26322, r26321, r26312, MPFR_RNDN);
        mpfr_cos(r26323, r26311, MPFR_RNDN);
        mpfr_mul(r26324, r26323, r26319, MPFR_RNDN);
        mpfr_fma(r26325, r26310, r26322, r26324, MPFR_RNDN);
        mpfr_fma(r26326, r26320, r26325, r26323, MPFR_RNDN);
        mpfr_atan2(r26327, r26318, r26326, MPFR_RNDN);
        mpfr_add(r26328, r26308, r26327, MPFR_RNDN);
        return mpfr_get_d(r26328, MPFR_RNDN);
}

static mpfr_t r26329, r26330, r26331, r26332, r26333, r26334, r26335, r26336, r26337, r26338, r26339, r26340, r26341, r26342, r26343, r26344, r26345, r26346, r26347, r26348, r26349;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26329);
        mpfr_init(r26330);
        mpfr_init(r26331);
        mpfr_init(r26332);
        mpfr_init(r26333);
        mpfr_init(r26334);
        mpfr_init(r26335);
        mpfr_init(r26336);
        mpfr_init(r26337);
        mpfr_init(r26338);
        mpfr_init(r26339);
        mpfr_init(r26340);
        mpfr_init(r26341);
        mpfr_init(r26342);
        mpfr_init(r26343);
        mpfr_init(r26344);
        mpfr_init(r26345);
        mpfr_init(r26346);
        mpfr_init(r26347);
        mpfr_init(r26348);
        mpfr_init(r26349);
}

double f_dm(double lambda1, double phi1, double __attribute__((unused)) phi2, double delta, double theta) {
        mpfr_set_d(r26329, lambda1, MPFR_RNDN);
        mpfr_set_d(r26330, phi1, MPFR_RNDN);
        mpfr_cos(r26331, r26330, MPFR_RNDN);
        mpfr_set_d(r26332, delta, MPFR_RNDN);
        mpfr_sin(r26333, r26332, MPFR_RNDN);
        mpfr_mul(r26334, r26331, r26333, MPFR_RNDN);
        mpfr_expm1(r26335, r26334, MPFR_RNDN);
        mpfr_log1p(r26336, r26335, MPFR_RNDN);
        mpfr_set_d(r26337, theta, MPFR_RNDN);
        mpfr_sin(r26338, r26337, MPFR_RNDN);
        mpfr_mul(r26339, r26336, r26338, MPFR_RNDN);
        mpfr_sin(r26340, r26330, MPFR_RNDN);
        mpfr_neg(r26341, r26340, MPFR_RNDN);
        mpfr_cos(r26342, r26337, MPFR_RNDN);
        mpfr_mul(r26343, r26342, r26333, MPFR_RNDN);
        mpfr_cos(r26344, r26332, MPFR_RNDN);
        mpfr_mul(r26345, r26344, r26340, MPFR_RNDN);
        mpfr_fma(r26346, r26331, r26343, r26345, MPFR_RNDN);
        mpfr_fma(r26347, r26341, r26346, r26344, MPFR_RNDN);
        mpfr_atan2(r26348, r26339, r26347, MPFR_RNDN);
        mpfr_add(r26349, r26329, r26348, MPFR_RNDN);
        return mpfr_get_d(r26349, MPFR_RNDN);
}