#include #include #include #include #include char *name = "Harley's example"; double f_if(float c_p, float c_n, float t, float s) { float r25312 = 1; float r25313 = s; float r25314 = -r25313; float r25315 = exp(r25314); float r25316 = r25312 + r25315; float r25317 = r25312 / r25316; float r25318 = c_p; float r25319 = pow(r25317, r25318); float r25320 = r25312 - r25317; float r25321 = c_n; float r25322 = pow(r25320, r25321); float r25323 = r25319 * r25322; float r25324 = t; float r25325 = -r25324; float r25326 = exp(r25325); float r25327 = r25312 + r25326; float r25328 = r25312 / r25327; float r25329 = pow(r25328, r25318); float r25330 = r25312 - r25328; float r25331 = pow(r25330, r25321); float r25332 = r25329 * r25331; float r25333 = r25323 / r25332; return r25333; } double f_id(double c_p, double c_n, double t, double s) { double r25334 = 1; double r25335 = s; double r25336 = -r25335; double r25337 = exp(r25336); double r25338 = r25334 + r25337; double r25339 = r25334 / r25338; double r25340 = c_p; double r25341 = pow(r25339, r25340); double r25342 = r25334 - r25339; double r25343 = c_n; double r25344 = pow(r25342, r25343); double r25345 = r25341 * r25344; double r25346 = t; double r25347 = -r25346; double r25348 = exp(r25347); double r25349 = r25334 + r25348; double r25350 = r25334 / r25349; double r25351 = pow(r25350, r25340); double r25352 = r25334 - r25350; double r25353 = pow(r25352, r25343); double r25354 = r25351 * r25353; double r25355 = r25345 / r25354; return r25355; } double f_of(float c_p, float c_n, float t, float s) { float r25356 = 1; float r25357 = s; float r25358 = -r25357; float r25359 = exp(r25358); float r25360 = r25356 + r25359; float r25361 = r25356 / r25360; float r25362 = r25356 - r25361; float r25363 = log(r25362); float r25364 = t; float r25365 = -r25364; float r25366 = exp(r25365); float r25367 = r25366 + r25356; float r25368 = r25356 / r25367; float r25369 = r25356 - r25368; float r25370 = log(r25369); float r25371 = r25363 - r25370; float r25372 = c_n; float r25373 = r25371 * r25372; float r25374 = log(r25360); float r25375 = log(r25367); float r25376 = r25374 - r25375; float r25377 = c_p; float r25378 = r25376 * r25377; float r25379 = r25373 - r25378; float r25380 = exp(r25379); return r25380; } double f_od(double c_p, double c_n, double t, double s) { double r25381 = 1; double r25382 = s; double r25383 = -r25382; double r25384 = exp(r25383); double r25385 = r25381 + r25384; double r25386 = r25381 / r25385; double r25387 = r25381 - r25386; double r25388 = log(r25387); double r25389 = t; double r25390 = -r25389; double r25391 = exp(r25390); double r25392 = r25391 + r25381; double r25393 = r25381 / r25392; double r25394 = r25381 - r25393; double r25395 = log(r25394); double r25396 = r25388 - r25395; double r25397 = c_n; double r25398 = r25396 * r25397; double r25399 = log(r25385); double r25400 = log(r25392); double r25401 = r25399 - r25400; double r25402 = c_p; double r25403 = r25401 * r25402; double r25404 = r25398 - r25403; double r25405 = exp(r25404); return r25405; } 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 r25406, r25407, r25408, r25409, r25410, r25411, r25412, r25413, r25414, r25415, r25416, r25417, r25418, r25419, r25420, r25421, r25422, r25423, r25424, r25425, r25426, r25427; void setup_mpfr_f_im() { mpfr_set_default_prec(2192); mpfr_init_set_str(r25406, "1", 10, MPFR_RNDN); mpfr_init(r25407); mpfr_init(r25408); mpfr_init(r25409); mpfr_init(r25410); mpfr_init(r25411); mpfr_init(r25412); mpfr_init(r25413); mpfr_init(r25414); mpfr_init(r25415); mpfr_init(r25416); mpfr_init(r25417); mpfr_init(r25418); mpfr_init(r25419); mpfr_init(r25420); mpfr_init(r25421); mpfr_init(r25422); mpfr_init(r25423); mpfr_init(r25424); mpfr_init(r25425); mpfr_init(r25426); mpfr_init(r25427); } double f_im(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r25407, s, MPFR_RNDN); mpfr_neg(r25408, r25407, MPFR_RNDN); mpfr_exp(r25409, r25408, MPFR_RNDN); mpfr_add(r25410, r25406, r25409, MPFR_RNDN); mpfr_div(r25411, r25406, r25410, MPFR_RNDN); mpfr_set_d(r25412, c_p, MPFR_RNDN); mpfr_pow(r25413, r25411, r25412, MPFR_RNDN); mpfr_sub(r25414, r25406, r25411, MPFR_RNDN); mpfr_set_d(r25415, c_n, MPFR_RNDN); mpfr_pow(r25416, r25414, r25415, MPFR_RNDN); mpfr_mul(r25417, r25413, r25416, MPFR_RNDN); mpfr_set_d(r25418, t, MPFR_RNDN); mpfr_neg(r25419, r25418, MPFR_RNDN); mpfr_exp(r25420, r25419, MPFR_RNDN); mpfr_add(r25421, r25406, r25420, MPFR_RNDN); mpfr_div(r25422, r25406, r25421, MPFR_RNDN); mpfr_pow(r25423, r25422, r25412, MPFR_RNDN); mpfr_sub(r25424, r25406, r25422, MPFR_RNDN); mpfr_pow(r25425, r25424, r25415, MPFR_RNDN); mpfr_mul(r25426, r25423, r25425, MPFR_RNDN); mpfr_div(r25427, r25417, r25426, MPFR_RNDN); return mpfr_get_d(r25427, MPFR_RNDN); } static mpfr_t r25428, r25429, r25430, r25431, r25432, r25433, r25434, r25435, r25436, r25437, r25438, r25439, r25440, r25441, r25442, r25443, r25444, r25445, r25446, r25447, r25448, r25449, r25450, r25451, r25452; void setup_mpfr_f_fm() { mpfr_set_default_prec(2192); mpfr_init_set_str(r25428, "1", 10, MPFR_RNDN); mpfr_init(r25429); mpfr_init(r25430); mpfr_init(r25431); mpfr_init(r25432); mpfr_init(r25433); mpfr_init(r25434); mpfr_init(r25435); mpfr_init(r25436); mpfr_init(r25437); mpfr_init(r25438); mpfr_init(r25439); mpfr_init(r25440); mpfr_init(r25441); mpfr_init(r25442); mpfr_init(r25443); mpfr_init(r25444); mpfr_init(r25445); mpfr_init(r25446); mpfr_init(r25447); mpfr_init(r25448); mpfr_init(r25449); mpfr_init(r25450); mpfr_init(r25451); mpfr_init(r25452); } double f_fm(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r25429, s, MPFR_RNDN); mpfr_neg(r25430, r25429, MPFR_RNDN); mpfr_exp(r25431, r25430, MPFR_RNDN); mpfr_add(r25432, r25428, r25431, MPFR_RNDN); mpfr_div(r25433, r25428, r25432, MPFR_RNDN); mpfr_sub(r25434, r25428, r25433, MPFR_RNDN); mpfr_log(r25435, r25434, MPFR_RNDN); mpfr_set_d(r25436, t, MPFR_RNDN); mpfr_neg(r25437, r25436, MPFR_RNDN); mpfr_exp(r25438, r25437, MPFR_RNDN); mpfr_add(r25439, r25438, r25428, MPFR_RNDN); mpfr_div(r25440, r25428, r25439, MPFR_RNDN); mpfr_sub(r25441, r25428, r25440, MPFR_RNDN); mpfr_log(r25442, r25441, MPFR_RNDN); mpfr_sub(r25443, r25435, r25442, MPFR_RNDN); mpfr_set_d(r25444, c_n, MPFR_RNDN); mpfr_mul(r25445, r25443, r25444, MPFR_RNDN); mpfr_log(r25446, r25432, MPFR_RNDN); mpfr_log(r25447, r25439, MPFR_RNDN); mpfr_sub(r25448, r25446, r25447, MPFR_RNDN); mpfr_set_d(r25449, c_p, MPFR_RNDN); mpfr_mul(r25450, r25448, r25449, MPFR_RNDN); mpfr_sub(r25451, r25445, r25450, MPFR_RNDN); mpfr_exp(r25452, r25451, MPFR_RNDN); return mpfr_get_d(r25452, MPFR_RNDN); } static mpfr_t r25453, r25454, r25455, r25456, r25457, r25458, r25459, r25460, r25461, r25462, r25463, r25464, r25465, r25466, r25467, r25468, r25469, r25470, r25471, r25472, r25473, r25474, r25475, r25476, r25477; void setup_mpfr_f_dm() { mpfr_set_default_prec(2192); mpfr_init_set_str(r25453, "1", 10, MPFR_RNDN); mpfr_init(r25454); mpfr_init(r25455); mpfr_init(r25456); mpfr_init(r25457); mpfr_init(r25458); mpfr_init(r25459); mpfr_init(r25460); mpfr_init(r25461); mpfr_init(r25462); mpfr_init(r25463); mpfr_init(r25464); mpfr_init(r25465); mpfr_init(r25466); mpfr_init(r25467); mpfr_init(r25468); mpfr_init(r25469); mpfr_init(r25470); mpfr_init(r25471); mpfr_init(r25472); mpfr_init(r25473); mpfr_init(r25474); mpfr_init(r25475); mpfr_init(r25476); mpfr_init(r25477); } double f_dm(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r25454, s, MPFR_RNDN); mpfr_neg(r25455, r25454, MPFR_RNDN); mpfr_exp(r25456, r25455, MPFR_RNDN); mpfr_add(r25457, r25453, r25456, MPFR_RNDN); mpfr_div(r25458, r25453, r25457, MPFR_RNDN); mpfr_sub(r25459, r25453, r25458, MPFR_RNDN); mpfr_log(r25460, r25459, MPFR_RNDN); mpfr_set_d(r25461, t, MPFR_RNDN); mpfr_neg(r25462, r25461, MPFR_RNDN); mpfr_exp(r25463, r25462, MPFR_RNDN); mpfr_add(r25464, r25463, r25453, MPFR_RNDN); mpfr_div(r25465, r25453, r25464, MPFR_RNDN); mpfr_sub(r25466, r25453, r25465, MPFR_RNDN); mpfr_log(r25467, r25466, MPFR_RNDN); mpfr_sub(r25468, r25460, r25467, MPFR_RNDN); mpfr_set_d(r25469, c_n, MPFR_RNDN); mpfr_mul(r25470, r25468, r25469, MPFR_RNDN); mpfr_log(r25471, r25457, MPFR_RNDN); mpfr_log(r25472, r25464, MPFR_RNDN); mpfr_sub(r25473, r25471, r25472, MPFR_RNDN); mpfr_set_d(r25474, c_p, MPFR_RNDN); mpfr_mul(r25475, r25473, r25474, MPFR_RNDN); mpfr_sub(r25476, r25470, r25475, MPFR_RNDN); mpfr_exp(r25477, r25476, MPFR_RNDN); return mpfr_get_d(r25477, MPFR_RNDN); }