#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 r26860 = a;
        float r26861 = r26860 * r26860;
        float r26862 = b;
        float r26863 = r26862 * r26862;
        float r26864 = r26861 + r26863;
        float r26865 = 2;
        float r26866 = pow(r26864, r26865);
        float r26867 = 4;
        float r26868 = 1;
        float r26869 = r26868 + r26860;
        float r26870 = r26861 * r26869;
        float r26871 = 3;
        float r26872 = r26871 * r26860;
        float r26873 = r26868 - r26872;
        float r26874 = r26863 * r26873;
        float r26875 = r26870 + r26874;
        float r26876 = r26867 * r26875;
        float r26877 = r26866 + r26876;
        float r26878 = r26877 - r26868;
        return r26878;
}

double f_id(double a, double b) {
        double r26879 = a;
        double r26880 = r26879 * r26879;
        double r26881 = b;
        double r26882 = r26881 * r26881;
        double r26883 = r26880 + r26882;
        double r26884 = 2;
        double r26885 = pow(r26883, r26884);
        double r26886 = 4;
        double r26887 = 1;
        double r26888 = r26887 + r26879;
        double r26889 = r26880 * r26888;
        double r26890 = 3;
        double r26891 = r26890 * r26879;
        double r26892 = r26887 - r26891;
        double r26893 = r26882 * r26892;
        double r26894 = r26889 + r26893;
        double r26895 = r26886 * r26894;
        double r26896 = r26885 + r26895;
        double r26897 = r26896 - r26887;
        return r26897;
}


double f_of(float a, float b) {
        float r26898 = a;
        float r26899 = r26898 * r26898;
        float r26900 = b;
        float r26901 = r26900 * r26900;
        float r26902 = r26899 + r26901;
        float r26903 = 2;
        float r26904 = pow(r26902, r26903);
        float r26905 = 4;
        float r26906 = 1;
        float r26907 = r26906 + r26898;
        float r26908 = r26899 * r26907;
        float r26909 = 3;
        float r26910 = r26909 * r26898;
        float r26911 = r26906 - r26910;
        float r26912 = r26901 * r26911;
        float r26913 = r26908 + r26912;
        float r26914 = r26905 * r26913;
        float r26915 = r26904 + r26914;
        float r26916 = sqrt(r26915);
        float r26917 = r26916 * r26916;
        float r26918 = r26917 - r26906;
        return r26918;
}

double f_od(double a, double b) {
        double r26919 = a;
        double r26920 = r26919 * r26919;
        double r26921 = b;
        double r26922 = r26921 * r26921;
        double r26923 = r26920 + r26922;
        double r26924 = 2;
        double r26925 = pow(r26923, r26924);
        double r26926 = 4;
        double r26927 = 1;
        double r26928 = r26927 + r26919;
        double r26929 = r26920 * r26928;
        double r26930 = 3;
        double r26931 = r26930 * r26919;
        double r26932 = r26927 - r26931;
        double r26933 = r26922 * r26932;
        double r26934 = r26929 + r26933;
        double r26935 = r26926 * r26934;
        double r26936 = r26925 + r26935;
        double r26937 = sqrt(r26936);
        double r26938 = r26937 * r26937;
        double r26939 = r26938 - r26927;
        return r26939;
}

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 r26940, r26941, r26942, r26943, r26944, r26945, r26946, r26947, r26948, r26949, r26950, r26951, r26952, r26953, r26954, r26955, r26956, r26957, r26958;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26940);
        mpfr_init(r26941);
        mpfr_init(r26942);
        mpfr_init(r26943);
        mpfr_init(r26944);
        mpfr_init_set_str(r26945, "2", 10, MPFR_RNDN);
        mpfr_init(r26946);
        mpfr_init_set_str(r26947, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r26948, "1", 10, MPFR_RNDN);
        mpfr_init(r26949);
        mpfr_init(r26950);
        mpfr_init_set_str(r26951, "3", 10, MPFR_RNDN);
        mpfr_init(r26952);
        mpfr_init(r26953);
        mpfr_init(r26954);
        mpfr_init(r26955);
        mpfr_init(r26956);
        mpfr_init(r26957);
        mpfr_init(r26958);
}

double f_im(double a, double b) {
        mpfr_set_d(r26940, a, MPFR_RNDN);
        mpfr_mul(r26941, r26940, r26940, MPFR_RNDN);
        mpfr_set_d(r26942, b, MPFR_RNDN);
        mpfr_mul(r26943, r26942, r26942, MPFR_RNDN);
        mpfr_add(r26944, r26941, r26943, MPFR_RNDN);
        ;
        mpfr_pow(r26946, r26944, r26945, MPFR_RNDN);
        ;
        ;
        mpfr_add(r26949, r26948, r26940, MPFR_RNDN);
        mpfr_mul(r26950, r26941, r26949, MPFR_RNDN);
        ;
        mpfr_mul(r26952, r26951, r26940, MPFR_RNDN);
        mpfr_sub(r26953, r26948, r26952, MPFR_RNDN);
        mpfr_mul(r26954, r26943, r26953, MPFR_RNDN);
        mpfr_add(r26955, r26950, r26954, MPFR_RNDN);
        mpfr_mul(r26956, r26947, r26955, MPFR_RNDN);
        mpfr_add(r26957, r26946, r26956, MPFR_RNDN);
        mpfr_sub(r26958, r26957, r26948, MPFR_RNDN);
        return mpfr_get_d(r26958, MPFR_RNDN);
}

static mpfr_t r26959, r26960, r26961, r26962, r26963, r26964, r26965, r26966, r26967, r26968, r26969, r26970, r26971, r26972, r26973, r26974, r26975, r26976, r26977, r26978, r26979;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26959);
        mpfr_init(r26960);
        mpfr_init(r26961);
        mpfr_init(r26962);
        mpfr_init(r26963);
        mpfr_init_set_str(r26964, "2", 10, MPFR_RNDN);
        mpfr_init(r26965);
        mpfr_init_set_str(r26966, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r26967, "1", 10, MPFR_RNDN);
        mpfr_init(r26968);
        mpfr_init(r26969);
        mpfr_init_set_str(r26970, "3", 10, MPFR_RNDN);
        mpfr_init(r26971);
        mpfr_init(r26972);
        mpfr_init(r26973);
        mpfr_init(r26974);
        mpfr_init(r26975);
        mpfr_init(r26976);
        mpfr_init(r26977);
        mpfr_init(r26978);
        mpfr_init(r26979);
}

double f_fm(double a, double b) {
        mpfr_set_d(r26959, a, MPFR_RNDN);
        mpfr_mul(r26960, r26959, r26959, MPFR_RNDN);
        mpfr_set_d(r26961, b, MPFR_RNDN);
        mpfr_mul(r26962, r26961, r26961, MPFR_RNDN);
        mpfr_add(r26963, r26960, r26962, MPFR_RNDN);
        ;
        mpfr_pow(r26965, r26963, r26964, MPFR_RNDN);
        ;
        ;
        mpfr_add(r26968, r26967, r26959, MPFR_RNDN);
        mpfr_mul(r26969, r26960, r26968, MPFR_RNDN);
        ;
        mpfr_mul(r26971, r26970, r26959, MPFR_RNDN);
        mpfr_sub(r26972, r26967, r26971, MPFR_RNDN);
        mpfr_mul(r26973, r26962, r26972, MPFR_RNDN);
        mpfr_add(r26974, r26969, r26973, MPFR_RNDN);
        mpfr_mul(r26975, r26966, r26974, MPFR_RNDN);
        mpfr_add(r26976, r26965, r26975, MPFR_RNDN);
        mpfr_sqrt(r26977, r26976, MPFR_RNDN);
        mpfr_mul(r26978, r26977, r26977, MPFR_RNDN);
        mpfr_sub(r26979, r26978, r26967, MPFR_RNDN);
        return mpfr_get_d(r26979, MPFR_RNDN);
}

static mpfr_t r26980, r26981, r26982, r26983, r26984, r26985, r26986, r26987, r26988, r26989, r26990, r26991, r26992, r26993, r26994, r26995, r26996, r26997, r26998, r26999, r27000;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26980);
        mpfr_init(r26981);
        mpfr_init(r26982);
        mpfr_init(r26983);
        mpfr_init(r26984);
        mpfr_init_set_str(r26985, "2", 10, MPFR_RNDN);
        mpfr_init(r26986);
        mpfr_init_set_str(r26987, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r26988, "1", 10, MPFR_RNDN);
        mpfr_init(r26989);
        mpfr_init(r26990);
        mpfr_init_set_str(r26991, "3", 10, MPFR_RNDN);
        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(r27000);
}

double f_dm(double a, double b) {
        mpfr_set_d(r26980, a, MPFR_RNDN);
        mpfr_mul(r26981, r26980, r26980, MPFR_RNDN);
        mpfr_set_d(r26982, b, MPFR_RNDN);
        mpfr_mul(r26983, r26982, r26982, MPFR_RNDN);
        mpfr_add(r26984, r26981, r26983, MPFR_RNDN);
        ;
        mpfr_pow(r26986, r26984, r26985, MPFR_RNDN);
        ;
        ;
        mpfr_add(r26989, r26988, r26980, MPFR_RNDN);
        mpfr_mul(r26990, r26981, r26989, MPFR_RNDN);
        ;
        mpfr_mul(r26992, r26991, r26980, MPFR_RNDN);
        mpfr_sub(r26993, r26988, r26992, MPFR_RNDN);
        mpfr_mul(r26994, r26983, r26993, MPFR_RNDN);
        mpfr_add(r26995, r26990, r26994, MPFR_RNDN);
        mpfr_mul(r26996, r26987, r26995, MPFR_RNDN);
        mpfr_add(r26997, r26986, r26996, MPFR_RNDN);
        mpfr_sqrt(r26998, r26997, MPFR_RNDN);
        mpfr_mul(r26999, r26998, r26998, MPFR_RNDN);
        mpfr_sub(r27000, r26999, r26988, MPFR_RNDN);
        return mpfr_get_d(r27000, MPFR_RNDN);
}