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

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

double f_if(float a, float b) {
        float r27405 = a;
        float r27406 = r27405 * r27405;
        float r27407 = b;
        float r27408 = r27407 * r27407;
        float r27409 = r27406 + r27408;
        float r27410 = 2;
        float r27411 = pow(r27409, r27410);
        float r27412 = 4;
        float r27413 = 1;
        float r27414 = r27413 + r27405;
        float r27415 = r27406 * r27414;
        float r27416 = 3;
        float r27417 = r27416 * r27405;
        float r27418 = r27413 - r27417;
        float r27419 = r27408 * r27418;
        float r27420 = r27415 + r27419;
        float r27421 = r27412 * r27420;
        float r27422 = r27411 + r27421;
        float r27423 = r27422 - r27413;
        return r27423;
}

double f_id(double a, double b) {
        double r27424 = a;
        double r27425 = r27424 * r27424;
        double r27426 = b;
        double r27427 = r27426 * r27426;
        double r27428 = r27425 + r27427;
        double r27429 = 2;
        double r27430 = pow(r27428, r27429);
        double r27431 = 4;
        double r27432 = 1;
        double r27433 = r27432 + r27424;
        double r27434 = r27425 * r27433;
        double r27435 = 3;
        double r27436 = r27435 * r27424;
        double r27437 = r27432 - r27436;
        double r27438 = r27427 * r27437;
        double r27439 = r27434 + r27438;
        double r27440 = r27431 * r27439;
        double r27441 = r27430 + r27440;
        double r27442 = r27441 - r27432;
        return r27442;
}


double f_of(float a, float b) {
        float r27443 = a;
        float r27444 = r27443 * r27443;
        float r27445 = b;
        float r27446 = r27445 * r27445;
        float r27447 = r27444 + r27446;
        float r27448 = 2;
        float r27449 = pow(r27447, r27448);
        float r27450 = 4;
        float r27451 = 1;
        float r27452 = r27451 + r27443;
        float r27453 = r27444 * r27452;
        float r27454 = 3;
        float r27455 = r27454 * r27443;
        float r27456 = r27451 - r27455;
        float r27457 = r27446 * r27456;
        float r27458 = r27453 + r27457;
        float r27459 = r27450 * r27458;
        float r27460 = r27449 + r27459;
        float r27461 = sqrt(r27460);
        float r27462 = r27461 * r27461;
        float r27463 = r27462 - r27451;
        return r27463;
}

double f_od(double a, double b) {
        double r27464 = a;
        double r27465 = r27464 * r27464;
        double r27466 = b;
        double r27467 = r27466 * r27466;
        double r27468 = r27465 + r27467;
        double r27469 = 2;
        double r27470 = pow(r27468, r27469);
        double r27471 = 4;
        double r27472 = 1;
        double r27473 = r27472 + r27464;
        double r27474 = r27465 * r27473;
        double r27475 = 3;
        double r27476 = r27475 * r27464;
        double r27477 = r27472 - r27476;
        double r27478 = r27467 * r27477;
        double r27479 = r27474 + r27478;
        double r27480 = r27471 * r27479;
        double r27481 = r27470 + r27480;
        double r27482 = sqrt(r27481);
        double r27483 = r27482 * r27482;
        double r27484 = r27483 - r27472;
        return r27484;
}

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 r27485, r27486, r27487, r27488, r27489, r27490, r27491, r27492, r27493, r27494, r27495, r27496, r27497, r27498, r27499, r27500, r27501, r27502, r27503;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r27485);
        mpfr_init(r27486);
        mpfr_init(r27487);
        mpfr_init(r27488);
        mpfr_init(r27489);
        mpfr_init_set_str(r27490, "2", 10, MPFR_RNDN);
        mpfr_init(r27491);
        mpfr_init_set_str(r27492, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27493, "1", 10, MPFR_RNDN);
        mpfr_init(r27494);
        mpfr_init(r27495);
        mpfr_init_set_str(r27496, "3", 10, MPFR_RNDN);
        mpfr_init(r27497);
        mpfr_init(r27498);
        mpfr_init(r27499);
        mpfr_init(r27500);
        mpfr_init(r27501);
        mpfr_init(r27502);
        mpfr_init(r27503);
}

double f_im(double a, double b) {
        mpfr_set_d(r27485, a, MPFR_RNDN);
        mpfr_mul(r27486, r27485, r27485, MPFR_RNDN);
        mpfr_set_d(r27487, b, MPFR_RNDN);
        mpfr_mul(r27488, r27487, r27487, MPFR_RNDN);
        mpfr_add(r27489, r27486, r27488, MPFR_RNDN);
        ;
        mpfr_pow(r27491, r27489, r27490, MPFR_RNDN);
        ;
        ;
        mpfr_add(r27494, r27493, r27485, MPFR_RNDN);
        mpfr_mul(r27495, r27486, r27494, MPFR_RNDN);
        ;
        mpfr_mul(r27497, r27496, r27485, MPFR_RNDN);
        mpfr_sub(r27498, r27493, r27497, MPFR_RNDN);
        mpfr_mul(r27499, r27488, r27498, MPFR_RNDN);
        mpfr_add(r27500, r27495, r27499, MPFR_RNDN);
        mpfr_mul(r27501, r27492, r27500, MPFR_RNDN);
        mpfr_add(r27502, r27491, r27501, MPFR_RNDN);
        mpfr_sub(r27503, r27502, r27493, MPFR_RNDN);
        return mpfr_get_d(r27503, MPFR_RNDN);
}

static mpfr_t r27504, r27505, r27506, r27507, r27508, r27509, r27510, r27511, r27512, r27513, r27514, r27515, r27516, r27517, r27518, r27519, r27520, r27521, r27522, r27523, r27524;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27504);
        mpfr_init(r27505);
        mpfr_init(r27506);
        mpfr_init(r27507);
        mpfr_init(r27508);
        mpfr_init_set_str(r27509, "2", 10, MPFR_RNDN);
        mpfr_init(r27510);
        mpfr_init_set_str(r27511, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27512, "1", 10, MPFR_RNDN);
        mpfr_init(r27513);
        mpfr_init(r27514);
        mpfr_init_set_str(r27515, "3", 10, MPFR_RNDN);
        mpfr_init(r27516);
        mpfr_init(r27517);
        mpfr_init(r27518);
        mpfr_init(r27519);
        mpfr_init(r27520);
        mpfr_init(r27521);
        mpfr_init(r27522);
        mpfr_init(r27523);
        mpfr_init(r27524);
}

double f_fm(double a, double b) {
        mpfr_set_d(r27504, a, MPFR_RNDN);
        mpfr_mul(r27505, r27504, r27504, MPFR_RNDN);
        mpfr_set_d(r27506, b, MPFR_RNDN);
        mpfr_mul(r27507, r27506, r27506, MPFR_RNDN);
        mpfr_add(r27508, r27505, r27507, MPFR_RNDN);
        ;
        mpfr_pow(r27510, r27508, r27509, MPFR_RNDN);
        ;
        ;
        mpfr_add(r27513, r27512, r27504, MPFR_RNDN);
        mpfr_mul(r27514, r27505, r27513, MPFR_RNDN);
        ;
        mpfr_mul(r27516, r27515, r27504, MPFR_RNDN);
        mpfr_sub(r27517, r27512, r27516, MPFR_RNDN);
        mpfr_mul(r27518, r27507, r27517, MPFR_RNDN);
        mpfr_add(r27519, r27514, r27518, MPFR_RNDN);
        mpfr_mul(r27520, r27511, r27519, MPFR_RNDN);
        mpfr_add(r27521, r27510, r27520, MPFR_RNDN);
        mpfr_sqrt(r27522, r27521, MPFR_RNDN);
        mpfr_mul(r27523, r27522, r27522, MPFR_RNDN);
        mpfr_sub(r27524, r27523, r27512, MPFR_RNDN);
        return mpfr_get_d(r27524, MPFR_RNDN);
}

static mpfr_t r27525, r27526, r27527, r27528, r27529, r27530, r27531, r27532, r27533, r27534, r27535, r27536, r27537, r27538, r27539, r27540, r27541, r27542, r27543, r27544, r27545;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27525);
        mpfr_init(r27526);
        mpfr_init(r27527);
        mpfr_init(r27528);
        mpfr_init(r27529);
        mpfr_init_set_str(r27530, "2", 10, MPFR_RNDN);
        mpfr_init(r27531);
        mpfr_init_set_str(r27532, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27533, "1", 10, MPFR_RNDN);
        mpfr_init(r27534);
        mpfr_init(r27535);
        mpfr_init_set_str(r27536, "3", 10, MPFR_RNDN);
        mpfr_init(r27537);
        mpfr_init(r27538);
        mpfr_init(r27539);
        mpfr_init(r27540);
        mpfr_init(r27541);
        mpfr_init(r27542);
        mpfr_init(r27543);
        mpfr_init(r27544);
        mpfr_init(r27545);
}

double f_dm(double a, double b) {
        mpfr_set_d(r27525, a, MPFR_RNDN);
        mpfr_mul(r27526, r27525, r27525, MPFR_RNDN);
        mpfr_set_d(r27527, b, MPFR_RNDN);
        mpfr_mul(r27528, r27527, r27527, MPFR_RNDN);
        mpfr_add(r27529, r27526, r27528, MPFR_RNDN);
        ;
        mpfr_pow(r27531, r27529, r27530, MPFR_RNDN);
        ;
        ;
        mpfr_add(r27534, r27533, r27525, MPFR_RNDN);
        mpfr_mul(r27535, r27526, r27534, MPFR_RNDN);
        ;
        mpfr_mul(r27537, r27536, r27525, MPFR_RNDN);
        mpfr_sub(r27538, r27533, r27537, MPFR_RNDN);
        mpfr_mul(r27539, r27528, r27538, MPFR_RNDN);
        mpfr_add(r27540, r27535, r27539, MPFR_RNDN);
        mpfr_mul(r27541, r27532, r27540, MPFR_RNDN);
        mpfr_add(r27542, r27531, r27541, MPFR_RNDN);
        mpfr_sqrt(r27543, r27542, MPFR_RNDN);
        mpfr_mul(r27544, r27543, r27543, MPFR_RNDN);
        mpfr_sub(r27545, r27544, r27533, MPFR_RNDN);
        return mpfr_get_d(r27545, MPFR_RNDN);
}