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

char *name = "Numeric.SpecFunctions:invIncompleteGamma from math-functions-0.1.5.2, C";

double f_if(float x) {
        float r43497 = 2.30753;
        float r43498 = x;
        float r43499 = 0.27061;
        float r43500 = r43498 * r43499;
        float r43501 = r43497 + r43500;
        float r43502 = 1.0;
        float r43503 = 0.99229;
        float r43504 = 0.04481;
        float r43505 = r43498 * r43504;
        float r43506 = r43503 + r43505;
        float r43507 = r43498 * r43506;
        float r43508 = r43502 + r43507;
        float r43509 = r43501 / r43508;
        float r43510 = r43509 - r43498;
        return r43510;
}

double f_id(double x) {
        double r43511 = 2.30753;
        double r43512 = x;
        double r43513 = 0.27061;
        double r43514 = r43512 * r43513;
        double r43515 = r43511 + r43514;
        double r43516 = 1.0;
        double r43517 = 0.99229;
        double r43518 = 0.04481;
        double r43519 = r43512 * r43518;
        double r43520 = r43517 + r43519;
        double r43521 = r43512 * r43520;
        double r43522 = r43516 + r43521;
        double r43523 = r43515 / r43522;
        double r43524 = r43523 - r43512;
        return r43524;
}


double f_of(float x) {
        float r43525 = 2.30753;
        float r43526 = x;
        float r43527 = 0.27061;
        float r43528 = r43526 * r43527;
        float r43529 = r43525 + r43528;
        float r43530 = 1.0;
        float r43531 = 0.99229;
        float r43532 = 0.04481;
        float r43533 = r43526 * r43532;
        float r43534 = r43531 + r43533;
        float r43535 = r43526 * r43534;
        float r43536 = r43530 + r43535;
        float r43537 = r43529 / r43536;
        float r43538 = r43537 - r43526;
        return r43538;
}

double f_od(double x) {
        double r43539 = 2.30753;
        double r43540 = x;
        double r43541 = 0.27061;
        double r43542 = r43540 * r43541;
        double r43543 = r43539 + r43542;
        double r43544 = 1.0;
        double r43545 = 0.99229;
        double r43546 = 0.04481;
        double r43547 = r43540 * r43546;
        double r43548 = r43545 + r43547;
        double r43549 = r43540 * r43548;
        double r43550 = r43544 + r43549;
        double r43551 = r43543 / r43550;
        double r43552 = r43551 - r43540;
        return r43552;
}

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 r43553, r43554, r43555, r43556, r43557, r43558, r43559, r43560, r43561, r43562, r43563, r43564, r43565, r43566;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r43553, "2.30753", 10, MPFR_RNDN);
        mpfr_init(r43554);
        mpfr_init_set_str(r43555, "0.27061", 10, MPFR_RNDN);
        mpfr_init(r43556);
        mpfr_init(r43557);
        mpfr_init_set_str(r43558, "1.0", 10, MPFR_RNDN);
        mpfr_init_set_str(r43559, "0.99229", 10, MPFR_RNDN);
        mpfr_init_set_str(r43560, "0.04481", 10, MPFR_RNDN);
        mpfr_init(r43561);
        mpfr_init(r43562);
        mpfr_init(r43563);
        mpfr_init(r43564);
        mpfr_init(r43565);
        mpfr_init(r43566);
}

double f_im(double x) {
        ;
        mpfr_set_d(r43554, x, MPFR_RNDN);
        ;
        mpfr_mul(r43556, r43554, r43555, MPFR_RNDN);
        mpfr_add(r43557, r43553, r43556, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_mul(r43561, r43554, r43560, MPFR_RNDN);
        mpfr_add(r43562, r43559, r43561, MPFR_RNDN);
        mpfr_mul(r43563, r43554, r43562, MPFR_RNDN);
        mpfr_add(r43564, r43558, r43563, MPFR_RNDN);
        mpfr_div(r43565, r43557, r43564, MPFR_RNDN);
        mpfr_sub(r43566, r43565, r43554, MPFR_RNDN);
        return mpfr_get_d(r43566, MPFR_RNDN);
}

static mpfr_t r43567, r43568, r43569, r43570, r43571, r43572, r43573, r43574, r43575, r43576, r43577, r43578, r43579, r43580;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r43567, "2.30753", 10, MPFR_RNDN);
        mpfr_init(r43568);
        mpfr_init_set_str(r43569, "0.27061", 10, MPFR_RNDN);
        mpfr_init(r43570);
        mpfr_init(r43571);
        mpfr_init_set_str(r43572, "1.0", 10, MPFR_RNDN);
        mpfr_init_set_str(r43573, "0.99229", 10, MPFR_RNDN);
        mpfr_init_set_str(r43574, "0.04481", 10, MPFR_RNDN);
        mpfr_init(r43575);
        mpfr_init(r43576);
        mpfr_init(r43577);
        mpfr_init(r43578);
        mpfr_init(r43579);
        mpfr_init(r43580);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r43568, x, MPFR_RNDN);
        ;
        mpfr_mul(r43570, r43568, r43569, MPFR_RNDN);
        mpfr_add(r43571, r43567, r43570, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_mul(r43575, r43568, r43574, MPFR_RNDN);
        mpfr_add(r43576, r43573, r43575, MPFR_RNDN);
        mpfr_mul(r43577, r43568, r43576, MPFR_RNDN);
        mpfr_add(r43578, r43572, r43577, MPFR_RNDN);
        mpfr_div(r43579, r43571, r43578, MPFR_RNDN);
        mpfr_sub(r43580, r43579, r43568, MPFR_RNDN);
        return mpfr_get_d(r43580, MPFR_RNDN);
}

static mpfr_t r43581, r43582, r43583, r43584, r43585, r43586, r43587, r43588, r43589, r43590, r43591, r43592, r43593, r43594;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r43581, "2.30753", 10, MPFR_RNDN);
        mpfr_init(r43582);
        mpfr_init_set_str(r43583, "0.27061", 10, MPFR_RNDN);
        mpfr_init(r43584);
        mpfr_init(r43585);
        mpfr_init_set_str(r43586, "1.0", 10, MPFR_RNDN);
        mpfr_init_set_str(r43587, "0.99229", 10, MPFR_RNDN);
        mpfr_init_set_str(r43588, "0.04481", 10, MPFR_RNDN);
        mpfr_init(r43589);
        mpfr_init(r43590);
        mpfr_init(r43591);
        mpfr_init(r43592);
        mpfr_init(r43593);
        mpfr_init(r43594);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r43582, x, MPFR_RNDN);
        ;
        mpfr_mul(r43584, r43582, r43583, MPFR_RNDN);
        mpfr_add(r43585, r43581, r43584, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_mul(r43589, r43582, r43588, MPFR_RNDN);
        mpfr_add(r43590, r43587, r43589, MPFR_RNDN);
        mpfr_mul(r43591, r43582, r43590, MPFR_RNDN);
        mpfr_add(r43592, r43586, r43591, MPFR_RNDN);
        mpfr_div(r43593, r43585, r43592, MPFR_RNDN);
        mpfr_sub(r43594, r43593, r43582, MPFR_RNDN);
        return mpfr_get_d(r43594, MPFR_RNDN);
}