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

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

double f_if(float x) {
        float r38509 = x;
        float r38510 = r38509 * r38509;
        float r38511 = 2.0f;
        float r38512 = r38510 * r38511;
        float r38513 = 1.0f;
        float r38514 = r38512 - r38513;
        return r38514;
}

double f_id(double x) {
        double r38515 = x;
        double r38516 = r38515 * r38515;
        double r38517 = 2.0;
        double r38518 = r38516 * r38517;
        double r38519 = 1.0;
        double r38520 = r38518 - r38519;
        return r38520;
}


double f_of(float x) {
        float r38521 = x;
        float r38522 = r38521 * r38521;
        float r38523 = 2.0f;
        float r38524 = r38522 * r38523;
        float r38525 = 1.0f;
        float r38526 = r38524 - r38525;
        return r38526;
}

double f_od(double x) {
        double r38527 = x;
        double r38528 = r38527 * r38527;
        double r38529 = 2.0;
        double r38530 = r38528 * r38529;
        double r38531 = 1.0;
        double r38532 = r38530 - r38531;
        return r38532;
}

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 r38533, r38534, r38535, r38536, r38537, r38538;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r38533);
        mpfr_init(r38534);
        mpfr_init_set_str(r38535, "2.0", 10, MPFR_RNDN);
        mpfr_init(r38536);
        mpfr_init_set_str(r38537, "1.0", 10, MPFR_RNDN);
        mpfr_init(r38538);
}

double f_im(double x) {
        mpfr_set_d(r38533, x, MPFR_RNDN);
        mpfr_mul(r38534, r38533, r38533, MPFR_RNDN);
        ;
        mpfr_mul(r38536, r38534, r38535, MPFR_RNDN);
        ;
        mpfr_sub(r38538, r38536, r38537, MPFR_RNDN);
        return mpfr_get_d(r38538, MPFR_RNDN);
}

static mpfr_t r38539, r38540, r38541, r38542, r38543, r38544;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r38539);
        mpfr_init(r38540);
        mpfr_init_set_str(r38541, "2.0", 10, MPFR_RNDN);
        mpfr_init(r38542);
        mpfr_init_set_str(r38543, "1.0", 10, MPFR_RNDN);
        mpfr_init(r38544);
}

double f_fm(double x) {
        mpfr_set_d(r38539, x, MPFR_RNDN);
        mpfr_mul(r38540, r38539, r38539, MPFR_RNDN);
        ;
        mpfr_mul(r38542, r38540, r38541, MPFR_RNDN);
        ;
        mpfr_sub(r38544, r38542, r38543, MPFR_RNDN);
        return mpfr_get_d(r38544, MPFR_RNDN);
}

static mpfr_t r38545, r38546, r38547, r38548, r38549, r38550;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r38545);
        mpfr_init(r38546);
        mpfr_init_set_str(r38547, "2.0", 10, MPFR_RNDN);
        mpfr_init(r38548);
        mpfr_init_set_str(r38549, "1.0", 10, MPFR_RNDN);
        mpfr_init(r38550);
}

double f_dm(double x) {
        mpfr_set_d(r38545, x, MPFR_RNDN);
        mpfr_mul(r38546, r38545, r38545, MPFR_RNDN);
        ;
        mpfr_mul(r38548, r38546, r38547, MPFR_RNDN);
        ;
        mpfr_sub(r38550, r38548, r38549, MPFR_RNDN);
        return mpfr_get_d(r38550, MPFR_RNDN);
}