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

char *name = "Octave 3.8, jcobi/4, as called";

double f_if(float i) {
        float r16967 = i;
        float r16968 = r16967 * r16967;
        float r16969 = r16968 * r16968;
        float r16970 = 2.0f;
        float r16971 = r16970 * r16967;
        float r16972 = r16971 * r16971;
        float r16973 = r16969 / r16972;
        float r16974 = 1.0f;
        float r16975 = r16972 - r16974;
        float r16976 = r16973 / r16975;
        return r16976;
}

double f_id(double i) {
        double r16977 = i;
        double r16978 = r16977 * r16977;
        double r16979 = r16978 * r16978;
        double r16980 = 2.0;
        double r16981 = r16980 * r16977;
        double r16982 = r16981 * r16981;
        double r16983 = r16979 / r16982;
        double r16984 = 1.0;
        double r16985 = r16982 - r16984;
        double r16986 = r16983 / r16985;
        return r16986;
}


double f_of(float i) {
        float r16987 = i;
        float r16988 = 49809.97373367609f;
        bool r16989 = r16987 <= r16988;
        float r16990 = 2.0f;
        float r16991 = r16987 / r16990;
        float r16992 = r16991 * r16991;
        float r16993 = r16987 * r16990;
        float r16994 = r16993 * r16993;
        float r16995 = 1.0f;
        float r16996 = r16994 - r16995;
        float r16997 = r16992 / r16996;
        float r16998 = 0.0625f;
        float r16999 = 0.00390625f;
        float r17000 = 4.0f;
        float r17001 = pow(r16987, r17000);
        float r17002 = r16999 / r17001;
        float r17003 = r16998 + r17002;
        float r17004 = 0.015625f;
        float r17005 = r16987 * r16987;
        float r17006 = r17004 / r17005;
        float r17007 = r17003 + r17006;
        float r17008 = r16989 ? r16997 : r17007;
        return r17008;
}

double f_od(double i) {
        double r17009 = i;
        double r17010 = 49809.97373367609;
        bool r17011 = r17009 <= r17010;
        double r17012 = 2.0;
        double r17013 = r17009 / r17012;
        double r17014 = r17013 * r17013;
        double r17015 = r17009 * r17012;
        double r17016 = r17015 * r17015;
        double r17017 = 1.0;
        double r17018 = r17016 - r17017;
        double r17019 = r17014 / r17018;
        double r17020 = 0.0625;
        double r17021 = 0.00390625;
        double r17022 = 4.0;
        double r17023 = pow(r17009, r17022);
        double r17024 = r17021 / r17023;
        double r17025 = r17020 + r17024;
        double r17026 = 0.015625;
        double r17027 = r17009 * r17009;
        double r17028 = r17026 / r17027;
        double r17029 = r17025 + r17028;
        double r17030 = r17011 ? r17019 : r17029;
        return r17030;
}

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 r17031, r17032, r17033, r17034, r17035, r17036, r17037, r17038, r17039, r17040;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17031);
        mpfr_init(r17032);
        mpfr_init(r17033);
        mpfr_init_set_str(r17034, "2", 10, MPFR_RNDN);
        mpfr_init(r17035);
        mpfr_init(r17036);
        mpfr_init(r17037);
        mpfr_init_set_str(r17038, "1.0", 10, MPFR_RNDN);
        mpfr_init(r17039);
        mpfr_init(r17040);
}

double f_im(double i) {
        mpfr_set_d(r17031, i, MPFR_RNDN);
        mpfr_mul(r17032, r17031, r17031, MPFR_RNDN);
        mpfr_mul(r17033, r17032, r17032, MPFR_RNDN);
        ;
        mpfr_mul(r17035, r17034, r17031, MPFR_RNDN);
        mpfr_mul(r17036, r17035, r17035, MPFR_RNDN);
        mpfr_div(r17037, r17033, r17036, MPFR_RNDN);
        ;
        mpfr_sub(r17039, r17036, r17038, MPFR_RNDN);
        mpfr_div(r17040, r17037, r17039, MPFR_RNDN);
        return mpfr_get_d(r17040, MPFR_RNDN);
}

static mpfr_t r17041, r17042, r17043, r17044, r17045, r17046, r17047, r17048, r17049, r17050, r17051, r17052, r17053, r17054, r17055, r17056, r17057, r17058, r17059, r17060, r17061, r17062;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17041);
        mpfr_init_set_str(r17042, "49809.97373367609", 10, MPFR_RNDN);
        mpfr_init(r17043);
        mpfr_init_set_str(r17044, "2", 10, MPFR_RNDN);
        mpfr_init(r17045);
        mpfr_init(r17046);
        mpfr_init(r17047);
        mpfr_init(r17048);
        mpfr_init_set_str(r17049, "1.0", 10, MPFR_RNDN);
        mpfr_init(r17050);
        mpfr_init(r17051);
        mpfr_init_set_str(r17052, "1/16", 10, MPFR_RNDN);
        mpfr_init_set_str(r17053, "0.00390625", 10, MPFR_RNDN);
        mpfr_init_set_str(r17054, "4", 10, MPFR_RNDN);
        mpfr_init(r17055);
        mpfr_init(r17056);
        mpfr_init(r17057);
        mpfr_init_set_str(r17058, "0.015625", 10, MPFR_RNDN);
        mpfr_init(r17059);
        mpfr_init(r17060);
        mpfr_init(r17061);
        mpfr_init(r17062);
}

double f_fm(double i) {
        mpfr_set_d(r17041, i, MPFR_RNDN);
        ;
        mpfr_set_si(r17043, mpfr_cmp(r17041, r17042) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r17045, r17041, r17044, MPFR_RNDN);
        mpfr_sqr(r17046, r17045, MPFR_RNDN);
        mpfr_mul(r17047, r17041, r17044, MPFR_RNDN);
        mpfr_mul(r17048, r17047, r17047, MPFR_RNDN);
        ;
        mpfr_sub(r17050, r17048, r17049, MPFR_RNDN);
        mpfr_div(r17051, r17046, r17050, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_pow(r17055, r17041, r17054, MPFR_RNDN);
        mpfr_div(r17056, r17053, r17055, MPFR_RNDN);
        mpfr_add(r17057, r17052, r17056, MPFR_RNDN);
        ;
        mpfr_mul(r17059, r17041, r17041, MPFR_RNDN);
        mpfr_div(r17060, r17058, r17059, MPFR_RNDN);
        mpfr_add(r17061, r17057, r17060, MPFR_RNDN);
        if (mpfr_get_si(r17043, MPFR_RNDN)) { mpfr_set(r17062, r17051, MPFR_RNDN); } else { mpfr_set(r17062, r17061, MPFR_RNDN); };
        return mpfr_get_d(r17062, MPFR_RNDN);
}

static mpfr_t r17063, r17064, r17065, r17066, r17067, r17068, r17069, r17070, r17071, r17072, r17073, r17074, r17075, r17076, r17077, r17078, r17079, r17080, r17081, r17082, r17083, r17084;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17063);
        mpfr_init_set_str(r17064, "49809.97373367609", 10, MPFR_RNDN);
        mpfr_init(r17065);
        mpfr_init_set_str(r17066, "2", 10, MPFR_RNDN);
        mpfr_init(r17067);
        mpfr_init(r17068);
        mpfr_init(r17069);
        mpfr_init(r17070);
        mpfr_init_set_str(r17071, "1.0", 10, MPFR_RNDN);
        mpfr_init(r17072);
        mpfr_init(r17073);
        mpfr_init_set_str(r17074, "1/16", 10, MPFR_RNDN);
        mpfr_init_set_str(r17075, "0.00390625", 10, MPFR_RNDN);
        mpfr_init_set_str(r17076, "4", 10, MPFR_RNDN);
        mpfr_init(r17077);
        mpfr_init(r17078);
        mpfr_init(r17079);
        mpfr_init_set_str(r17080, "0.015625", 10, MPFR_RNDN);
        mpfr_init(r17081);
        mpfr_init(r17082);
        mpfr_init(r17083);
        mpfr_init(r17084);
}

double f_dm(double i) {
        mpfr_set_d(r17063, i, MPFR_RNDN);
        ;
        mpfr_set_si(r17065, mpfr_cmp(r17063, r17064) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r17067, r17063, r17066, MPFR_RNDN);
        mpfr_sqr(r17068, r17067, MPFR_RNDN);
        mpfr_mul(r17069, r17063, r17066, MPFR_RNDN);
        mpfr_mul(r17070, r17069, r17069, MPFR_RNDN);
        ;
        mpfr_sub(r17072, r17070, r17071, MPFR_RNDN);
        mpfr_div(r17073, r17068, r17072, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_pow(r17077, r17063, r17076, MPFR_RNDN);
        mpfr_div(r17078, r17075, r17077, MPFR_RNDN);
        mpfr_add(r17079, r17074, r17078, MPFR_RNDN);
        ;
        mpfr_mul(r17081, r17063, r17063, MPFR_RNDN);
        mpfr_div(r17082, r17080, r17081, MPFR_RNDN);
        mpfr_add(r17083, r17079, r17082, MPFR_RNDN);
        if (mpfr_get_si(r17065, MPFR_RNDN)) { mpfr_set(r17084, r17073, MPFR_RNDN); } else { mpfr_set(r17084, r17083, MPFR_RNDN); };
        return mpfr_get_d(r17084, MPFR_RNDN);
}