#include #include #include #include #include char *name = "Asymptote C"; double f_if(float x) { float r24395 = x; float r24396 = 1; float r24397 = r24395 + r24396; float r24398 = r24395 / r24397; float r24399 = r24395 - r24396; float r24400 = r24397 / r24399; float r24401 = r24398 - r24400; return r24401; } double f_id(double x) { double r24402 = x; double r24403 = 1; double r24404 = r24402 + r24403; double r24405 = r24402 / r24404; double r24406 = r24402 - r24403; double r24407 = r24404 / r24406; double r24408 = r24405 - r24407; return r24408; } double f_of(float x) { float r24409 = x; float r24410 = 1; float r24411 = r24409 + r24410; float r24412 = r24409 / r24411; float r24413 = r24409 - r24410; float r24414 = r24411 / r24413; float r24415 = r24412 - r24414; float r24416 = 2.18952811964851e-47; bool r24417 = r24415 <= r24416; float r24418 = 3; float r24419 = r24418 / r24409; float r24420 = r24410 + r24419; float r24421 = -r24410; float r24422 = r24409 * r24409; float r24423 = r24421 / r24422; float r24424 = -r24419; float r24425 = fma(r24420, r24423, r24424); float r24426 = r24410 / r24411; float r24427 = -r24414; float r24428 = fma(r24409, r24426, r24427); float r24429 = r24417 ? r24425 : r24428; return r24429; } double f_od(double x) { double r24430 = x; double r24431 = 1; double r24432 = r24430 + r24431; double r24433 = r24430 / r24432; double r24434 = r24430 - r24431; double r24435 = r24432 / r24434; double r24436 = r24433 - r24435; double r24437 = 2.18952811964851e-47; bool r24438 = r24436 <= r24437; double r24439 = 3; double r24440 = r24439 / r24430; double r24441 = r24431 + r24440; double r24442 = -r24431; double r24443 = r24430 * r24430; double r24444 = r24442 / r24443; double r24445 = -r24440; double r24446 = fma(r24441, r24444, r24445); double r24447 = r24431 / r24432; double r24448 = -r24435; double r24449 = fma(r24430, r24447, r24448); double r24450 = r24438 ? r24446 : r24449; return r24450; } 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 r24451, r24452, r24453, r24454, r24455, r24456, r24457; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r24451); mpfr_init_set_str(r24452, "1", 10, MPFR_RNDN); mpfr_init(r24453); mpfr_init(r24454); mpfr_init(r24455); mpfr_init(r24456); mpfr_init(r24457); } double f_im(double x) { mpfr_set_d(r24451, x, MPFR_RNDN); ; mpfr_add(r24453, r24451, r24452, MPFR_RNDN); mpfr_div(r24454, r24451, r24453, MPFR_RNDN); mpfr_sub(r24455, r24451, r24452, MPFR_RNDN); mpfr_div(r24456, r24453, r24455, MPFR_RNDN); mpfr_sub(r24457, r24454, r24456, MPFR_RNDN); return mpfr_get_d(r24457, MPFR_RNDN); } static mpfr_t r24458, r24459, r24460, r24461, r24462, r24463, r24464, r24465, r24466, r24467, r24468, r24469, r24470, r24471, r24472, r24473, r24474, r24475, r24476, r24477, r24478; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24458); mpfr_init_set_str(r24459, "1", 10, MPFR_RNDN); mpfr_init(r24460); mpfr_init(r24461); mpfr_init(r24462); mpfr_init(r24463); mpfr_init(r24464); mpfr_init_set_str(r24465, "2.18952811964851e-47", 10, MPFR_RNDN); mpfr_init(r24466); mpfr_init_set_str(r24467, "3", 10, MPFR_RNDN); mpfr_init(r24468); mpfr_init(r24469); mpfr_init(r24470); mpfr_init(r24471); mpfr_init(r24472); mpfr_init(r24473); mpfr_init(r24474); mpfr_init(r24475); mpfr_init(r24476); mpfr_init(r24477); mpfr_init(r24478); } double f_fm(double x) { mpfr_set_d(r24458, x, MPFR_RNDN); ; mpfr_add(r24460, r24458, r24459, MPFR_RNDN); mpfr_div(r24461, r24458, r24460, MPFR_RNDN); mpfr_sub(r24462, r24458, r24459, MPFR_RNDN); mpfr_div(r24463, r24460, r24462, MPFR_RNDN); mpfr_sub(r24464, r24461, r24463, MPFR_RNDN); ; mpfr_set_si(r24466, mpfr_cmp(r24464, r24465) <= 0, MPFR_RNDN); ; mpfr_div(r24468, r24467, r24458, MPFR_RNDN); mpfr_add(r24469, r24459, r24468, MPFR_RNDN); mpfr_neg(r24470, r24459, MPFR_RNDN); mpfr_mul(r24471, r24458, r24458, MPFR_RNDN); mpfr_div(r24472, r24470, r24471, MPFR_RNDN); mpfr_neg(r24473, r24468, MPFR_RNDN); mpfr_fma(r24474, r24469, r24472, r24473, MPFR_RNDN); mpfr_div(r24475, r24459, r24460, MPFR_RNDN); mpfr_neg(r24476, r24463, MPFR_RNDN); mpfr_fma(r24477, r24458, r24475, r24476, MPFR_RNDN); if (mpfr_get_si(r24466, MPFR_RNDN)) { mpfr_set(r24478, r24474, MPFR_RNDN); } else { mpfr_set(r24478, r24477, MPFR_RNDN); }; return mpfr_get_d(r24478, MPFR_RNDN); } static mpfr_t r24479, r24480, r24481, r24482, r24483, r24484, r24485, r24486, r24487, r24488, r24489, r24490, r24491, r24492, r24493, r24494, r24495, r24496, r24497, r24498, r24499; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24479); mpfr_init_set_str(r24480, "1", 10, MPFR_RNDN); mpfr_init(r24481); mpfr_init(r24482); mpfr_init(r24483); mpfr_init(r24484); mpfr_init(r24485); mpfr_init_set_str(r24486, "2.18952811964851e-47", 10, MPFR_RNDN); mpfr_init(r24487); mpfr_init_set_str(r24488, "3", 10, MPFR_RNDN); mpfr_init(r24489); mpfr_init(r24490); mpfr_init(r24491); mpfr_init(r24492); mpfr_init(r24493); mpfr_init(r24494); mpfr_init(r24495); mpfr_init(r24496); mpfr_init(r24497); mpfr_init(r24498); mpfr_init(r24499); } double f_dm(double x) { mpfr_set_d(r24479, x, MPFR_RNDN); ; mpfr_add(r24481, r24479, r24480, MPFR_RNDN); mpfr_div(r24482, r24479, r24481, MPFR_RNDN); mpfr_sub(r24483, r24479, r24480, MPFR_RNDN); mpfr_div(r24484, r24481, r24483, MPFR_RNDN); mpfr_sub(r24485, r24482, r24484, MPFR_RNDN); ; mpfr_set_si(r24487, mpfr_cmp(r24485, r24486) <= 0, MPFR_RNDN); ; mpfr_div(r24489, r24488, r24479, MPFR_RNDN); mpfr_add(r24490, r24480, r24489, MPFR_RNDN); mpfr_neg(r24491, r24480, MPFR_RNDN); mpfr_mul(r24492, r24479, r24479, MPFR_RNDN); mpfr_div(r24493, r24491, r24492, MPFR_RNDN); mpfr_neg(r24494, r24489, MPFR_RNDN); mpfr_fma(r24495, r24490, r24493, r24494, MPFR_RNDN); mpfr_div(r24496, r24480, r24481, MPFR_RNDN); mpfr_neg(r24497, r24484, MPFR_RNDN); mpfr_fma(r24498, r24479, r24496, r24497, MPFR_RNDN); if (mpfr_get_si(r24487, MPFR_RNDN)) { mpfr_set(r24499, r24495, MPFR_RNDN); } else { mpfr_set(r24499, r24498, MPFR_RNDN); }; return mpfr_get_d(r24499, MPFR_RNDN); }