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

char *name = "Numeric.Integration.TanhSinh:everywhere from integration-0.2.1";

double f_if(float x, float y) {
        float r46365 = x;
        float r46366 = 1.0;
        float r46367 = y;
        float r46368 = r46367 * r46367;
        float r46369 = r46366 + r46368;
        float r46370 = r46365 * r46369;
        return r46370;
}

double f_id(double x, double y) {
        double r46371 = x;
        double r46372 = 1.0;
        double r46373 = y;
        double r46374 = r46373 * r46373;
        double r46375 = r46372 + r46374;
        double r46376 = r46371 * r46375;
        return r46376;
}


double f_of(float x, float y) {
        float r46377 = x;
        float r46378 = 1.0;
        float r46379 = y;
        float r46380 = r46379 * r46379;
        float r46381 = r46378 + r46380;
        float r46382 = r46377 * r46381;
        return r46382;
}

double f_od(double x, double y) {
        double r46383 = x;
        double r46384 = 1.0;
        double r46385 = y;
        double r46386 = r46385 * r46385;
        double r46387 = r46384 + r46386;
        double r46388 = r46383 * r46387;
        return r46388;
}

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 r46389, r46390, r46391, r46392, r46393, r46394;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r46389);
        mpfr_init_set_str(r46390, "1.0", 10, MPFR_RNDN);
        mpfr_init(r46391);
        mpfr_init(r46392);
        mpfr_init(r46393);
        mpfr_init(r46394);
}

double f_im(double x, double y) {
        mpfr_set_d(r46389, x, MPFR_RNDN);
        ;
        mpfr_set_d(r46391, y, MPFR_RNDN);
        mpfr_mul(r46392, r46391, r46391, MPFR_RNDN);
        mpfr_add(r46393, r46390, r46392, MPFR_RNDN);
        mpfr_mul(r46394, r46389, r46393, MPFR_RNDN);
        return mpfr_get_d(r46394, MPFR_RNDN);
}

static mpfr_t r46395, r46396, r46397, r46398, r46399, r46400;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r46395);
        mpfr_init_set_str(r46396, "1.0", 10, MPFR_RNDN);
        mpfr_init(r46397);
        mpfr_init(r46398);
        mpfr_init(r46399);
        mpfr_init(r46400);
}

double f_fm(double x, double y) {
        mpfr_set_d(r46395, x, MPFR_RNDN);
        ;
        mpfr_set_d(r46397, y, MPFR_RNDN);
        mpfr_mul(r46398, r46397, r46397, MPFR_RNDN);
        mpfr_add(r46399, r46396, r46398, MPFR_RNDN);
        mpfr_mul(r46400, r46395, r46399, MPFR_RNDN);
        return mpfr_get_d(r46400, MPFR_RNDN);
}

static mpfr_t r46401, r46402, r46403, r46404, r46405, r46406;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r46401);
        mpfr_init_set_str(r46402, "1.0", 10, MPFR_RNDN);
        mpfr_init(r46403);
        mpfr_init(r46404);
        mpfr_init(r46405);
        mpfr_init(r46406);
}

double f_dm(double x, double y) {
        mpfr_set_d(r46401, x, MPFR_RNDN);
        ;
        mpfr_set_d(r46403, y, MPFR_RNDN);
        mpfr_mul(r46404, r46403, r46403, MPFR_RNDN);
        mpfr_add(r46405, r46402, r46404, MPFR_RNDN);
        mpfr_mul(r46406, r46401, r46405, MPFR_RNDN);
        return mpfr_get_d(r46406, MPFR_RNDN);
}