#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 r26710 = 1;
        float r26711 = 5;
        float r26712 = v;
        float r26713 = r26712 * r26712;
        float r26714 = r26711 * r26713;
        float r26715 = r26710 - r26714;
        float r26716 = atan2(1.0, 0.0);
        float r26717 = t;
        float r26718 = r26716 * r26717;
        float r26719 = 2;
        float r26720 = 3;
        float r26721 = r26720 * r26713;
        float r26722 = r26710 - r26721;
        float r26723 = r26719 * r26722;
        float r26724 = sqrt(r26723);
        float r26725 = r26718 * r26724;
        float r26726 = r26710 - r26713;
        float r26727 = r26725 * r26726;
        float r26728 = r26715 / r26727;
        return r26728;
}

double f_id(double v, double t) {
        double r26729 = 1;
        double r26730 = 5;
        double r26731 = v;
        double r26732 = r26731 * r26731;
        double r26733 = r26730 * r26732;
        double r26734 = r26729 - r26733;
        double r26735 = atan2(1.0, 0.0);
        double r26736 = t;
        double r26737 = r26735 * r26736;
        double r26738 = 2;
        double r26739 = 3;
        double r26740 = r26739 * r26732;
        double r26741 = r26729 - r26740;
        double r26742 = r26738 * r26741;
        double r26743 = sqrt(r26742);
        double r26744 = r26737 * r26743;
        double r26745 = r26729 - r26732;
        double r26746 = r26744 * r26745;
        double r26747 = r26734 / r26746;
        return r26747;
}


double f_of(float v, float t) {
        float r26748 = 1;
        float r26749 = v;
        float r26750 = 5;
        float r26751 = r26749 * r26750;
        float r26752 = r26749 * r26751;
        float r26753 = r26748 - r26752;
        float r26754 = t;
        float r26755 = r26753 / r26754;
        float r26756 = atan2(1.0, 0.0);
        float r26757 = r26749 * r26749;
        float r26758 = r26757 * r26757;
        float r26759 = r26748 - r26758;
        float r26760 = r26756 * r26759;
        float r26761 = 2;
        float r26762 = 3;
        float r26763 = r26762 * r26757;
        float r26764 = r26748 - r26763;
        float r26765 = r26761 * r26764;
        float r26766 = sqrt(r26765);
        float r26767 = r26760 * r26766;
        float r26768 = r26755 / r26767;
        float r26769 = r26748 + r26757;
        float r26770 = r26768 * r26769;
        return r26770;
}

double f_od(double v, double t) {
        double r26771 = 1;
        double r26772 = v;
        double r26773 = 5;
        double r26774 = r26772 * r26773;
        double r26775 = r26772 * r26774;
        double r26776 = r26771 - r26775;
        double r26777 = t;
        double r26778 = r26776 / r26777;
        double r26779 = atan2(1.0, 0.0);
        double r26780 = r26772 * r26772;
        double r26781 = r26780 * r26780;
        double r26782 = r26771 - r26781;
        double r26783 = r26779 * r26782;
        double r26784 = 2;
        double r26785 = 3;
        double r26786 = r26785 * r26780;
        double r26787 = r26771 - r26786;
        double r26788 = r26784 * r26787;
        double r26789 = sqrt(r26788);
        double r26790 = r26783 * r26789;
        double r26791 = r26778 / r26790;
        double r26792 = r26771 + r26780;
        double r26793 = r26791 * r26792;
        return r26793;
}

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 r26794, r26795, r26796, r26797, r26798, r26799, r26800, r26801, r26802, r26803, r26804, r26805, r26806, r26807, r26808, r26809, r26810, r26811, r26812;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r26794, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r26795, "5", 10, MPFR_RNDN);
        mpfr_init(r26796);
        mpfr_init(r26797);
        mpfr_init(r26798);
        mpfr_init(r26799);
        mpfr_init(r26800);
        mpfr_init(r26801);
        mpfr_init(r26802);
        mpfr_init_set_str(r26803, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r26804, "3", 10, MPFR_RNDN);
        mpfr_init(r26805);
        mpfr_init(r26806);
        mpfr_init(r26807);
        mpfr_init(r26808);
        mpfr_init(r26809);
        mpfr_init(r26810);
        mpfr_init(r26811);
        mpfr_init(r26812);
}

double f_im(double v, double t) {
        ;
        ;
        mpfr_set_d(r26796, v, MPFR_RNDN);
        mpfr_mul(r26797, r26796, r26796, MPFR_RNDN);
        mpfr_mul(r26798, r26795, r26797, MPFR_RNDN);
        mpfr_sub(r26799, r26794, r26798, MPFR_RNDN);
        mpfr_const_pi(r26800, MPFR_RNDN);
        mpfr_set_d(r26801, t, MPFR_RNDN);
        mpfr_mul(r26802, r26800, r26801, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r26805, r26804, r26797, MPFR_RNDN);
        mpfr_sub(r26806, r26794, r26805, MPFR_RNDN);
        mpfr_mul(r26807, r26803, r26806, MPFR_RNDN);
        mpfr_sqrt(r26808, r26807, MPFR_RNDN);
        mpfr_mul(r26809, r26802, r26808, MPFR_RNDN);
        mpfr_sub(r26810, r26794, r26797, MPFR_RNDN);
        mpfr_mul(r26811, r26809, r26810, MPFR_RNDN);
        mpfr_div(r26812, r26799, r26811, MPFR_RNDN);
        return mpfr_get_d(r26812, MPFR_RNDN);
}

static mpfr_t r26813, r26814, r26815, r26816, r26817, r26818, r26819, r26820, r26821, r26822, r26823, r26824, r26825, r26826, r26827, r26828, r26829, r26830, r26831, r26832, r26833, r26834, r26835;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r26813, "1", 10, MPFR_RNDN);
        mpfr_init(r26814);
        mpfr_init_set_str(r26815, "5", 10, MPFR_RNDN);
        mpfr_init(r26816);
        mpfr_init(r26817);
        mpfr_init(r26818);
        mpfr_init(r26819);
        mpfr_init(r26820);
        mpfr_init(r26821);
        mpfr_init(r26822);
        mpfr_init(r26823);
        mpfr_init(r26824);
        mpfr_init(r26825);
        mpfr_init_set_str(r26826, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r26827, "3", 10, MPFR_RNDN);
        mpfr_init(r26828);
        mpfr_init(r26829);
        mpfr_init(r26830);
        mpfr_init(r26831);
        mpfr_init(r26832);
        mpfr_init(r26833);
        mpfr_init(r26834);
        mpfr_init(r26835);
}

double f_fm(double v, double t) {
        ;
        mpfr_set_d(r26814, v, MPFR_RNDN);
        ;
        mpfr_mul(r26816, r26814, r26815, MPFR_RNDN);
        mpfr_mul(r26817, r26814, r26816, MPFR_RNDN);
        mpfr_sub(r26818, r26813, r26817, MPFR_RNDN);
        mpfr_set_d(r26819, t, MPFR_RNDN);
        mpfr_div(r26820, r26818, r26819, MPFR_RNDN);
        mpfr_const_pi(r26821, MPFR_RNDN);
        mpfr_mul(r26822, r26814, r26814, MPFR_RNDN);
        mpfr_mul(r26823, r26822, r26822, MPFR_RNDN);
        mpfr_sub(r26824, r26813, r26823, MPFR_RNDN);
        mpfr_mul(r26825, r26821, r26824, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r26828, r26827, r26822, MPFR_RNDN);
        mpfr_sub(r26829, r26813, r26828, MPFR_RNDN);
        mpfr_mul(r26830, r26826, r26829, MPFR_RNDN);
        mpfr_sqrt(r26831, r26830, MPFR_RNDN);
        mpfr_mul(r26832, r26825, r26831, MPFR_RNDN);
        mpfr_div(r26833, r26820, r26832, MPFR_RNDN);
        mpfr_add(r26834, r26813, r26822, MPFR_RNDN);
        mpfr_mul(r26835, r26833, r26834, MPFR_RNDN);
        return mpfr_get_d(r26835, MPFR_RNDN);
}

static mpfr_t r26836, r26837, r26838, r26839, r26840, r26841, r26842, r26843, r26844, r26845, r26846, r26847, r26848, r26849, r26850, r26851, r26852, r26853, r26854, r26855, r26856, r26857, r26858;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r26836, "1", 10, MPFR_RNDN);
        mpfr_init(r26837);
        mpfr_init_set_str(r26838, "5", 10, MPFR_RNDN);
        mpfr_init(r26839);
        mpfr_init(r26840);
        mpfr_init(r26841);
        mpfr_init(r26842);
        mpfr_init(r26843);
        mpfr_init(r26844);
        mpfr_init(r26845);
        mpfr_init(r26846);
        mpfr_init(r26847);
        mpfr_init(r26848);
        mpfr_init_set_str(r26849, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r26850, "3", 10, MPFR_RNDN);
        mpfr_init(r26851);
        mpfr_init(r26852);
        mpfr_init(r26853);
        mpfr_init(r26854);
        mpfr_init(r26855);
        mpfr_init(r26856);
        mpfr_init(r26857);
        mpfr_init(r26858);
}

double f_dm(double v, double t) {
        ;
        mpfr_set_d(r26837, v, MPFR_RNDN);
        ;
        mpfr_mul(r26839, r26837, r26838, MPFR_RNDN);
        mpfr_mul(r26840, r26837, r26839, MPFR_RNDN);
        mpfr_sub(r26841, r26836, r26840, MPFR_RNDN);
        mpfr_set_d(r26842, t, MPFR_RNDN);
        mpfr_div(r26843, r26841, r26842, MPFR_RNDN);
        mpfr_const_pi(r26844, MPFR_RNDN);
        mpfr_mul(r26845, r26837, r26837, MPFR_RNDN);
        mpfr_mul(r26846, r26845, r26845, MPFR_RNDN);
        mpfr_sub(r26847, r26836, r26846, MPFR_RNDN);
        mpfr_mul(r26848, r26844, r26847, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r26851, r26850, r26845, MPFR_RNDN);
        mpfr_sub(r26852, r26836, r26851, MPFR_RNDN);
        mpfr_mul(r26853, r26849, r26852, MPFR_RNDN);
        mpfr_sqrt(r26854, r26853, MPFR_RNDN);
        mpfr_mul(r26855, r26848, r26854, MPFR_RNDN);
        mpfr_div(r26856, r26843, r26855, MPFR_RNDN);
        mpfr_add(r26857, r26836, r26845, MPFR_RNDN);
        mpfr_mul(r26858, r26856, r26857, MPFR_RNDN);
        return mpfr_get_d(r26858, MPFR_RNDN);
}