#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 r25127 = a;
        float r25128 = r25127 * r25127;
        float r25129 = b;
        float r25130 = r25129 * r25129;
        float r25131 = r25128 + r25130;
        float r25132 = 2;
        float r25133 = pow(r25131, r25132);
        float r25134 = 4;
        float r25135 = 1;
        float r25136 = r25135 + r25127;
        float r25137 = r25128 * r25136;
        float r25138 = 3;
        float r25139 = r25138 * r25127;
        float r25140 = r25135 - r25139;
        float r25141 = r25130 * r25140;
        float r25142 = r25137 + r25141;
        float r25143 = r25134 * r25142;
        float r25144 = r25133 + r25143;
        float r25145 = r25144 - r25135;
        return r25145;
}

double f_id(double a, double b) {
        double r25146 = a;
        double r25147 = r25146 * r25146;
        double r25148 = b;
        double r25149 = r25148 * r25148;
        double r25150 = r25147 + r25149;
        double r25151 = 2;
        double r25152 = pow(r25150, r25151);
        double r25153 = 4;
        double r25154 = 1;
        double r25155 = r25154 + r25146;
        double r25156 = r25147 * r25155;
        double r25157 = 3;
        double r25158 = r25157 * r25146;
        double r25159 = r25154 - r25158;
        double r25160 = r25149 * r25159;
        double r25161 = r25156 + r25160;
        double r25162 = r25153 * r25161;
        double r25163 = r25152 + r25162;
        double r25164 = r25163 - r25154;
        return r25164;
}


double f_of(float a, float b) {
        float r25165 = a;
        float r25166 = r25165 * r25165;
        float r25167 = b;
        float r25168 = r25167 * r25167;
        float r25169 = r25166 + r25168;
        float r25170 = 2;
        float r25171 = pow(r25169, r25170);
        float r25172 = 4;
        float r25173 = 1;
        float r25174 = r25173 + r25165;
        float r25175 = r25166 * r25174;
        float r25176 = 3;
        float r25177 = r25176 * r25165;
        float r25178 = r25173 - r25177;
        float r25179 = r25168 * r25178;
        float r25180 = r25175 + r25179;
        float r25181 = r25172 * r25180;
        float r25182 = r25171 + r25181;
        float r25183 = sqrt(r25182);
        float r25184 = r25183 * r25183;
        float r25185 = r25184 - r25173;
        return r25185;
}

double f_od(double a, double b) {
        double r25186 = a;
        double r25187 = r25186 * r25186;
        double r25188 = b;
        double r25189 = r25188 * r25188;
        double r25190 = r25187 + r25189;
        double r25191 = 2;
        double r25192 = pow(r25190, r25191);
        double r25193 = 4;
        double r25194 = 1;
        double r25195 = r25194 + r25186;
        double r25196 = r25187 * r25195;
        double r25197 = 3;
        double r25198 = r25197 * r25186;
        double r25199 = r25194 - r25198;
        double r25200 = r25189 * r25199;
        double r25201 = r25196 + r25200;
        double r25202 = r25193 * r25201;
        double r25203 = r25192 + r25202;
        double r25204 = sqrt(r25203);
        double r25205 = r25204 * r25204;
        double r25206 = r25205 - r25194;
        return r25206;
}

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 r25207, r25208, r25209, r25210, r25211, r25212, r25213, r25214, r25215, r25216, r25217, r25218, r25219, r25220, r25221, r25222, r25223, r25224, r25225;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r25207);
        mpfr_init(r25208);
        mpfr_init(r25209);
        mpfr_init(r25210);
        mpfr_init(r25211);
        mpfr_init_set_str(r25212, "2", 10, MPFR_RNDN);
        mpfr_init(r25213);
        mpfr_init_set_str(r25214, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25215, "1", 10, MPFR_RNDN);
        mpfr_init(r25216);
        mpfr_init(r25217);
        mpfr_init_set_str(r25218, "3", 10, MPFR_RNDN);
        mpfr_init(r25219);
        mpfr_init(r25220);
        mpfr_init(r25221);
        mpfr_init(r25222);
        mpfr_init(r25223);
        mpfr_init(r25224);
        mpfr_init(r25225);
}

double f_im(double a, double b) {
        mpfr_set_d(r25207, a, MPFR_RNDN);
        mpfr_mul(r25208, r25207, r25207, MPFR_RNDN);
        mpfr_set_d(r25209, b, MPFR_RNDN);
        mpfr_mul(r25210, r25209, r25209, MPFR_RNDN);
        mpfr_add(r25211, r25208, r25210, MPFR_RNDN);
        ;
        mpfr_pow(r25213, r25211, r25212, MPFR_RNDN);
        ;
        ;
        mpfr_add(r25216, r25215, r25207, MPFR_RNDN);
        mpfr_mul(r25217, r25208, r25216, MPFR_RNDN);
        ;
        mpfr_mul(r25219, r25218, r25207, MPFR_RNDN);
        mpfr_sub(r25220, r25215, r25219, MPFR_RNDN);
        mpfr_mul(r25221, r25210, r25220, MPFR_RNDN);
        mpfr_add(r25222, r25217, r25221, MPFR_RNDN);
        mpfr_mul(r25223, r25214, r25222, MPFR_RNDN);
        mpfr_add(r25224, r25213, r25223, MPFR_RNDN);
        mpfr_sub(r25225, r25224, r25215, MPFR_RNDN);
        return mpfr_get_d(r25225, MPFR_RNDN);
}

static mpfr_t r25226, r25227, r25228, r25229, r25230, r25231, r25232, r25233, r25234, r25235, r25236, r25237, r25238, r25239, r25240, r25241, r25242, r25243, r25244, r25245, r25246;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r25226);
        mpfr_init(r25227);
        mpfr_init(r25228);
        mpfr_init(r25229);
        mpfr_init(r25230);
        mpfr_init_set_str(r25231, "2", 10, MPFR_RNDN);
        mpfr_init(r25232);
        mpfr_init_set_str(r25233, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25234, "1", 10, MPFR_RNDN);
        mpfr_init(r25235);
        mpfr_init(r25236);
        mpfr_init_set_str(r25237, "3", 10, MPFR_RNDN);
        mpfr_init(r25238);
        mpfr_init(r25239);
        mpfr_init(r25240);
        mpfr_init(r25241);
        mpfr_init(r25242);
        mpfr_init(r25243);
        mpfr_init(r25244);
        mpfr_init(r25245);
        mpfr_init(r25246);
}

double f_fm(double a, double b) {
        mpfr_set_d(r25226, a, MPFR_RNDN);
        mpfr_mul(r25227, r25226, r25226, MPFR_RNDN);
        mpfr_set_d(r25228, b, MPFR_RNDN);
        mpfr_mul(r25229, r25228, r25228, MPFR_RNDN);
        mpfr_add(r25230, r25227, r25229, MPFR_RNDN);
        ;
        mpfr_pow(r25232, r25230, r25231, MPFR_RNDN);
        ;
        ;
        mpfr_add(r25235, r25234, r25226, MPFR_RNDN);
        mpfr_mul(r25236, r25227, r25235, MPFR_RNDN);
        ;
        mpfr_mul(r25238, r25237, r25226, MPFR_RNDN);
        mpfr_sub(r25239, r25234, r25238, MPFR_RNDN);
        mpfr_mul(r25240, r25229, r25239, MPFR_RNDN);
        mpfr_add(r25241, r25236, r25240, MPFR_RNDN);
        mpfr_mul(r25242, r25233, r25241, MPFR_RNDN);
        mpfr_add(r25243, r25232, r25242, MPFR_RNDN);
        mpfr_sqrt(r25244, r25243, MPFR_RNDN);
        mpfr_mul(r25245, r25244, r25244, MPFR_RNDN);
        mpfr_sub(r25246, r25245, r25234, MPFR_RNDN);
        return mpfr_get_d(r25246, MPFR_RNDN);
}

static mpfr_t r25247, r25248, r25249, r25250, r25251, r25252, r25253, r25254, r25255, r25256, r25257, r25258, r25259, r25260, r25261, r25262, r25263, r25264, r25265, r25266, r25267;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r25247);
        mpfr_init(r25248);
        mpfr_init(r25249);
        mpfr_init(r25250);
        mpfr_init(r25251);
        mpfr_init_set_str(r25252, "2", 10, MPFR_RNDN);
        mpfr_init(r25253);
        mpfr_init_set_str(r25254, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r25255, "1", 10, MPFR_RNDN);
        mpfr_init(r25256);
        mpfr_init(r25257);
        mpfr_init_set_str(r25258, "3", 10, MPFR_RNDN);
        mpfr_init(r25259);
        mpfr_init(r25260);
        mpfr_init(r25261);
        mpfr_init(r25262);
        mpfr_init(r25263);
        mpfr_init(r25264);
        mpfr_init(r25265);
        mpfr_init(r25266);
        mpfr_init(r25267);
}

double f_dm(double a, double b) {
        mpfr_set_d(r25247, a, MPFR_RNDN);
        mpfr_mul(r25248, r25247, r25247, MPFR_RNDN);
        mpfr_set_d(r25249, b, MPFR_RNDN);
        mpfr_mul(r25250, r25249, r25249, MPFR_RNDN);
        mpfr_add(r25251, r25248, r25250, MPFR_RNDN);
        ;
        mpfr_pow(r25253, r25251, r25252, MPFR_RNDN);
        ;
        ;
        mpfr_add(r25256, r25255, r25247, MPFR_RNDN);
        mpfr_mul(r25257, r25248, r25256, MPFR_RNDN);
        ;
        mpfr_mul(r25259, r25258, r25247, MPFR_RNDN);
        mpfr_sub(r25260, r25255, r25259, MPFR_RNDN);
        mpfr_mul(r25261, r25250, r25260, MPFR_RNDN);
        mpfr_add(r25262, r25257, r25261, MPFR_RNDN);
        mpfr_mul(r25263, r25254, r25262, MPFR_RNDN);
        mpfr_add(r25264, r25253, r25263, MPFR_RNDN);
        mpfr_sqrt(r25265, r25264, MPFR_RNDN);
        mpfr_mul(r25266, r25265, r25265, MPFR_RNDN);
        mpfr_sub(r25267, r25266, r25255, MPFR_RNDN);
        return mpfr_get_d(r25267, MPFR_RNDN);
}