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