#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 r26867 = a;
        float r26868 = r26867 * r26867;
        float r26869 = b;
        float r26870 = r26869 * r26869;
        float r26871 = r26868 + r26870;
        float r26872 = 2;
        float r26873 = pow(r26871, r26872);
        float r26874 = 4;
        float r26875 = 1;
        float r26876 = r26875 + r26867;
        float r26877 = r26868 * r26876;
        float r26878 = 3;
        float r26879 = r26878 * r26867;
        float r26880 = r26875 - r26879;
        float r26881 = r26870 * r26880;
        float r26882 = r26877 + r26881;
        float r26883 = r26874 * r26882;
        float r26884 = r26873 + r26883;
        float r26885 = r26884 - r26875;
        return r26885;
}

double f_id(double a, double b) {
        double r26886 = a;
        double r26887 = r26886 * r26886;
        double r26888 = b;
        double r26889 = r26888 * r26888;
        double r26890 = r26887 + r26889;
        double r26891 = 2;
        double r26892 = pow(r26890, r26891);
        double r26893 = 4;
        double r26894 = 1;
        double r26895 = r26894 + r26886;
        double r26896 = r26887 * r26895;
        double r26897 = 3;
        double r26898 = r26897 * r26886;
        double r26899 = r26894 - r26898;
        double r26900 = r26889 * r26899;
        double r26901 = r26896 + r26900;
        double r26902 = r26893 * r26901;
        double r26903 = r26892 + r26902;
        double r26904 = r26903 - r26894;
        return r26904;
}


double f_of(float a, float b) {
        float r26905 = a;
        float r26906 = b;
        float r26907 = r26906 * r26906;
        float r26908 = fma(r26905, r26905, r26907);
        float r26909 = hypot(r26906, r26905);
        float r26910 = r26908 * r26909;
        float r26911 = r26905 * r26905;
        float r26912 = fma(r26906, r26906, r26911);
        float r26913 = sqrt(r26912);
        float r26914 = r26910 * r26913;
        float r26915 = 4;
        float r26916 = 1;
        float r26917 = r26916 + r26905;
        float r26918 = r26911 * r26917;
        float r26919 = 3;
        float r26920 = r26919 * r26905;
        float r26921 = r26916 - r26920;
        float r26922 = r26907 * r26921;
        float r26923 = r26918 + r26922;
        float r26924 = r26915 * r26923;
        float r26925 = r26914 + r26924;
        float r26926 = r26925 - r26916;
        return r26926;
}

double f_od(double a, double b) {
        double r26927 = a;
        double r26928 = b;
        double r26929 = r26928 * r26928;
        double r26930 = fma(r26927, r26927, r26929);
        double r26931 = hypot(r26928, r26927);
        double r26932 = r26930 * r26931;
        double r26933 = r26927 * r26927;
        double r26934 = fma(r26928, r26928, r26933);
        double r26935 = sqrt(r26934);
        double r26936 = r26932 * r26935;
        double r26937 = 4;
        double r26938 = 1;
        double r26939 = r26938 + r26927;
        double r26940 = r26933 * r26939;
        double r26941 = 3;
        double r26942 = r26941 * r26927;
        double r26943 = r26938 - r26942;
        double r26944 = r26929 * r26943;
        double r26945 = r26940 + r26944;
        double r26946 = r26937 * r26945;
        double r26947 = r26936 + r26946;
        double r26948 = r26947 - r26938;
        return r26948;
}

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 r26949, r26950, r26951, r26952, r26953, r26954, r26955, r26956, r26957, r26958, r26959, r26960, r26961, r26962, r26963, r26964, r26965, r26966, r26967;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r26949);
        mpfr_init(r26950);
        mpfr_init(r26951);
        mpfr_init(r26952);
        mpfr_init(r26953);
        mpfr_init_set_str(r26954, "2", 10, MPFR_RNDN);
        mpfr_init(r26955);
        mpfr_init_set_str(r26956, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r26957, "1", 10, MPFR_RNDN);
        mpfr_init(r26958);
        mpfr_init(r26959);
        mpfr_init_set_str(r26960, "3", 10, MPFR_RNDN);
        mpfr_init(r26961);
        mpfr_init(r26962);
        mpfr_init(r26963);
        mpfr_init(r26964);
        mpfr_init(r26965);
        mpfr_init(r26966);
        mpfr_init(r26967);
}

double f_im(double a, double b) {
        mpfr_set_d(r26949, a, MPFR_RNDN);
        mpfr_mul(r26950, r26949, r26949, MPFR_RNDN);
        mpfr_set_d(r26951, b, MPFR_RNDN);
        mpfr_mul(r26952, r26951, r26951, MPFR_RNDN);
        mpfr_add(r26953, r26950, r26952, MPFR_RNDN);
        ;
        mpfr_pow(r26955, r26953, r26954, MPFR_RNDN);
        ;
        ;
        mpfr_add(r26958, r26957, r26949, MPFR_RNDN);
        mpfr_mul(r26959, r26950, r26958, MPFR_RNDN);
        ;
        mpfr_mul(r26961, r26960, r26949, MPFR_RNDN);
        mpfr_sub(r26962, r26957, r26961, MPFR_RNDN);
        mpfr_mul(r26963, r26952, r26962, MPFR_RNDN);
        mpfr_add(r26964, r26959, r26963, MPFR_RNDN);
        mpfr_mul(r26965, r26956, r26964, MPFR_RNDN);
        mpfr_add(r26966, r26955, r26965, MPFR_RNDN);
        mpfr_sub(r26967, r26966, r26957, MPFR_RNDN);
        return mpfr_get_d(r26967, MPFR_RNDN);
}

static mpfr_t r26968, r26969, r26970, r26971, r26972, r26973, r26974, r26975, r26976, r26977, r26978, r26979, r26980, r26981, r26982, r26983, r26984, r26985, r26986, r26987, r26988, r26989;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r26968);
        mpfr_init(r26969);
        mpfr_init(r26970);
        mpfr_init(r26971);
        mpfr_init(r26972);
        mpfr_init(r26973);
        mpfr_init(r26974);
        mpfr_init(r26975);
        mpfr_init(r26976);
        mpfr_init(r26977);
        mpfr_init_set_str(r26978, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r26979, "1", 10, MPFR_RNDN);
        mpfr_init(r26980);
        mpfr_init(r26981);
        mpfr_init_set_str(r26982, "3", 10, MPFR_RNDN);
        mpfr_init(r26983);
        mpfr_init(r26984);
        mpfr_init(r26985);
        mpfr_init(r26986);
        mpfr_init(r26987);
        mpfr_init(r26988);
        mpfr_init(r26989);
}

double f_fm(double a, double b) {
        mpfr_set_d(r26968, a, MPFR_RNDN);
        mpfr_set_d(r26969, b, MPFR_RNDN);
        mpfr_mul(r26970, r26969, r26969, MPFR_RNDN);
        mpfr_fma(r26971, r26968, r26968, r26970, MPFR_RNDN);
        mpfr_hypot(r26972, r26969, r26968, MPFR_RNDN);
        mpfr_mul(r26973, r26971, r26972, MPFR_RNDN);
        mpfr_mul(r26974, r26968, r26968, MPFR_RNDN);
        mpfr_fma(r26975, r26969, r26969, r26974, MPFR_RNDN);
        mpfr_sqrt(r26976, r26975, MPFR_RNDN);
        mpfr_mul(r26977, r26973, r26976, MPFR_RNDN);
        ;
        ;
        mpfr_add(r26980, r26979, r26968, MPFR_RNDN);
        mpfr_mul(r26981, r26974, r26980, MPFR_RNDN);
        ;
        mpfr_mul(r26983, r26982, r26968, MPFR_RNDN);
        mpfr_sub(r26984, r26979, r26983, MPFR_RNDN);
        mpfr_mul(r26985, r26970, r26984, MPFR_RNDN);
        mpfr_add(r26986, r26981, r26985, MPFR_RNDN);
        mpfr_mul(r26987, r26978, r26986, MPFR_RNDN);
        mpfr_add(r26988, r26977, r26987, MPFR_RNDN);
        mpfr_sub(r26989, r26988, r26979, MPFR_RNDN);
        return mpfr_get_d(r26989, MPFR_RNDN);
}

static mpfr_t r26990, r26991, r26992, r26993, r26994, r26995, r26996, r26997, r26998, r26999, r27000, r27001, r27002, r27003, r27004, r27005, r27006, r27007, r27008, r27009, r27010, r27011;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r26990);
        mpfr_init(r26991);
        mpfr_init(r26992);
        mpfr_init(r26993);
        mpfr_init(r26994);
        mpfr_init(r26995);
        mpfr_init(r26996);
        mpfr_init(r26997);
        mpfr_init(r26998);
        mpfr_init(r26999);
        mpfr_init_set_str(r27000, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27001, "1", 10, MPFR_RNDN);
        mpfr_init(r27002);
        mpfr_init(r27003);
        mpfr_init_set_str(r27004, "3", 10, MPFR_RNDN);
        mpfr_init(r27005);
        mpfr_init(r27006);
        mpfr_init(r27007);
        mpfr_init(r27008);
        mpfr_init(r27009);
        mpfr_init(r27010);
        mpfr_init(r27011);
}

double f_dm(double a, double b) {
        mpfr_set_d(r26990, a, MPFR_RNDN);
        mpfr_set_d(r26991, b, MPFR_RNDN);
        mpfr_mul(r26992, r26991, r26991, MPFR_RNDN);
        mpfr_fma(r26993, r26990, r26990, r26992, MPFR_RNDN);
        mpfr_hypot(r26994, r26991, r26990, MPFR_RNDN);
        mpfr_mul(r26995, r26993, r26994, MPFR_RNDN);
        mpfr_mul(r26996, r26990, r26990, MPFR_RNDN);
        mpfr_fma(r26997, r26991, r26991, r26996, MPFR_RNDN);
        mpfr_sqrt(r26998, r26997, MPFR_RNDN);
        mpfr_mul(r26999, r26995, r26998, MPFR_RNDN);
        ;
        ;
        mpfr_add(r27002, r27001, r26990, MPFR_RNDN);
        mpfr_mul(r27003, r26996, r27002, MPFR_RNDN);
        ;
        mpfr_mul(r27005, r27004, r26990, MPFR_RNDN);
        mpfr_sub(r27006, r27001, r27005, MPFR_RNDN);
        mpfr_mul(r27007, r26992, r27006, MPFR_RNDN);
        mpfr_add(r27008, r27003, r27007, MPFR_RNDN);
        mpfr_mul(r27009, r27000, r27008, MPFR_RNDN);
        mpfr_add(r27010, r26999, r27009, MPFR_RNDN);
        mpfr_sub(r27011, r27010, r27001, MPFR_RNDN);
        return mpfr_get_d(r27011, MPFR_RNDN);
}