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

char *name = "Jmat.Real.erfi, branch x less than or equal to 0.5";

double f_if(float x) {
        float r18000 = 1.0f;
        float r18001 = atan2(1.0, 0.0);
        float r18002 = sqrt(r18001);
        float r18003 = r18000 / r18002;
        float r18004 = 2.0f;
        float r18005 = x;
        float r18006 = fabs(r18005);
        float r18007 = r18004 * r18006;
        float r18008 = 3.0f;
        float r18009 = r18004 / r18008;
        float r18010 = r18006 * r18006;
        float r18011 = r18010 * r18006;
        float r18012 = r18009 * r18011;
        float r18013 = r18007 + r18012;
        float r18014 = 5.0f;
        float r18015 = r18000 / r18014;
        float r18016 = r18011 * r18006;
        float r18017 = r18016 * r18006;
        float r18018 = r18015 * r18017;
        float r18019 = r18013 + r18018;
        float r18020 = 21.0f;
        float r18021 = r18000 / r18020;
        float r18022 = r18017 * r18006;
        float r18023 = r18022 * r18006;
        float r18024 = r18021 * r18023;
        float r18025 = r18019 + r18024;
        float r18026 = r18003 * r18025;
        float r18027 = fabs(r18026);
        return r18027;
}

double f_id(double x) {
        double r18028 = 1.0;
        double r18029 = atan2(1.0, 0.0);
        double r18030 = sqrt(r18029);
        double r18031 = r18028 / r18030;
        double r18032 = 2.0;
        double r18033 = x;
        double r18034 = fabs(r18033);
        double r18035 = r18032 * r18034;
        double r18036 = 3.0;
        double r18037 = r18032 / r18036;
        double r18038 = r18034 * r18034;
        double r18039 = r18038 * r18034;
        double r18040 = r18037 * r18039;
        double r18041 = r18035 + r18040;
        double r18042 = 5.0;
        double r18043 = r18028 / r18042;
        double r18044 = r18039 * r18034;
        double r18045 = r18044 * r18034;
        double r18046 = r18043 * r18045;
        double r18047 = r18041 + r18046;
        double r18048 = 21.0;
        double r18049 = r18028 / r18048;
        double r18050 = r18045 * r18034;
        double r18051 = r18050 * r18034;
        double r18052 = r18049 * r18051;
        double r18053 = r18047 + r18052;
        double r18054 = r18031 * r18053;
        double r18055 = fabs(r18054);
        return r18055;
}


double f_of(float x) {
        float r18056 = 1.0f;
        float r18057 = atan2(1.0, 0.0);
        float r18058 = r18056 / r18057;
        float r18059 = sqrt(r18058);
        float r18060 = 0.2f;
        float r18061 = x;
        float r18062 = fabs(r18061);
        float r18063 = r18060 * r18062;
        float r18064 = r18062 * (r18062 * r18062);
        float r18065 = r18063 * r18064;
        float r18066 = 0.6666666666666666f;
        float r18067 = 2.0f;
        float r18068 = r18062 * r18067;
        float r18069 = fma(r18066, r18064, r18068);
        float r18070 = fma(r18065, r18062, r18069);
        float r18071 = 0.047619047619047616f;
        float r18072 = r18062 * r18071;
        float r18073 = r18059 * r18072;
        float r18074 = r18062 * r18062;
        float r18075 = r18074 * (r18074 * r18074);
        float r18076 = r18073 * r18075;
        float r18077 = fma(r18059, r18070, r18076);
        float r18078 = fabs(r18077);
        return r18078;
}

double f_od(double x) {
        double r18079 = 1.0;
        double r18080 = atan2(1.0, 0.0);
        double r18081 = r18079 / r18080;
        double r18082 = sqrt(r18081);
        double r18083 = 0.2;
        double r18084 = x;
        double r18085 = fabs(r18084);
        double r18086 = r18083 * r18085;
        double r18087 = r18085 * (r18085 * r18085);
        double r18088 = r18086 * r18087;
        double r18089 = 0.6666666666666666;
        double r18090 = 2.0;
        double r18091 = r18085 * r18090;
        double r18092 = fma(r18089, r18087, r18091);
        double r18093 = fma(r18088, r18085, r18092);
        double r18094 = 0.047619047619047616;
        double r18095 = r18085 * r18094;
        double r18096 = r18082 * r18095;
        double r18097 = r18085 * r18085;
        double r18098 = r18097 * (r18097 * r18097);
        double r18099 = r18096 * r18098;
        double r18100 = fma(r18082, r18093, r18099);
        double r18101 = fabs(r18100);
        return r18101;
}

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 r18102, r18103, r18104, r18105, r18106, r18107, r18108, r18109, r18110, r18111, r18112, r18113, r18114, r18115, r18116, r18117, r18118, r18119, r18120, r18121, r18122, r18123, r18124, r18125, r18126, r18127, r18128, r18129;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18102, "1", 10, MPFR_RNDN);
        mpfr_init(r18103);
        mpfr_init(r18104);
        mpfr_init(r18105);
        mpfr_init_set_str(r18106, "2", 10, MPFR_RNDN);
        mpfr_init(r18107);
        mpfr_init(r18108);
        mpfr_init(r18109);
        mpfr_init_set_str(r18110, "3", 10, MPFR_RNDN);
        mpfr_init(r18111);
        mpfr_init(r18112);
        mpfr_init(r18113);
        mpfr_init(r18114);
        mpfr_init(r18115);
        mpfr_init_set_str(r18116, "5", 10, MPFR_RNDN);
        mpfr_init(r18117);
        mpfr_init(r18118);
        mpfr_init(r18119);
        mpfr_init(r18120);
        mpfr_init(r18121);
        mpfr_init_set_str(r18122, "21", 10, MPFR_RNDN);
        mpfr_init(r18123);
        mpfr_init(r18124);
        mpfr_init(r18125);
        mpfr_init(r18126);
        mpfr_init(r18127);
        mpfr_init(r18128);
        mpfr_init(r18129);
}

double f_im(double x) {
        ;
        mpfr_const_pi(r18103, MPFR_RNDN);
        mpfr_sqrt(r18104, r18103, MPFR_RNDN);
        mpfr_div(r18105, r18102, r18104, MPFR_RNDN);
        ;
        mpfr_set_d(r18107, x, MPFR_RNDN);
        mpfr_abs(r18108, r18107, MPFR_RNDN);
        mpfr_mul(r18109, r18106, r18108, MPFR_RNDN);
        ;
        mpfr_div(r18111, r18106, r18110, MPFR_RNDN);
        mpfr_mul(r18112, r18108, r18108, MPFR_RNDN);
        mpfr_mul(r18113, r18112, r18108, MPFR_RNDN);
        mpfr_mul(r18114, r18111, r18113, MPFR_RNDN);
        mpfr_add(r18115, r18109, r18114, MPFR_RNDN);
        ;
        mpfr_div(r18117, r18102, r18116, MPFR_RNDN);
        mpfr_mul(r18118, r18113, r18108, MPFR_RNDN);
        mpfr_mul(r18119, r18118, r18108, MPFR_RNDN);
        mpfr_mul(r18120, r18117, r18119, MPFR_RNDN);
        mpfr_add(r18121, r18115, r18120, MPFR_RNDN);
        ;
        mpfr_div(r18123, r18102, r18122, MPFR_RNDN);
        mpfr_mul(r18124, r18119, r18108, MPFR_RNDN);
        mpfr_mul(r18125, r18124, r18108, MPFR_RNDN);
        mpfr_mul(r18126, r18123, r18125, MPFR_RNDN);
        mpfr_add(r18127, r18121, r18126, MPFR_RNDN);
        mpfr_mul(r18128, r18105, r18127, MPFR_RNDN);
        mpfr_abs(r18129, r18128, MPFR_RNDN);
        return mpfr_get_d(r18129, MPFR_RNDN);
}

static mpfr_t r18130, r18131, r18132, r18133, r18134, r18135, r18136, r18137, r18138, r18139, r18140, r18141, r18142, r18143, r18144, r18145, r18146, r18147, r18148, r18149, r18150, r18151, r18152;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18130, "1", 10, MPFR_RNDN);
        mpfr_init(r18131);
        mpfr_init(r18132);
        mpfr_init(r18133);
        mpfr_init_set_str(r18134, "1/5", 10, MPFR_RNDN);
        mpfr_init(r18135);
        mpfr_init(r18136);
        mpfr_init(r18137);
        mpfr_init(r18138);
        mpfr_init(r18139);
        mpfr_init_set_str(r18140, "2/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r18141, "2", 10, MPFR_RNDN);
        mpfr_init(r18142);
        mpfr_init(r18143);
        mpfr_init(r18144);
        mpfr_init_set_str(r18145, "1/21", 10, MPFR_RNDN);
        mpfr_init(r18146);
        mpfr_init(r18147);
        mpfr_init(r18148);
        mpfr_init(r18149);
        mpfr_init(r18150);
        mpfr_init(r18151);
        mpfr_init(r18152);
}

double f_fm(double x) {
        ;
        mpfr_const_pi(r18131, MPFR_RNDN);
        mpfr_div(r18132, r18130, r18131, MPFR_RNDN);
        mpfr_sqrt(r18133, r18132, MPFR_RNDN);
        ;
        mpfr_set_d(r18135, x, MPFR_RNDN);
        mpfr_abs(r18136, r18135, MPFR_RNDN);
        mpfr_mul(r18137, r18134, r18136, MPFR_RNDN);
        mpfr_mul(r18138, r18136, r18136, MPFR_RNDN); mpfr_mul(r18138, r18138, r18136, MPFR_RNDN);
        mpfr_mul(r18139, r18137, r18138, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r18142, r18136, r18141, MPFR_RNDN);
        mpfr_fma(r18143, r18140, r18138, r18142, MPFR_RNDN);
        mpfr_fma(r18144, r18139, r18136, r18143, MPFR_RNDN);
        ;
        mpfr_mul(r18146, r18136, r18145, MPFR_RNDN);
        mpfr_mul(r18147, r18133, r18146, MPFR_RNDN);
        mpfr_mul(r18148, r18136, r18136, MPFR_RNDN);
        mpfr_mul(r18149, r18148, r18148, MPFR_RNDN); mpfr_mul(r18149, r18149, r18148, MPFR_RNDN);
        mpfr_mul(r18150, r18147, r18149, MPFR_RNDN);
        mpfr_fma(r18151, r18133, r18144, r18150, MPFR_RNDN);
        mpfr_abs(r18152, r18151, MPFR_RNDN);
        return mpfr_get_d(r18152, MPFR_RNDN);
}

static mpfr_t r18153, r18154, r18155, r18156, r18157, r18158, r18159, r18160, r18161, r18162, r18163, r18164, r18165, r18166, r18167, r18168, r18169, r18170, r18171, r18172, r18173, r18174, r18175;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18153, "1", 10, MPFR_RNDN);
        mpfr_init(r18154);
        mpfr_init(r18155);
        mpfr_init(r18156);
        mpfr_init_set_str(r18157, "1/5", 10, MPFR_RNDN);
        mpfr_init(r18158);
        mpfr_init(r18159);
        mpfr_init(r18160);
        mpfr_init(r18161);
        mpfr_init(r18162);
        mpfr_init_set_str(r18163, "2/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r18164, "2", 10, MPFR_RNDN);
        mpfr_init(r18165);
        mpfr_init(r18166);
        mpfr_init(r18167);
        mpfr_init_set_str(r18168, "1/21", 10, MPFR_RNDN);
        mpfr_init(r18169);
        mpfr_init(r18170);
        mpfr_init(r18171);
        mpfr_init(r18172);
        mpfr_init(r18173);
        mpfr_init(r18174);
        mpfr_init(r18175);
}

double f_dm(double x) {
        ;
        mpfr_const_pi(r18154, MPFR_RNDN);
        mpfr_div(r18155, r18153, r18154, MPFR_RNDN);
        mpfr_sqrt(r18156, r18155, MPFR_RNDN);
        ;
        mpfr_set_d(r18158, x, MPFR_RNDN);
        mpfr_abs(r18159, r18158, MPFR_RNDN);
        mpfr_mul(r18160, r18157, r18159, MPFR_RNDN);
        mpfr_mul(r18161, r18159, r18159, MPFR_RNDN); mpfr_mul(r18161, r18161, r18159, MPFR_RNDN);
        mpfr_mul(r18162, r18160, r18161, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r18165, r18159, r18164, MPFR_RNDN);
        mpfr_fma(r18166, r18163, r18161, r18165, MPFR_RNDN);
        mpfr_fma(r18167, r18162, r18159, r18166, MPFR_RNDN);
        ;
        mpfr_mul(r18169, r18159, r18168, MPFR_RNDN);
        mpfr_mul(r18170, r18156, r18169, MPFR_RNDN);
        mpfr_mul(r18171, r18159, r18159, MPFR_RNDN);
        mpfr_mul(r18172, r18171, r18171, MPFR_RNDN); mpfr_mul(r18172, r18172, r18171, MPFR_RNDN);
        mpfr_mul(r18173, r18170, r18172, MPFR_RNDN);
        mpfr_fma(r18174, r18156, r18167, r18173, MPFR_RNDN);
        mpfr_abs(r18175, r18174, MPFR_RNDN);
        return mpfr_get_d(r18175, MPFR_RNDN);
}