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

char *name = "Numeric.SpecFunctions:logGammaL from math-functions-0.1.5.2";

double f_if(float x, float y, float z, float t, float a) {
        float r27307 = x;
        float r27308 = y;
        float r27309 = r27307 + r27308;
        float r27310 = log(r27309);
        float r27311 = z;
        float r27312 = log(r27311);
        float r27313 = r27310 + r27312;
        float r27314 = t;
        float r27315 = r27313 - r27314;
        float r27316 = a;
        float r27317 = 0.5;
        float r27318 = r27316 - r27317;
        float r27319 = log(r27314);
        float r27320 = r27318 * r27319;
        float r27321 = r27315 + r27320;
        return r27321;
}

double f_id(double x, double y, double z, double t, double a) {
        double r27322 = x;
        double r27323 = y;
        double r27324 = r27322 + r27323;
        double r27325 = log(r27324);
        double r27326 = z;
        double r27327 = log(r27326);
        double r27328 = r27325 + r27327;
        double r27329 = t;
        double r27330 = r27328 - r27329;
        double r27331 = a;
        double r27332 = 0.5;
        double r27333 = r27331 - r27332;
        double r27334 = log(r27329);
        double r27335 = r27333 * r27334;
        double r27336 = r27330 + r27335;
        return r27336;
}


double f_of(float x, float y, float z, float t, float a) {
        float r27337 = x;
        float r27338 = y;
        float r27339 = r27337 + r27338;
        float r27340 = cbrt(r27339);
        float r27341 = r27340 * r27340;
        float r27342 = log(r27341);
        float r27343 = log(r27340);
        float r27344 = z;
        float r27345 = log(r27344);
        float r27346 = r27343 + r27345;
        float r27347 = r27342 + r27346;
        float r27348 = t;
        float r27349 = r27347 - r27348;
        float r27350 = a;
        float r27351 = 0.5;
        float r27352 = r27350 - r27351;
        float r27353 = log(r27348);
        float r27354 = r27352 * r27353;
        float r27355 = r27349 + r27354;
        return r27355;
}

double f_od(double x, double y, double z, double t, double a) {
        double r27356 = x;
        double r27357 = y;
        double r27358 = r27356 + r27357;
        double r27359 = cbrt(r27358);
        double r27360 = r27359 * r27359;
        double r27361 = log(r27360);
        double r27362 = log(r27359);
        double r27363 = z;
        double r27364 = log(r27363);
        double r27365 = r27362 + r27364;
        double r27366 = r27361 + r27365;
        double r27367 = t;
        double r27368 = r27366 - r27367;
        double r27369 = a;
        double r27370 = 0.5;
        double r27371 = r27369 - r27370;
        double r27372 = log(r27367);
        double r27373 = r27371 * r27372;
        double r27374 = r27368 + r27373;
        return r27374;
}

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 r27375, r27376, r27377, r27378, r27379, r27380, r27381, r27382, r27383, r27384, r27385, r27386, r27387, r27388, r27389;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r27375);
        mpfr_init(r27376);
        mpfr_init(r27377);
        mpfr_init(r27378);
        mpfr_init(r27379);
        mpfr_init(r27380);
        mpfr_init(r27381);
        mpfr_init(r27382);
        mpfr_init(r27383);
        mpfr_init(r27384);
        mpfr_init_set_str(r27385, "0.5", 10, MPFR_RNDN);
        mpfr_init(r27386);
        mpfr_init(r27387);
        mpfr_init(r27388);
        mpfr_init(r27389);
}

double f_im(double x, double y, double z, double t, double a) {
        mpfr_set_d(r27375, x, MPFR_RNDN);
        mpfr_set_d(r27376, y, MPFR_RNDN);
        mpfr_add(r27377, r27375, r27376, MPFR_RNDN);
        mpfr_log(r27378, r27377, MPFR_RNDN);
        mpfr_set_d(r27379, z, MPFR_RNDN);
        mpfr_log(r27380, r27379, MPFR_RNDN);
        mpfr_add(r27381, r27378, r27380, MPFR_RNDN);
        mpfr_set_d(r27382, t, MPFR_RNDN);
        mpfr_sub(r27383, r27381, r27382, MPFR_RNDN);
        mpfr_set_d(r27384, a, MPFR_RNDN);
        ;
        mpfr_sub(r27386, r27384, r27385, MPFR_RNDN);
        mpfr_log(r27387, r27382, MPFR_RNDN);
        mpfr_mul(r27388, r27386, r27387, MPFR_RNDN);
        mpfr_add(r27389, r27383, r27388, MPFR_RNDN);
        return mpfr_get_d(r27389, MPFR_RNDN);
}

static mpfr_t r27390, r27391, r27392, r27393, r27394, r27395, r27396, r27397, r27398, r27399, r27400, r27401, r27402, r27403, r27404, r27405, r27406, r27407, r27408;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r27390);
        mpfr_init(r27391);
        mpfr_init(r27392);
        mpfr_init(r27393);
        mpfr_init(r27394);
        mpfr_init(r27395);
        mpfr_init(r27396);
        mpfr_init(r27397);
        mpfr_init(r27398);
        mpfr_init(r27399);
        mpfr_init(r27400);
        mpfr_init(r27401);
        mpfr_init(r27402);
        mpfr_init(r27403);
        mpfr_init_set_str(r27404, "0.5", 10, MPFR_RNDN);
        mpfr_init(r27405);
        mpfr_init(r27406);
        mpfr_init(r27407);
        mpfr_init(r27408);
}

double f_fm(double x, double y, double z, double t, double a) {
        mpfr_set_d(r27390, x, MPFR_RNDN);
        mpfr_set_d(r27391, y, MPFR_RNDN);
        mpfr_add(r27392, r27390, r27391, MPFR_RNDN);
        mpfr_cbrt(r27393, r27392, MPFR_RNDN);
        mpfr_mul(r27394, r27393, r27393, MPFR_RNDN);
        mpfr_log(r27395, r27394, MPFR_RNDN);
        mpfr_log(r27396, r27393, MPFR_RNDN);
        mpfr_set_d(r27397, z, MPFR_RNDN);
        mpfr_log(r27398, r27397, MPFR_RNDN);
        mpfr_add(r27399, r27396, r27398, MPFR_RNDN);
        mpfr_add(r27400, r27395, r27399, MPFR_RNDN);
        mpfr_set_d(r27401, t, MPFR_RNDN);
        mpfr_sub(r27402, r27400, r27401, MPFR_RNDN);
        mpfr_set_d(r27403, a, MPFR_RNDN);
        ;
        mpfr_sub(r27405, r27403, r27404, MPFR_RNDN);
        mpfr_log(r27406, r27401, MPFR_RNDN);
        mpfr_mul(r27407, r27405, r27406, MPFR_RNDN);
        mpfr_add(r27408, r27402, r27407, MPFR_RNDN);
        return mpfr_get_d(r27408, MPFR_RNDN);
}

static mpfr_t r27409, r27410, r27411, r27412, r27413, r27414, r27415, r27416, r27417, r27418, r27419, r27420, r27421, r27422, r27423, r27424, r27425, r27426, r27427;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r27409);
        mpfr_init(r27410);
        mpfr_init(r27411);
        mpfr_init(r27412);
        mpfr_init(r27413);
        mpfr_init(r27414);
        mpfr_init(r27415);
        mpfr_init(r27416);
        mpfr_init(r27417);
        mpfr_init(r27418);
        mpfr_init(r27419);
        mpfr_init(r27420);
        mpfr_init(r27421);
        mpfr_init(r27422);
        mpfr_init_set_str(r27423, "0.5", 10, MPFR_RNDN);
        mpfr_init(r27424);
        mpfr_init(r27425);
        mpfr_init(r27426);
        mpfr_init(r27427);
}

double f_dm(double x, double y, double z, double t, double a) {
        mpfr_set_d(r27409, x, MPFR_RNDN);
        mpfr_set_d(r27410, y, MPFR_RNDN);
        mpfr_add(r27411, r27409, r27410, MPFR_RNDN);
        mpfr_cbrt(r27412, r27411, MPFR_RNDN);
        mpfr_mul(r27413, r27412, r27412, MPFR_RNDN);
        mpfr_log(r27414, r27413, MPFR_RNDN);
        mpfr_log(r27415, r27412, MPFR_RNDN);
        mpfr_set_d(r27416, z, MPFR_RNDN);
        mpfr_log(r27417, r27416, MPFR_RNDN);
        mpfr_add(r27418, r27415, r27417, MPFR_RNDN);
        mpfr_add(r27419, r27414, r27418, MPFR_RNDN);
        mpfr_set_d(r27420, t, MPFR_RNDN);
        mpfr_sub(r27421, r27419, r27420, MPFR_RNDN);
        mpfr_set_d(r27422, a, MPFR_RNDN);
        ;
        mpfr_sub(r27424, r27422, r27423, MPFR_RNDN);
        mpfr_log(r27425, r27420, MPFR_RNDN);
        mpfr_mul(r27426, r27424, r27425, MPFR_RNDN);
        mpfr_add(r27427, r27421, r27426, MPFR_RNDN);
        return mpfr_get_d(r27427, MPFR_RNDN);
}