#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 r27302 = a;
        float r27303 = r27302 * r27302;
        float r27304 = b;
        float r27305 = r27304 * r27304;
        float r27306 = r27303 + r27305;
        float r27307 = 2;
        float r27308 = pow(r27306, r27307);
        float r27309 = 4;
        float r27310 = 1;
        float r27311 = r27310 + r27302;
        float r27312 = r27303 * r27311;
        float r27313 = 3;
        float r27314 = r27313 * r27302;
        float r27315 = r27310 - r27314;
        float r27316 = r27305 * r27315;
        float r27317 = r27312 + r27316;
        float r27318 = r27309 * r27317;
        float r27319 = r27308 + r27318;
        float r27320 = r27319 - r27310;
        return r27320;
}

double f_id(double a, double b) {
        double r27321 = a;
        double r27322 = r27321 * r27321;
        double r27323 = b;
        double r27324 = r27323 * r27323;
        double r27325 = r27322 + r27324;
        double r27326 = 2;
        double r27327 = pow(r27325, r27326);
        double r27328 = 4;
        double r27329 = 1;
        double r27330 = r27329 + r27321;
        double r27331 = r27322 * r27330;
        double r27332 = 3;
        double r27333 = r27332 * r27321;
        double r27334 = r27329 - r27333;
        double r27335 = r27324 * r27334;
        double r27336 = r27331 + r27335;
        double r27337 = r27328 * r27336;
        double r27338 = r27327 + r27337;
        double r27339 = r27338 - r27329;
        return r27339;
}


double f_of(float a, float b) {
        float r27340 = a;
        float r27341 = r27340 * r27340;
        float r27342 = b;
        float r27343 = r27342 * r27342;
        float r27344 = r27341 + r27343;
        float r27345 = 2;
        float r27346 = pow(r27344, r27345);
        float r27347 = 4;
        float r27348 = 1;
        float r27349 = r27348 + r27340;
        float r27350 = r27341 * r27349;
        float r27351 = 3;
        float r27352 = r27351 * r27340;
        float r27353 = r27348 - r27352;
        float r27354 = r27343 * r27353;
        float r27355 = r27350 + r27354;
        float r27356 = r27347 * r27355;
        float r27357 = r27346 + r27356;
        float r27358 = r27357 - r27348;
        return r27358;
}

double f_od(double a, double b) {
        double r27359 = a;
        double r27360 = r27359 * r27359;
        double r27361 = b;
        double r27362 = r27361 * r27361;
        double r27363 = r27360 + r27362;
        double r27364 = 2;
        double r27365 = pow(r27363, r27364);
        double r27366 = 4;
        double r27367 = 1;
        double r27368 = r27367 + r27359;
        double r27369 = r27360 * r27368;
        double r27370 = 3;
        double r27371 = r27370 * r27359;
        double r27372 = r27367 - r27371;
        double r27373 = r27362 * r27372;
        double r27374 = r27369 + r27373;
        double r27375 = r27366 * r27374;
        double r27376 = r27365 + r27375;
        double r27377 = r27376 - r27367;
        return r27377;
}

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 r27378, r27379, r27380, r27381, r27382, r27383, r27384, r27385, r27386, r27387, r27388, r27389, r27390, r27391, r27392, r27393, r27394, r27395, r27396;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r27378);
        mpfr_init(r27379);
        mpfr_init(r27380);
        mpfr_init(r27381);
        mpfr_init(r27382);
        mpfr_init_set_str(r27383, "2", 10, MPFR_RNDN);
        mpfr_init(r27384);
        mpfr_init_set_str(r27385, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27386, "1", 10, MPFR_RNDN);
        mpfr_init(r27387);
        mpfr_init(r27388);
        mpfr_init_set_str(r27389, "3", 10, MPFR_RNDN);
        mpfr_init(r27390);
        mpfr_init(r27391);
        mpfr_init(r27392);
        mpfr_init(r27393);
        mpfr_init(r27394);
        mpfr_init(r27395);
        mpfr_init(r27396);
}

double f_im(double a, double b) {
        mpfr_set_d(r27378, a, MPFR_RNDN);
        mpfr_mul(r27379, r27378, r27378, MPFR_RNDN);
        mpfr_set_d(r27380, b, MPFR_RNDN);
        mpfr_mul(r27381, r27380, r27380, MPFR_RNDN);
        mpfr_add(r27382, r27379, r27381, MPFR_RNDN);
        ;
        mpfr_pow(r27384, r27382, r27383, MPFR_RNDN);
        ;
        ;
        mpfr_add(r27387, r27386, r27378, MPFR_RNDN);
        mpfr_mul(r27388, r27379, r27387, MPFR_RNDN);
        ;
        mpfr_mul(r27390, r27389, r27378, MPFR_RNDN);
        mpfr_sub(r27391, r27386, r27390, MPFR_RNDN);
        mpfr_mul(r27392, r27381, r27391, MPFR_RNDN);
        mpfr_add(r27393, r27388, r27392, MPFR_RNDN);
        mpfr_mul(r27394, r27385, r27393, MPFR_RNDN);
        mpfr_add(r27395, r27384, r27394, MPFR_RNDN);
        mpfr_sub(r27396, r27395, r27386, MPFR_RNDN);
        return mpfr_get_d(r27396, MPFR_RNDN);
}

static mpfr_t r27397, r27398, r27399, r27400, r27401, r27402, r27403, r27404, r27405, r27406, r27407, r27408, r27409, r27410, r27411, r27412, r27413, r27414, r27415;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27397);
        mpfr_init(r27398);
        mpfr_init(r27399);
        mpfr_init(r27400);
        mpfr_init(r27401);
        mpfr_init_set_str(r27402, "2", 10, MPFR_RNDN);
        mpfr_init(r27403);
        mpfr_init_set_str(r27404, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27405, "1", 10, MPFR_RNDN);
        mpfr_init(r27406);
        mpfr_init(r27407);
        mpfr_init_set_str(r27408, "3", 10, MPFR_RNDN);
        mpfr_init(r27409);
        mpfr_init(r27410);
        mpfr_init(r27411);
        mpfr_init(r27412);
        mpfr_init(r27413);
        mpfr_init(r27414);
        mpfr_init(r27415);
}

double f_fm(double a, double b) {
        mpfr_set_d(r27397, a, MPFR_RNDN);
        mpfr_mul(r27398, r27397, r27397, MPFR_RNDN);
        mpfr_set_d(r27399, b, MPFR_RNDN);
        mpfr_mul(r27400, r27399, r27399, MPFR_RNDN);
        mpfr_add(r27401, r27398, r27400, MPFR_RNDN);
        ;
        mpfr_pow(r27403, r27401, r27402, MPFR_RNDN);
        ;
        ;
        mpfr_add(r27406, r27405, r27397, MPFR_RNDN);
        mpfr_mul(r27407, r27398, r27406, MPFR_RNDN);
        ;
        mpfr_mul(r27409, r27408, r27397, MPFR_RNDN);
        mpfr_sub(r27410, r27405, r27409, MPFR_RNDN);
        mpfr_mul(r27411, r27400, r27410, MPFR_RNDN);
        mpfr_add(r27412, r27407, r27411, MPFR_RNDN);
        mpfr_mul(r27413, r27404, r27412, MPFR_RNDN);
        mpfr_add(r27414, r27403, r27413, MPFR_RNDN);
        mpfr_sub(r27415, r27414, r27405, MPFR_RNDN);
        return mpfr_get_d(r27415, MPFR_RNDN);
}

static mpfr_t r27416, r27417, r27418, r27419, r27420, r27421, r27422, r27423, r27424, r27425, r27426, r27427, r27428, r27429, r27430, r27431, r27432, r27433, r27434;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27416);
        mpfr_init(r27417);
        mpfr_init(r27418);
        mpfr_init(r27419);
        mpfr_init(r27420);
        mpfr_init_set_str(r27421, "2", 10, MPFR_RNDN);
        mpfr_init(r27422);
        mpfr_init_set_str(r27423, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27424, "1", 10, MPFR_RNDN);
        mpfr_init(r27425);
        mpfr_init(r27426);
        mpfr_init_set_str(r27427, "3", 10, MPFR_RNDN);
        mpfr_init(r27428);
        mpfr_init(r27429);
        mpfr_init(r27430);
        mpfr_init(r27431);
        mpfr_init(r27432);
        mpfr_init(r27433);
        mpfr_init(r27434);
}

double f_dm(double a, double b) {
        mpfr_set_d(r27416, a, MPFR_RNDN);
        mpfr_mul(r27417, r27416, r27416, MPFR_RNDN);
        mpfr_set_d(r27418, b, MPFR_RNDN);
        mpfr_mul(r27419, r27418, r27418, MPFR_RNDN);
        mpfr_add(r27420, r27417, r27419, MPFR_RNDN);
        ;
        mpfr_pow(r27422, r27420, r27421, MPFR_RNDN);
        ;
        ;
        mpfr_add(r27425, r27424, r27416, MPFR_RNDN);
        mpfr_mul(r27426, r27417, r27425, MPFR_RNDN);
        ;
        mpfr_mul(r27428, r27427, r27416, MPFR_RNDN);
        mpfr_sub(r27429, r27424, r27428, MPFR_RNDN);
        mpfr_mul(r27430, r27419, r27429, MPFR_RNDN);
        mpfr_add(r27431, r27426, r27430, MPFR_RNDN);
        mpfr_mul(r27432, r27423, r27431, MPFR_RNDN);
        mpfr_add(r27433, r27422, r27432, MPFR_RNDN);
        mpfr_sub(r27434, r27433, r27424, MPFR_RNDN);
        return mpfr_get_d(r27434, MPFR_RNDN);
}