#include #include #include #include #include char *name = "Harley's example"; double f_if(float c_p, float c_n, float t, float s) { float r24313 = 1; float r24314 = s; float r24315 = -r24314; float r24316 = exp(r24315); float r24317 = r24313 + r24316; float r24318 = r24313 / r24317; float r24319 = c_p; float r24320 = pow(r24318, r24319); float r24321 = r24313 - r24318; float r24322 = c_n; float r24323 = pow(r24321, r24322); float r24324 = r24320 * r24323; float r24325 = t; float r24326 = -r24325; float r24327 = exp(r24326); float r24328 = r24313 + r24327; float r24329 = r24313 / r24328; float r24330 = pow(r24329, r24319); float r24331 = r24313 - r24329; float r24332 = pow(r24331, r24322); float r24333 = r24330 * r24332; float r24334 = r24324 / r24333; return r24334; } double f_id(double c_p, double c_n, double t, double s) { double r24335 = 1; double r24336 = s; double r24337 = -r24336; double r24338 = exp(r24337); double r24339 = r24335 + r24338; double r24340 = r24335 / r24339; double r24341 = c_p; double r24342 = pow(r24340, r24341); double r24343 = r24335 - r24340; double r24344 = c_n; double r24345 = pow(r24343, r24344); double r24346 = r24342 * r24345; double r24347 = t; double r24348 = -r24347; double r24349 = exp(r24348); double r24350 = r24335 + r24349; double r24351 = r24335 / r24350; double r24352 = pow(r24351, r24341); double r24353 = r24335 - r24351; double r24354 = pow(r24353, r24344); double r24355 = r24352 * r24354; double r24356 = r24346 / r24355; return r24356; } double f_of(float c_p, float c_n, float t, float s) { float r24357 = 1; float r24358 = s; float r24359 = -r24358; float r24360 = exp(r24359); float r24361 = r24357 + r24360; float r24362 = log(r24361); float r24363 = t; float r24364 = -r24363; float r24365 = exp(r24364); float r24366 = r24365 + r24357; float r24367 = log(r24366); float r24368 = r24362 - r24367; float r24369 = c_p; float r24370 = -r24369; float r24371 = r24368 * r24370; float r24372 = r24357 / r24361; float r24373 = r24357 - r24372; float r24374 = log(r24373); float r24375 = r24357 / r24366; float r24376 = r24357 - r24375; float r24377 = log(r24376); float r24378 = r24374 - r24377; float r24379 = c_n; float r24380 = r24378 * r24379; float r24381 = r24371 + r24380; float r24382 = exp(r24381); return r24382; } double f_od(double c_p, double c_n, double t, double s) { double r24383 = 1; double r24384 = s; double r24385 = -r24384; double r24386 = exp(r24385); double r24387 = r24383 + r24386; double r24388 = log(r24387); double r24389 = t; double r24390 = -r24389; double r24391 = exp(r24390); double r24392 = r24391 + r24383; double r24393 = log(r24392); double r24394 = r24388 - r24393; double r24395 = c_p; double r24396 = -r24395; double r24397 = r24394 * r24396; double r24398 = r24383 / r24387; double r24399 = r24383 - r24398; double r24400 = log(r24399); double r24401 = r24383 / r24392; double r24402 = r24383 - r24401; double r24403 = log(r24402); double r24404 = r24400 - r24403; double r24405 = c_n; double r24406 = r24404 * r24405; double r24407 = r24397 + r24406; double r24408 = exp(r24407); return r24408; } 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 r24409, r24410, r24411, r24412, r24413, r24414, r24415, r24416, r24417, r24418, r24419, r24420, r24421, r24422, r24423, r24424, r24425, r24426, r24427, r24428, r24429, r24430; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init_set_str(r24409, "1", 10, MPFR_RNDN); mpfr_init(r24410); mpfr_init(r24411); mpfr_init(r24412); mpfr_init(r24413); mpfr_init(r24414); mpfr_init(r24415); mpfr_init(r24416); mpfr_init(r24417); mpfr_init(r24418); mpfr_init(r24419); mpfr_init(r24420); mpfr_init(r24421); mpfr_init(r24422); mpfr_init(r24423); mpfr_init(r24424); mpfr_init(r24425); mpfr_init(r24426); mpfr_init(r24427); mpfr_init(r24428); mpfr_init(r24429); mpfr_init(r24430); } double f_im(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r24410, s, MPFR_RNDN); mpfr_neg(r24411, r24410, MPFR_RNDN); mpfr_exp(r24412, r24411, MPFR_RNDN); mpfr_add(r24413, r24409, r24412, MPFR_RNDN); mpfr_div(r24414, r24409, r24413, MPFR_RNDN); mpfr_set_d(r24415, c_p, MPFR_RNDN); mpfr_pow(r24416, r24414, r24415, MPFR_RNDN); mpfr_sub(r24417, r24409, r24414, MPFR_RNDN); mpfr_set_d(r24418, c_n, MPFR_RNDN); mpfr_pow(r24419, r24417, r24418, MPFR_RNDN); mpfr_mul(r24420, r24416, r24419, MPFR_RNDN); mpfr_set_d(r24421, t, MPFR_RNDN); mpfr_neg(r24422, r24421, MPFR_RNDN); mpfr_exp(r24423, r24422, MPFR_RNDN); mpfr_add(r24424, r24409, r24423, MPFR_RNDN); mpfr_div(r24425, r24409, r24424, MPFR_RNDN); mpfr_pow(r24426, r24425, r24415, MPFR_RNDN); mpfr_sub(r24427, r24409, r24425, MPFR_RNDN); mpfr_pow(r24428, r24427, r24418, MPFR_RNDN); mpfr_mul(r24429, r24426, r24428, MPFR_RNDN); mpfr_div(r24430, r24420, r24429, MPFR_RNDN); return mpfr_get_d(r24430, MPFR_RNDN); } static mpfr_t r24431, r24432, r24433, r24434, r24435, r24436, r24437, r24438, r24439, r24440, r24441, r24442, r24443, r24444, r24445, r24446, r24447, r24448, r24449, r24450, r24451, r24452, r24453, r24454, r24455, r24456; void setup_mpfr_f_fm() { mpfr_set_default_prec(2448); mpfr_init_set_str(r24431, "1", 10, MPFR_RNDN); mpfr_init(r24432); mpfr_init(r24433); mpfr_init(r24434); mpfr_init(r24435); mpfr_init(r24436); mpfr_init(r24437); mpfr_init(r24438); mpfr_init(r24439); mpfr_init(r24440); mpfr_init(r24441); mpfr_init(r24442); mpfr_init(r24443); mpfr_init(r24444); mpfr_init(r24445); mpfr_init(r24446); mpfr_init(r24447); mpfr_init(r24448); mpfr_init(r24449); mpfr_init(r24450); mpfr_init(r24451); mpfr_init(r24452); mpfr_init(r24453); mpfr_init(r24454); mpfr_init(r24455); mpfr_init(r24456); } double f_fm(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r24432, s, MPFR_RNDN); mpfr_neg(r24433, r24432, MPFR_RNDN); mpfr_exp(r24434, r24433, MPFR_RNDN); mpfr_add(r24435, r24431, r24434, MPFR_RNDN); mpfr_log(r24436, r24435, MPFR_RNDN); mpfr_set_d(r24437, t, MPFR_RNDN); mpfr_neg(r24438, r24437, MPFR_RNDN); mpfr_exp(r24439, r24438, MPFR_RNDN); mpfr_add(r24440, r24439, r24431, MPFR_RNDN); mpfr_log(r24441, r24440, MPFR_RNDN); mpfr_sub(r24442, r24436, r24441, MPFR_RNDN); mpfr_set_d(r24443, c_p, MPFR_RNDN); mpfr_neg(r24444, r24443, MPFR_RNDN); mpfr_mul(r24445, r24442, r24444, MPFR_RNDN); mpfr_div(r24446, r24431, r24435, MPFR_RNDN); mpfr_sub(r24447, r24431, r24446, MPFR_RNDN); mpfr_log(r24448, r24447, MPFR_RNDN); mpfr_div(r24449, r24431, r24440, MPFR_RNDN); mpfr_sub(r24450, r24431, r24449, MPFR_RNDN); mpfr_log(r24451, r24450, MPFR_RNDN); mpfr_sub(r24452, r24448, r24451, MPFR_RNDN); mpfr_set_d(r24453, c_n, MPFR_RNDN); mpfr_mul(r24454, r24452, r24453, MPFR_RNDN); mpfr_add(r24455, r24445, r24454, MPFR_RNDN); mpfr_exp(r24456, r24455, MPFR_RNDN); return mpfr_get_d(r24456, MPFR_RNDN); } static mpfr_t r24457, r24458, r24459, r24460, r24461, r24462, r24463, r24464, r24465, r24466, r24467, r24468, r24469, r24470, r24471, r24472, r24473, r24474, r24475, r24476, r24477, r24478, r24479, r24480, r24481, r24482; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init_set_str(r24457, "1", 10, MPFR_RNDN); mpfr_init(r24458); mpfr_init(r24459); mpfr_init(r24460); mpfr_init(r24461); mpfr_init(r24462); mpfr_init(r24463); mpfr_init(r24464); mpfr_init(r24465); mpfr_init(r24466); mpfr_init(r24467); 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); mpfr_init(r24479); mpfr_init(r24480); mpfr_init(r24481); mpfr_init(r24482); } double f_dm(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r24458, s, MPFR_RNDN); mpfr_neg(r24459, r24458, MPFR_RNDN); mpfr_exp(r24460, r24459, MPFR_RNDN); mpfr_add(r24461, r24457, r24460, MPFR_RNDN); mpfr_log(r24462, r24461, MPFR_RNDN); mpfr_set_d(r24463, t, MPFR_RNDN); mpfr_neg(r24464, r24463, MPFR_RNDN); mpfr_exp(r24465, r24464, MPFR_RNDN); mpfr_add(r24466, r24465, r24457, MPFR_RNDN); mpfr_log(r24467, r24466, MPFR_RNDN); mpfr_sub(r24468, r24462, r24467, MPFR_RNDN); mpfr_set_d(r24469, c_p, MPFR_RNDN); mpfr_neg(r24470, r24469, MPFR_RNDN); mpfr_mul(r24471, r24468, r24470, MPFR_RNDN); mpfr_div(r24472, r24457, r24461, MPFR_RNDN); mpfr_sub(r24473, r24457, r24472, MPFR_RNDN); mpfr_log(r24474, r24473, MPFR_RNDN); mpfr_div(r24475, r24457, r24466, MPFR_RNDN); mpfr_sub(r24476, r24457, r24475, MPFR_RNDN); mpfr_log(r24477, r24476, MPFR_RNDN); mpfr_sub(r24478, r24474, r24477, MPFR_RNDN); mpfr_set_d(r24479, c_n, MPFR_RNDN); mpfr_mul(r24480, r24478, r24479, MPFR_RNDN); mpfr_add(r24481, r24471, r24480, MPFR_RNDN); mpfr_exp(r24482, r24481, MPFR_RNDN); return mpfr_get_d(r24482, MPFR_RNDN); }