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

char *name = "Falkner and Boettcher, Appendix A";

double f_if(float a, float k, float m) {
        float r10855 = a;
        float r10856 = k;
        float r10857 = m;
        float r10858 = pow(r10856, r10857);
        float r10859 = r10855 * r10858;
        float r10860 = 1.0f;
        float r10861 = 10.0f;
        float r10862 = r10861 * r10856;
        float r10863 = r10860 + r10862;
        float r10864 = r10856 * r10856;
        float r10865 = r10863 + r10864;
        float r10866 = r10859 / r10865;
        return r10866;
}

double f_id(double a, double k, double m) {
        double r10867 = a;
        double r10868 = k;
        double r10869 = m;
        double r10870 = pow(r10868, r10869);
        double r10871 = r10867 * r10870;
        double r10872 = 1.0;
        double r10873 = 10.0;
        double r10874 = r10873 * r10868;
        double r10875 = r10872 + r10874;
        double r10876 = r10868 * r10868;
        double r10877 = r10875 + r10876;
        double r10878 = r10871 / r10877;
        return r10878;
}


double f_of(float a, float k, float m) {
        float r10879 = k;
        float r10880 = 3.9827950205668443e+118f;
        bool r10881 = r10879 <= r10880;
        float r10882 = a;
        float r10883 = 10.0f;
        float r10884 = r10879 + r10883;
        float r10885 = 1.0f;
        float r10886 = fma(r10879, r10884, r10885);
        float r10887 = sqrt(r10886);
        float r10888 = r10882 / r10887;
        float r10889 = m;
        float r10890 = pow(r10879, r10889);
        float r10891 = r10888 * r10890;
        float r10892 = r10891 / r10887;
        float r10893 = exp(r10889);
        float r10894 = log(r10879);
        float r10895 = pow(r10893, r10894);
        float r10896 = 4.0f;
        float r10897 = pow(r10879, r10896);
        float r10898 = r10882 / r10897;
        float r10899 = r10895 * r10898;
        float r10900 = 99.0f;
        float r10901 = r10882 / r10879;
        float r10902 = r10883 / r10879;
        float r10903 = r10902 * r10901;
        float r10904 = r10901 - r10903;
        float r10905 = r10895 / r10879;
        float r10906 = r10904 * r10905;
        float r10907 = fma(r10899, r10900, r10906);
        float r10908 = r10881 ? r10892 : r10907;
        return r10908;
}

double f_od(double a, double k, double m) {
        double r10909 = k;
        double r10910 = 3.9827950205668443e+118;
        bool r10911 = r10909 <= r10910;
        double r10912 = a;
        double r10913 = 10.0;
        double r10914 = r10909 + r10913;
        double r10915 = 1.0;
        double r10916 = fma(r10909, r10914, r10915);
        double r10917 = sqrt(r10916);
        double r10918 = r10912 / r10917;
        double r10919 = m;
        double r10920 = pow(r10909, r10919);
        double r10921 = r10918 * r10920;
        double r10922 = r10921 / r10917;
        double r10923 = exp(r10919);
        double r10924 = log(r10909);
        double r10925 = pow(r10923, r10924);
        double r10926 = 4.0;
        double r10927 = pow(r10909, r10926);
        double r10928 = r10912 / r10927;
        double r10929 = r10925 * r10928;
        double r10930 = 99.0;
        double r10931 = r10912 / r10909;
        double r10932 = r10913 / r10909;
        double r10933 = r10932 * r10931;
        double r10934 = r10931 - r10933;
        double r10935 = r10925 / r10909;
        double r10936 = r10934 * r10935;
        double r10937 = fma(r10929, r10930, r10936);
        double r10938 = r10911 ? r10922 : r10937;
        return r10938;
}

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 r10939, r10940, r10941, r10942, r10943, r10944, r10945, r10946, r10947, r10948, r10949, r10950;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r10939);
        mpfr_init(r10940);
        mpfr_init(r10941);
        mpfr_init(r10942);
        mpfr_init(r10943);
        mpfr_init_set_str(r10944, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r10945, "10", 10, MPFR_RNDN);
        mpfr_init(r10946);
        mpfr_init(r10947);
        mpfr_init(r10948);
        mpfr_init(r10949);
        mpfr_init(r10950);
}

double f_im(double a, double k, double m) {
        mpfr_set_d(r10939, a, MPFR_RNDN);
        mpfr_set_d(r10940, k, MPFR_RNDN);
        mpfr_set_d(r10941, m, MPFR_RNDN);
        mpfr_pow(r10942, r10940, r10941, MPFR_RNDN);
        mpfr_mul(r10943, r10939, r10942, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r10946, r10945, r10940, MPFR_RNDN);
        mpfr_add(r10947, r10944, r10946, MPFR_RNDN);
        mpfr_mul(r10948, r10940, r10940, MPFR_RNDN);
        mpfr_add(r10949, r10947, r10948, MPFR_RNDN);
        mpfr_div(r10950, r10943, r10949, MPFR_RNDN);
        return mpfr_get_d(r10950, MPFR_RNDN);
}

static mpfr_t r10951, r10952, r10953, r10954, r10955, r10956, r10957, r10958, r10959, r10960, r10961, r10962, r10963, r10964, r10965, r10966, r10967, r10968, r10969, r10970, r10971, r10972, r10973, r10974, r10975, r10976, r10977, r10978, r10979, r10980;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r10951);
        mpfr_init_set_str(r10952, "3.9827950205668443e+118", 10, MPFR_RNDN);
        mpfr_init(r10953);
        mpfr_init(r10954);
        mpfr_init_set_str(r10955, "10", 10, MPFR_RNDN);
        mpfr_init(r10956);
        mpfr_init_set_str(r10957, "1", 10, MPFR_RNDN);
        mpfr_init(r10958);
        mpfr_init(r10959);
        mpfr_init(r10960);
        mpfr_init(r10961);
        mpfr_init(r10962);
        mpfr_init(r10963);
        mpfr_init(r10964);
        mpfr_init(r10965);
        mpfr_init(r10966);
        mpfr_init(r10967);
        mpfr_init_set_str(r10968, "4", 10, MPFR_RNDN);
        mpfr_init(r10969);
        mpfr_init(r10970);
        mpfr_init(r10971);
        mpfr_init_set_str(r10972, "99", 10, MPFR_RNDN);
        mpfr_init(r10973);
        mpfr_init(r10974);
        mpfr_init(r10975);
        mpfr_init(r10976);
        mpfr_init(r10977);
        mpfr_init(r10978);
        mpfr_init(r10979);
        mpfr_init(r10980);
}

double f_fm(double a, double k, double m) {
        mpfr_set_d(r10951, k, MPFR_RNDN);
        ;
        mpfr_set_si(r10953, mpfr_cmp(r10951, r10952) <= 0, MPFR_RNDN);
        mpfr_set_d(r10954, a, MPFR_RNDN);
        ;
        mpfr_add(r10956, r10951, r10955, MPFR_RNDN);
        ;
        mpfr_fma(r10958, r10951, r10956, r10957, MPFR_RNDN);
        mpfr_sqrt(r10959, r10958, MPFR_RNDN);
        mpfr_div(r10960, r10954, r10959, MPFR_RNDN);
        mpfr_set_d(r10961, m, MPFR_RNDN);
        mpfr_pow(r10962, r10951, r10961, MPFR_RNDN);
        mpfr_mul(r10963, r10960, r10962, MPFR_RNDN);
        mpfr_div(r10964, r10963, r10959, MPFR_RNDN);
        mpfr_exp(r10965, r10961, MPFR_RNDN);
        mpfr_log(r10966, r10951, MPFR_RNDN);
        mpfr_pow(r10967, r10965, r10966, MPFR_RNDN);
        ;
        mpfr_pow(r10969, r10951, r10968, MPFR_RNDN);
        mpfr_div(r10970, r10954, r10969, MPFR_RNDN);
        mpfr_mul(r10971, r10967, r10970, MPFR_RNDN);
        ;
        mpfr_div(r10973, r10954, r10951, MPFR_RNDN);
        mpfr_div(r10974, r10955, r10951, MPFR_RNDN);
        mpfr_mul(r10975, r10974, r10973, MPFR_RNDN);
        mpfr_sub(r10976, r10973, r10975, MPFR_RNDN);
        mpfr_div(r10977, r10967, r10951, MPFR_RNDN);
        mpfr_mul(r10978, r10976, r10977, MPFR_RNDN);
        mpfr_fma(r10979, r10971, r10972, r10978, MPFR_RNDN);
        if (mpfr_get_si(r10953, MPFR_RNDN)) { mpfr_set(r10980, r10964, MPFR_RNDN); } else { mpfr_set(r10980, r10979, MPFR_RNDN); };
        return mpfr_get_d(r10980, MPFR_RNDN);
}

static mpfr_t r10981, r10982, r10983, r10984, r10985, r10986, r10987, r10988, r10989, r10990, r10991, r10992, r10993, r10994, r10995, r10996, r10997, r10998, r10999, r11000, r11001, r11002, r11003, r11004, r11005, r11006, r11007, r11008, r11009, r11010;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r10981);
        mpfr_init_set_str(r10982, "3.9827950205668443e+118", 10, MPFR_RNDN);
        mpfr_init(r10983);
        mpfr_init(r10984);
        mpfr_init_set_str(r10985, "10", 10, MPFR_RNDN);
        mpfr_init(r10986);
        mpfr_init_set_str(r10987, "1", 10, MPFR_RNDN);
        mpfr_init(r10988);
        mpfr_init(r10989);
        mpfr_init(r10990);
        mpfr_init(r10991);
        mpfr_init(r10992);
        mpfr_init(r10993);
        mpfr_init(r10994);
        mpfr_init(r10995);
        mpfr_init(r10996);
        mpfr_init(r10997);
        mpfr_init_set_str(r10998, "4", 10, MPFR_RNDN);
        mpfr_init(r10999);
        mpfr_init(r11000);
        mpfr_init(r11001);
        mpfr_init_set_str(r11002, "99", 10, MPFR_RNDN);
        mpfr_init(r11003);
        mpfr_init(r11004);
        mpfr_init(r11005);
        mpfr_init(r11006);
        mpfr_init(r11007);
        mpfr_init(r11008);
        mpfr_init(r11009);
        mpfr_init(r11010);
}

double f_dm(double a, double k, double m) {
        mpfr_set_d(r10981, k, MPFR_RNDN);
        ;
        mpfr_set_si(r10983, mpfr_cmp(r10981, r10982) <= 0, MPFR_RNDN);
        mpfr_set_d(r10984, a, MPFR_RNDN);
        ;
        mpfr_add(r10986, r10981, r10985, MPFR_RNDN);
        ;
        mpfr_fma(r10988, r10981, r10986, r10987, MPFR_RNDN);
        mpfr_sqrt(r10989, r10988, MPFR_RNDN);
        mpfr_div(r10990, r10984, r10989, MPFR_RNDN);
        mpfr_set_d(r10991, m, MPFR_RNDN);
        mpfr_pow(r10992, r10981, r10991, MPFR_RNDN);
        mpfr_mul(r10993, r10990, r10992, MPFR_RNDN);
        mpfr_div(r10994, r10993, r10989, MPFR_RNDN);
        mpfr_exp(r10995, r10991, MPFR_RNDN);
        mpfr_log(r10996, r10981, MPFR_RNDN);
        mpfr_pow(r10997, r10995, r10996, MPFR_RNDN);
        ;
        mpfr_pow(r10999, r10981, r10998, MPFR_RNDN);
        mpfr_div(r11000, r10984, r10999, MPFR_RNDN);
        mpfr_mul(r11001, r10997, r11000, MPFR_RNDN);
        ;
        mpfr_div(r11003, r10984, r10981, MPFR_RNDN);
        mpfr_div(r11004, r10985, r10981, MPFR_RNDN);
        mpfr_mul(r11005, r11004, r11003, MPFR_RNDN);
        mpfr_sub(r11006, r11003, r11005, MPFR_RNDN);
        mpfr_div(r11007, r10997, r10981, MPFR_RNDN);
        mpfr_mul(r11008, r11006, r11007, MPFR_RNDN);
        mpfr_fma(r11009, r11001, r11002, r11008, MPFR_RNDN);
        if (mpfr_get_si(r10983, MPFR_RNDN)) { mpfr_set(r11010, r10994, MPFR_RNDN); } else { mpfr_set(r11010, r11009, MPFR_RNDN); };
        return mpfr_get_d(r11010, MPFR_RNDN);
}