#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 r9266 = x;
        float r9267 = y;
        float r9268 = r9266 * r9267;
        float r9269 = z;
        float r9270 = r9268 + r9269;
        float r9271 = r9270 * r9267;
        float r9272 = 27464.7644705;
        float r9273 = r9271 + r9272;
        float r9274 = r9273 * r9267;
        float r9275 = 230661.510616;
        float r9276 = r9274 + r9275;
        float r9277 = r9276 * r9267;
        float r9278 = t;
        float r9279 = r9277 + r9278;
        float r9280 = a;
        float r9281 = r9267 + r9280;
        float r9282 = r9281 * r9267;
        float r9283 = b;
        float r9284 = r9282 + r9283;
        float r9285 = r9284 * r9267;
        float r9286 = c;
        float r9287 = r9285 + r9286;
        float r9288 = r9287 * r9267;
        float r9289 = i;
        float r9290 = r9288 + r9289;
        float r9291 = r9279 / r9290;
        return r9291;
}

double f_id(double x, double y, double z, double t, double a, double b, double c, double i) {
        double r9292 = x;
        double r9293 = y;
        double r9294 = r9292 * r9293;
        double r9295 = z;
        double r9296 = r9294 + r9295;
        double r9297 = r9296 * r9293;
        double r9298 = 27464.7644705;
        double r9299 = r9297 + r9298;
        double r9300 = r9299 * r9293;
        double r9301 = 230661.510616;
        double r9302 = r9300 + r9301;
        double r9303 = r9302 * r9293;
        double r9304 = t;
        double r9305 = r9303 + r9304;
        double r9306 = a;
        double r9307 = r9293 + r9306;
        double r9308 = r9307 * r9293;
        double r9309 = b;
        double r9310 = r9308 + r9309;
        double r9311 = r9310 * r9293;
        double r9312 = c;
        double r9313 = r9311 + r9312;
        double r9314 = r9313 * r9293;
        double r9315 = i;
        double r9316 = r9314 + r9315;
        double r9317 = r9305 / r9316;
        return r9317;
}


double f_of(float x, float y, float z, float t, float a, float b, float c, float i) {
        float r9318 = x;
        float r9319 = y;
        float r9320 = r9318 * r9319;
        float r9321 = z;
        float r9322 = r9320 + r9321;
        float r9323 = r9322 * r9319;
        float r9324 = 27464.7644705;
        float r9325 = r9323 + r9324;
        float r9326 = r9325 * r9319;
        float r9327 = 230661.510616;
        float r9328 = r9326 + r9327;
        float r9329 = r9328 * r9319;
        float r9330 = t;
        float r9331 = r9329 + r9330;
        float r9332 = a;
        float r9333 = r9319 + r9332;
        float r9334 = r9333 * r9319;
        float r9335 = b;
        float r9336 = r9334 + r9335;
        float r9337 = r9336 * r9319;
        float r9338 = c;
        float r9339 = r9337 + r9338;
        float r9340 = r9339 * r9319;
        float r9341 = i;
        float r9342 = r9340 + r9341;
        float r9343 = r9331 / r9342;
        return r9343;
}

double f_od(double x, double y, double z, double t, double a, double b, double c, double i) {
        double r9344 = x;
        double r9345 = y;
        double r9346 = r9344 * r9345;
        double r9347 = z;
        double r9348 = r9346 + r9347;
        double r9349 = r9348 * r9345;
        double r9350 = 27464.7644705;
        double r9351 = r9349 + r9350;
        double r9352 = r9351 * r9345;
        double r9353 = 230661.510616;
        double r9354 = r9352 + r9353;
        double r9355 = r9354 * r9345;
        double r9356 = t;
        double r9357 = r9355 + r9356;
        double r9358 = a;
        double r9359 = r9345 + r9358;
        double r9360 = r9359 * r9345;
        double r9361 = b;
        double r9362 = r9360 + r9361;
        double r9363 = r9362 * r9345;
        double r9364 = c;
        double r9365 = r9363 + r9364;
        double r9366 = r9365 * r9345;
        double r9367 = i;
        double r9368 = r9366 + r9367;
        double r9369 = r9357 / r9368;
        return r9369;
}

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 r9370, r9371, r9372, r9373, r9374, r9375, r9376, r9377, r9378, r9379, r9380, r9381, r9382, r9383, r9384, r9385, r9386, r9387, r9388, r9389, r9390, r9391, r9392, r9393, r9394, r9395;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9370);
        mpfr_init(r9371);
        mpfr_init(r9372);
        mpfr_init(r9373);
        mpfr_init(r9374);
        mpfr_init(r9375);
        mpfr_init_set_str(r9376, "27464.7644705", 10, MPFR_RNDN);
        mpfr_init(r9377);
        mpfr_init(r9378);
        mpfr_init_set_str(r9379, "230661.510616", 10, MPFR_RNDN);
        mpfr_init(r9380);
        mpfr_init(r9381);
        mpfr_init(r9382);
        mpfr_init(r9383);
        mpfr_init(r9384);
        mpfr_init(r9385);
        mpfr_init(r9386);
        mpfr_init(r9387);
        mpfr_init(r9388);
        mpfr_init(r9389);
        mpfr_init(r9390);
        mpfr_init(r9391);
        mpfr_init(r9392);
        mpfr_init(r9393);
        mpfr_init(r9394);
        mpfr_init(r9395);
}

double f_im(double x, double y, double z, double t, double a, double b, double c, double i) {
        mpfr_set_d(r9370, x, MPFR_RNDN);
        mpfr_set_d(r9371, y, MPFR_RNDN);
        mpfr_mul(r9372, r9370, r9371, MPFR_RNDN);
        mpfr_set_d(r9373, z, MPFR_RNDN);
        mpfr_add(r9374, r9372, r9373, MPFR_RNDN);
        mpfr_mul(r9375, r9374, r9371, MPFR_RNDN);
        ;
        mpfr_add(r9377, r9375, r9376, MPFR_RNDN);
        mpfr_mul(r9378, r9377, r9371, MPFR_RNDN);
        ;
        mpfr_add(r9380, r9378, r9379, MPFR_RNDN);
        mpfr_mul(r9381, r9380, r9371, MPFR_RNDN);
        mpfr_set_d(r9382, t, MPFR_RNDN);
        mpfr_add(r9383, r9381, r9382, MPFR_RNDN);
        mpfr_set_d(r9384, a, MPFR_RNDN);
        mpfr_add(r9385, r9371, r9384, MPFR_RNDN);
        mpfr_mul(r9386, r9385, r9371, MPFR_RNDN);
        mpfr_set_d(r9387, b, MPFR_RNDN);
        mpfr_add(r9388, r9386, r9387, MPFR_RNDN);
        mpfr_mul(r9389, r9388, r9371, MPFR_RNDN);
        mpfr_set_d(r9390, c, MPFR_RNDN);
        mpfr_add(r9391, r9389, r9390, MPFR_RNDN);
        mpfr_mul(r9392, r9391, r9371, MPFR_RNDN);
        mpfr_set_d(r9393, i, MPFR_RNDN);
        mpfr_add(r9394, r9392, r9393, MPFR_RNDN);
        mpfr_div(r9395, r9383, r9394, MPFR_RNDN);
        return mpfr_get_d(r9395, MPFR_RNDN);
}

static mpfr_t r9396, r9397, r9398, r9399, r9400, r9401, r9402, r9403, r9404, r9405, r9406, r9407, r9408, r9409, r9410, r9411, r9412, r9413, r9414, r9415, r9416, r9417, r9418, r9419, r9420, r9421;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9396);
        mpfr_init(r9397);
        mpfr_init(r9398);
        mpfr_init(r9399);
        mpfr_init(r9400);
        mpfr_init(r9401);
        mpfr_init_set_str(r9402, "27464.7644705", 10, MPFR_RNDN);
        mpfr_init(r9403);
        mpfr_init(r9404);
        mpfr_init_set_str(r9405, "230661.510616", 10, MPFR_RNDN);
        mpfr_init(r9406);
        mpfr_init(r9407);
        mpfr_init(r9408);
        mpfr_init(r9409);
        mpfr_init(r9410);
        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);
}

double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) {
        mpfr_set_d(r9396, x, MPFR_RNDN);
        mpfr_set_d(r9397, y, MPFR_RNDN);
        mpfr_mul(r9398, r9396, r9397, MPFR_RNDN);
        mpfr_set_d(r9399, z, MPFR_RNDN);
        mpfr_add(r9400, r9398, r9399, MPFR_RNDN);
        mpfr_mul(r9401, r9400, r9397, MPFR_RNDN);
        ;
        mpfr_add(r9403, r9401, r9402, MPFR_RNDN);
        mpfr_mul(r9404, r9403, r9397, MPFR_RNDN);
        ;
        mpfr_add(r9406, r9404, r9405, MPFR_RNDN);
        mpfr_mul(r9407, r9406, r9397, MPFR_RNDN);
        mpfr_set_d(r9408, t, MPFR_RNDN);
        mpfr_add(r9409, r9407, r9408, MPFR_RNDN);
        mpfr_set_d(r9410, a, MPFR_RNDN);
        mpfr_add(r9411, r9397, r9410, MPFR_RNDN);
        mpfr_mul(r9412, r9411, r9397, MPFR_RNDN);
        mpfr_set_d(r9413, b, MPFR_RNDN);
        mpfr_add(r9414, r9412, r9413, MPFR_RNDN);
        mpfr_mul(r9415, r9414, r9397, MPFR_RNDN);
        mpfr_set_d(r9416, c, MPFR_RNDN);
        mpfr_add(r9417, r9415, r9416, MPFR_RNDN);
        mpfr_mul(r9418, r9417, r9397, MPFR_RNDN);
        mpfr_set_d(r9419, i, MPFR_RNDN);
        mpfr_add(r9420, r9418, r9419, MPFR_RNDN);
        mpfr_div(r9421, r9409, r9420, MPFR_RNDN);
        return mpfr_get_d(r9421, MPFR_RNDN);
}

static mpfr_t r9422, r9423, r9424, r9425, r9426, r9427, r9428, r9429, r9430, r9431, r9432, r9433, r9434, r9435, r9436, r9437, r9438, r9439, r9440, r9441, r9442, r9443, r9444, r9445, r9446, r9447;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9422);
        mpfr_init(r9423);
        mpfr_init(r9424);
        mpfr_init(r9425);
        mpfr_init(r9426);
        mpfr_init(r9427);
        mpfr_init_set_str(r9428, "27464.7644705", 10, MPFR_RNDN);
        mpfr_init(r9429);
        mpfr_init(r9430);
        mpfr_init_set_str(r9431, "230661.510616", 10, MPFR_RNDN);
        mpfr_init(r9432);
        mpfr_init(r9433);
        mpfr_init(r9434);
        mpfr_init(r9435);
        mpfr_init(r9436);
        mpfr_init(r9437);
        mpfr_init(r9438);
        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);
}

double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) {
        mpfr_set_d(r9422, x, MPFR_RNDN);
        mpfr_set_d(r9423, y, MPFR_RNDN);
        mpfr_mul(r9424, r9422, r9423, MPFR_RNDN);
        mpfr_set_d(r9425, z, MPFR_RNDN);
        mpfr_add(r9426, r9424, r9425, MPFR_RNDN);
        mpfr_mul(r9427, r9426, r9423, MPFR_RNDN);
        ;
        mpfr_add(r9429, r9427, r9428, MPFR_RNDN);
        mpfr_mul(r9430, r9429, r9423, MPFR_RNDN);
        ;
        mpfr_add(r9432, r9430, r9431, MPFR_RNDN);
        mpfr_mul(r9433, r9432, r9423, MPFR_RNDN);
        mpfr_set_d(r9434, t, MPFR_RNDN);
        mpfr_add(r9435, r9433, r9434, MPFR_RNDN);
        mpfr_set_d(r9436, a, MPFR_RNDN);
        mpfr_add(r9437, r9423, r9436, MPFR_RNDN);
        mpfr_mul(r9438, r9437, r9423, MPFR_RNDN);
        mpfr_set_d(r9439, b, MPFR_RNDN);
        mpfr_add(r9440, r9438, r9439, MPFR_RNDN);
        mpfr_mul(r9441, r9440, r9423, MPFR_RNDN);
        mpfr_set_d(r9442, c, MPFR_RNDN);
        mpfr_add(r9443, r9441, r9442, MPFR_RNDN);
        mpfr_mul(r9444, r9443, r9423, MPFR_RNDN);
        mpfr_set_d(r9445, i, MPFR_RNDN);
        mpfr_add(r9446, r9444, r9445, MPFR_RNDN);
        mpfr_div(r9447, r9435, r9446, MPFR_RNDN);
        return mpfr_get_d(r9447, MPFR_RNDN);
}