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

char *name = "Asymptote C";

double f_if(float x) {
        float r24176 = x;
        float r24177 = 1;
        float r24178 = r24176 + r24177;
        float r24179 = r24176 / r24178;
        float r24180 = r24176 - r24177;
        float r24181 = r24178 / r24180;
        float r24182 = r24179 - r24181;
        return r24182;
}

double f_id(double x) {
        double r24183 = x;
        double r24184 = 1;
        double r24185 = r24183 + r24184;
        double r24186 = r24183 / r24185;
        double r24187 = r24183 - r24184;
        double r24188 = r24185 / r24187;
        double r24189 = r24186 - r24188;
        return r24189;
}


double f_of(float x) {
        float r24190 = x;
        float r24191 = -12542.503822746461;
        bool r24192 = r24190 <= r24191;
        float r24193 = 3;
        float r24194 = r24193 / r24190;
        float r24195 = 1;
        float r24196 = r24194 + r24195;
        float r24197 = -r24195;
        float r24198 = r24190 * r24190;
        float r24199 = r24197 / r24198;
        float r24200 = -r24194;
        float r24201 = fma(r24196, r24199, r24200);
        float r24202 = 14709.455020816105;
        bool r24203 = r24190 <= r24202;
        float r24204 = pow(r24190, r24193);
        float r24205 = pow(r24195, r24193);
        float r24206 = r24204 + r24205;
        float r24207 = r24190 / r24206;
        float r24208 = r24195 - r24190;
        float r24209 = r24198 + r24208;
        float r24210 = r24190 + r24195;
        float r24211 = r24190 - r24195;
        float r24212 = r24210 / r24211;
        float r24213 = -r24212;
        float r24214 = fma(r24207, r24209, r24213);
        float r24215 = r24203 ? r24214 : r24201;
        float r24216 = r24192 ? r24201 : r24215;
        return r24216;
}

double f_od(double x) {
        double r24217 = x;
        double r24218 = -12542.503822746461;
        bool r24219 = r24217 <= r24218;
        double r24220 = 3;
        double r24221 = r24220 / r24217;
        double r24222 = 1;
        double r24223 = r24221 + r24222;
        double r24224 = -r24222;
        double r24225 = r24217 * r24217;
        double r24226 = r24224 / r24225;
        double r24227 = -r24221;
        double r24228 = fma(r24223, r24226, r24227);
        double r24229 = 14709.455020816105;
        bool r24230 = r24217 <= r24229;
        double r24231 = pow(r24217, r24220);
        double r24232 = pow(r24222, r24220);
        double r24233 = r24231 + r24232;
        double r24234 = r24217 / r24233;
        double r24235 = r24222 - r24217;
        double r24236 = r24225 + r24235;
        double r24237 = r24217 + r24222;
        double r24238 = r24217 - r24222;
        double r24239 = r24237 / r24238;
        double r24240 = -r24239;
        double r24241 = fma(r24234, r24236, r24240);
        double r24242 = r24230 ? r24241 : r24228;
        double r24243 = r24219 ? r24228 : r24242;
        return r24243;
}

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 r24244, r24245, r24246, r24247, r24248, r24249, r24250;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r24244);
        mpfr_init_set_str(r24245, "1", 10, MPFR_RNDN);
        mpfr_init(r24246);
        mpfr_init(r24247);
        mpfr_init(r24248);
        mpfr_init(r24249);
        mpfr_init(r24250);
}

double f_im(double x) {
        mpfr_set_d(r24244, x, MPFR_RNDN);
        ;
        mpfr_add(r24246, r24244, r24245, MPFR_RNDN);
        mpfr_div(r24247, r24244, r24246, MPFR_RNDN);
        mpfr_sub(r24248, r24244, r24245, MPFR_RNDN);
        mpfr_div(r24249, r24246, r24248, MPFR_RNDN);
        mpfr_sub(r24250, r24247, r24249, MPFR_RNDN);
        return mpfr_get_d(r24250, MPFR_RNDN);
}

static mpfr_t r24251, r24252, r24253, r24254, r24255, r24256, r24257, r24258, r24259, r24260, r24261, r24262, r24263, r24264, r24265, r24266, r24267, r24268, r24269, r24270, r24271, r24272, r24273, r24274, r24275, r24276, r24277;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r24251);
        mpfr_init_set_str(r24252, "-12542.503822746461", 10, MPFR_RNDN);
        mpfr_init(r24253);
        mpfr_init_set_str(r24254, "3", 10, MPFR_RNDN);
        mpfr_init(r24255);
        mpfr_init_set_str(r24256, "1", 10, MPFR_RNDN);
        mpfr_init(r24257);
        mpfr_init(r24258);
        mpfr_init(r24259);
        mpfr_init(r24260);
        mpfr_init(r24261);
        mpfr_init(r24262);
        mpfr_init_set_str(r24263, "14709.455020816105", 10, MPFR_RNDN);
        mpfr_init(r24264);
        mpfr_init(r24265);
        mpfr_init(r24266);
        mpfr_init(r24267);
        mpfr_init(r24268);
        mpfr_init(r24269);
        mpfr_init(r24270);
        mpfr_init(r24271);
        mpfr_init(r24272);
        mpfr_init(r24273);
        mpfr_init(r24274);
        mpfr_init(r24275);
        mpfr_init(r24276);
        mpfr_init(r24277);
}

double f_fm(double x) {
        mpfr_set_d(r24251, x, MPFR_RNDN);
        ;
        mpfr_set_si(r24253, mpfr_cmp(r24251, r24252) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r24255, r24254, r24251, MPFR_RNDN);
        ;
        mpfr_add(r24257, r24255, r24256, MPFR_RNDN);
        mpfr_neg(r24258, r24256, MPFR_RNDN);
        mpfr_mul(r24259, r24251, r24251, MPFR_RNDN);
        mpfr_div(r24260, r24258, r24259, MPFR_RNDN);
        mpfr_neg(r24261, r24255, MPFR_RNDN);
        mpfr_fma(r24262, r24257, r24260, r24261, MPFR_RNDN);
        ;
        mpfr_set_si(r24264, mpfr_cmp(r24251, r24263) <= 0, MPFR_RNDN);
        mpfr_pow(r24265, r24251, r24254, MPFR_RNDN);
        mpfr_pow(r24266, r24256, r24254, MPFR_RNDN);
        mpfr_add(r24267, r24265, r24266, MPFR_RNDN);
        mpfr_div(r24268, r24251, r24267, MPFR_RNDN);
        mpfr_sub(r24269, r24256, r24251, MPFR_RNDN);
        mpfr_add(r24270, r24259, r24269, MPFR_RNDN);
        mpfr_add(r24271, r24251, r24256, MPFR_RNDN);
        mpfr_sub(r24272, r24251, r24256, MPFR_RNDN);
        mpfr_div(r24273, r24271, r24272, MPFR_RNDN);
        mpfr_neg(r24274, r24273, MPFR_RNDN);
        mpfr_fma(r24275, r24268, r24270, r24274, MPFR_RNDN);
        if (mpfr_get_si(r24264, MPFR_RNDN)) { mpfr_set(r24276, r24275, MPFR_RNDN); } else { mpfr_set(r24276, r24262, MPFR_RNDN); };
        if (mpfr_get_si(r24253, MPFR_RNDN)) { mpfr_set(r24277, r24262, MPFR_RNDN); } else { mpfr_set(r24277, r24276, MPFR_RNDN); };
        return mpfr_get_d(r24277, MPFR_RNDN);
}

static mpfr_t r24278, r24279, r24280, r24281, r24282, r24283, r24284, r24285, r24286, r24287, r24288, r24289, r24290, r24291, r24292, r24293, r24294, r24295, r24296, r24297, r24298, r24299, r24300, r24301, r24302, r24303, r24304;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r24278);
        mpfr_init_set_str(r24279, "-12542.503822746461", 10, MPFR_RNDN);
        mpfr_init(r24280);
        mpfr_init_set_str(r24281, "3", 10, MPFR_RNDN);
        mpfr_init(r24282);
        mpfr_init_set_str(r24283, "1", 10, MPFR_RNDN);
        mpfr_init(r24284);
        mpfr_init(r24285);
        mpfr_init(r24286);
        mpfr_init(r24287);
        mpfr_init(r24288);
        mpfr_init(r24289);
        mpfr_init_set_str(r24290, "14709.455020816105", 10, MPFR_RNDN);
        mpfr_init(r24291);
        mpfr_init(r24292);
        mpfr_init(r24293);
        mpfr_init(r24294);
        mpfr_init(r24295);
        mpfr_init(r24296);
        mpfr_init(r24297);
        mpfr_init(r24298);
        mpfr_init(r24299);
        mpfr_init(r24300);
        mpfr_init(r24301);
        mpfr_init(r24302);
        mpfr_init(r24303);
        mpfr_init(r24304);
}

double f_dm(double x) {
        mpfr_set_d(r24278, x, MPFR_RNDN);
        ;
        mpfr_set_si(r24280, mpfr_cmp(r24278, r24279) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r24282, r24281, r24278, MPFR_RNDN);
        ;
        mpfr_add(r24284, r24282, r24283, MPFR_RNDN);
        mpfr_neg(r24285, r24283, MPFR_RNDN);
        mpfr_mul(r24286, r24278, r24278, MPFR_RNDN);
        mpfr_div(r24287, r24285, r24286, MPFR_RNDN);
        mpfr_neg(r24288, r24282, MPFR_RNDN);
        mpfr_fma(r24289, r24284, r24287, r24288, MPFR_RNDN);
        ;
        mpfr_set_si(r24291, mpfr_cmp(r24278, r24290) <= 0, MPFR_RNDN);
        mpfr_pow(r24292, r24278, r24281, MPFR_RNDN);
        mpfr_pow(r24293, r24283, r24281, MPFR_RNDN);
        mpfr_add(r24294, r24292, r24293, MPFR_RNDN);
        mpfr_div(r24295, r24278, r24294, MPFR_RNDN);
        mpfr_sub(r24296, r24283, r24278, MPFR_RNDN);
        mpfr_add(r24297, r24286, r24296, MPFR_RNDN);
        mpfr_add(r24298, r24278, r24283, MPFR_RNDN);
        mpfr_sub(r24299, r24278, r24283, MPFR_RNDN);
        mpfr_div(r24300, r24298, r24299, MPFR_RNDN);
        mpfr_neg(r24301, r24300, MPFR_RNDN);
        mpfr_fma(r24302, r24295, r24297, r24301, MPFR_RNDN);
        if (mpfr_get_si(r24291, MPFR_RNDN)) { mpfr_set(r24303, r24302, MPFR_RNDN); } else { mpfr_set(r24303, r24289, MPFR_RNDN); };
        if (mpfr_get_si(r24280, MPFR_RNDN)) { mpfr_set(r24304, r24289, MPFR_RNDN); } else { mpfr_set(r24304, r24303, MPFR_RNDN); };
        return mpfr_get_d(r24304, MPFR_RNDN);
}