#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 r25251 = a;
        float r25252 = r25251 * r25251;
        float r25253 = b;
        float r25254 = r25253 * r25253;
        float r25255 = r25252 + r25254;
        float r25256 = 2;
        float r25257 = pow(r25255, r25256);
        float r25258 = 4;
        float r25259 = 1;
        float r25260 = r25259 - r25251;
        float r25261 = r25252 * r25260;
        float r25262 = 3;
        float r25263 = r25262 + r25251;
        float r25264 = r25254 * r25263;
        float r25265 = r25261 + r25264;
        float r25266 = r25258 * r25265;
        float r25267 = r25257 + r25266;
        float r25268 = r25267 - r25259;
        return r25268;
}

double f_id(double a, double b) {
        double r25269 = a;
        double r25270 = r25269 * r25269;
        double r25271 = b;
        double r25272 = r25271 * r25271;
        double r25273 = r25270 + r25272;
        double r25274 = 2;
        double r25275 = pow(r25273, r25274);
        double r25276 = 4;
        double r25277 = 1;
        double r25278 = r25277 - r25269;
        double r25279 = r25270 * r25278;
        double r25280 = 3;
        double r25281 = r25280 + r25269;
        double r25282 = r25272 * r25281;
        double r25283 = r25279 + r25282;
        double r25284 = r25276 * r25283;
        double r25285 = r25275 + r25284;
        double r25286 = r25285 - r25277;
        return r25286;
}


double f_of(float a, float b) {
        float r25287 = a;
        float r25288 = r25287 * r25287;
        float r25289 = b;
        float r25290 = r25289 * r25289;
        float r25291 = r25288 + r25290;
        float r25292 = 2;
        float r25293 = pow(r25291, r25292);
        float r25294 = cbrt(r25293);
        float r25295 = r25294 * r25294;
        float r25296 = fma(r25289, r25289, r25288);
        float r25297 = cbrt(r25296);
        float r25298 = r25297 * r25296;
        float r25299 = cbrt(r25298);
        float r25300 = cbrt(r25291);
        float r25301 = pow(r25300, r25292);
        float r25302 = cbrt(r25301);
        float r25303 = r25299 * r25302;
        float r25304 = r25295 * r25303;
        float r25305 = 4;
        float r25306 = 1;
        float r25307 = r25306 - r25287;
        float r25308 = r25288 * r25307;
        float r25309 = 3;
        float r25310 = r25309 + r25287;
        float r25311 = r25290 * r25310;
        float r25312 = r25308 + r25311;
        float r25313 = r25305 * r25312;
        float r25314 = r25304 + r25313;
        float r25315 = r25314 - r25306;
        return r25315;
}

double f_od(double a, double b) {
        double r25316 = a;
        double r25317 = r25316 * r25316;
        double r25318 = b;
        double r25319 = r25318 * r25318;
        double r25320 = r25317 + r25319;
        double r25321 = 2;
        double r25322 = pow(r25320, r25321);
        double r25323 = cbrt(r25322);
        double r25324 = r25323 * r25323;
        double r25325 = fma(r25318, r25318, r25317);
        double r25326 = cbrt(r25325);
        double r25327 = r25326 * r25325;
        double r25328 = cbrt(r25327);
        double r25329 = cbrt(r25320);
        double r25330 = pow(r25329, r25321);
        double r25331 = cbrt(r25330);
        double r25332 = r25328 * r25331;
        double r25333 = r25324 * r25332;
        double r25334 = 4;
        double r25335 = 1;
        double r25336 = r25335 - r25316;
        double r25337 = r25317 * r25336;
        double r25338 = 3;
        double r25339 = r25338 + r25316;
        double r25340 = r25319 * r25339;
        double r25341 = r25337 + r25340;
        double r25342 = r25334 * r25341;
        double r25343 = r25333 + r25342;
        double r25344 = r25343 - r25335;
        return r25344;
}

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 r25345, r25346, r25347, r25348, r25349, r25350, r25351, r25352, r25353, r25354, r25355, r25356, r25357, r25358, r25359, r25360, r25361, r25362;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r25345);
        mpfr_init(r25346);
        mpfr_init(r25347);
        mpfr_init(r25348);
        mpfr_init(r25349);
        mpfr_init_set_str(r25350, "2", 10, MPFR_RNDN);
        mpfr_init(r25351);
        mpfr_init_set_str(r25352, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25353, "1", 10, MPFR_RNDN);
        mpfr_init(r25354);
        mpfr_init(r25355);
        mpfr_init_set_str(r25356, "3", 10, MPFR_RNDN);
        mpfr_init(r25357);
        mpfr_init(r25358);
        mpfr_init(r25359);
        mpfr_init(r25360);
        mpfr_init(r25361);
        mpfr_init(r25362);
}

double f_im(double a, double b) {
        mpfr_set_d(r25345, a, MPFR_RNDN);
        mpfr_mul(r25346, r25345, r25345, MPFR_RNDN);
        mpfr_set_d(r25347, b, MPFR_RNDN);
        mpfr_mul(r25348, r25347, r25347, MPFR_RNDN);
        mpfr_add(r25349, r25346, r25348, MPFR_RNDN);
        ;
        mpfr_pow(r25351, r25349, r25350, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25354, r25353, r25345, MPFR_RNDN);
        mpfr_mul(r25355, r25346, r25354, MPFR_RNDN);
        ;
        mpfr_add(r25357, r25356, r25345, MPFR_RNDN);
        mpfr_mul(r25358, r25348, r25357, MPFR_RNDN);
        mpfr_add(r25359, r25355, r25358, MPFR_RNDN);
        mpfr_mul(r25360, r25352, r25359, MPFR_RNDN);
        mpfr_add(r25361, r25351, r25360, MPFR_RNDN);
        mpfr_sub(r25362, r25361, r25353, MPFR_RNDN);
        return mpfr_get_d(r25362, MPFR_RNDN);
}

static mpfr_t r25363, r25364, r25365, r25366, r25367, r25368, r25369, r25370, r25371, r25372, r25373, r25374, r25375, r25376, r25377, r25378, r25379, r25380, r25381, r25382, r25383, r25384, r25385, r25386, r25387, r25388, r25389, r25390, r25391;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r25363);
        mpfr_init(r25364);
        mpfr_init(r25365);
        mpfr_init(r25366);
        mpfr_init(r25367);
        mpfr_init_set_str(r25368, "2", 10, MPFR_RNDN);
        mpfr_init(r25369);
        mpfr_init(r25370);
        mpfr_init(r25371);
        mpfr_init(r25372);
        mpfr_init(r25373);
        mpfr_init(r25374);
        mpfr_init(r25375);
        mpfr_init(r25376);
        mpfr_init(r25377);
        mpfr_init(r25378);
        mpfr_init(r25379);
        mpfr_init(r25380);
        mpfr_init_set_str(r25381, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25382, "1", 10, MPFR_RNDN);
        mpfr_init(r25383);
        mpfr_init(r25384);
        mpfr_init_set_str(r25385, "3", 10, MPFR_RNDN);
        mpfr_init(r25386);
        mpfr_init(r25387);
        mpfr_init(r25388);
        mpfr_init(r25389);
        mpfr_init(r25390);
        mpfr_init(r25391);
}

double f_fm(double a, double b) {
        mpfr_set_d(r25363, a, MPFR_RNDN);
        mpfr_mul(r25364, r25363, r25363, MPFR_RNDN);
        mpfr_set_d(r25365, b, MPFR_RNDN);
        mpfr_mul(r25366, r25365, r25365, MPFR_RNDN);
        mpfr_add(r25367, r25364, r25366, MPFR_RNDN);
        ;
        mpfr_pow(r25369, r25367, r25368, MPFR_RNDN);
        mpfr_cbrt(r25370, r25369, MPFR_RNDN);
        mpfr_mul(r25371, r25370, r25370, MPFR_RNDN);
        mpfr_fma(r25372, r25365, r25365, r25364, MPFR_RNDN);
        mpfr_cbrt(r25373, r25372, MPFR_RNDN);
        mpfr_mul(r25374, r25373, r25372, MPFR_RNDN);
        mpfr_cbrt(r25375, r25374, MPFR_RNDN);
        mpfr_cbrt(r25376, r25367, MPFR_RNDN);
        mpfr_pow(r25377, r25376, r25368, MPFR_RNDN);
        mpfr_cbrt(r25378, r25377, MPFR_RNDN);
        mpfr_mul(r25379, r25375, r25378, MPFR_RNDN);
        mpfr_mul(r25380, r25371, r25379, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25383, r25382, r25363, MPFR_RNDN);
        mpfr_mul(r25384, r25364, r25383, MPFR_RNDN);
        ;
        mpfr_add(r25386, r25385, r25363, MPFR_RNDN);
        mpfr_mul(r25387, r25366, r25386, MPFR_RNDN);
        mpfr_add(r25388, r25384, r25387, MPFR_RNDN);
        mpfr_mul(r25389, r25381, r25388, MPFR_RNDN);
        mpfr_add(r25390, r25380, r25389, MPFR_RNDN);
        mpfr_sub(r25391, r25390, r25382, MPFR_RNDN);
        return mpfr_get_d(r25391, MPFR_RNDN);
}

static mpfr_t r25392, r25393, r25394, r25395, r25396, r25397, r25398, r25399, r25400, r25401, r25402, r25403, r25404, r25405, r25406, r25407, r25408, r25409, r25410, r25411, r25412, r25413, r25414, r25415, r25416, r25417, r25418, r25419, r25420;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r25392);
        mpfr_init(r25393);
        mpfr_init(r25394);
        mpfr_init(r25395);
        mpfr_init(r25396);
        mpfr_init_set_str(r25397, "2", 10, MPFR_RNDN);
        mpfr_init(r25398);
        mpfr_init(r25399);
        mpfr_init(r25400);
        mpfr_init(r25401);
        mpfr_init(r25402);
        mpfr_init(r25403);
        mpfr_init(r25404);
        mpfr_init(r25405);
        mpfr_init(r25406);
        mpfr_init(r25407);
        mpfr_init(r25408);
        mpfr_init(r25409);
        mpfr_init_set_str(r25410, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25411, "1", 10, MPFR_RNDN);
        mpfr_init(r25412);
        mpfr_init(r25413);
        mpfr_init_set_str(r25414, "3", 10, MPFR_RNDN);
        mpfr_init(r25415);
        mpfr_init(r25416);
        mpfr_init(r25417);
        mpfr_init(r25418);
        mpfr_init(r25419);
        mpfr_init(r25420);
}

double f_dm(double a, double b) {
        mpfr_set_d(r25392, a, MPFR_RNDN);
        mpfr_mul(r25393, r25392, r25392, MPFR_RNDN);
        mpfr_set_d(r25394, b, MPFR_RNDN);
        mpfr_mul(r25395, r25394, r25394, MPFR_RNDN);
        mpfr_add(r25396, r25393, r25395, MPFR_RNDN);
        ;
        mpfr_pow(r25398, r25396, r25397, MPFR_RNDN);
        mpfr_cbrt(r25399, r25398, MPFR_RNDN);
        mpfr_mul(r25400, r25399, r25399, MPFR_RNDN);
        mpfr_fma(r25401, r25394, r25394, r25393, MPFR_RNDN);
        mpfr_cbrt(r25402, r25401, MPFR_RNDN);
        mpfr_mul(r25403, r25402, r25401, MPFR_RNDN);
        mpfr_cbrt(r25404, r25403, MPFR_RNDN);
        mpfr_cbrt(r25405, r25396, MPFR_RNDN);
        mpfr_pow(r25406, r25405, r25397, MPFR_RNDN);
        mpfr_cbrt(r25407, r25406, MPFR_RNDN);
        mpfr_mul(r25408, r25404, r25407, MPFR_RNDN);
        mpfr_mul(r25409, r25400, r25408, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25412, r25411, r25392, MPFR_RNDN);
        mpfr_mul(r25413, r25393, r25412, MPFR_RNDN);
        ;
        mpfr_add(r25415, r25414, r25392, MPFR_RNDN);
        mpfr_mul(r25416, r25395, r25415, MPFR_RNDN);
        mpfr_add(r25417, r25413, r25416, MPFR_RNDN);
        mpfr_mul(r25418, r25410, r25417, MPFR_RNDN);
        mpfr_add(r25419, r25409, r25418, MPFR_RNDN);
        mpfr_sub(r25420, r25419, r25411, MPFR_RNDN);
        return mpfr_get_d(r25420, MPFR_RNDN);
}