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

char *name = "NMSE Section 6.1 mentioned, A";

double f_if(float x, float eps) {
        float r18811 = 1.0f;
        float r18812 = eps;
        float r18813 = r18811 / r18812;
        float r18814 = r18811 + r18813;
        float r18815 = r18811 - r18812;
        float r18816 = x;
        float r18817 = r18815 * r18816;
        float r18818 = -r18817;
        float r18819 = exp(r18818);
        float r18820 = r18814 * r18819;
        float r18821 = r18813 - r18811;
        float r18822 = r18811 + r18812;
        float r18823 = r18822 * r18816;
        float r18824 = -r18823;
        float r18825 = exp(r18824);
        float r18826 = r18821 * r18825;
        float r18827 = r18820 - r18826;
        float r18828 = 2.0f;
        float r18829 = r18827 / r18828;
        return r18829;
}

double f_id(double x, double eps) {
        double r18830 = 1.0;
        double r18831 = eps;
        double r18832 = r18830 / r18831;
        double r18833 = r18830 + r18832;
        double r18834 = r18830 - r18831;
        double r18835 = x;
        double r18836 = r18834 * r18835;
        double r18837 = -r18836;
        double r18838 = exp(r18837);
        double r18839 = r18833 * r18838;
        double r18840 = r18832 - r18830;
        double r18841 = r18830 + r18831;
        double r18842 = r18841 * r18835;
        double r18843 = -r18842;
        double r18844 = exp(r18843);
        double r18845 = r18840 * r18844;
        double r18846 = r18839 - r18845;
        double r18847 = 2.0;
        double r18848 = r18846 / r18847;
        return r18848;
}


double f_of(float x, float eps) {
        float r18849 = x;
        float r18850 = 6252088.196986908f;
        bool r18851 = r18849 <= r18850;
        float r18852 = 1.0f;
        float r18853 = eps;
        float r18854 = r18852 / r18853;
        float r18855 = r18852 + r18854;
        float r18856 = r18852 - r18853;
        float r18857 = r18856 * r18849;
        float r18858 = -r18857;
        float r18859 = exp(r18858);
        float r18860 = r18855 * r18859;
        float r18861 = r18854 - r18852;
        float r18862 = r18852 + r18853;
        float r18863 = r18862 * r18849;
        float r18864 = -r18863;
        float r18865 = exp(r18864);
        float r18866 = r18861 * r18865;
        float r18867 = r18860 - r18866;
        float r18868 = 2.0f;
        float r18869 = r18867 / r18868;
        float r18870 = r18851 ? r18852 : r18869;
        return r18870;
}

double f_od(double x, double eps) {
        double r18871 = x;
        double r18872 = 6252088.196986908;
        bool r18873 = r18871 <= r18872;
        double r18874 = 1.0;
        double r18875 = eps;
        double r18876 = r18874 / r18875;
        double r18877 = r18874 + r18876;
        double r18878 = r18874 - r18875;
        double r18879 = r18878 * r18871;
        double r18880 = -r18879;
        double r18881 = exp(r18880);
        double r18882 = r18877 * r18881;
        double r18883 = r18876 - r18874;
        double r18884 = r18874 + r18875;
        double r18885 = r18884 * r18871;
        double r18886 = -r18885;
        double r18887 = exp(r18886);
        double r18888 = r18883 * r18887;
        double r18889 = r18882 - r18888;
        double r18890 = 2.0;
        double r18891 = r18889 / r18890;
        double r18892 = r18873 ? r18874 : r18891;
        return r18892;
}

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 r18893, r18894, r18895, r18896, r18897, r18898, r18899, r18900, r18901, r18902, r18903, r18904, r18905, r18906, r18907, r18908, r18909, r18910, r18911;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r18893, "1", 10, MPFR_RNDN);
        mpfr_init(r18894);
        mpfr_init(r18895);
        mpfr_init(r18896);
        mpfr_init(r18897);
        mpfr_init(r18898);
        mpfr_init(r18899);
        mpfr_init(r18900);
        mpfr_init(r18901);
        mpfr_init(r18902);
        mpfr_init(r18903);
        mpfr_init(r18904);
        mpfr_init(r18905);
        mpfr_init(r18906);
        mpfr_init(r18907);
        mpfr_init(r18908);
        mpfr_init(r18909);
        mpfr_init_set_str(r18910, "2", 10, MPFR_RNDN);
        mpfr_init(r18911);
}

double f_im(double x, double eps) {
        ;
        mpfr_set_d(r18894, eps, MPFR_RNDN);
        mpfr_div(r18895, r18893, r18894, MPFR_RNDN);
        mpfr_add(r18896, r18893, r18895, MPFR_RNDN);
        mpfr_sub(r18897, r18893, r18894, MPFR_RNDN);
        mpfr_set_d(r18898, x, MPFR_RNDN);
        mpfr_mul(r18899, r18897, r18898, MPFR_RNDN);
        mpfr_neg(r18900, r18899, MPFR_RNDN);
        mpfr_exp(r18901, r18900, MPFR_RNDN);
        mpfr_mul(r18902, r18896, r18901, MPFR_RNDN);
        mpfr_sub(r18903, r18895, r18893, MPFR_RNDN);
        mpfr_add(r18904, r18893, r18894, MPFR_RNDN);
        mpfr_mul(r18905, r18904, r18898, MPFR_RNDN);
        mpfr_neg(r18906, r18905, MPFR_RNDN);
        mpfr_exp(r18907, r18906, MPFR_RNDN);
        mpfr_mul(r18908, r18903, r18907, MPFR_RNDN);
        mpfr_sub(r18909, r18902, r18908, MPFR_RNDN);
        ;
        mpfr_div(r18911, r18909, r18910, MPFR_RNDN);
        return mpfr_get_d(r18911, MPFR_RNDN);
}

static mpfr_t r18912, r18913, r18914, r18915, r18916, r18917, r18918, r18919, r18920, r18921, r18922, r18923, r18924, r18925, r18926, r18927, r18928, r18929, r18930, r18931, r18932, r18933;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r18912);
        mpfr_init_set_str(r18913, "6252088.196986908", 10, MPFR_RNDN);
        mpfr_init(r18914);
        mpfr_init_set_str(r18915, "1", 10, MPFR_RNDN);
        mpfr_init(r18916);
        mpfr_init(r18917);
        mpfr_init(r18918);
        mpfr_init(r18919);
        mpfr_init(r18920);
        mpfr_init(r18921);
        mpfr_init(r18922);
        mpfr_init(r18923);
        mpfr_init(r18924);
        mpfr_init(r18925);
        mpfr_init(r18926);
        mpfr_init(r18927);
        mpfr_init(r18928);
        mpfr_init(r18929);
        mpfr_init(r18930);
        mpfr_init_set_str(r18931, "2", 10, MPFR_RNDN);
        mpfr_init(r18932);
        mpfr_init(r18933);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r18912, x, MPFR_RNDN);
        ;
        mpfr_set_si(r18914, mpfr_cmp(r18912, r18913) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r18916, eps, MPFR_RNDN);
        mpfr_div(r18917, r18915, r18916, MPFR_RNDN);
        mpfr_add(r18918, r18915, r18917, MPFR_RNDN);
        mpfr_sub(r18919, r18915, r18916, MPFR_RNDN);
        mpfr_mul(r18920, r18919, r18912, MPFR_RNDN);
        mpfr_neg(r18921, r18920, MPFR_RNDN);
        mpfr_exp(r18922, r18921, MPFR_RNDN);
        mpfr_mul(r18923, r18918, r18922, MPFR_RNDN);
        mpfr_sub(r18924, r18917, r18915, MPFR_RNDN);
        mpfr_add(r18925, r18915, r18916, MPFR_RNDN);
        mpfr_mul(r18926, r18925, r18912, MPFR_RNDN);
        mpfr_neg(r18927, r18926, MPFR_RNDN);
        mpfr_exp(r18928, r18927, MPFR_RNDN);
        mpfr_mul(r18929, r18924, r18928, MPFR_RNDN);
        mpfr_sub(r18930, r18923, r18929, MPFR_RNDN);
        ;
        mpfr_div(r18932, r18930, r18931, MPFR_RNDN);
        if (mpfr_get_si(r18914, MPFR_RNDN)) { mpfr_set(r18933, r18915, MPFR_RNDN); } else { mpfr_set(r18933, r18932, MPFR_RNDN); };
        return mpfr_get_d(r18933, MPFR_RNDN);
}

static mpfr_t r18934, r18935, r18936, r18937, r18938, r18939, r18940, r18941, r18942, r18943, r18944, r18945, r18946, r18947, r18948, r18949, r18950, r18951, r18952, r18953, r18954, r18955;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r18934);
        mpfr_init_set_str(r18935, "6252088.196986908", 10, MPFR_RNDN);
        mpfr_init(r18936);
        mpfr_init_set_str(r18937, "1", 10, MPFR_RNDN);
        mpfr_init(r18938);
        mpfr_init(r18939);
        mpfr_init(r18940);
        mpfr_init(r18941);
        mpfr_init(r18942);
        mpfr_init(r18943);
        mpfr_init(r18944);
        mpfr_init(r18945);
        mpfr_init(r18946);
        mpfr_init(r18947);
        mpfr_init(r18948);
        mpfr_init(r18949);
        mpfr_init(r18950);
        mpfr_init(r18951);
        mpfr_init(r18952);
        mpfr_init_set_str(r18953, "2", 10, MPFR_RNDN);
        mpfr_init(r18954);
        mpfr_init(r18955);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r18934, x, MPFR_RNDN);
        ;
        mpfr_set_si(r18936, mpfr_cmp(r18934, r18935) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r18938, eps, MPFR_RNDN);
        mpfr_div(r18939, r18937, r18938, MPFR_RNDN);
        mpfr_add(r18940, r18937, r18939, MPFR_RNDN);
        mpfr_sub(r18941, r18937, r18938, MPFR_RNDN);
        mpfr_mul(r18942, r18941, r18934, MPFR_RNDN);
        mpfr_neg(r18943, r18942, MPFR_RNDN);
        mpfr_exp(r18944, r18943, MPFR_RNDN);
        mpfr_mul(r18945, r18940, r18944, MPFR_RNDN);
        mpfr_sub(r18946, r18939, r18937, MPFR_RNDN);
        mpfr_add(r18947, r18937, r18938, MPFR_RNDN);
        mpfr_mul(r18948, r18947, r18934, MPFR_RNDN);
        mpfr_neg(r18949, r18948, MPFR_RNDN);
        mpfr_exp(r18950, r18949, MPFR_RNDN);
        mpfr_mul(r18951, r18946, r18950, MPFR_RNDN);
        mpfr_sub(r18952, r18945, r18951, MPFR_RNDN);
        ;
        mpfr_div(r18954, r18952, r18953, MPFR_RNDN);
        if (mpfr_get_si(r18936, MPFR_RNDN)) { mpfr_set(r18955, r18937, MPFR_RNDN); } else { mpfr_set(r18955, r18954, MPFR_RNDN); };
        return mpfr_get_d(r18955, MPFR_RNDN);
}