#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 r18040 = 1.0f;
        float r18041 = atan2(1.0, 0.0);
        float r18042 = sqrt(r18041);
        float r18043 = r18040 / r18042;
        float r18044 = 2.0f;
        float r18045 = x;
        float r18046 = fabs(r18045);
        float r18047 = r18044 * r18046;
        float r18048 = 3.0f;
        float r18049 = r18044 / r18048;
        float r18050 = r18046 * r18046;
        float r18051 = r18050 * r18046;
        float r18052 = r18049 * r18051;
        float r18053 = r18047 + r18052;
        float r18054 = 5.0f;
        float r18055 = r18040 / r18054;
        float r18056 = r18051 * r18046;
        float r18057 = r18056 * r18046;
        float r18058 = r18055 * r18057;
        float r18059 = r18053 + r18058;
        float r18060 = 21.0f;
        float r18061 = r18040 / r18060;
        float r18062 = r18057 * r18046;
        float r18063 = r18062 * r18046;
        float r18064 = r18061 * r18063;
        float r18065 = r18059 + r18064;
        float r18066 = r18043 * r18065;
        float r18067 = fabs(r18066);
        return r18067;
}

double f_id(double x) {
        double r18068 = 1.0;
        double r18069 = atan2(1.0, 0.0);
        double r18070 = sqrt(r18069);
        double r18071 = r18068 / r18070;
        double r18072 = 2.0;
        double r18073 = x;
        double r18074 = fabs(r18073);
        double r18075 = r18072 * r18074;
        double r18076 = 3.0;
        double r18077 = r18072 / r18076;
        double r18078 = r18074 * r18074;
        double r18079 = r18078 * r18074;
        double r18080 = r18077 * r18079;
        double r18081 = r18075 + r18080;
        double r18082 = 5.0;
        double r18083 = r18068 / r18082;
        double r18084 = r18079 * r18074;
        double r18085 = r18084 * r18074;
        double r18086 = r18083 * r18085;
        double r18087 = r18081 + r18086;
        double r18088 = 21.0;
        double r18089 = r18068 / r18088;
        double r18090 = r18085 * r18074;
        double r18091 = r18090 * r18074;
        double r18092 = r18089 * r18091;
        double r18093 = r18087 + r18092;
        double r18094 = r18071 * r18093;
        double r18095 = fabs(r18094);
        return r18095;
}


double f_of(float x) {
        float r18096 = 1.0f;
        float r18097 = atan2(1.0, 0.0);
        float r18098 = r18096 / r18097;
        float r18099 = sqrt(r18098);
        float r18100 = 0.2f;
        float r18101 = x;
        float r18102 = fabs(r18101);
        float r18103 = r18102 * (r18102 * r18102);
        float r18104 = r18102 * r18103;
        float r18105 = r18100 * r18104;
        float r18106 = 0.6666666666666666f;
        float r18107 = 2.0f;
        float r18108 = r18107 * r18102;
        float r18109 = fma(r18106, r18103, r18108);
        float r18110 = fma(r18105, r18102, r18109);
        float r18111 = 0.047619047619047616f;
        float r18112 = 3.0f;
        float r18113 = r18112 + r18112;
        float r18114 = pow(r18102, r18113);
        float r18115 = r18114 * r18102;
        float r18116 = r18111 * r18115;
        float r18117 = r18110 + r18116;
        float r18118 = r18099 * r18117;
        float r18119 = fabs(r18118);
        return r18119;
}

double f_od(double x) {
        double r18120 = 1.0;
        double r18121 = atan2(1.0, 0.0);
        double r18122 = r18120 / r18121;
        double r18123 = sqrt(r18122);
        double r18124 = 0.2;
        double r18125 = x;
        double r18126 = fabs(r18125);
        double r18127 = r18126 * (r18126 * r18126);
        double r18128 = r18126 * r18127;
        double r18129 = r18124 * r18128;
        double r18130 = 0.6666666666666666;
        double r18131 = 2.0;
        double r18132 = r18131 * r18126;
        double r18133 = fma(r18130, r18127, r18132);
        double r18134 = fma(r18129, r18126, r18133);
        double r18135 = 0.047619047619047616;
        double r18136 = 3.0;
        double r18137 = r18136 + r18136;
        double r18138 = pow(r18126, r18137);
        double r18139 = r18138 * r18126;
        double r18140 = r18135 * r18139;
        double r18141 = r18134 + r18140;
        double r18142 = r18123 * r18141;
        double r18143 = fabs(r18142);
        return r18143;
}

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 r18144, r18145, r18146, r18147, r18148, r18149, r18150, r18151, r18152, r18153, r18154, r18155, r18156, r18157, r18158, r18159, r18160, r18161, r18162, r18163, r18164, r18165, r18166, r18167, r18168, r18169, r18170, r18171;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18144, "1", 10, MPFR_RNDN);
        mpfr_init(r18145);
        mpfr_init(r18146);
        mpfr_init(r18147);
        mpfr_init_set_str(r18148, "2", 10, MPFR_RNDN);
        mpfr_init(r18149);
        mpfr_init(r18150);
        mpfr_init(r18151);
        mpfr_init_set_str(r18152, "3", 10, MPFR_RNDN);
        mpfr_init(r18153);
        mpfr_init(r18154);
        mpfr_init(r18155);
        mpfr_init(r18156);
        mpfr_init(r18157);
        mpfr_init_set_str(r18158, "5", 10, MPFR_RNDN);
        mpfr_init(r18159);
        mpfr_init(r18160);
        mpfr_init(r18161);
        mpfr_init(r18162);
        mpfr_init(r18163);
        mpfr_init_set_str(r18164, "21", 10, MPFR_RNDN);
        mpfr_init(r18165);
        mpfr_init(r18166);
        mpfr_init(r18167);
        mpfr_init(r18168);
        mpfr_init(r18169);
        mpfr_init(r18170);
        mpfr_init(r18171);
}

double f_im(double x) {
        ;
        mpfr_const_pi(r18145, MPFR_RNDN);
        mpfr_sqrt(r18146, r18145, MPFR_RNDN);
        mpfr_div(r18147, r18144, r18146, MPFR_RNDN);
        ;
        mpfr_set_d(r18149, x, MPFR_RNDN);
        mpfr_abs(r18150, r18149, MPFR_RNDN);
        mpfr_mul(r18151, r18148, r18150, MPFR_RNDN);
        ;
        mpfr_div(r18153, r18148, r18152, MPFR_RNDN);
        mpfr_mul(r18154, r18150, r18150, MPFR_RNDN);
        mpfr_mul(r18155, r18154, r18150, MPFR_RNDN);
        mpfr_mul(r18156, r18153, r18155, MPFR_RNDN);
        mpfr_add(r18157, r18151, r18156, MPFR_RNDN);
        ;
        mpfr_div(r18159, r18144, r18158, MPFR_RNDN);
        mpfr_mul(r18160, r18155, r18150, MPFR_RNDN);
        mpfr_mul(r18161, r18160, r18150, MPFR_RNDN);
        mpfr_mul(r18162, r18159, r18161, MPFR_RNDN);
        mpfr_add(r18163, r18157, r18162, MPFR_RNDN);
        ;
        mpfr_div(r18165, r18144, r18164, MPFR_RNDN);
        mpfr_mul(r18166, r18161, r18150, MPFR_RNDN);
        mpfr_mul(r18167, r18166, r18150, MPFR_RNDN);
        mpfr_mul(r18168, r18165, r18167, MPFR_RNDN);
        mpfr_add(r18169, r18163, r18168, MPFR_RNDN);
        mpfr_mul(r18170, r18147, r18169, MPFR_RNDN);
        mpfr_abs(r18171, r18170, MPFR_RNDN);
        return mpfr_get_d(r18171, MPFR_RNDN);
}

static mpfr_t r18172, r18173, r18174, r18175, r18176, r18177, r18178, r18179, r18180, r18181, r18182, r18183, r18184, r18185, r18186, r18187, r18188, r18189, r18190, r18191, r18192, r18193, r18194, r18195;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18172, "1", 10, MPFR_RNDN);
        mpfr_init(r18173);
        mpfr_init(r18174);
        mpfr_init(r18175);
        mpfr_init_set_str(r18176, "1/5", 10, MPFR_RNDN);
        mpfr_init(r18177);
        mpfr_init(r18178);
        mpfr_init(r18179);
        mpfr_init(r18180);
        mpfr_init(r18181);
        mpfr_init_set_str(r18182, "2/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r18183, "2", 10, MPFR_RNDN);
        mpfr_init(r18184);
        mpfr_init(r18185);
        mpfr_init(r18186);
        mpfr_init_set_str(r18187, "1/21", 10, MPFR_RNDN);
        mpfr_init_set_str(r18188, "3", 10, MPFR_RNDN);
        mpfr_init(r18189);
        mpfr_init(r18190);
        mpfr_init(r18191);
        mpfr_init(r18192);
        mpfr_init(r18193);
        mpfr_init(r18194);
        mpfr_init(r18195);
}

double f_fm(double x) {
        ;
        mpfr_const_pi(r18173, MPFR_RNDN);
        mpfr_div(r18174, r18172, r18173, MPFR_RNDN);
        mpfr_sqrt(r18175, r18174, MPFR_RNDN);
        ;
        mpfr_set_d(r18177, x, MPFR_RNDN);
        mpfr_abs(r18178, r18177, MPFR_RNDN);
        mpfr_mul(r18179, r18178, r18178, MPFR_RNDN); mpfr_mul(r18179, r18179, r18178, MPFR_RNDN);
        mpfr_mul(r18180, r18178, r18179, MPFR_RNDN);
        mpfr_mul(r18181, r18176, r18180, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r18184, r18183, r18178, MPFR_RNDN);
        mpfr_fma(r18185, r18182, r18179, r18184, MPFR_RNDN);
        mpfr_fma(r18186, r18181, r18178, r18185, MPFR_RNDN);
        ;
        ;
        mpfr_add(r18189, r18188, r18188, MPFR_RNDN);
        mpfr_pow(r18190, r18178, r18189, MPFR_RNDN);
        mpfr_mul(r18191, r18190, r18178, MPFR_RNDN);
        mpfr_mul(r18192, r18187, r18191, MPFR_RNDN);
        mpfr_add(r18193, r18186, r18192, MPFR_RNDN);
        mpfr_mul(r18194, r18175, r18193, MPFR_RNDN);
        mpfr_abs(r18195, r18194, MPFR_RNDN);
        return mpfr_get_d(r18195, MPFR_RNDN);
}

static mpfr_t r18196, r18197, r18198, r18199, r18200, r18201, r18202, r18203, r18204, r18205, r18206, r18207, r18208, r18209, r18210, r18211, r18212, r18213, r18214, r18215, r18216, r18217, r18218, r18219;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18196, "1", 10, MPFR_RNDN);
        mpfr_init(r18197);
        mpfr_init(r18198);
        mpfr_init(r18199);
        mpfr_init_set_str(r18200, "1/5", 10, MPFR_RNDN);
        mpfr_init(r18201);
        mpfr_init(r18202);
        mpfr_init(r18203);
        mpfr_init(r18204);
        mpfr_init(r18205);
        mpfr_init_set_str(r18206, "2/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r18207, "2", 10, MPFR_RNDN);
        mpfr_init(r18208);
        mpfr_init(r18209);
        mpfr_init(r18210);
        mpfr_init_set_str(r18211, "1/21", 10, MPFR_RNDN);
        mpfr_init_set_str(r18212, "3", 10, MPFR_RNDN);
        mpfr_init(r18213);
        mpfr_init(r18214);
        mpfr_init(r18215);
        mpfr_init(r18216);
        mpfr_init(r18217);
        mpfr_init(r18218);
        mpfr_init(r18219);
}

double f_dm(double x) {
        ;
        mpfr_const_pi(r18197, MPFR_RNDN);
        mpfr_div(r18198, r18196, r18197, MPFR_RNDN);
        mpfr_sqrt(r18199, r18198, MPFR_RNDN);
        ;
        mpfr_set_d(r18201, x, MPFR_RNDN);
        mpfr_abs(r18202, r18201, MPFR_RNDN);
        mpfr_mul(r18203, r18202, r18202, MPFR_RNDN); mpfr_mul(r18203, r18203, r18202, MPFR_RNDN);
        mpfr_mul(r18204, r18202, r18203, MPFR_RNDN);
        mpfr_mul(r18205, r18200, r18204, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r18208, r18207, r18202, MPFR_RNDN);
        mpfr_fma(r18209, r18206, r18203, r18208, MPFR_RNDN);
        mpfr_fma(r18210, r18205, r18202, r18209, MPFR_RNDN);
        ;
        ;
        mpfr_add(r18213, r18212, r18212, MPFR_RNDN);
        mpfr_pow(r18214, r18202, r18213, MPFR_RNDN);
        mpfr_mul(r18215, r18214, r18202, MPFR_RNDN);
        mpfr_mul(r18216, r18211, r18215, MPFR_RNDN);
        mpfr_add(r18217, r18210, r18216, MPFR_RNDN);
        mpfr_mul(r18218, r18199, r18217, MPFR_RNDN);
        mpfr_abs(r18219, r18218, MPFR_RNDN);
        return mpfr_get_d(r18219, MPFR_RNDN);
}