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