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

char *name = "VandenBroeck and Keller, Equation (20)";

double f_if(float f) {
        float r20783 = 1;
        float r20784 = atan2(1.0, 0.0);
        float r20785 = 4;
        float r20786 = r20784 / r20785;
        float r20787 = r20783 / r20786;
        float r20788 = f;
        float r20789 = r20786 * r20788;
        float r20790 = exp(r20789);
        float r20791 = -r20789;
        float r20792 = exp(r20791);
        float r20793 = r20790 + r20792;
        float r20794 = r20790 - r20792;
        float r20795 = r20793 / r20794;
        float r20796 = log(r20795);
        float r20797 = r20787 * r20796;
        float r20798 = -r20797;
        return r20798;
}

double f_id(double f) {
        double r20799 = 1;
        double r20800 = atan2(1.0, 0.0);
        double r20801 = 4;
        double r20802 = r20800 / r20801;
        double r20803 = r20799 / r20802;
        double r20804 = f;
        double r20805 = r20802 * r20804;
        double r20806 = exp(r20805);
        double r20807 = -r20805;
        double r20808 = exp(r20807);
        double r20809 = r20806 + r20808;
        double r20810 = r20806 - r20808;
        double r20811 = r20809 / r20810;
        double r20812 = log(r20811);
        double r20813 = r20803 * r20812;
        double r20814 = -r20813;
        return r20814;
}


double f_of(float f) {
        float r20815 = 4;
        float r20816 = atan2(1.0, 0.0);
        float r20817 = r20815 / r20816;
        float r20818 = log(r20817);
        float r20819 = f;
        float r20820 = log(r20819);
        float r20821 = r20818 - r20820;
        float r20822 = r20817 * r20821;
        float r20823 = r20819 * r20819;
        float r20824 = 1/12;
        float r20825 = r20824 * r20816;
        float r20826 = r20823 * r20825;
        float r20827 = 7/5760;
        float r20828 = r20827 * r20816;
        float r20829 = pow(r20819, r20815);
        float r20830 = r20828 * r20829;
        float r20831 = r20816 * r20816;
        float r20832 = r20830 * r20831;
        float r20833 = r20826 - r20832;
        float r20834 = r20822 + r20833;
        float r20835 = -r20834;
        return r20835;
}

double f_od(double f) {
        double r20836 = 4;
        double r20837 = atan2(1.0, 0.0);
        double r20838 = r20836 / r20837;
        double r20839 = log(r20838);
        double r20840 = f;
        double r20841 = log(r20840);
        double r20842 = r20839 - r20841;
        double r20843 = r20838 * r20842;
        double r20844 = r20840 * r20840;
        double r20845 = 1/12;
        double r20846 = r20845 * r20837;
        double r20847 = r20844 * r20846;
        double r20848 = 7/5760;
        double r20849 = r20848 * r20837;
        double r20850 = pow(r20840, r20836);
        double r20851 = r20849 * r20850;
        double r20852 = r20837 * r20837;
        double r20853 = r20851 * r20852;
        double r20854 = r20847 - r20853;
        double r20855 = r20843 + r20854;
        double r20856 = -r20855;
        return r20856;
}

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 r20857, r20858, r20859, r20860, r20861, r20862, r20863, r20864, r20865, r20866, r20867, r20868, r20869, r20870, r20871, r20872;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r20857, "1", 10, MPFR_RNDN);
        mpfr_init(r20858);
        mpfr_init_set_str(r20859, "4", 10, MPFR_RNDN);
        mpfr_init(r20860);
        mpfr_init(r20861);
        mpfr_init(r20862);
        mpfr_init(r20863);
        mpfr_init(r20864);
        mpfr_init(r20865);
        mpfr_init(r20866);
        mpfr_init(r20867);
        mpfr_init(r20868);
        mpfr_init(r20869);
        mpfr_init(r20870);
        mpfr_init(r20871);
        mpfr_init(r20872);
}

double f_im(double f) {
        ;
        mpfr_const_pi(r20858, MPFR_RNDN);
        ;
        mpfr_div(r20860, r20858, r20859, MPFR_RNDN);
        mpfr_div(r20861, r20857, r20860, MPFR_RNDN);
        mpfr_set_d(r20862, f, MPFR_RNDN);
        mpfr_mul(r20863, r20860, r20862, MPFR_RNDN);
        mpfr_exp(r20864, r20863, MPFR_RNDN);
        mpfr_neg(r20865, r20863, MPFR_RNDN);
        mpfr_exp(r20866, r20865, MPFR_RNDN);
        mpfr_add(r20867, r20864, r20866, MPFR_RNDN);
        mpfr_sub(r20868, r20864, r20866, MPFR_RNDN);
        mpfr_div(r20869, r20867, r20868, MPFR_RNDN);
        mpfr_log(r20870, r20869, MPFR_RNDN);
        mpfr_mul(r20871, r20861, r20870, MPFR_RNDN);
        mpfr_neg(r20872, r20871, MPFR_RNDN);
        return mpfr_get_d(r20872, MPFR_RNDN);
}

static mpfr_t r20873, r20874, r20875, r20876, r20877, r20878, r20879, r20880, r20881, r20882, r20883, r20884, r20885, r20886, r20887, r20888, r20889, r20890, r20891, r20892, r20893;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r20873, "4", 10, MPFR_RNDN);
        mpfr_init(r20874);
        mpfr_init(r20875);
        mpfr_init(r20876);
        mpfr_init(r20877);
        mpfr_init(r20878);
        mpfr_init(r20879);
        mpfr_init(r20880);
        mpfr_init(r20881);
        mpfr_init_set_str(r20882, "1/12", 10, MPFR_RNDN);
        mpfr_init(r20883);
        mpfr_init(r20884);
        mpfr_init_set_str(r20885, "7/5760", 10, MPFR_RNDN);
        mpfr_init(r20886);
        mpfr_init(r20887);
        mpfr_init(r20888);
        mpfr_init(r20889);
        mpfr_init(r20890);
        mpfr_init(r20891);
        mpfr_init(r20892);
        mpfr_init(r20893);
}

double f_fm(double f) {
        ;
        mpfr_const_pi(r20874, MPFR_RNDN);
        mpfr_div(r20875, r20873, r20874, MPFR_RNDN);
        mpfr_log(r20876, r20875, MPFR_RNDN);
        mpfr_set_d(r20877, f, MPFR_RNDN);
        mpfr_log(r20878, r20877, MPFR_RNDN);
        mpfr_sub(r20879, r20876, r20878, MPFR_RNDN);
        mpfr_mul(r20880, r20875, r20879, MPFR_RNDN);
        mpfr_mul(r20881, r20877, r20877, MPFR_RNDN);
        ;
        mpfr_mul(r20883, r20882, r20874, MPFR_RNDN);
        mpfr_mul(r20884, r20881, r20883, MPFR_RNDN);
        ;
        mpfr_mul(r20886, r20885, r20874, MPFR_RNDN);
        mpfr_pow(r20887, r20877, r20873, MPFR_RNDN);
        mpfr_mul(r20888, r20886, r20887, MPFR_RNDN);
        mpfr_mul(r20889, r20874, r20874, MPFR_RNDN);
        mpfr_mul(r20890, r20888, r20889, MPFR_RNDN);
        mpfr_sub(r20891, r20884, r20890, MPFR_RNDN);
        mpfr_add(r20892, r20880, r20891, MPFR_RNDN);
        mpfr_neg(r20893, r20892, MPFR_RNDN);
        return mpfr_get_d(r20893, MPFR_RNDN);
}

static mpfr_t r20894, r20895, r20896, r20897, r20898, r20899, r20900, r20901, r20902, r20903, r20904, r20905, r20906, r20907, r20908, r20909, r20910, r20911, r20912, r20913, r20914;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r20894, "4", 10, MPFR_RNDN);
        mpfr_init(r20895);
        mpfr_init(r20896);
        mpfr_init(r20897);
        mpfr_init(r20898);
        mpfr_init(r20899);
        mpfr_init(r20900);
        mpfr_init(r20901);
        mpfr_init(r20902);
        mpfr_init_set_str(r20903, "1/12", 10, MPFR_RNDN);
        mpfr_init(r20904);
        mpfr_init(r20905);
        mpfr_init_set_str(r20906, "7/5760", 10, MPFR_RNDN);
        mpfr_init(r20907);
        mpfr_init(r20908);
        mpfr_init(r20909);
        mpfr_init(r20910);
        mpfr_init(r20911);
        mpfr_init(r20912);
        mpfr_init(r20913);
        mpfr_init(r20914);
}

double f_dm(double f) {
        ;
        mpfr_const_pi(r20895, MPFR_RNDN);
        mpfr_div(r20896, r20894, r20895, MPFR_RNDN);
        mpfr_log(r20897, r20896, MPFR_RNDN);
        mpfr_set_d(r20898, f, MPFR_RNDN);
        mpfr_log(r20899, r20898, MPFR_RNDN);
        mpfr_sub(r20900, r20897, r20899, MPFR_RNDN);
        mpfr_mul(r20901, r20896, r20900, MPFR_RNDN);
        mpfr_mul(r20902, r20898, r20898, MPFR_RNDN);
        ;
        mpfr_mul(r20904, r20903, r20895, MPFR_RNDN);
        mpfr_mul(r20905, r20902, r20904, MPFR_RNDN);
        ;
        mpfr_mul(r20907, r20906, r20895, MPFR_RNDN);
        mpfr_pow(r20908, r20898, r20894, MPFR_RNDN);
        mpfr_mul(r20909, r20907, r20908, MPFR_RNDN);
        mpfr_mul(r20910, r20895, r20895, MPFR_RNDN);
        mpfr_mul(r20911, r20909, r20910, MPFR_RNDN);
        mpfr_sub(r20912, r20905, r20911, MPFR_RNDN);
        mpfr_add(r20913, r20901, r20912, MPFR_RNDN);
        mpfr_neg(r20914, r20913, MPFR_RNDN);
        return mpfr_get_d(r20914, MPFR_RNDN);
}