#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 r25442 = a;
        float r25443 = r25442 * r25442;
        float r25444 = b;
        float r25445 = r25444 * r25444;
        float r25446 = r25443 + r25445;
        float r25447 = 2;
        float r25448 = pow(r25446, r25447);
        float r25449 = 4;
        float r25450 = 1;
        float r25451 = r25450 - r25442;
        float r25452 = r25443 * r25451;
        float r25453 = 3;
        float r25454 = r25453 + r25442;
        float r25455 = r25445 * r25454;
        float r25456 = r25452 + r25455;
        float r25457 = r25449 * r25456;
        float r25458 = r25448 + r25457;
        float r25459 = r25458 - r25450;
        return r25459;
}

double f_id(double a, double b) {
        double r25460 = a;
        double r25461 = r25460 * r25460;
        double r25462 = b;
        double r25463 = r25462 * r25462;
        double r25464 = r25461 + r25463;
        double r25465 = 2;
        double r25466 = pow(r25464, r25465);
        double r25467 = 4;
        double r25468 = 1;
        double r25469 = r25468 - r25460;
        double r25470 = r25461 * r25469;
        double r25471 = 3;
        double r25472 = r25471 + r25460;
        double r25473 = r25463 * r25472;
        double r25474 = r25470 + r25473;
        double r25475 = r25467 * r25474;
        double r25476 = r25466 + r25475;
        double r25477 = r25476 - r25468;
        return r25477;
}


double f_of(float a, float b) {
        float r25478 = a;
        float r25479 = r25478 * r25478;
        float r25480 = b;
        float r25481 = r25480 * r25480;
        float r25482 = r25479 + r25481;
        float r25483 = 2;
        float r25484 = pow(r25482, r25483);
        float r25485 = 4;
        float r25486 = 1;
        float r25487 = r25486 - r25478;
        float r25488 = r25479 * r25487;
        float r25489 = 3;
        float r25490 = r25489 + r25478;
        float r25491 = r25481 * r25490;
        float r25492 = r25488 + r25491;
        float r25493 = r25485 * r25492;
        float r25494 = r25484 + r25493;
        float r25495 = r25494 - r25486;
        return r25495;
}

double f_od(double a, double b) {
        double r25496 = a;
        double r25497 = r25496 * r25496;
        double r25498 = b;
        double r25499 = r25498 * r25498;
        double r25500 = r25497 + r25499;
        double r25501 = 2;
        double r25502 = pow(r25500, r25501);
        double r25503 = 4;
        double r25504 = 1;
        double r25505 = r25504 - r25496;
        double r25506 = r25497 * r25505;
        double r25507 = 3;
        double r25508 = r25507 + r25496;
        double r25509 = r25499 * r25508;
        double r25510 = r25506 + r25509;
        double r25511 = r25503 * r25510;
        double r25512 = r25502 + r25511;
        double r25513 = r25512 - r25504;
        return r25513;
}

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 r25514, r25515, r25516, r25517, r25518, r25519, r25520, r25521, r25522, r25523, r25524, r25525, r25526, r25527, r25528, r25529, r25530, r25531;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r25514);
        mpfr_init(r25515);
        mpfr_init(r25516);
        mpfr_init(r25517);
        mpfr_init(r25518);
        mpfr_init_set_str(r25519, "2", 10, MPFR_RNDN);
        mpfr_init(r25520);
        mpfr_init_set_str(r25521, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25522, "1", 10, MPFR_RNDN);
        mpfr_init(r25523);
        mpfr_init(r25524);
        mpfr_init_set_str(r25525, "3", 10, MPFR_RNDN);
        mpfr_init(r25526);
        mpfr_init(r25527);
        mpfr_init(r25528);
        mpfr_init(r25529);
        mpfr_init(r25530);
        mpfr_init(r25531);
}

double f_im(double a, double b) {
        mpfr_set_d(r25514, a, MPFR_RNDN);
        mpfr_mul(r25515, r25514, r25514, MPFR_RNDN);
        mpfr_set_d(r25516, b, MPFR_RNDN);
        mpfr_mul(r25517, r25516, r25516, MPFR_RNDN);
        mpfr_add(r25518, r25515, r25517, MPFR_RNDN);
        ;
        mpfr_pow(r25520, r25518, r25519, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25523, r25522, r25514, MPFR_RNDN);
        mpfr_mul(r25524, r25515, r25523, MPFR_RNDN);
        ;
        mpfr_add(r25526, r25525, r25514, MPFR_RNDN);
        mpfr_mul(r25527, r25517, r25526, MPFR_RNDN);
        mpfr_add(r25528, r25524, r25527, MPFR_RNDN);
        mpfr_mul(r25529, r25521, r25528, MPFR_RNDN);
        mpfr_add(r25530, r25520, r25529, MPFR_RNDN);
        mpfr_sub(r25531, r25530, r25522, MPFR_RNDN);
        return mpfr_get_d(r25531, MPFR_RNDN);
}

static mpfr_t r25532, r25533, r25534, r25535, r25536, r25537, r25538, r25539, r25540, r25541, r25542, r25543, r25544, r25545, r25546, r25547, r25548, r25549;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r25532);
        mpfr_init(r25533);
        mpfr_init(r25534);
        mpfr_init(r25535);
        mpfr_init(r25536);
        mpfr_init_set_str(r25537, "2", 10, MPFR_RNDN);
        mpfr_init(r25538);
        mpfr_init_set_str(r25539, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25540, "1", 10, MPFR_RNDN);
        mpfr_init(r25541);
        mpfr_init(r25542);
        mpfr_init_set_str(r25543, "3", 10, MPFR_RNDN);
        mpfr_init(r25544);
        mpfr_init(r25545);
        mpfr_init(r25546);
        mpfr_init(r25547);
        mpfr_init(r25548);
        mpfr_init(r25549);
}

double f_fm(double a, double b) {
        mpfr_set_d(r25532, a, MPFR_RNDN);
        mpfr_mul(r25533, r25532, r25532, MPFR_RNDN);
        mpfr_set_d(r25534, b, MPFR_RNDN);
        mpfr_mul(r25535, r25534, r25534, MPFR_RNDN);
        mpfr_add(r25536, r25533, r25535, MPFR_RNDN);
        ;
        mpfr_pow(r25538, r25536, r25537, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25541, r25540, r25532, MPFR_RNDN);
        mpfr_mul(r25542, r25533, r25541, MPFR_RNDN);
        ;
        mpfr_add(r25544, r25543, r25532, MPFR_RNDN);
        mpfr_mul(r25545, r25535, r25544, MPFR_RNDN);
        mpfr_add(r25546, r25542, r25545, MPFR_RNDN);
        mpfr_mul(r25547, r25539, r25546, MPFR_RNDN);
        mpfr_add(r25548, r25538, r25547, MPFR_RNDN);
        mpfr_sub(r25549, r25548, r25540, MPFR_RNDN);
        return mpfr_get_d(r25549, MPFR_RNDN);
}

static mpfr_t r25550, r25551, r25552, r25553, r25554, r25555, r25556, r25557, r25558, r25559, r25560, r25561, r25562, r25563, r25564, r25565, r25566, r25567;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r25550);
        mpfr_init(r25551);
        mpfr_init(r25552);
        mpfr_init(r25553);
        mpfr_init(r25554);
        mpfr_init_set_str(r25555, "2", 10, MPFR_RNDN);
        mpfr_init(r25556);
        mpfr_init_set_str(r25557, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25558, "1", 10, MPFR_RNDN);
        mpfr_init(r25559);
        mpfr_init(r25560);
        mpfr_init_set_str(r25561, "3", 10, MPFR_RNDN);
        mpfr_init(r25562);
        mpfr_init(r25563);
        mpfr_init(r25564);
        mpfr_init(r25565);
        mpfr_init(r25566);
        mpfr_init(r25567);
}

double f_dm(double a, double b) {
        mpfr_set_d(r25550, a, MPFR_RNDN);
        mpfr_mul(r25551, r25550, r25550, MPFR_RNDN);
        mpfr_set_d(r25552, b, MPFR_RNDN);
        mpfr_mul(r25553, r25552, r25552, MPFR_RNDN);
        mpfr_add(r25554, r25551, r25553, MPFR_RNDN);
        ;
        mpfr_pow(r25556, r25554, r25555, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25559, r25558, r25550, MPFR_RNDN);
        mpfr_mul(r25560, r25551, r25559, MPFR_RNDN);
        ;
        mpfr_add(r25562, r25561, r25550, MPFR_RNDN);
        mpfr_mul(r25563, r25553, r25562, MPFR_RNDN);
        mpfr_add(r25564, r25560, r25563, MPFR_RNDN);
        mpfr_mul(r25565, r25557, r25564, MPFR_RNDN);
        mpfr_add(r25566, r25556, r25565, MPFR_RNDN);
        mpfr_sub(r25567, r25566, r25558, MPFR_RNDN);
        return mpfr_get_d(r25567, MPFR_RNDN);
}