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

char *name = "Toniolo and Linder, Equation (7)";

double f_if(float x, float l, float t) {
        float r21684 = 2;
        float r21685 = sqrt(r21684);
        float r21686 = t;
        float r21687 = r21685 * r21686;
        float r21688 = x;
        float r21689 = 1;
        float r21690 = r21688 + r21689;
        float r21691 = r21688 - r21689;
        float r21692 = r21690 / r21691;
        float r21693 = l;
        float r21694 = r21693 * r21693;
        float r21695 = r21686 * r21686;
        float r21696 = r21684 * r21695;
        float r21697 = r21694 + r21696;
        float r21698 = r21692 * r21697;
        float r21699 = r21698 - r21694;
        float r21700 = sqrt(r21699);
        float r21701 = r21687 / r21700;
        return r21701;
}

double f_id(double x, double l, double t) {
        double r21702 = 2;
        double r21703 = sqrt(r21702);
        double r21704 = t;
        double r21705 = r21703 * r21704;
        double r21706 = x;
        double r21707 = 1;
        double r21708 = r21706 + r21707;
        double r21709 = r21706 - r21707;
        double r21710 = r21708 / r21709;
        double r21711 = l;
        double r21712 = r21711 * r21711;
        double r21713 = r21704 * r21704;
        double r21714 = r21702 * r21713;
        double r21715 = r21712 + r21714;
        double r21716 = r21710 * r21715;
        double r21717 = r21716 - r21712;
        double r21718 = sqrt(r21717);
        double r21719 = r21705 / r21718;
        return r21719;
}


double f_of(float x, float l, float t) {
        float r21720 = t;
        float r21721 = -6.447314359762627e+120;
        bool r21722 = r21720 <= r21721;
        float r21723 = 2;
        float r21724 = sqrt(r21723);
        float r21725 = r21720 * r21724;
        float r21726 = r21720 / r21724;
        float r21727 = x;
        float r21728 = r21727 * r21727;
        float r21729 = r21726 / r21728;
        float r21730 = 1;
        float r21731 = r21730 - r21723;
        float r21732 = r21729 * r21731;
        float r21733 = r21723 / r21727;
        float r21734 = r21733 / r21724;
        float r21735 = r21724 + r21734;
        float r21736 = r21720 * r21735;
        float r21737 = r21732 - r21736;
        float r21738 = r21725 / r21737;
        float r21739 = 2.2327750609834536e-213;
        bool r21740 = r21720 <= r21739;
        float r21741 = sqrt(r21724);
        float r21742 = r21720 * r21741;
        float r21743 = r21742 * r21741;
        float r21744 = 4;
        float r21745 = r21744 / r21727;
        float r21746 = r21745 + r21723;
        float r21747 = r21720 * r21720;
        float r21748 = r21746 * r21747;
        float r21749 = l;
        float r21750 = r21723 * r21749;
        float r21751 = r21727 / r21749;
        float r21752 = r21750 / r21751;
        float r21753 = r21748 + r21752;
        float r21754 = sqrt(r21753);
        float r21755 = r21743 / r21754;
        float r21756 = 1.7113149206191654e-132;
        bool r21757 = r21720 <= r21756;
        float r21758 = r21728 * r21724;
        float r21759 = r21720 / r21758;
        float r21760 = r21723 - r21730;
        float r21761 = r21759 * r21760;
        float r21762 = r21736 + r21761;
        float r21763 = r21725 / r21762;
        float r21764 = 4.776380960244624e+22;
        bool r21765 = r21720 <= r21764;
        float r21766 = sqrt(r21754);
        float r21767 = r21766 * r21766;
        float r21768 = r21725 / r21767;
        float r21769 = r21765 ? r21768 : r21763;
        float r21770 = r21757 ? r21763 : r21769;
        float r21771 = r21740 ? r21755 : r21770;
        float r21772 = r21722 ? r21738 : r21771;
        return r21772;
}

double f_od(double x, double l, double t) {
        double r21773 = t;
        double r21774 = -6.447314359762627e+120;
        bool r21775 = r21773 <= r21774;
        double r21776 = 2;
        double r21777 = sqrt(r21776);
        double r21778 = r21773 * r21777;
        double r21779 = r21773 / r21777;
        double r21780 = x;
        double r21781 = r21780 * r21780;
        double r21782 = r21779 / r21781;
        double r21783 = 1;
        double r21784 = r21783 - r21776;
        double r21785 = r21782 * r21784;
        double r21786 = r21776 / r21780;
        double r21787 = r21786 / r21777;
        double r21788 = r21777 + r21787;
        double r21789 = r21773 * r21788;
        double r21790 = r21785 - r21789;
        double r21791 = r21778 / r21790;
        double r21792 = 2.2327750609834536e-213;
        bool r21793 = r21773 <= r21792;
        double r21794 = sqrt(r21777);
        double r21795 = r21773 * r21794;
        double r21796 = r21795 * r21794;
        double r21797 = 4;
        double r21798 = r21797 / r21780;
        double r21799 = r21798 + r21776;
        double r21800 = r21773 * r21773;
        double r21801 = r21799 * r21800;
        double r21802 = l;
        double r21803 = r21776 * r21802;
        double r21804 = r21780 / r21802;
        double r21805 = r21803 / r21804;
        double r21806 = r21801 + r21805;
        double r21807 = sqrt(r21806);
        double r21808 = r21796 / r21807;
        double r21809 = 1.7113149206191654e-132;
        bool r21810 = r21773 <= r21809;
        double r21811 = r21781 * r21777;
        double r21812 = r21773 / r21811;
        double r21813 = r21776 - r21783;
        double r21814 = r21812 * r21813;
        double r21815 = r21789 + r21814;
        double r21816 = r21778 / r21815;
        double r21817 = 4.776380960244624e+22;
        bool r21818 = r21773 <= r21817;
        double r21819 = sqrt(r21807);
        double r21820 = r21819 * r21819;
        double r21821 = r21778 / r21820;
        double r21822 = r21818 ? r21821 : r21816;
        double r21823 = r21810 ? r21816 : r21822;
        double r21824 = r21793 ? r21808 : r21823;
        double r21825 = r21775 ? r21791 : r21824;
        return r21825;
}

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 r21826, r21827, r21828, r21829, r21830, r21831, r21832, r21833, r21834, r21835, r21836, r21837, r21838, r21839, r21840, r21841, r21842, r21843;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r21826, "2", 10, MPFR_RNDN);
        mpfr_init(r21827);
        mpfr_init(r21828);
        mpfr_init(r21829);
        mpfr_init(r21830);
        mpfr_init_set_str(r21831, "1", 10, MPFR_RNDN);
        mpfr_init(r21832);
        mpfr_init(r21833);
        mpfr_init(r21834);
        mpfr_init(r21835);
        mpfr_init(r21836);
        mpfr_init(r21837);
        mpfr_init(r21838);
        mpfr_init(r21839);
        mpfr_init(r21840);
        mpfr_init(r21841);
        mpfr_init(r21842);
        mpfr_init(r21843);
}

double f_im(double x, double l, double t) {
        ;
        mpfr_sqrt(r21827, r21826, MPFR_RNDN);
        mpfr_set_d(r21828, t, MPFR_RNDN);
        mpfr_mul(r21829, r21827, r21828, MPFR_RNDN);
        mpfr_set_d(r21830, x, MPFR_RNDN);
        ;
        mpfr_add(r21832, r21830, r21831, MPFR_RNDN);
        mpfr_sub(r21833, r21830, r21831, MPFR_RNDN);
        mpfr_div(r21834, r21832, r21833, MPFR_RNDN);
        mpfr_set_d(r21835, l, MPFR_RNDN);
        mpfr_mul(r21836, r21835, r21835, MPFR_RNDN);
        mpfr_mul(r21837, r21828, r21828, MPFR_RNDN);
        mpfr_mul(r21838, r21826, r21837, MPFR_RNDN);
        mpfr_add(r21839, r21836, r21838, MPFR_RNDN);
        mpfr_mul(r21840, r21834, r21839, MPFR_RNDN);
        mpfr_sub(r21841, r21840, r21836, MPFR_RNDN);
        mpfr_sqrt(r21842, r21841, MPFR_RNDN);
        mpfr_div(r21843, r21829, r21842, MPFR_RNDN);
        return mpfr_get_d(r21843, MPFR_RNDN);
}

static mpfr_t r21844, r21845, r21846, r21847, r21848, r21849, r21850, r21851, r21852, r21853, r21854, r21855, r21856, r21857, r21858, r21859, r21860, r21861, r21862, r21863, r21864, r21865, r21866, r21867, r21868, r21869, r21870, r21871, r21872, r21873, r21874, r21875, r21876, r21877, r21878, r21879, r21880, r21881, r21882, r21883, r21884, r21885, r21886, r21887, r21888, r21889, r21890, r21891, r21892, r21893, r21894, r21895, r21896;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r21844);
        mpfr_init_set_str(r21845, "-6.447314359762627e+120", 10, MPFR_RNDN);
        mpfr_init(r21846);
        mpfr_init_set_str(r21847, "2", 10, MPFR_RNDN);
        mpfr_init(r21848);
        mpfr_init(r21849);
        mpfr_init(r21850);
        mpfr_init(r21851);
        mpfr_init(r21852);
        mpfr_init(r21853);
        mpfr_init_set_str(r21854, "1", 10, MPFR_RNDN);
        mpfr_init(r21855);
        mpfr_init(r21856);
        mpfr_init(r21857);
        mpfr_init(r21858);
        mpfr_init(r21859);
        mpfr_init(r21860);
        mpfr_init(r21861);
        mpfr_init(r21862);
        mpfr_init_set_str(r21863, "2.2327750609834536e-213", 10, MPFR_RNDN);
        mpfr_init(r21864);
        mpfr_init(r21865);
        mpfr_init(r21866);
        mpfr_init(r21867);
        mpfr_init_set_str(r21868, "4", 10, MPFR_RNDN);
        mpfr_init(r21869);
        mpfr_init(r21870);
        mpfr_init(r21871);
        mpfr_init(r21872);
        mpfr_init(r21873);
        mpfr_init(r21874);
        mpfr_init(r21875);
        mpfr_init(r21876);
        mpfr_init(r21877);
        mpfr_init(r21878);
        mpfr_init(r21879);
        mpfr_init_set_str(r21880, "1.7113149206191654e-132", 10, MPFR_RNDN);
        mpfr_init(r21881);
        mpfr_init(r21882);
        mpfr_init(r21883);
        mpfr_init(r21884);
        mpfr_init(r21885);
        mpfr_init(r21886);
        mpfr_init(r21887);
        mpfr_init_set_str(r21888, "4.776380960244624e+22", 10, MPFR_RNDN);
        mpfr_init(r21889);
        mpfr_init(r21890);
        mpfr_init(r21891);
        mpfr_init(r21892);
        mpfr_init(r21893);
        mpfr_init(r21894);
        mpfr_init(r21895);
        mpfr_init(r21896);
}

double f_fm(double x, double l, double t) {
        mpfr_set_d(r21844, t, MPFR_RNDN);
        ;
        mpfr_set_si(r21846, mpfr_cmp(r21844, r21845) <= 0, MPFR_RNDN);
        ;
        mpfr_sqrt(r21848, r21847, MPFR_RNDN);
        mpfr_mul(r21849, r21844, r21848, MPFR_RNDN);
        mpfr_div(r21850, r21844, r21848, MPFR_RNDN);
        mpfr_set_d(r21851, x, MPFR_RNDN);
        mpfr_mul(r21852, r21851, r21851, MPFR_RNDN);
        mpfr_div(r21853, r21850, r21852, MPFR_RNDN);
        ;
        mpfr_sub(r21855, r21854, r21847, MPFR_RNDN);
        mpfr_mul(r21856, r21853, r21855, MPFR_RNDN);
        mpfr_div(r21857, r21847, r21851, MPFR_RNDN);
        mpfr_div(r21858, r21857, r21848, MPFR_RNDN);
        mpfr_add(r21859, r21848, r21858, MPFR_RNDN);
        mpfr_mul(r21860, r21844, r21859, MPFR_RNDN);
        mpfr_sub(r21861, r21856, r21860, MPFR_RNDN);
        mpfr_div(r21862, r21849, r21861, MPFR_RNDN);
        ;
        mpfr_set_si(r21864, mpfr_cmp(r21844, r21863) <= 0, MPFR_RNDN);
        mpfr_sqrt(r21865, r21848, MPFR_RNDN);
        mpfr_mul(r21866, r21844, r21865, MPFR_RNDN);
        mpfr_mul(r21867, r21866, r21865, MPFR_RNDN);
        ;
        mpfr_div(r21869, r21868, r21851, MPFR_RNDN);
        mpfr_add(r21870, r21869, r21847, MPFR_RNDN);
        mpfr_mul(r21871, r21844, r21844, MPFR_RNDN);
        mpfr_mul(r21872, r21870, r21871, MPFR_RNDN);
        mpfr_set_d(r21873, l, MPFR_RNDN);
        mpfr_mul(r21874, r21847, r21873, MPFR_RNDN);
        mpfr_div(r21875, r21851, r21873, MPFR_RNDN);
        mpfr_div(r21876, r21874, r21875, MPFR_RNDN);
        mpfr_add(r21877, r21872, r21876, MPFR_RNDN);
        mpfr_sqrt(r21878, r21877, MPFR_RNDN);
        mpfr_div(r21879, r21867, r21878, MPFR_RNDN);
        ;
        mpfr_set_si(r21881, mpfr_cmp(r21844, r21880) <= 0, MPFR_RNDN);
        mpfr_mul(r21882, r21852, r21848, MPFR_RNDN);
        mpfr_div(r21883, r21844, r21882, MPFR_RNDN);
        mpfr_sub(r21884, r21847, r21854, MPFR_RNDN);
        mpfr_mul(r21885, r21883, r21884, MPFR_RNDN);
        mpfr_add(r21886, r21860, r21885, MPFR_RNDN);
        mpfr_div(r21887, r21849, r21886, MPFR_RNDN);
        ;
        mpfr_set_si(r21889, mpfr_cmp(r21844, r21888) <= 0, MPFR_RNDN);
        mpfr_sqrt(r21890, r21878, MPFR_RNDN);
        mpfr_mul(r21891, r21890, r21890, MPFR_RNDN);
        mpfr_div(r21892, r21849, r21891, MPFR_RNDN);
        if (mpfr_get_si(r21889, MPFR_RNDN)) { mpfr_set(r21893, r21892, MPFR_RNDN); } else { mpfr_set(r21893, r21887, MPFR_RNDN); };
        if (mpfr_get_si(r21881, MPFR_RNDN)) { mpfr_set(r21894, r21887, MPFR_RNDN); } else { mpfr_set(r21894, r21893, MPFR_RNDN); };
        if (mpfr_get_si(r21864, MPFR_RNDN)) { mpfr_set(r21895, r21879, MPFR_RNDN); } else { mpfr_set(r21895, r21894, MPFR_RNDN); };
        if (mpfr_get_si(r21846, MPFR_RNDN)) { mpfr_set(r21896, r21862, MPFR_RNDN); } else { mpfr_set(r21896, r21895, MPFR_RNDN); };
        return mpfr_get_d(r21896, MPFR_RNDN);
}

static mpfr_t r21897, r21898, r21899, r21900, r21901, r21902, r21903, r21904, r21905, r21906, r21907, r21908, r21909, r21910, r21911, r21912, r21913, r21914, r21915, r21916, r21917, r21918, r21919, r21920, r21921, r21922, r21923, r21924, r21925, r21926, r21927, r21928, r21929, r21930, r21931, r21932, r21933, r21934, r21935, r21936, r21937, r21938, r21939, r21940, r21941, r21942, r21943, r21944, r21945, r21946, r21947, r21948, r21949;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r21897);
        mpfr_init_set_str(r21898, "-6.447314359762627e+120", 10, MPFR_RNDN);
        mpfr_init(r21899);
        mpfr_init_set_str(r21900, "2", 10, MPFR_RNDN);
        mpfr_init(r21901);
        mpfr_init(r21902);
        mpfr_init(r21903);
        mpfr_init(r21904);
        mpfr_init(r21905);
        mpfr_init(r21906);
        mpfr_init_set_str(r21907, "1", 10, MPFR_RNDN);
        mpfr_init(r21908);
        mpfr_init(r21909);
        mpfr_init(r21910);
        mpfr_init(r21911);
        mpfr_init(r21912);
        mpfr_init(r21913);
        mpfr_init(r21914);
        mpfr_init(r21915);
        mpfr_init_set_str(r21916, "2.2327750609834536e-213", 10, MPFR_RNDN);
        mpfr_init(r21917);
        mpfr_init(r21918);
        mpfr_init(r21919);
        mpfr_init(r21920);
        mpfr_init_set_str(r21921, "4", 10, MPFR_RNDN);
        mpfr_init(r21922);
        mpfr_init(r21923);
        mpfr_init(r21924);
        mpfr_init(r21925);
        mpfr_init(r21926);
        mpfr_init(r21927);
        mpfr_init(r21928);
        mpfr_init(r21929);
        mpfr_init(r21930);
        mpfr_init(r21931);
        mpfr_init(r21932);
        mpfr_init_set_str(r21933, "1.7113149206191654e-132", 10, MPFR_RNDN);
        mpfr_init(r21934);
        mpfr_init(r21935);
        mpfr_init(r21936);
        mpfr_init(r21937);
        mpfr_init(r21938);
        mpfr_init(r21939);
        mpfr_init(r21940);
        mpfr_init_set_str(r21941, "4.776380960244624e+22", 10, MPFR_RNDN);
        mpfr_init(r21942);
        mpfr_init(r21943);
        mpfr_init(r21944);
        mpfr_init(r21945);
        mpfr_init(r21946);
        mpfr_init(r21947);
        mpfr_init(r21948);
        mpfr_init(r21949);
}

double f_dm(double x, double l, double t) {
        mpfr_set_d(r21897, t, MPFR_RNDN);
        ;
        mpfr_set_si(r21899, mpfr_cmp(r21897, r21898) <= 0, MPFR_RNDN);
        ;
        mpfr_sqrt(r21901, r21900, MPFR_RNDN);
        mpfr_mul(r21902, r21897, r21901, MPFR_RNDN);
        mpfr_div(r21903, r21897, r21901, MPFR_RNDN);
        mpfr_set_d(r21904, x, MPFR_RNDN);
        mpfr_mul(r21905, r21904, r21904, MPFR_RNDN);
        mpfr_div(r21906, r21903, r21905, MPFR_RNDN);
        ;
        mpfr_sub(r21908, r21907, r21900, MPFR_RNDN);
        mpfr_mul(r21909, r21906, r21908, MPFR_RNDN);
        mpfr_div(r21910, r21900, r21904, MPFR_RNDN);
        mpfr_div(r21911, r21910, r21901, MPFR_RNDN);
        mpfr_add(r21912, r21901, r21911, MPFR_RNDN);
        mpfr_mul(r21913, r21897, r21912, MPFR_RNDN);
        mpfr_sub(r21914, r21909, r21913, MPFR_RNDN);
        mpfr_div(r21915, r21902, r21914, MPFR_RNDN);
        ;
        mpfr_set_si(r21917, mpfr_cmp(r21897, r21916) <= 0, MPFR_RNDN);
        mpfr_sqrt(r21918, r21901, MPFR_RNDN);
        mpfr_mul(r21919, r21897, r21918, MPFR_RNDN);
        mpfr_mul(r21920, r21919, r21918, MPFR_RNDN);
        ;
        mpfr_div(r21922, r21921, r21904, MPFR_RNDN);
        mpfr_add(r21923, r21922, r21900, MPFR_RNDN);
        mpfr_mul(r21924, r21897, r21897, MPFR_RNDN);
        mpfr_mul(r21925, r21923, r21924, MPFR_RNDN);
        mpfr_set_d(r21926, l, MPFR_RNDN);
        mpfr_mul(r21927, r21900, r21926, MPFR_RNDN);
        mpfr_div(r21928, r21904, r21926, MPFR_RNDN);
        mpfr_div(r21929, r21927, r21928, MPFR_RNDN);
        mpfr_add(r21930, r21925, r21929, MPFR_RNDN);
        mpfr_sqrt(r21931, r21930, MPFR_RNDN);
        mpfr_div(r21932, r21920, r21931, MPFR_RNDN);
        ;
        mpfr_set_si(r21934, mpfr_cmp(r21897, r21933) <= 0, MPFR_RNDN);
        mpfr_mul(r21935, r21905, r21901, MPFR_RNDN);
        mpfr_div(r21936, r21897, r21935, MPFR_RNDN);
        mpfr_sub(r21937, r21900, r21907, MPFR_RNDN);
        mpfr_mul(r21938, r21936, r21937, MPFR_RNDN);
        mpfr_add(r21939, r21913, r21938, MPFR_RNDN);
        mpfr_div(r21940, r21902, r21939, MPFR_RNDN);
        ;
        mpfr_set_si(r21942, mpfr_cmp(r21897, r21941) <= 0, MPFR_RNDN);
        mpfr_sqrt(r21943, r21931, MPFR_RNDN);
        mpfr_mul(r21944, r21943, r21943, MPFR_RNDN);
        mpfr_div(r21945, r21902, r21944, MPFR_RNDN);
        if (mpfr_get_si(r21942, MPFR_RNDN)) { mpfr_set(r21946, r21945, MPFR_RNDN); } else { mpfr_set(r21946, r21940, MPFR_RNDN); };
        if (mpfr_get_si(r21934, MPFR_RNDN)) { mpfr_set(r21947, r21940, MPFR_RNDN); } else { mpfr_set(r21947, r21946, MPFR_RNDN); };
        if (mpfr_get_si(r21917, MPFR_RNDN)) { mpfr_set(r21948, r21932, MPFR_RNDN); } else { mpfr_set(r21948, r21947, MPFR_RNDN); };
        if (mpfr_get_si(r21899, MPFR_RNDN)) { mpfr_set(r21949, r21915, MPFR_RNDN); } else { mpfr_set(r21949, r21948, MPFR_RNDN); };
        return mpfr_get_d(r21949, MPFR_RNDN);
}