#include #include #include #include #include char *name = "simple fma test"; double f_if(float x, float y, float z) { float r19292 = x; float r19293 = y; float r19294 = z; float r19295 = fma(r19292, r19293, r19294); float r19296 = 1.0f; float r19297 = r19292 * r19293; float r19298 = r19297 + r19294; float r19299 = r19296 + r19298; float r19300 = r19295 - r19299; return r19300; } double f_id(double x, double y, double z) { double r19301 = x; double r19302 = y; double r19303 = z; double r19304 = fma(r19301, r19302, r19303); double r19305 = 1.0; double r19306 = r19301 * r19302; double r19307 = r19306 + r19303; double r19308 = r19305 + r19307; double r19309 = r19304 - r19308; return r19309; } double f_of(float x, float y, float z) { float r19310 = z; float r19311 = -6.793618503523319e-17f; bool r19312 = r19310 <= r19311; float r19313 = 1.0f; float r19314 = x; float r19315 = r19313 / r19314; float r19316 = y; float r19317 = r19313 / r19316; float r19318 = r19313 / r19310; float r19319 = fma(r19315, r19317, r19318); float r19320 = r19316 * r19314; float r19321 = r19313 / r19320; float r19322 = r19319 - r19321; float r19323 = r19322 - r19313; float r19324 = exp(r19323); float r19325 = log(r19324); float r19326 = 1.6241870003327005e+17f; bool r19327 = r19310 <= r19326; float r19328 = fma(r19314, r19316, r19310); float r19329 = r19314 * r19316; float r19330 = r19329 + r19310; float r19331 = cbrt(r19330); float r19332 = r19330 * r19330; float r19333 = cbrt(r19332); float r19334 = r19331 * r19333; float r19335 = r19313 + r19334; float r19336 = r19328 - r19335; float r19337 = r19327 ? r19336 : r19325; float r19338 = r19312 ? r19325 : r19337; return r19338; } double f_od(double x, double y, double z) { double r19339 = z; double r19340 = -6.793618503523319e-17; bool r19341 = r19339 <= r19340; double r19342 = 1.0; double r19343 = x; double r19344 = r19342 / r19343; double r19345 = y; double r19346 = r19342 / r19345; double r19347 = r19342 / r19339; double r19348 = fma(r19344, r19346, r19347); double r19349 = r19345 * r19343; double r19350 = r19342 / r19349; double r19351 = r19348 - r19350; double r19352 = r19351 - r19342; double r19353 = exp(r19352); double r19354 = log(r19353); double r19355 = 1.6241870003327005e+17; bool r19356 = r19339 <= r19355; double r19357 = fma(r19343, r19345, r19339); double r19358 = r19343 * r19345; double r19359 = r19358 + r19339; double r19360 = cbrt(r19359); double r19361 = r19359 * r19359; double r19362 = cbrt(r19361); double r19363 = r19360 * r19362; double r19364 = r19342 + r19363; double r19365 = r19357 - r19364; double r19366 = r19356 ? r19365 : r19354; double r19367 = r19341 ? r19354 : r19366; return r19367; } 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 r19368, r19369, r19370, r19371, r19372, r19373, r19374, r19375, r19376; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init(r19368); mpfr_init(r19369); mpfr_init(r19370); mpfr_init(r19371); mpfr_init_set_str(r19372, "1", 10, MPFR_RNDN); mpfr_init(r19373); mpfr_init(r19374); mpfr_init(r19375); mpfr_init(r19376); } double f_im(double x, double y, double z) { mpfr_set_d(r19368, x, MPFR_RNDN); mpfr_set_d(r19369, y, MPFR_RNDN); mpfr_set_d(r19370, z, MPFR_RNDN); mpfr_fma(r19371, r19368, r19369, r19370, MPFR_RNDN); ; mpfr_mul(r19373, r19368, r19369, MPFR_RNDN); mpfr_add(r19374, r19373, r19370, MPFR_RNDN); mpfr_add(r19375, r19372, r19374, MPFR_RNDN); mpfr_sub(r19376, r19371, r19375, MPFR_RNDN); return mpfr_get_d(r19376, MPFR_RNDN); } static mpfr_t r19377, r19378, r19379, r19380, r19381, r19382, r19383, r19384, r19385, r19386, r19387, r19388, r19389, r19390, r19391, r19392, r19393, r19394, r19395, r19396, r19397, r19398, r19399, r19400, r19401, r19402, r19403, r19404, r19405; void setup_mpfr_f_fm() { mpfr_set_default_prec(2448); mpfr_init(r19377); mpfr_init_set_str(r19378, "-6.793618503523319e-17", 10, MPFR_RNDN); mpfr_init(r19379); mpfr_init_set_str(r19380, "1", 10, MPFR_RNDN); mpfr_init(r19381); mpfr_init(r19382); mpfr_init(r19383); mpfr_init(r19384); mpfr_init(r19385); mpfr_init(r19386); mpfr_init(r19387); mpfr_init(r19388); mpfr_init(r19389); mpfr_init(r19390); mpfr_init(r19391); mpfr_init(r19392); mpfr_init_set_str(r19393, "1.6241870003327005e+17", 10, MPFR_RNDN); mpfr_init(r19394); mpfr_init(r19395); mpfr_init(r19396); mpfr_init(r19397); mpfr_init(r19398); mpfr_init(r19399); mpfr_init(r19400); mpfr_init(r19401); mpfr_init(r19402); mpfr_init(r19403); mpfr_init(r19404); mpfr_init(r19405); } double f_fm(double x, double y, double z) { mpfr_set_d(r19377, z, MPFR_RNDN); ; mpfr_set_si(r19379, mpfr_cmp(r19377, r19378) <= 0, MPFR_RNDN); ; mpfr_set_d(r19381, x, MPFR_RNDN); mpfr_div(r19382, r19380, r19381, MPFR_RNDN); mpfr_set_d(r19383, y, MPFR_RNDN); mpfr_div(r19384, r19380, r19383, MPFR_RNDN); mpfr_div(r19385, r19380, r19377, MPFR_RNDN); mpfr_fma(r19386, r19382, r19384, r19385, MPFR_RNDN); mpfr_mul(r19387, r19383, r19381, MPFR_RNDN); mpfr_div(r19388, r19380, r19387, MPFR_RNDN); mpfr_sub(r19389, r19386, r19388, MPFR_RNDN); mpfr_sub(r19390, r19389, r19380, MPFR_RNDN); mpfr_exp(r19391, r19390, MPFR_RNDN); mpfr_log(r19392, r19391, MPFR_RNDN); ; mpfr_set_si(r19394, mpfr_cmp(r19377, r19393) <= 0, MPFR_RNDN); mpfr_fma(r19395, r19381, r19383, r19377, MPFR_RNDN); mpfr_mul(r19396, r19381, r19383, MPFR_RNDN); mpfr_add(r19397, r19396, r19377, MPFR_RNDN); mpfr_cbrt(r19398, r19397, MPFR_RNDN); mpfr_mul(r19399, r19397, r19397, MPFR_RNDN); mpfr_cbrt(r19400, r19399, MPFR_RNDN); mpfr_mul(r19401, r19398, r19400, MPFR_RNDN); mpfr_add(r19402, r19380, r19401, MPFR_RNDN); mpfr_sub(r19403, r19395, r19402, MPFR_RNDN); if (mpfr_get_si(r19394, MPFR_RNDN)) { mpfr_set(r19404, r19403, MPFR_RNDN); } else { mpfr_set(r19404, r19392, MPFR_RNDN); }; if (mpfr_get_si(r19379, MPFR_RNDN)) { mpfr_set(r19405, r19392, MPFR_RNDN); } else { mpfr_set(r19405, r19404, MPFR_RNDN); }; return mpfr_get_d(r19405, MPFR_RNDN); } static mpfr_t r19406, r19407, r19408, r19409, r19410, r19411, r19412, r19413, r19414, r19415, r19416, r19417, r19418, r19419, r19420, r19421, r19422, r19423, r19424, r19425, r19426, r19427, r19428, r19429, r19430, r19431, r19432, r19433, r19434; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init(r19406); mpfr_init_set_str(r19407, "-6.793618503523319e-17", 10, MPFR_RNDN); mpfr_init(r19408); mpfr_init_set_str(r19409, "1", 10, MPFR_RNDN); mpfr_init(r19410); mpfr_init(r19411); mpfr_init(r19412); mpfr_init(r19413); mpfr_init(r19414); mpfr_init(r19415); mpfr_init(r19416); mpfr_init(r19417); mpfr_init(r19418); mpfr_init(r19419); mpfr_init(r19420); mpfr_init(r19421); mpfr_init_set_str(r19422, "1.6241870003327005e+17", 10, MPFR_RNDN); mpfr_init(r19423); mpfr_init(r19424); mpfr_init(r19425); mpfr_init(r19426); mpfr_init(r19427); mpfr_init(r19428); mpfr_init(r19429); mpfr_init(r19430); mpfr_init(r19431); mpfr_init(r19432); mpfr_init(r19433); mpfr_init(r19434); } double f_dm(double x, double y, double z) { mpfr_set_d(r19406, z, MPFR_RNDN); ; mpfr_set_si(r19408, mpfr_cmp(r19406, r19407) <= 0, MPFR_RNDN); ; mpfr_set_d(r19410, x, MPFR_RNDN); mpfr_div(r19411, r19409, r19410, MPFR_RNDN); mpfr_set_d(r19412, y, MPFR_RNDN); mpfr_div(r19413, r19409, r19412, MPFR_RNDN); mpfr_div(r19414, r19409, r19406, MPFR_RNDN); mpfr_fma(r19415, r19411, r19413, r19414, MPFR_RNDN); mpfr_mul(r19416, r19412, r19410, MPFR_RNDN); mpfr_div(r19417, r19409, r19416, MPFR_RNDN); mpfr_sub(r19418, r19415, r19417, MPFR_RNDN); mpfr_sub(r19419, r19418, r19409, MPFR_RNDN); mpfr_exp(r19420, r19419, MPFR_RNDN); mpfr_log(r19421, r19420, MPFR_RNDN); ; mpfr_set_si(r19423, mpfr_cmp(r19406, r19422) <= 0, MPFR_RNDN); mpfr_fma(r19424, r19410, r19412, r19406, MPFR_RNDN); mpfr_mul(r19425, r19410, r19412, MPFR_RNDN); mpfr_add(r19426, r19425, r19406, MPFR_RNDN); mpfr_cbrt(r19427, r19426, MPFR_RNDN); mpfr_mul(r19428, r19426, r19426, MPFR_RNDN); mpfr_cbrt(r19429, r19428, MPFR_RNDN); mpfr_mul(r19430, r19427, r19429, MPFR_RNDN); mpfr_add(r19431, r19409, r19430, MPFR_RNDN); mpfr_sub(r19432, r19424, r19431, MPFR_RNDN); if (mpfr_get_si(r19423, MPFR_RNDN)) { mpfr_set(r19433, r19432, MPFR_RNDN); } else { mpfr_set(r19433, r19421, MPFR_RNDN); }; if (mpfr_get_si(r19408, MPFR_RNDN)) { mpfr_set(r19434, r19421, MPFR_RNDN); } else { mpfr_set(r19434, r19433, MPFR_RNDN); }; return mpfr_get_d(r19434, MPFR_RNDN); }