#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Falkner and Boettcher, Equation (20:1,3)";

double f_if(float v, float t) {
        float r26857 = 1;
        float r26858 = 5;
        float r26859 = v;
        float r26860 = r26859 * r26859;
        float r26861 = r26858 * r26860;
        float r26862 = r26857 - r26861;
        float r26863 = atan2(1.0, 0.0);
        float r26864 = t;
        float r26865 = r26863 * r26864;
        float r26866 = 2;
        float r26867 = 3;
        float r26868 = r26867 * r26860;
        float r26869 = r26857 - r26868;
        float r26870 = r26866 * r26869;
        float r26871 = sqrt(r26870);
        float r26872 = r26865 * r26871;
        float r26873 = r26857 - r26860;
        float r26874 = r26872 * r26873;
        float r26875 = r26862 / r26874;
        return r26875;
}

double f_id(double v, double t) {
        double r26876 = 1;
        double r26877 = 5;
        double r26878 = v;
        double r26879 = r26878 * r26878;
        double r26880 = r26877 * r26879;
        double r26881 = r26876 - r26880;
        double r26882 = atan2(1.0, 0.0);
        double r26883 = t;
        double r26884 = r26882 * r26883;
        double r26885 = 2;
        double r26886 = 3;
        double r26887 = r26886 * r26879;
        double r26888 = r26876 - r26887;
        double r26889 = r26885 * r26888;
        double r26890 = sqrt(r26889);
        double r26891 = r26884 * r26890;
        double r26892 = r26876 - r26879;
        double r26893 = r26891 * r26892;
        double r26894 = r26881 / r26893;
        return r26894;
}


double f_of(float v, float t) {
        float r26895 = 5;
        float r26896 = -r26895;
        float r26897 = v;
        float r26898 = r26897 * r26897;
        float r26899 = 1;
        float r26900 = fma(r26896, r26898, r26899);
        float r26901 = 3;
        float r26902 = -r26901;
        float r26903 = r26898 * r26902;
        float r26904 = r26898 * r26901;
        float r26905 = fma(r26903, r26904, r26899);
        float r26906 = 2;
        float r26907 = r26905 * r26906;
        float r26908 = sqrt(r26907);
        float r26909 = r26900 / r26908;
        float r26910 = pow(r26897, r26901);
        float r26911 = -r26897;
        float r26912 = fma(r26910, r26911, r26899);
        float r26913 = t;
        float r26914 = atan2(1.0, 0.0);
        float r26915 = r26913 * r26914;
        float r26916 = r26912 * r26915;
        float r26917 = r26909 / r26916;
        float r26918 = r26901 * r26898;
        float r26919 = r26899 + r26918;
        float r26920 = sqrt(r26919);
        float r26921 = r26899 + r26898;
        float r26922 = r26920 * r26921;
        float r26923 = r26917 * r26922;
        return r26923;
}

double f_od(double v, double t) {
        double r26924 = 5;
        double r26925 = -r26924;
        double r26926 = v;
        double r26927 = r26926 * r26926;
        double r26928 = 1;
        double r26929 = fma(r26925, r26927, r26928);
        double r26930 = 3;
        double r26931 = -r26930;
        double r26932 = r26927 * r26931;
        double r26933 = r26927 * r26930;
        double r26934 = fma(r26932, r26933, r26928);
        double r26935 = 2;
        double r26936 = r26934 * r26935;
        double r26937 = sqrt(r26936);
        double r26938 = r26929 / r26937;
        double r26939 = pow(r26926, r26930);
        double r26940 = -r26926;
        double r26941 = fma(r26939, r26940, r26928);
        double r26942 = t;
        double r26943 = atan2(1.0, 0.0);
        double r26944 = r26942 * r26943;
        double r26945 = r26941 * r26944;
        double r26946 = r26938 / r26945;
        double r26947 = r26930 * r26927;
        double r26948 = r26928 + r26947;
        double r26949 = sqrt(r26948);
        double r26950 = r26928 + r26927;
        double r26951 = r26949 * r26950;
        double r26952 = r26946 * r26951;
        return r26952;
}

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 r26953, r26954, r26955, r26956, r26957, r26958, r26959, r26960, r26961, r26962, r26963, r26964, r26965, r26966, r26967, r26968, r26969, r26970, r26971;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r26953, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r26954, "5", 10, MPFR_RNDN);
        mpfr_init(r26955);
        mpfr_init(r26956);
        mpfr_init(r26957);
        mpfr_init(r26958);
        mpfr_init(r26959);
        mpfr_init(r26960);
        mpfr_init(r26961);
        mpfr_init_set_str(r26962, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r26963, "3", 10, MPFR_RNDN);
        mpfr_init(r26964);
        mpfr_init(r26965);
        mpfr_init(r26966);
        mpfr_init(r26967);
        mpfr_init(r26968);
        mpfr_init(r26969);
        mpfr_init(r26970);
        mpfr_init(r26971);
}

double f_im(double v, double t) {
        ;
        ;
        mpfr_set_d(r26955, v, MPFR_RNDN);
        mpfr_mul(r26956, r26955, r26955, MPFR_RNDN);
        mpfr_mul(r26957, r26954, r26956, MPFR_RNDN);
        mpfr_sub(r26958, r26953, r26957, MPFR_RNDN);
        mpfr_const_pi(r26959, MPFR_RNDN);
        mpfr_set_d(r26960, t, MPFR_RNDN);
        mpfr_mul(r26961, r26959, r26960, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r26964, r26963, r26956, MPFR_RNDN);
        mpfr_sub(r26965, r26953, r26964, MPFR_RNDN);
        mpfr_mul(r26966, r26962, r26965, MPFR_RNDN);
        mpfr_sqrt(r26967, r26966, MPFR_RNDN);
        mpfr_mul(r26968, r26961, r26967, MPFR_RNDN);
        mpfr_sub(r26969, r26953, r26956, MPFR_RNDN);
        mpfr_mul(r26970, r26968, r26969, MPFR_RNDN);
        mpfr_div(r26971, r26958, r26970, MPFR_RNDN);
        return mpfr_get_d(r26971, MPFR_RNDN);
}

static mpfr_t r26972, r26973, r26974, r26975, r26976, r26977, r26978, r26979, 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_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r26972, "5", 10, MPFR_RNDN);
        mpfr_init(r26973);
        mpfr_init(r26974);
        mpfr_init(r26975);
        mpfr_init_set_str(r26976, "1", 10, MPFR_RNDN);
        mpfr_init(r26977);
        mpfr_init_set_str(r26978, "3", 10, MPFR_RNDN);
        mpfr_init(r26979);
        mpfr_init(r26980);
        mpfr_init(r26981);
        mpfr_init(r26982);
        mpfr_init_set_str(r26983, "2", 10, MPFR_RNDN);
        mpfr_init(r26984);
        mpfr_init(r26985);
        mpfr_init(r26986);
        mpfr_init(r26987);
        mpfr_init(r26988);
        mpfr_init(r26989);
        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(r27000);
}

double f_fm(double v, double t) {
        ;
        mpfr_neg(r26973, r26972, MPFR_RNDN);
        mpfr_set_d(r26974, v, MPFR_RNDN);
        mpfr_mul(r26975, r26974, r26974, MPFR_RNDN);
        ;
        mpfr_fma(r26977, r26973, r26975, r26976, MPFR_RNDN);
        ;
        mpfr_neg(r26979, r26978, MPFR_RNDN);
        mpfr_mul(r26980, r26975, r26979, MPFR_RNDN);
        mpfr_mul(r26981, r26975, r26978, MPFR_RNDN);
        mpfr_fma(r26982, r26980, r26981, r26976, MPFR_RNDN);
        ;
        mpfr_mul(r26984, r26982, r26983, MPFR_RNDN);
        mpfr_sqrt(r26985, r26984, MPFR_RNDN);
        mpfr_div(r26986, r26977, r26985, MPFR_RNDN);
        mpfr_pow(r26987, r26974, r26978, MPFR_RNDN);
        mpfr_neg(r26988, r26974, MPFR_RNDN);
        mpfr_fma(r26989, r26987, r26988, r26976, MPFR_RNDN);
        mpfr_set_d(r26990, t, MPFR_RNDN);
        mpfr_const_pi(r26991, MPFR_RNDN);
        mpfr_mul(r26992, r26990, r26991, MPFR_RNDN);
        mpfr_mul(r26993, r26989, r26992, MPFR_RNDN);
        mpfr_div(r26994, r26986, r26993, MPFR_RNDN);
        mpfr_mul(r26995, r26978, r26975, MPFR_RNDN);
        mpfr_add(r26996, r26976, r26995, MPFR_RNDN);
        mpfr_sqrt(r26997, r26996, MPFR_RNDN);
        mpfr_add(r26998, r26976, r26975, MPFR_RNDN);
        mpfr_mul(r26999, r26997, r26998, MPFR_RNDN);
        mpfr_mul(r27000, r26994, r26999, MPFR_RNDN);
        return mpfr_get_d(r27000, MPFR_RNDN);
}

static mpfr_t r27001, r27002, r27003, r27004, r27005, r27006, r27007, r27008, r27009, r27010, r27011, r27012, r27013, r27014, r27015, r27016, r27017, r27018, r27019, r27020, r27021, r27022, r27023, r27024, r27025, r27026, r27027, r27028, r27029;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27001, "5", 10, MPFR_RNDN);
        mpfr_init(r27002);
        mpfr_init(r27003);
        mpfr_init(r27004);
        mpfr_init_set_str(r27005, "1", 10, MPFR_RNDN);
        mpfr_init(r27006);
        mpfr_init_set_str(r27007, "3", 10, MPFR_RNDN);
        mpfr_init(r27008);
        mpfr_init(r27009);
        mpfr_init(r27010);
        mpfr_init(r27011);
        mpfr_init_set_str(r27012, "2", 10, MPFR_RNDN);
        mpfr_init(r27013);
        mpfr_init(r27014);
        mpfr_init(r27015);
        mpfr_init(r27016);
        mpfr_init(r27017);
        mpfr_init(r27018);
        mpfr_init(r27019);
        mpfr_init(r27020);
        mpfr_init(r27021);
        mpfr_init(r27022);
        mpfr_init(r27023);
        mpfr_init(r27024);
        mpfr_init(r27025);
        mpfr_init(r27026);
        mpfr_init(r27027);
        mpfr_init(r27028);
        mpfr_init(r27029);
}

double f_dm(double v, double t) {
        ;
        mpfr_neg(r27002, r27001, MPFR_RNDN);
        mpfr_set_d(r27003, v, MPFR_RNDN);
        mpfr_mul(r27004, r27003, r27003, MPFR_RNDN);
        ;
        mpfr_fma(r27006, r27002, r27004, r27005, MPFR_RNDN);
        ;
        mpfr_neg(r27008, r27007, MPFR_RNDN);
        mpfr_mul(r27009, r27004, r27008, MPFR_RNDN);
        mpfr_mul(r27010, r27004, r27007, MPFR_RNDN);
        mpfr_fma(r27011, r27009, r27010, r27005, MPFR_RNDN);
        ;
        mpfr_mul(r27013, r27011, r27012, MPFR_RNDN);
        mpfr_sqrt(r27014, r27013, MPFR_RNDN);
        mpfr_div(r27015, r27006, r27014, MPFR_RNDN);
        mpfr_pow(r27016, r27003, r27007, MPFR_RNDN);
        mpfr_neg(r27017, r27003, MPFR_RNDN);
        mpfr_fma(r27018, r27016, r27017, r27005, MPFR_RNDN);
        mpfr_set_d(r27019, t, MPFR_RNDN);
        mpfr_const_pi(r27020, MPFR_RNDN);
        mpfr_mul(r27021, r27019, r27020, MPFR_RNDN);
        mpfr_mul(r27022, r27018, r27021, MPFR_RNDN);
        mpfr_div(r27023, r27015, r27022, MPFR_RNDN);
        mpfr_mul(r27024, r27007, r27004, MPFR_RNDN);
        mpfr_add(r27025, r27005, r27024, MPFR_RNDN);
        mpfr_sqrt(r27026, r27025, MPFR_RNDN);
        mpfr_add(r27027, r27005, r27004, MPFR_RNDN);
        mpfr_mul(r27028, r27026, r27027, MPFR_RNDN);
        mpfr_mul(r27029, r27023, r27028, MPFR_RNDN);
        return mpfr_get_d(r27029, MPFR_RNDN);
}