#include #include #include #include #include char *name = "Asymptote C"; double f_if(float x) { float r24147 = x; float r24148 = 1; float r24149 = r24147 + r24148; float r24150 = r24147 / r24149; float r24151 = r24147 - r24148; float r24152 = r24149 / r24151; float r24153 = r24150 - r24152; return r24153; } double f_id(double x) { double r24154 = x; double r24155 = 1; double r24156 = r24154 + r24155; double r24157 = r24154 / r24156; double r24158 = r24154 - r24155; double r24159 = r24156 / r24158; double r24160 = r24157 - r24159; return r24160; } double f_of(float x) { float r24161 = x; float r24162 = 1; float r24163 = r24161 + r24162; float r24164 = r24161 / r24163; float r24165 = r24161 - r24162; float r24166 = r24163 / r24165; float r24167 = r24164 - r24166; float r24168 = 6.613776447506048e-06; bool r24169 = r24167 <= r24168; float r24170 = 3; float r24171 = r24170 / r24161; float r24172 = r24162 + r24171; float r24173 = -r24162; float r24174 = r24161 * r24161; float r24175 = r24173 / r24174; float r24176 = -r24171; float r24177 = fma(r24172, r24175, r24176); float r24178 = r24170 + r24161; float r24179 = fma(r24161, r24178, r24162); float r24180 = r24169 ? r24177 : r24179; return r24180; } double f_od(double x) { double r24181 = x; double r24182 = 1; double r24183 = r24181 + r24182; double r24184 = r24181 / r24183; double r24185 = r24181 - r24182; double r24186 = r24183 / r24185; double r24187 = r24184 - r24186; double r24188 = 6.613776447506048e-06; bool r24189 = r24187 <= r24188; double r24190 = 3; double r24191 = r24190 / r24181; double r24192 = r24182 + r24191; double r24193 = -r24182; double r24194 = r24181 * r24181; double r24195 = r24193 / r24194; double r24196 = -r24191; double r24197 = fma(r24192, r24195, r24196); double r24198 = r24190 + r24181; double r24199 = fma(r24181, r24198, r24182); double r24200 = r24189 ? r24197 : r24199; return r24200; } 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 r24201, r24202, r24203, r24204, r24205, r24206, r24207; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r24201); mpfr_init_set_str(r24202, "1", 10, MPFR_RNDN); mpfr_init(r24203); mpfr_init(r24204); mpfr_init(r24205); mpfr_init(r24206); mpfr_init(r24207); } double f_im(double x) { mpfr_set_d(r24201, x, MPFR_RNDN); ; mpfr_add(r24203, r24201, r24202, MPFR_RNDN); mpfr_div(r24204, r24201, r24203, MPFR_RNDN); mpfr_sub(r24205, r24201, r24202, MPFR_RNDN); mpfr_div(r24206, r24203, r24205, MPFR_RNDN); mpfr_sub(r24207, r24204, r24206, MPFR_RNDN); return mpfr_get_d(r24207, MPFR_RNDN); } static mpfr_t r24208, r24209, r24210, r24211, r24212, r24213, r24214, r24215, r24216, r24217, r24218, r24219, r24220, r24221, r24222, r24223, r24224, r24225, r24226, r24227; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24208); mpfr_init_set_str(r24209, "1", 10, MPFR_RNDN); mpfr_init(r24210); mpfr_init(r24211); mpfr_init(r24212); mpfr_init(r24213); mpfr_init(r24214); mpfr_init_set_str(r24215, "6.613776447506048e-06", 10, MPFR_RNDN); mpfr_init(r24216); mpfr_init_set_str(r24217, "3", 10, MPFR_RNDN); mpfr_init(r24218); mpfr_init(r24219); mpfr_init(r24220); mpfr_init(r24221); mpfr_init(r24222); mpfr_init(r24223); mpfr_init(r24224); mpfr_init(r24225); mpfr_init(r24226); mpfr_init(r24227); } double f_fm(double x) { mpfr_set_d(r24208, x, MPFR_RNDN); ; mpfr_add(r24210, r24208, r24209, MPFR_RNDN); mpfr_div(r24211, r24208, r24210, MPFR_RNDN); mpfr_sub(r24212, r24208, r24209, MPFR_RNDN); mpfr_div(r24213, r24210, r24212, MPFR_RNDN); mpfr_sub(r24214, r24211, r24213, MPFR_RNDN); ; mpfr_set_si(r24216, mpfr_cmp(r24214, r24215) <= 0, MPFR_RNDN); ; mpfr_div(r24218, r24217, r24208, MPFR_RNDN); mpfr_add(r24219, r24209, r24218, MPFR_RNDN); mpfr_neg(r24220, r24209, MPFR_RNDN); mpfr_mul(r24221, r24208, r24208, MPFR_RNDN); mpfr_div(r24222, r24220, r24221, MPFR_RNDN); mpfr_neg(r24223, r24218, MPFR_RNDN); mpfr_fma(r24224, r24219, r24222, r24223, MPFR_RNDN); mpfr_add(r24225, r24217, r24208, MPFR_RNDN); mpfr_fma(r24226, r24208, r24225, r24209, MPFR_RNDN); if (mpfr_get_si(r24216, MPFR_RNDN)) { mpfr_set(r24227, r24224, MPFR_RNDN); } else { mpfr_set(r24227, r24226, MPFR_RNDN); }; return mpfr_get_d(r24227, MPFR_RNDN); } static mpfr_t r24228, r24229, r24230, r24231, r24232, r24233, r24234, r24235, r24236, r24237, r24238, r24239, r24240, r24241, r24242, r24243, r24244, r24245, r24246, r24247; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24228); mpfr_init_set_str(r24229, "1", 10, MPFR_RNDN); mpfr_init(r24230); mpfr_init(r24231); mpfr_init(r24232); mpfr_init(r24233); mpfr_init(r24234); mpfr_init_set_str(r24235, "6.613776447506048e-06", 10, MPFR_RNDN); mpfr_init(r24236); mpfr_init_set_str(r24237, "3", 10, MPFR_RNDN); mpfr_init(r24238); mpfr_init(r24239); mpfr_init(r24240); mpfr_init(r24241); mpfr_init(r24242); mpfr_init(r24243); mpfr_init(r24244); mpfr_init(r24245); mpfr_init(r24246); mpfr_init(r24247); } double f_dm(double x) { mpfr_set_d(r24228, x, MPFR_RNDN); ; mpfr_add(r24230, r24228, r24229, MPFR_RNDN); mpfr_div(r24231, r24228, r24230, MPFR_RNDN); mpfr_sub(r24232, r24228, r24229, MPFR_RNDN); mpfr_div(r24233, r24230, r24232, MPFR_RNDN); mpfr_sub(r24234, r24231, r24233, MPFR_RNDN); ; mpfr_set_si(r24236, mpfr_cmp(r24234, r24235) <= 0, MPFR_RNDN); ; mpfr_div(r24238, r24237, r24228, MPFR_RNDN); mpfr_add(r24239, r24229, r24238, MPFR_RNDN); mpfr_neg(r24240, r24229, MPFR_RNDN); mpfr_mul(r24241, r24228, r24228, MPFR_RNDN); mpfr_div(r24242, r24240, r24241, MPFR_RNDN); mpfr_neg(r24243, r24238, MPFR_RNDN); mpfr_fma(r24244, r24239, r24242, r24243, MPFR_RNDN); mpfr_add(r24245, r24237, r24228, MPFR_RNDN); mpfr_fma(r24246, r24228, r24245, r24229, MPFR_RNDN); if (mpfr_get_si(r24236, MPFR_RNDN)) { mpfr_set(r24247, r24244, MPFR_RNDN); } else { mpfr_set(r24247, r24246, MPFR_RNDN); }; return mpfr_get_d(r24247, MPFR_RNDN); }