#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 r27319 = x;
        float r27320 = y;
        float r27321 = r27319 + r27320;
        float r27322 = log(r27321);
        float r27323 = z;
        float r27324 = log(r27323);
        float r27325 = r27322 + r27324;
        float r27326 = t;
        float r27327 = r27325 - r27326;
        float r27328 = a;
        float r27329 = 0.5;
        float r27330 = r27328 - r27329;
        float r27331 = log(r27326);
        float r27332 = r27330 * r27331;
        float r27333 = r27327 + r27332;
        return r27333;
}

double f_id(double x, double y, double z, double t, double a) {
        double r27334 = x;
        double r27335 = y;
        double r27336 = r27334 + r27335;
        double r27337 = log(r27336);
        double r27338 = z;
        double r27339 = log(r27338);
        double r27340 = r27337 + r27339;
        double r27341 = t;
        double r27342 = r27340 - r27341;
        double r27343 = a;
        double r27344 = 0.5;
        double r27345 = r27343 - r27344;
        double r27346 = log(r27341);
        double r27347 = r27345 * r27346;
        double r27348 = r27342 + r27347;
        return r27348;
}


double f_of(float x, float y, float z, float t, float a) {
        float r27349 = x;
        float r27350 = y;
        float r27351 = r27349 + r27350;
        float r27352 = cbrt(r27351);
        float r27353 = r27352 * r27352;
        float r27354 = log(r27353);
        float r27355 = log(r27352);
        float r27356 = z;
        float r27357 = log(r27356);
        float r27358 = r27355 + r27357;
        float r27359 = r27354 + r27358;
        float r27360 = t;
        float r27361 = r27359 - r27360;
        float r27362 = a;
        float r27363 = 0.5;
        float r27364 = r27362 - r27363;
        float r27365 = log(r27360);
        float r27366 = r27364 * r27365;
        float r27367 = r27361 + r27366;
        return r27367;
}

double f_od(double x, double y, double z, double t, double a) {
        double r27368 = x;
        double r27369 = y;
        double r27370 = r27368 + r27369;
        double r27371 = cbrt(r27370);
        double r27372 = r27371 * r27371;
        double r27373 = log(r27372);
        double r27374 = log(r27371);
        double r27375 = z;
        double r27376 = log(r27375);
        double r27377 = r27374 + r27376;
        double r27378 = r27373 + r27377;
        double r27379 = t;
        double r27380 = r27378 - r27379;
        double r27381 = a;
        double r27382 = 0.5;
        double r27383 = r27381 - r27382;
        double r27384 = log(r27379);
        double r27385 = r27383 * r27384;
        double r27386 = r27380 + r27385;
        return r27386;
}

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 r27387, r27388, r27389, r27390, r27391, r27392, r27393, r27394, r27395, r27396, r27397, r27398, r27399, r27400, r27401;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r27387);
        mpfr_init(r27388);
        mpfr_init(r27389);
        mpfr_init(r27390);
        mpfr_init(r27391);
        mpfr_init(r27392);
        mpfr_init(r27393);
        mpfr_init(r27394);
        mpfr_init(r27395);
        mpfr_init(r27396);
        mpfr_init_set_str(r27397, "0.5", 10, MPFR_RNDN);
        mpfr_init(r27398);
        mpfr_init(r27399);
        mpfr_init(r27400);
        mpfr_init(r27401);
}

double f_im(double x, double y, double z, double t, double a) {
        mpfr_set_d(r27387, x, MPFR_RNDN);
        mpfr_set_d(r27388, y, MPFR_RNDN);
        mpfr_add(r27389, r27387, r27388, MPFR_RNDN);
        mpfr_log(r27390, r27389, MPFR_RNDN);
        mpfr_set_d(r27391, z, MPFR_RNDN);
        mpfr_log(r27392, r27391, MPFR_RNDN);
        mpfr_add(r27393, r27390, r27392, MPFR_RNDN);
        mpfr_set_d(r27394, t, MPFR_RNDN);
        mpfr_sub(r27395, r27393, r27394, MPFR_RNDN);
        mpfr_set_d(r27396, a, MPFR_RNDN);
        ;
        mpfr_sub(r27398, r27396, r27397, MPFR_RNDN);
        mpfr_log(r27399, r27394, MPFR_RNDN);
        mpfr_mul(r27400, r27398, r27399, MPFR_RNDN);
        mpfr_add(r27401, r27395, r27400, MPFR_RNDN);
        return mpfr_get_d(r27401, MPFR_RNDN);
}

static mpfr_t r27402, r27403, r27404, r27405, r27406, r27407, r27408, r27409, r27410, r27411, r27412, r27413, r27414, r27415, r27416, r27417, r27418, r27419, r27420;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r27402);
        mpfr_init(r27403);
        mpfr_init(r27404);
        mpfr_init(r27405);
        mpfr_init(r27406);
        mpfr_init(r27407);
        mpfr_init(r27408);
        mpfr_init(r27409);
        mpfr_init(r27410);
        mpfr_init(r27411);
        mpfr_init(r27412);
        mpfr_init(r27413);
        mpfr_init(r27414);
        mpfr_init(r27415);
        mpfr_init_set_str(r27416, "0.5", 10, MPFR_RNDN);
        mpfr_init(r27417);
        mpfr_init(r27418);
        mpfr_init(r27419);
        mpfr_init(r27420);
}

double f_fm(double x, double y, double z, double t, double a) {
        mpfr_set_d(r27402, x, MPFR_RNDN);
        mpfr_set_d(r27403, y, MPFR_RNDN);
        mpfr_add(r27404, r27402, r27403, MPFR_RNDN);
        mpfr_cbrt(r27405, r27404, MPFR_RNDN);
        mpfr_mul(r27406, r27405, r27405, MPFR_RNDN);
        mpfr_log(r27407, r27406, MPFR_RNDN);
        mpfr_log(r27408, r27405, MPFR_RNDN);
        mpfr_set_d(r27409, z, MPFR_RNDN);
        mpfr_log(r27410, r27409, MPFR_RNDN);
        mpfr_add(r27411, r27408, r27410, MPFR_RNDN);
        mpfr_add(r27412, r27407, r27411, MPFR_RNDN);
        mpfr_set_d(r27413, t, MPFR_RNDN);
        mpfr_sub(r27414, r27412, r27413, MPFR_RNDN);
        mpfr_set_d(r27415, a, MPFR_RNDN);
        ;
        mpfr_sub(r27417, r27415, r27416, MPFR_RNDN);
        mpfr_log(r27418, r27413, MPFR_RNDN);
        mpfr_mul(r27419, r27417, r27418, MPFR_RNDN);
        mpfr_add(r27420, r27414, r27419, MPFR_RNDN);
        return mpfr_get_d(r27420, MPFR_RNDN);
}

static mpfr_t r27421, r27422, r27423, r27424, r27425, r27426, r27427, r27428, r27429, r27430, r27431, r27432, r27433, r27434, r27435, r27436, r27437, r27438, r27439;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r27421);
        mpfr_init(r27422);
        mpfr_init(r27423);
        mpfr_init(r27424);
        mpfr_init(r27425);
        mpfr_init(r27426);
        mpfr_init(r27427);
        mpfr_init(r27428);
        mpfr_init(r27429);
        mpfr_init(r27430);
        mpfr_init(r27431);
        mpfr_init(r27432);
        mpfr_init(r27433);
        mpfr_init(r27434);
        mpfr_init_set_str(r27435, "0.5", 10, MPFR_RNDN);
        mpfr_init(r27436);
        mpfr_init(r27437);
        mpfr_init(r27438);
        mpfr_init(r27439);
}

double f_dm(double x, double y, double z, double t, double a) {
        mpfr_set_d(r27421, x, MPFR_RNDN);
        mpfr_set_d(r27422, y, MPFR_RNDN);
        mpfr_add(r27423, r27421, r27422, MPFR_RNDN);
        mpfr_cbrt(r27424, r27423, MPFR_RNDN);
        mpfr_mul(r27425, r27424, r27424, MPFR_RNDN);
        mpfr_log(r27426, r27425, MPFR_RNDN);
        mpfr_log(r27427, r27424, MPFR_RNDN);
        mpfr_set_d(r27428, z, MPFR_RNDN);
        mpfr_log(r27429, r27428, MPFR_RNDN);
        mpfr_add(r27430, r27427, r27429, MPFR_RNDN);
        mpfr_add(r27431, r27426, r27430, MPFR_RNDN);
        mpfr_set_d(r27432, t, MPFR_RNDN);
        mpfr_sub(r27433, r27431, r27432, MPFR_RNDN);
        mpfr_set_d(r27434, a, MPFR_RNDN);
        ;
        mpfr_sub(r27436, r27434, r27435, MPFR_RNDN);
        mpfr_log(r27437, r27432, MPFR_RNDN);
        mpfr_mul(r27438, r27436, r27437, MPFR_RNDN);
        mpfr_add(r27439, r27433, r27438, MPFR_RNDN);
        return mpfr_get_d(r27439, MPFR_RNDN);
}