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

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

double f_if(float x, float y, float z, float t, float a, float b, float c, float i) {
        float r25250 = x;
        float r25251 = y;
        float r25252 = r25250 * r25251;
        float r25253 = z;
        float r25254 = r25252 + r25253;
        float r25255 = r25254 * r25251;
        float r25256 = 27464.7644705;
        float r25257 = r25255 + r25256;
        float r25258 = r25257 * r25251;
        float r25259 = 230661.510616;
        float r25260 = r25258 + r25259;
        float r25261 = r25260 * r25251;
        float r25262 = t;
        float r25263 = r25261 + r25262;
        float r25264 = a;
        float r25265 = r25251 + r25264;
        float r25266 = r25265 * r25251;
        float r25267 = b;
        float r25268 = r25266 + r25267;
        float r25269 = r25268 * r25251;
        float r25270 = c;
        float r25271 = r25269 + r25270;
        float r25272 = r25271 * r25251;
        float r25273 = i;
        float r25274 = r25272 + r25273;
        float r25275 = r25263 / r25274;
        return r25275;
}

double f_id(double x, double y, double z, double t, double a, double b, double c, double i) {
        double r25276 = x;
        double r25277 = y;
        double r25278 = r25276 * r25277;
        double r25279 = z;
        double r25280 = r25278 + r25279;
        double r25281 = r25280 * r25277;
        double r25282 = 27464.7644705;
        double r25283 = r25281 + r25282;
        double r25284 = r25283 * r25277;
        double r25285 = 230661.510616;
        double r25286 = r25284 + r25285;
        double r25287 = r25286 * r25277;
        double r25288 = t;
        double r25289 = r25287 + r25288;
        double r25290 = a;
        double r25291 = r25277 + r25290;
        double r25292 = r25291 * r25277;
        double r25293 = b;
        double r25294 = r25292 + r25293;
        double r25295 = r25294 * r25277;
        double r25296 = c;
        double r25297 = r25295 + r25296;
        double r25298 = r25297 * r25277;
        double r25299 = i;
        double r25300 = r25298 + r25299;
        double r25301 = r25289 / r25300;
        return r25301;
}


double f_of(float x, float y, float z, float t, float a, float b, float c, float i) {
        float r25302 = x;
        float r25303 = y;
        float r25304 = r25302 * r25303;
        float r25305 = z;
        float r25306 = r25304 + r25305;
        float r25307 = r25306 * r25303;
        float r25308 = 27464.7644705;
        float r25309 = r25307 + r25308;
        float r25310 = r25309 * r25303;
        float r25311 = 230661.510616;
        float r25312 = r25310 + r25311;
        float r25313 = r25312 * r25303;
        float r25314 = t;
        float r25315 = r25313 + r25314;
        float r25316 = -1.7831438797702915e+308;
        bool r25317 = r25315 <= r25316;
        float r25318 = r25305 / r25303;
        float r25319 = r25318 + r25302;
        float r25320 = 1.779498167574997e+308;
        bool r25321 = r25315 <= r25320;
        float r25322 = a;
        float r25323 = r25303 + r25322;
        float r25324 = r25323 * r25303;
        float r25325 = b;
        float r25326 = r25324 + r25325;
        float r25327 = r25326 * r25303;
        float r25328 = c;
        float r25329 = r25327 + r25328;
        float r25330 = r25329 * r25303;
        float r25331 = i;
        float r25332 = r25330 + r25331;
        float r25333 = r25315 / r25332;
        float r25334 = r25321 ? r25333 : r25319;
        float r25335 = r25317 ? r25319 : r25334;
        return r25335;
}

double f_od(double x, double y, double z, double t, double a, double b, double c, double i) {
        double r25336 = x;
        double r25337 = y;
        double r25338 = r25336 * r25337;
        double r25339 = z;
        double r25340 = r25338 + r25339;
        double r25341 = r25340 * r25337;
        double r25342 = 27464.7644705;
        double r25343 = r25341 + r25342;
        double r25344 = r25343 * r25337;
        double r25345 = 230661.510616;
        double r25346 = r25344 + r25345;
        double r25347 = r25346 * r25337;
        double r25348 = t;
        double r25349 = r25347 + r25348;
        double r25350 = -1.7831438797702915e+308;
        bool r25351 = r25349 <= r25350;
        double r25352 = r25339 / r25337;
        double r25353 = r25352 + r25336;
        double r25354 = 1.779498167574997e+308;
        bool r25355 = r25349 <= r25354;
        double r25356 = a;
        double r25357 = r25337 + r25356;
        double r25358 = r25357 * r25337;
        double r25359 = b;
        double r25360 = r25358 + r25359;
        double r25361 = r25360 * r25337;
        double r25362 = c;
        double r25363 = r25361 + r25362;
        double r25364 = r25363 * r25337;
        double r25365 = i;
        double r25366 = r25364 + r25365;
        double r25367 = r25349 / r25366;
        double r25368 = r25355 ? r25367 : r25353;
        double r25369 = r25351 ? r25353 : r25368;
        return r25369;
}

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 r25370, r25371, r25372, r25373, r25374, r25375, r25376, r25377, r25378, r25379, r25380, r25381, r25382, r25383, r25384, r25385, r25386, r25387, r25388, r25389, r25390, r25391, r25392, r25393, r25394, r25395;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r25370);
        mpfr_init(r25371);
        mpfr_init(r25372);
        mpfr_init(r25373);
        mpfr_init(r25374);
        mpfr_init(r25375);
        mpfr_init_set_str(r25376, "27464.7644705", 10, MPFR_RNDN);
        mpfr_init(r25377);
        mpfr_init(r25378);
        mpfr_init_set_str(r25379, "230661.510616", 10, MPFR_RNDN);
        mpfr_init(r25380);
        mpfr_init(r25381);
        mpfr_init(r25382);
        mpfr_init(r25383);
        mpfr_init(r25384);
        mpfr_init(r25385);
        mpfr_init(r25386);
        mpfr_init(r25387);
        mpfr_init(r25388);
        mpfr_init(r25389);
        mpfr_init(r25390);
        mpfr_init(r25391);
        mpfr_init(r25392);
        mpfr_init(r25393);
        mpfr_init(r25394);
        mpfr_init(r25395);
}

double f_im(double x, double y, double z, double t, double a, double b, double c, double i) {
        mpfr_set_d(r25370, x, MPFR_RNDN);
        mpfr_set_d(r25371, y, MPFR_RNDN);
        mpfr_mul(r25372, r25370, r25371, MPFR_RNDN);
        mpfr_set_d(r25373, z, MPFR_RNDN);
        mpfr_add(r25374, r25372, r25373, MPFR_RNDN);
        mpfr_mul(r25375, r25374, r25371, MPFR_RNDN);
        ;
        mpfr_add(r25377, r25375, r25376, MPFR_RNDN);
        mpfr_mul(r25378, r25377, r25371, MPFR_RNDN);
        ;
        mpfr_add(r25380, r25378, r25379, MPFR_RNDN);
        mpfr_mul(r25381, r25380, r25371, MPFR_RNDN);
        mpfr_set_d(r25382, t, MPFR_RNDN);
        mpfr_add(r25383, r25381, r25382, MPFR_RNDN);
        mpfr_set_d(r25384, a, MPFR_RNDN);
        mpfr_add(r25385, r25371, r25384, MPFR_RNDN);
        mpfr_mul(r25386, r25385, r25371, MPFR_RNDN);
        mpfr_set_d(r25387, b, MPFR_RNDN);
        mpfr_add(r25388, r25386, r25387, MPFR_RNDN);
        mpfr_mul(r25389, r25388, r25371, MPFR_RNDN);
        mpfr_set_d(r25390, c, MPFR_RNDN);
        mpfr_add(r25391, r25389, r25390, MPFR_RNDN);
        mpfr_mul(r25392, r25391, r25371, MPFR_RNDN);
        mpfr_set_d(r25393, i, MPFR_RNDN);
        mpfr_add(r25394, r25392, r25393, MPFR_RNDN);
        mpfr_div(r25395, r25383, r25394, MPFR_RNDN);
        return mpfr_get_d(r25395, MPFR_RNDN);
}

static mpfr_t r25396, r25397, r25398, r25399, r25400, r25401, r25402, r25403, r25404, r25405, r25406, r25407, r25408, r25409, r25410, r25411, r25412, r25413, r25414, r25415, r25416, r25417, r25418, r25419, r25420, r25421, r25422, r25423, r25424, r25425, r25426, r25427, r25428, r25429;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25396);
        mpfr_init(r25397);
        mpfr_init(r25398);
        mpfr_init(r25399);
        mpfr_init(r25400);
        mpfr_init(r25401);
        mpfr_init_set_str(r25402, "27464.7644705", 10, MPFR_RNDN);
        mpfr_init(r25403);
        mpfr_init(r25404);
        mpfr_init_set_str(r25405, "230661.510616", 10, MPFR_RNDN);
        mpfr_init(r25406);
        mpfr_init(r25407);
        mpfr_init(r25408);
        mpfr_init(r25409);
        mpfr_init_set_str(r25410, "-1.7831438797702915e+308", 10, MPFR_RNDN);
        mpfr_init(r25411);
        mpfr_init(r25412);
        mpfr_init(r25413);
        mpfr_init_set_str(r25414, "1.779498167574997e+308", 10, MPFR_RNDN);
        mpfr_init(r25415);
        mpfr_init(r25416);
        mpfr_init(r25417);
        mpfr_init(r25418);
        mpfr_init(r25419);
        mpfr_init(r25420);
        mpfr_init(r25421);
        mpfr_init(r25422);
        mpfr_init(r25423);
        mpfr_init(r25424);
        mpfr_init(r25425);
        mpfr_init(r25426);
        mpfr_init(r25427);
        mpfr_init(r25428);
        mpfr_init(r25429);
}

double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) {
        mpfr_set_d(r25396, x, MPFR_RNDN);
        mpfr_set_d(r25397, y, MPFR_RNDN);
        mpfr_mul(r25398, r25396, r25397, MPFR_RNDN);
        mpfr_set_d(r25399, z, MPFR_RNDN);
        mpfr_add(r25400, r25398, r25399, MPFR_RNDN);
        mpfr_mul(r25401, r25400, r25397, MPFR_RNDN);
        ;
        mpfr_add(r25403, r25401, r25402, MPFR_RNDN);
        mpfr_mul(r25404, r25403, r25397, MPFR_RNDN);
        ;
        mpfr_add(r25406, r25404, r25405, MPFR_RNDN);
        mpfr_mul(r25407, r25406, r25397, MPFR_RNDN);
        mpfr_set_d(r25408, t, MPFR_RNDN);
        mpfr_add(r25409, r25407, r25408, MPFR_RNDN);
        ;
        mpfr_set_si(r25411, mpfr_cmp(r25409, r25410) <= 0, MPFR_RNDN);
        mpfr_div(r25412, r25399, r25397, MPFR_RNDN);
        mpfr_add(r25413, r25412, r25396, MPFR_RNDN);
        ;
        mpfr_set_si(r25415, mpfr_cmp(r25409, r25414) <= 0, MPFR_RNDN);
        mpfr_set_d(r25416, a, MPFR_RNDN);
        mpfr_add(r25417, r25397, r25416, MPFR_RNDN);
        mpfr_mul(r25418, r25417, r25397, MPFR_RNDN);
        mpfr_set_d(r25419, b, MPFR_RNDN);
        mpfr_add(r25420, r25418, r25419, MPFR_RNDN);
        mpfr_mul(r25421, r25420, r25397, MPFR_RNDN);
        mpfr_set_d(r25422, c, MPFR_RNDN);
        mpfr_add(r25423, r25421, r25422, MPFR_RNDN);
        mpfr_mul(r25424, r25423, r25397, MPFR_RNDN);
        mpfr_set_d(r25425, i, MPFR_RNDN);
        mpfr_add(r25426, r25424, r25425, MPFR_RNDN);
        mpfr_div(r25427, r25409, r25426, MPFR_RNDN);
        if (mpfr_get_si(r25415, MPFR_RNDN)) { mpfr_set(r25428, r25427, MPFR_RNDN); } else { mpfr_set(r25428, r25413, MPFR_RNDN); };
        if (mpfr_get_si(r25411, MPFR_RNDN)) { mpfr_set(r25429, r25413, MPFR_RNDN); } else { mpfr_set(r25429, r25428, MPFR_RNDN); };
        return mpfr_get_d(r25429, MPFR_RNDN);
}

static mpfr_t r25430, r25431, r25432, r25433, r25434, r25435, r25436, r25437, r25438, r25439, r25440, r25441, r25442, r25443, r25444, r25445, r25446, r25447, r25448, r25449, r25450, r25451, r25452, r25453, r25454, r25455, r25456, r25457, r25458, r25459, r25460, r25461, r25462, r25463;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25430);
        mpfr_init(r25431);
        mpfr_init(r25432);
        mpfr_init(r25433);
        mpfr_init(r25434);
        mpfr_init(r25435);
        mpfr_init_set_str(r25436, "27464.7644705", 10, MPFR_RNDN);
        mpfr_init(r25437);
        mpfr_init(r25438);
        mpfr_init_set_str(r25439, "230661.510616", 10, MPFR_RNDN);
        mpfr_init(r25440);
        mpfr_init(r25441);
        mpfr_init(r25442);
        mpfr_init(r25443);
        mpfr_init_set_str(r25444, "-1.7831438797702915e+308", 10, MPFR_RNDN);
        mpfr_init(r25445);
        mpfr_init(r25446);
        mpfr_init(r25447);
        mpfr_init_set_str(r25448, "1.779498167574997e+308", 10, MPFR_RNDN);
        mpfr_init(r25449);
        mpfr_init(r25450);
        mpfr_init(r25451);
        mpfr_init(r25452);
        mpfr_init(r25453);
        mpfr_init(r25454);
        mpfr_init(r25455);
        mpfr_init(r25456);
        mpfr_init(r25457);
        mpfr_init(r25458);
        mpfr_init(r25459);
        mpfr_init(r25460);
        mpfr_init(r25461);
        mpfr_init(r25462);
        mpfr_init(r25463);
}

double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) {
        mpfr_set_d(r25430, x, MPFR_RNDN);
        mpfr_set_d(r25431, y, MPFR_RNDN);
        mpfr_mul(r25432, r25430, r25431, MPFR_RNDN);
        mpfr_set_d(r25433, z, MPFR_RNDN);
        mpfr_add(r25434, r25432, r25433, MPFR_RNDN);
        mpfr_mul(r25435, r25434, r25431, MPFR_RNDN);
        ;
        mpfr_add(r25437, r25435, r25436, MPFR_RNDN);
        mpfr_mul(r25438, r25437, r25431, MPFR_RNDN);
        ;
        mpfr_add(r25440, r25438, r25439, MPFR_RNDN);
        mpfr_mul(r25441, r25440, r25431, MPFR_RNDN);
        mpfr_set_d(r25442, t, MPFR_RNDN);
        mpfr_add(r25443, r25441, r25442, MPFR_RNDN);
        ;
        mpfr_set_si(r25445, mpfr_cmp(r25443, r25444) <= 0, MPFR_RNDN);
        mpfr_div(r25446, r25433, r25431, MPFR_RNDN);
        mpfr_add(r25447, r25446, r25430, MPFR_RNDN);
        ;
        mpfr_set_si(r25449, mpfr_cmp(r25443, r25448) <= 0, MPFR_RNDN);
        mpfr_set_d(r25450, a, MPFR_RNDN);
        mpfr_add(r25451, r25431, r25450, MPFR_RNDN);
        mpfr_mul(r25452, r25451, r25431, MPFR_RNDN);
        mpfr_set_d(r25453, b, MPFR_RNDN);
        mpfr_add(r25454, r25452, r25453, MPFR_RNDN);
        mpfr_mul(r25455, r25454, r25431, MPFR_RNDN);
        mpfr_set_d(r25456, c, MPFR_RNDN);
        mpfr_add(r25457, r25455, r25456, MPFR_RNDN);
        mpfr_mul(r25458, r25457, r25431, MPFR_RNDN);
        mpfr_set_d(r25459, i, MPFR_RNDN);
        mpfr_add(r25460, r25458, r25459, MPFR_RNDN);
        mpfr_div(r25461, r25443, r25460, MPFR_RNDN);
        if (mpfr_get_si(r25449, MPFR_RNDN)) { mpfr_set(r25462, r25461, MPFR_RNDN); } else { mpfr_set(r25462, r25447, MPFR_RNDN); };
        if (mpfr_get_si(r25445, MPFR_RNDN)) { mpfr_set(r25463, r25447, MPFR_RNDN); } else { mpfr_set(r25463, r25462, MPFR_RNDN); };
        return mpfr_get_d(r25463, MPFR_RNDN);
}