#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 r24278 = 1;
        float r24279 = atan2(1.0, 0.0);
        float r24280 = sqrt(r24279);
        float r24281 = r24278 / r24280;
        float r24282 = 2;
        float r24283 = x;
        float r24284 = fabs(r24283);
        float r24285 = r24282 * r24284;
        float r24286 = 3;
        float r24287 = r24282 / r24286;
        float r24288 = r24284 * r24284;
        float r24289 = r24288 * r24284;
        float r24290 = r24287 * r24289;
        float r24291 = r24285 + r24290;
        float r24292 = 5;
        float r24293 = r24278 / r24292;
        float r24294 = r24289 * r24284;
        float r24295 = r24294 * r24284;
        float r24296 = r24293 * r24295;
        float r24297 = r24291 + r24296;
        float r24298 = 21;
        float r24299 = r24278 / r24298;
        float r24300 = r24295 * r24284;
        float r24301 = r24300 * r24284;
        float r24302 = r24299 * r24301;
        float r24303 = r24297 + r24302;
        float r24304 = r24281 * r24303;
        float r24305 = fabs(r24304);
        return r24305;
}

double f_id(double x) {
        double r24306 = 1;
        double r24307 = atan2(1.0, 0.0);
        double r24308 = sqrt(r24307);
        double r24309 = r24306 / r24308;
        double r24310 = 2;
        double r24311 = x;
        double r24312 = fabs(r24311);
        double r24313 = r24310 * r24312;
        double r24314 = 3;
        double r24315 = r24310 / r24314;
        double r24316 = r24312 * r24312;
        double r24317 = r24316 * r24312;
        double r24318 = r24315 * r24317;
        double r24319 = r24313 + r24318;
        double r24320 = 5;
        double r24321 = r24306 / r24320;
        double r24322 = r24317 * r24312;
        double r24323 = r24322 * r24312;
        double r24324 = r24321 * r24323;
        double r24325 = r24319 + r24324;
        double r24326 = 21;
        double r24327 = r24306 / r24326;
        double r24328 = r24323 * r24312;
        double r24329 = r24328 * r24312;
        double r24330 = r24327 * r24329;
        double r24331 = r24325 + r24330;
        double r24332 = r24309 * r24331;
        double r24333 = fabs(r24332);
        return r24333;
}


double f_of(float x) {
        float r24334 = 1;
        float r24335 = atan2(1.0, 0.0);
        float r24336 = r24334 / r24335;
        float r24337 = sqrt(r24336);
        float r24338 = 2/3;
        float r24339 = x;
        float r24340 = fabs(r24339);
        float r24341 = 3;
        float r24342 = pow(r24340, r24341);
        float r24343 = r24338 * r24342;
        float r24344 = 2;
        float r24345 = r24344 * r24340;
        float r24346 = 1/21;
        float r24347 = 7;
        float r24348 = pow(r24340, r24347);
        float r24349 = r24346 * r24348;
        float r24350 = 1/5;
        float r24351 = 5;
        float r24352 = pow(r24340, r24351);
        float r24353 = r24350 * r24352;
        float r24354 = r24349 + r24353;
        float r24355 = r24345 + r24354;
        float r24356 = r24343 + r24355;
        float r24357 = r24337 * r24356;
        float r24358 = fabs(r24357);
        return r24358;
}

double f_od(double x) {
        double r24359 = 1;
        double r24360 = atan2(1.0, 0.0);
        double r24361 = r24359 / r24360;
        double r24362 = sqrt(r24361);
        double r24363 = 2/3;
        double r24364 = x;
        double r24365 = fabs(r24364);
        double r24366 = 3;
        double r24367 = pow(r24365, r24366);
        double r24368 = r24363 * r24367;
        double r24369 = 2;
        double r24370 = r24369 * r24365;
        double r24371 = 1/21;
        double r24372 = 7;
        double r24373 = pow(r24365, r24372);
        double r24374 = r24371 * r24373;
        double r24375 = 1/5;
        double r24376 = 5;
        double r24377 = pow(r24365, r24376);
        double r24378 = r24375 * r24377;
        double r24379 = r24374 + r24378;
        double r24380 = r24370 + r24379;
        double r24381 = r24368 + r24380;
        double r24382 = r24362 * r24381;
        double r24383 = fabs(r24382);
        return r24383;
}

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 r24384, r24385, r24386, r24387, r24388, r24389, r24390, r24391, r24392, r24393, r24394, r24395, r24396, r24397, r24398, r24399, r24400, r24401, r24402, r24403, r24404, r24405, r24406, r24407, r24408, r24409, r24410, r24411;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r24384, "1", 10, MPFR_RNDN);
        mpfr_init(r24385);
        mpfr_init(r24386);
        mpfr_init(r24387);
        mpfr_init_set_str(r24388, "2", 10, MPFR_RNDN);
        mpfr_init(r24389);
        mpfr_init(r24390);
        mpfr_init(r24391);
        mpfr_init_set_str(r24392, "3", 10, MPFR_RNDN);
        mpfr_init(r24393);
        mpfr_init(r24394);
        mpfr_init(r24395);
        mpfr_init(r24396);
        mpfr_init(r24397);
        mpfr_init_set_str(r24398, "5", 10, MPFR_RNDN);
        mpfr_init(r24399);
        mpfr_init(r24400);
        mpfr_init(r24401);
        mpfr_init(r24402);
        mpfr_init(r24403);
        mpfr_init_set_str(r24404, "21", 10, MPFR_RNDN);
        mpfr_init(r24405);
        mpfr_init(r24406);
        mpfr_init(r24407);
        mpfr_init(r24408);
        mpfr_init(r24409);
        mpfr_init(r24410);
        mpfr_init(r24411);
}

double f_im(double x) {
        ;
        mpfr_const_pi(r24385, MPFR_RNDN);
        mpfr_sqrt(r24386, r24385, MPFR_RNDN);
        mpfr_div(r24387, r24384, r24386, MPFR_RNDN);
        ;
        mpfr_set_d(r24389, x, MPFR_RNDN);
        mpfr_abs(r24390, r24389, MPFR_RNDN);
        mpfr_mul(r24391, r24388, r24390, MPFR_RNDN);
        ;
        mpfr_div(r24393, r24388, r24392, MPFR_RNDN);
        mpfr_mul(r24394, r24390, r24390, MPFR_RNDN);
        mpfr_mul(r24395, r24394, r24390, MPFR_RNDN);
        mpfr_mul(r24396, r24393, r24395, MPFR_RNDN);
        mpfr_add(r24397, r24391, r24396, MPFR_RNDN);
        ;
        mpfr_div(r24399, r24384, r24398, MPFR_RNDN);
        mpfr_mul(r24400, r24395, r24390, MPFR_RNDN);
        mpfr_mul(r24401, r24400, r24390, MPFR_RNDN);
        mpfr_mul(r24402, r24399, r24401, MPFR_RNDN);
        mpfr_add(r24403, r24397, r24402, MPFR_RNDN);
        ;
        mpfr_div(r24405, r24384, r24404, MPFR_RNDN);
        mpfr_mul(r24406, r24401, r24390, MPFR_RNDN);
        mpfr_mul(r24407, r24406, r24390, MPFR_RNDN);
        mpfr_mul(r24408, r24405, r24407, MPFR_RNDN);
        mpfr_add(r24409, r24403, r24408, MPFR_RNDN);
        mpfr_mul(r24410, r24387, r24409, MPFR_RNDN);
        mpfr_abs(r24411, r24410, MPFR_RNDN);
        return mpfr_get_d(r24411, MPFR_RNDN);
}

static mpfr_t r24412, r24413, r24414, r24415, r24416, r24417, r24418, r24419, r24420, r24421, r24422, r24423, r24424, r24425, r24426, r24427, r24428, r24429, r24430, r24431, r24432, r24433, r24434, r24435, r24436;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r24412, "1", 10, MPFR_RNDN);
        mpfr_init(r24413);
        mpfr_init(r24414);
        mpfr_init(r24415);
        mpfr_init_set_str(r24416, "2/3", 10, MPFR_RNDN);
        mpfr_init(r24417);
        mpfr_init(r24418);
        mpfr_init_set_str(r24419, "3", 10, MPFR_RNDN);
        mpfr_init(r24420);
        mpfr_init(r24421);
        mpfr_init_set_str(r24422, "2", 10, MPFR_RNDN);
        mpfr_init(r24423);
        mpfr_init_set_str(r24424, "1/21", 10, MPFR_RNDN);
        mpfr_init_set_str(r24425, "7", 10, MPFR_RNDN);
        mpfr_init(r24426);
        mpfr_init(r24427);
        mpfr_init_set_str(r24428, "1/5", 10, MPFR_RNDN);
        mpfr_init_set_str(r24429, "5", 10, MPFR_RNDN);
        mpfr_init(r24430);
        mpfr_init(r24431);
        mpfr_init(r24432);
        mpfr_init(r24433);
        mpfr_init(r24434);
        mpfr_init(r24435);
        mpfr_init(r24436);
}

double f_fm(double x) {
        ;
        mpfr_const_pi(r24413, MPFR_RNDN);
        mpfr_div(r24414, r24412, r24413, MPFR_RNDN);
        mpfr_sqrt(r24415, r24414, MPFR_RNDN);
        ;
        mpfr_set_d(r24417, x, MPFR_RNDN);
        mpfr_abs(r24418, r24417, MPFR_RNDN);
        ;
        mpfr_pow(r24420, r24418, r24419, MPFR_RNDN);
        mpfr_mul(r24421, r24416, r24420, MPFR_RNDN);
        ;
        mpfr_mul(r24423, r24422, r24418, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r24426, r24418, r24425, MPFR_RNDN);
        mpfr_mul(r24427, r24424, r24426, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r24430, r24418, r24429, MPFR_RNDN);
        mpfr_mul(r24431, r24428, r24430, MPFR_RNDN);
        mpfr_add(r24432, r24427, r24431, MPFR_RNDN);
        mpfr_add(r24433, r24423, r24432, MPFR_RNDN);
        mpfr_add(r24434, r24421, r24433, MPFR_RNDN);
        mpfr_mul(r24435, r24415, r24434, MPFR_RNDN);
        mpfr_abs(r24436, r24435, MPFR_RNDN);
        return mpfr_get_d(r24436, MPFR_RNDN);
}

static mpfr_t r24437, r24438, r24439, r24440, r24441, r24442, r24443, r24444, r24445, r24446, r24447, r24448, r24449, r24450, r24451, r24452, r24453, r24454, r24455, r24456, r24457, r24458, r24459, r24460, r24461;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r24437, "1", 10, MPFR_RNDN);
        mpfr_init(r24438);
        mpfr_init(r24439);
        mpfr_init(r24440);
        mpfr_init_set_str(r24441, "2/3", 10, MPFR_RNDN);
        mpfr_init(r24442);
        mpfr_init(r24443);
        mpfr_init_set_str(r24444, "3", 10, MPFR_RNDN);
        mpfr_init(r24445);
        mpfr_init(r24446);
        mpfr_init_set_str(r24447, "2", 10, MPFR_RNDN);
        mpfr_init(r24448);
        mpfr_init_set_str(r24449, "1/21", 10, MPFR_RNDN);
        mpfr_init_set_str(r24450, "7", 10, MPFR_RNDN);
        mpfr_init(r24451);
        mpfr_init(r24452);
        mpfr_init_set_str(r24453, "1/5", 10, MPFR_RNDN);
        mpfr_init_set_str(r24454, "5", 10, MPFR_RNDN);
        mpfr_init(r24455);
        mpfr_init(r24456);
        mpfr_init(r24457);
        mpfr_init(r24458);
        mpfr_init(r24459);
        mpfr_init(r24460);
        mpfr_init(r24461);
}

double f_dm(double x) {
        ;
        mpfr_const_pi(r24438, MPFR_RNDN);
        mpfr_div(r24439, r24437, r24438, MPFR_RNDN);
        mpfr_sqrt(r24440, r24439, MPFR_RNDN);
        ;
        mpfr_set_d(r24442, x, MPFR_RNDN);
        mpfr_abs(r24443, r24442, MPFR_RNDN);
        ;
        mpfr_pow(r24445, r24443, r24444, MPFR_RNDN);
        mpfr_mul(r24446, r24441, r24445, MPFR_RNDN);
        ;
        mpfr_mul(r24448, r24447, r24443, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r24451, r24443, r24450, MPFR_RNDN);
        mpfr_mul(r24452, r24449, r24451, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r24455, r24443, r24454, MPFR_RNDN);
        mpfr_mul(r24456, r24453, r24455, MPFR_RNDN);
        mpfr_add(r24457, r24452, r24456, MPFR_RNDN);
        mpfr_add(r24458, r24448, r24457, MPFR_RNDN);
        mpfr_add(r24459, r24446, r24458, MPFR_RNDN);
        mpfr_mul(r24460, r24440, r24459, MPFR_RNDN);
        mpfr_abs(r24461, r24460, MPFR_RNDN);
        return mpfr_get_d(r24461, MPFR_RNDN);
}