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

char *name = "Bouland and Aaronson, Equation (24)";

double f_if(float a, float b) {
        float r25464 = a;
        float r25465 = r25464 * r25464;
        float r25466 = b;
        float r25467 = r25466 * r25466;
        float r25468 = r25465 + r25467;
        float r25469 = 2;
        float r25470 = pow(r25468, r25469);
        float r25471 = 4;
        float r25472 = 1;
        float r25473 = r25472 - r25464;
        float r25474 = r25465 * r25473;
        float r25475 = 3;
        float r25476 = r25475 + r25464;
        float r25477 = r25467 * r25476;
        float r25478 = r25474 + r25477;
        float r25479 = r25471 * r25478;
        float r25480 = r25470 + r25479;
        float r25481 = r25480 - r25472;
        return r25481;
}

double f_id(double a, double b) {
        double r25482 = a;
        double r25483 = r25482 * r25482;
        double r25484 = b;
        double r25485 = r25484 * r25484;
        double r25486 = r25483 + r25485;
        double r25487 = 2;
        double r25488 = pow(r25486, r25487);
        double r25489 = 4;
        double r25490 = 1;
        double r25491 = r25490 - r25482;
        double r25492 = r25483 * r25491;
        double r25493 = 3;
        double r25494 = r25493 + r25482;
        double r25495 = r25485 * r25494;
        double r25496 = r25492 + r25495;
        double r25497 = r25489 * r25496;
        double r25498 = r25488 + r25497;
        double r25499 = r25498 - r25490;
        return r25499;
}


double f_of(float a, float b) {
        float r25500 = b;
        float r25501 = r25500 * r25500;
        float r25502 = a;
        float r25503 = r25502 * r25502;
        float r25504 = r25501 + r25503;
        float r25505 = cbrt(r25504);
        float r25506 = r25504 * r25505;
        float r25507 = r25503 + r25501;
        float r25508 = cbrt(r25507);
        float r25509 = 2;
        float r25510 = pow(r25508, r25509);
        float r25511 = r25506 * r25510;
        float r25512 = 4;
        float r25513 = 1;
        float r25514 = r25513 - r25502;
        float r25515 = r25503 * r25514;
        float r25516 = 3;
        float r25517 = r25516 + r25502;
        float r25518 = r25501 * r25517;
        float r25519 = r25515 + r25518;
        float r25520 = r25512 * r25519;
        float r25521 = r25511 + r25520;
        float r25522 = r25521 - r25513;
        return r25522;
}

double f_od(double a, double b) {
        double r25523 = b;
        double r25524 = r25523 * r25523;
        double r25525 = a;
        double r25526 = r25525 * r25525;
        double r25527 = r25524 + r25526;
        double r25528 = cbrt(r25527);
        double r25529 = r25527 * r25528;
        double r25530 = r25526 + r25524;
        double r25531 = cbrt(r25530);
        double r25532 = 2;
        double r25533 = pow(r25531, r25532);
        double r25534 = r25529 * r25533;
        double r25535 = 4;
        double r25536 = 1;
        double r25537 = r25536 - r25525;
        double r25538 = r25526 * r25537;
        double r25539 = 3;
        double r25540 = r25539 + r25525;
        double r25541 = r25524 * r25540;
        double r25542 = r25538 + r25541;
        double r25543 = r25535 * r25542;
        double r25544 = r25534 + r25543;
        double r25545 = r25544 - r25536;
        return r25545;
}

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 r25546, r25547, r25548, r25549, r25550, r25551, r25552, r25553, r25554, r25555, r25556, r25557, r25558, r25559, r25560, r25561, r25562, r25563;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(4112);
        mpfr_init(r25546);
        mpfr_init(r25547);
        mpfr_init(r25548);
        mpfr_init(r25549);
        mpfr_init(r25550);
        mpfr_init_set_str(r25551, "2", 10, MPFR_RNDN);
        mpfr_init(r25552);
        mpfr_init_set_str(r25553, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25554, "1", 10, MPFR_RNDN);
        mpfr_init(r25555);
        mpfr_init(r25556);
        mpfr_init_set_str(r25557, "3", 10, MPFR_RNDN);
        mpfr_init(r25558);
        mpfr_init(r25559);
        mpfr_init(r25560);
        mpfr_init(r25561);
        mpfr_init(r25562);
        mpfr_init(r25563);
}

double f_im(double a, double b) {
        mpfr_set_d(r25546, a, MPFR_RNDN);
        mpfr_mul(r25547, r25546, r25546, MPFR_RNDN);
        mpfr_set_d(r25548, b, MPFR_RNDN);
        mpfr_mul(r25549, r25548, r25548, MPFR_RNDN);
        mpfr_add(r25550, r25547, r25549, MPFR_RNDN);
        ;
        mpfr_pow(r25552, r25550, r25551, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25555, r25554, r25546, MPFR_RNDN);
        mpfr_mul(r25556, r25547, r25555, MPFR_RNDN);
        ;
        mpfr_add(r25558, r25557, r25546, MPFR_RNDN);
        mpfr_mul(r25559, r25549, r25558, MPFR_RNDN);
        mpfr_add(r25560, r25556, r25559, MPFR_RNDN);
        mpfr_mul(r25561, r25553, r25560, MPFR_RNDN);
        mpfr_add(r25562, r25552, r25561, MPFR_RNDN);
        mpfr_sub(r25563, r25562, r25554, MPFR_RNDN);
        return mpfr_get_d(r25563, MPFR_RNDN);
}

static mpfr_t r25564, r25565, r25566, r25567, r25568, r25569, r25570, r25571, r25572, r25573, r25574, r25575, r25576, r25577, r25578, r25579, r25580, r25581, r25582, r25583, r25584, r25585, r25586;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(4112);
        mpfr_init(r25564);
        mpfr_init(r25565);
        mpfr_init(r25566);
        mpfr_init(r25567);
        mpfr_init(r25568);
        mpfr_init(r25569);
        mpfr_init(r25570);
        mpfr_init(r25571);
        mpfr_init(r25572);
        mpfr_init_set_str(r25573, "2", 10, MPFR_RNDN);
        mpfr_init(r25574);
        mpfr_init(r25575);
        mpfr_init_set_str(r25576, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25577, "1", 10, MPFR_RNDN);
        mpfr_init(r25578);
        mpfr_init(r25579);
        mpfr_init_set_str(r25580, "3", 10, MPFR_RNDN);
        mpfr_init(r25581);
        mpfr_init(r25582);
        mpfr_init(r25583);
        mpfr_init(r25584);
        mpfr_init(r25585);
        mpfr_init(r25586);
}

double f_fm(double a, double b) {
        mpfr_set_d(r25564, b, MPFR_RNDN);
        mpfr_mul(r25565, r25564, r25564, MPFR_RNDN);
        mpfr_set_d(r25566, a, MPFR_RNDN);
        mpfr_mul(r25567, r25566, r25566, MPFR_RNDN);
        mpfr_add(r25568, r25565, r25567, MPFR_RNDN);
        mpfr_cbrt(r25569, r25568, MPFR_RNDN);
        mpfr_mul(r25570, r25568, r25569, MPFR_RNDN);
        mpfr_add(r25571, r25567, r25565, MPFR_RNDN);
        mpfr_cbrt(r25572, r25571, MPFR_RNDN);
        ;
        mpfr_pow(r25574, r25572, r25573, MPFR_RNDN);
        mpfr_mul(r25575, r25570, r25574, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25578, r25577, r25566, MPFR_RNDN);
        mpfr_mul(r25579, r25567, r25578, MPFR_RNDN);
        ;
        mpfr_add(r25581, r25580, r25566, MPFR_RNDN);
        mpfr_mul(r25582, r25565, r25581, MPFR_RNDN);
        mpfr_add(r25583, r25579, r25582, MPFR_RNDN);
        mpfr_mul(r25584, r25576, r25583, MPFR_RNDN);
        mpfr_add(r25585, r25575, r25584, MPFR_RNDN);
        mpfr_sub(r25586, r25585, r25577, MPFR_RNDN);
        return mpfr_get_d(r25586, MPFR_RNDN);
}

static mpfr_t r25587, r25588, r25589, r25590, r25591, r25592, r25593, r25594, r25595, r25596, r25597, r25598, r25599, r25600, r25601, r25602, r25603, r25604, r25605, r25606, r25607, r25608, r25609;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(4112);
        mpfr_init(r25587);
        mpfr_init(r25588);
        mpfr_init(r25589);
        mpfr_init(r25590);
        mpfr_init(r25591);
        mpfr_init(r25592);
        mpfr_init(r25593);
        mpfr_init(r25594);
        mpfr_init(r25595);
        mpfr_init_set_str(r25596, "2", 10, MPFR_RNDN);
        mpfr_init(r25597);
        mpfr_init(r25598);
        mpfr_init_set_str(r25599, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25600, "1", 10, MPFR_RNDN);
        mpfr_init(r25601);
        mpfr_init(r25602);
        mpfr_init_set_str(r25603, "3", 10, MPFR_RNDN);
        mpfr_init(r25604);
        mpfr_init(r25605);
        mpfr_init(r25606);
        mpfr_init(r25607);
        mpfr_init(r25608);
        mpfr_init(r25609);
}

double f_dm(double a, double b) {
        mpfr_set_d(r25587, b, MPFR_RNDN);
        mpfr_mul(r25588, r25587, r25587, MPFR_RNDN);
        mpfr_set_d(r25589, a, MPFR_RNDN);
        mpfr_mul(r25590, r25589, r25589, MPFR_RNDN);
        mpfr_add(r25591, r25588, r25590, MPFR_RNDN);
        mpfr_cbrt(r25592, r25591, MPFR_RNDN);
        mpfr_mul(r25593, r25591, r25592, MPFR_RNDN);
        mpfr_add(r25594, r25590, r25588, MPFR_RNDN);
        mpfr_cbrt(r25595, r25594, MPFR_RNDN);
        ;
        mpfr_pow(r25597, r25595, r25596, MPFR_RNDN);
        mpfr_mul(r25598, r25593, r25597, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25601, r25600, r25589, MPFR_RNDN);
        mpfr_mul(r25602, r25590, r25601, MPFR_RNDN);
        ;
        mpfr_add(r25604, r25603, r25589, MPFR_RNDN);
        mpfr_mul(r25605, r25588, r25604, MPFR_RNDN);
        mpfr_add(r25606, r25602, r25605, MPFR_RNDN);
        mpfr_mul(r25607, r25599, r25606, MPFR_RNDN);
        mpfr_add(r25608, r25598, r25607, MPFR_RNDN);
        mpfr_sub(r25609, r25608, r25600, MPFR_RNDN);
        return mpfr_get_d(r25609, MPFR_RNDN);
}