#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 r25401 = a;
        float r25402 = r25401 * r25401;
        float r25403 = b;
        float r25404 = r25403 * r25403;
        float r25405 = r25402 + r25404;
        float r25406 = 2;
        float r25407 = pow(r25405, r25406);
        float r25408 = 4;
        float r25409 = 1;
        float r25410 = r25409 - r25401;
        float r25411 = r25402 * r25410;
        float r25412 = 3;
        float r25413 = r25412 + r25401;
        float r25414 = r25404 * r25413;
        float r25415 = r25411 + r25414;
        float r25416 = r25408 * r25415;
        float r25417 = r25407 + r25416;
        float r25418 = r25417 - r25409;
        return r25418;
}

double f_id(double a, double b) {
        double r25419 = a;
        double r25420 = r25419 * r25419;
        double r25421 = b;
        double r25422 = r25421 * r25421;
        double r25423 = r25420 + r25422;
        double r25424 = 2;
        double r25425 = pow(r25423, r25424);
        double r25426 = 4;
        double r25427 = 1;
        double r25428 = r25427 - r25419;
        double r25429 = r25420 * r25428;
        double r25430 = 3;
        double r25431 = r25430 + r25419;
        double r25432 = r25422 * r25431;
        double r25433 = r25429 + r25432;
        double r25434 = r25426 * r25433;
        double r25435 = r25425 + r25434;
        double r25436 = r25435 - r25427;
        return r25436;
}


double f_of(float a, float b) {
        float r25437 = a;
        float r25438 = r25437 * r25437;
        float r25439 = b;
        float r25440 = r25439 * r25439;
        float r25441 = r25438 + r25440;
        float r25442 = 2;
        float r25443 = pow(r25441, r25442);
        float r25444 = 4;
        float r25445 = 1;
        float r25446 = r25445 - r25437;
        float r25447 = r25438 * r25446;
        float r25448 = 3;
        float r25449 = r25448 + r25437;
        float r25450 = r25440 * r25449;
        float r25451 = r25447 + r25450;
        float r25452 = r25444 * r25451;
        float r25453 = r25443 + r25452;
        float r25454 = r25453 - r25445;
        return r25454;
}

double f_od(double a, double b) {
        double r25455 = a;
        double r25456 = r25455 * r25455;
        double r25457 = b;
        double r25458 = r25457 * r25457;
        double r25459 = r25456 + r25458;
        double r25460 = 2;
        double r25461 = pow(r25459, r25460);
        double r25462 = 4;
        double r25463 = 1;
        double r25464 = r25463 - r25455;
        double r25465 = r25456 * r25464;
        double r25466 = 3;
        double r25467 = r25466 + r25455;
        double r25468 = r25458 * r25467;
        double r25469 = r25465 + r25468;
        double r25470 = r25462 * r25469;
        double r25471 = r25461 + r25470;
        double r25472 = r25471 - r25463;
        return r25472;
}

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 r25473, r25474, r25475, r25476, r25477, r25478, r25479, r25480, r25481, r25482, r25483, r25484, r25485, r25486, r25487, r25488, r25489, r25490;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r25473);
        mpfr_init(r25474);
        mpfr_init(r25475);
        mpfr_init(r25476);
        mpfr_init(r25477);
        mpfr_init_set_str(r25478, "2", 10, MPFR_RNDN);
        mpfr_init(r25479);
        mpfr_init_set_str(r25480, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25481, "1", 10, MPFR_RNDN);
        mpfr_init(r25482);
        mpfr_init(r25483);
        mpfr_init_set_str(r25484, "3", 10, MPFR_RNDN);
        mpfr_init(r25485);
        mpfr_init(r25486);
        mpfr_init(r25487);
        mpfr_init(r25488);
        mpfr_init(r25489);
        mpfr_init(r25490);
}

double f_im(double a, double b) {
        mpfr_set_d(r25473, a, MPFR_RNDN);
        mpfr_mul(r25474, r25473, r25473, MPFR_RNDN);
        mpfr_set_d(r25475, b, MPFR_RNDN);
        mpfr_mul(r25476, r25475, r25475, MPFR_RNDN);
        mpfr_add(r25477, r25474, r25476, MPFR_RNDN);
        ;
        mpfr_pow(r25479, r25477, r25478, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25482, r25481, r25473, MPFR_RNDN);
        mpfr_mul(r25483, r25474, r25482, MPFR_RNDN);
        ;
        mpfr_add(r25485, r25484, r25473, MPFR_RNDN);
        mpfr_mul(r25486, r25476, r25485, MPFR_RNDN);
        mpfr_add(r25487, r25483, r25486, MPFR_RNDN);
        mpfr_mul(r25488, r25480, r25487, MPFR_RNDN);
        mpfr_add(r25489, r25479, r25488, MPFR_RNDN);
        mpfr_sub(r25490, r25489, r25481, MPFR_RNDN);
        return mpfr_get_d(r25490, MPFR_RNDN);
}

static mpfr_t r25491, r25492, r25493, r25494, r25495, r25496, r25497, r25498, r25499, r25500, r25501, r25502, r25503, r25504, r25505, r25506, r25507, r25508;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25491);
        mpfr_init(r25492);
        mpfr_init(r25493);
        mpfr_init(r25494);
        mpfr_init(r25495);
        mpfr_init_set_str(r25496, "2", 10, MPFR_RNDN);
        mpfr_init(r25497);
        mpfr_init_set_str(r25498, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25499, "1", 10, MPFR_RNDN);
        mpfr_init(r25500);
        mpfr_init(r25501);
        mpfr_init_set_str(r25502, "3", 10, MPFR_RNDN);
        mpfr_init(r25503);
        mpfr_init(r25504);
        mpfr_init(r25505);
        mpfr_init(r25506);
        mpfr_init(r25507);
        mpfr_init(r25508);
}

double f_fm(double a, double b) {
        mpfr_set_d(r25491, a, MPFR_RNDN);
        mpfr_mul(r25492, r25491, r25491, MPFR_RNDN);
        mpfr_set_d(r25493, b, MPFR_RNDN);
        mpfr_mul(r25494, r25493, r25493, MPFR_RNDN);
        mpfr_add(r25495, r25492, r25494, MPFR_RNDN);
        ;
        mpfr_pow(r25497, r25495, r25496, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25500, r25499, r25491, MPFR_RNDN);
        mpfr_mul(r25501, r25492, r25500, MPFR_RNDN);
        ;
        mpfr_add(r25503, r25502, r25491, MPFR_RNDN);
        mpfr_mul(r25504, r25494, r25503, MPFR_RNDN);
        mpfr_add(r25505, r25501, r25504, MPFR_RNDN);
        mpfr_mul(r25506, r25498, r25505, MPFR_RNDN);
        mpfr_add(r25507, r25497, r25506, MPFR_RNDN);
        mpfr_sub(r25508, r25507, r25499, MPFR_RNDN);
        return mpfr_get_d(r25508, MPFR_RNDN);
}

static mpfr_t r25509, r25510, r25511, r25512, r25513, r25514, r25515, r25516, r25517, r25518, r25519, r25520, r25521, r25522, r25523, r25524, r25525, r25526;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25509);
        mpfr_init(r25510);
        mpfr_init(r25511);
        mpfr_init(r25512);
        mpfr_init(r25513);
        mpfr_init_set_str(r25514, "2", 10, MPFR_RNDN);
        mpfr_init(r25515);
        mpfr_init_set_str(r25516, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25517, "1", 10, MPFR_RNDN);
        mpfr_init(r25518);
        mpfr_init(r25519);
        mpfr_init_set_str(r25520, "3", 10, MPFR_RNDN);
        mpfr_init(r25521);
        mpfr_init(r25522);
        mpfr_init(r25523);
        mpfr_init(r25524);
        mpfr_init(r25525);
        mpfr_init(r25526);
}

double f_dm(double a, double b) {
        mpfr_set_d(r25509, a, MPFR_RNDN);
        mpfr_mul(r25510, r25509, r25509, MPFR_RNDN);
        mpfr_set_d(r25511, b, MPFR_RNDN);
        mpfr_mul(r25512, r25511, r25511, MPFR_RNDN);
        mpfr_add(r25513, r25510, r25512, MPFR_RNDN);
        ;
        mpfr_pow(r25515, r25513, r25514, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25518, r25517, r25509, MPFR_RNDN);
        mpfr_mul(r25519, r25510, r25518, MPFR_RNDN);
        ;
        mpfr_add(r25521, r25520, r25509, MPFR_RNDN);
        mpfr_mul(r25522, r25512, r25521, MPFR_RNDN);
        mpfr_add(r25523, r25519, r25522, MPFR_RNDN);
        mpfr_mul(r25524, r25516, r25523, MPFR_RNDN);
        mpfr_add(r25525, r25515, r25524, MPFR_RNDN);
        mpfr_sub(r25526, r25525, r25517, MPFR_RNDN);
        return mpfr_get_d(r25526, MPFR_RNDN);
}