#include #include #include #include #include char *name = "simple fma test"; double f_if(float x, float y, float z) { float r23041 = x; float r23042 = y; float r23043 = z; float r23044 = fma(r23041, r23042, r23043); float r23045 = 1; float r23046 = r23041 * r23042; float r23047 = r23046 + r23043; float r23048 = r23045 + r23047; float r23049 = r23044 - r23048; return r23049; } double f_id(double x, double y, double z) { double r23050 = x; double r23051 = y; double r23052 = z; double r23053 = fma(r23050, r23051, r23052); double r23054 = 1; double r23055 = r23050 * r23051; double r23056 = r23055 + r23052; double r23057 = r23054 + r23056; double r23058 = r23053 - r23057; return r23058; } double f_of(float x, float y, float z) { float r23059 = x; float r23060 = y; float r23061 = z; float r23062 = fma(r23059, r23060, r23061); float r23063 = r23062 - r23061; float r23064 = r23060 * r23059; float r23065 = r23063 - r23064; float r23066 = 1; float r23067 = r23065 - r23066; float r23068 = exp(r23067); float r23069 = log(r23068); float r23070 = exp(r23069); float r23071 = log(r23070); float r23072 = -1.0000000000010343; bool r23073 = r23071 <= r23072; float r23074 = r23062 - r23064; float r23075 = r23074 - r23061; float r23076 = r23075 - r23066; float r23077 = exp(r23076); float r23078 = log(r23077); float r23079 = -1.0; bool r23080 = r23071 <= r23079; float r23081 = r23080 ? r23071 : r23078; float r23082 = r23073 ? r23078 : r23081; return r23082; } double f_od(double x, double y, double z) { double r23083 = x; double r23084 = y; double r23085 = z; double r23086 = fma(r23083, r23084, r23085); double r23087 = r23086 - r23085; double r23088 = r23084 * r23083; double r23089 = r23087 - r23088; double r23090 = 1; double r23091 = r23089 - r23090; double r23092 = exp(r23091); double r23093 = log(r23092); double r23094 = exp(r23093); double r23095 = log(r23094); double r23096 = -1.0000000000010343; bool r23097 = r23095 <= r23096; double r23098 = r23086 - r23088; double r23099 = r23098 - r23085; double r23100 = r23099 - r23090; double r23101 = exp(r23100); double r23102 = log(r23101); double r23103 = -1.0; bool r23104 = r23095 <= r23103; double r23105 = r23104 ? r23095 : r23102; double r23106 = r23097 ? r23102 : r23105; return r23106; } 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 r23107, r23108, r23109, r23110, r23111, r23112, r23113, r23114, r23115; void setup_mpfr_f_im() { mpfr_set_default_prec(2384); mpfr_init(r23107); mpfr_init(r23108); mpfr_init(r23109); mpfr_init(r23110); mpfr_init_set_str(r23111, "1", 10, MPFR_RNDN); mpfr_init(r23112); mpfr_init(r23113); mpfr_init(r23114); mpfr_init(r23115); } double f_im(double x, double y, double z) { mpfr_set_d(r23107, x, MPFR_RNDN); mpfr_set_d(r23108, y, MPFR_RNDN); mpfr_set_d(r23109, z, MPFR_RNDN); mpfr_fma(r23110, r23107, r23108, r23109, MPFR_RNDN); ; mpfr_mul(r23112, r23107, r23108, MPFR_RNDN); mpfr_add(r23113, r23112, r23109, MPFR_RNDN); mpfr_add(r23114, r23111, r23113, MPFR_RNDN); mpfr_sub(r23115, r23110, r23114, MPFR_RNDN); return mpfr_get_d(r23115, MPFR_RNDN); } static mpfr_t r23116, r23117, r23118, r23119, r23120, r23121, r23122, r23123, r23124, r23125, r23126, r23127, r23128, r23129, r23130, r23131, r23132, r23133, r23134, r23135, r23136, r23137, r23138, r23139; void setup_mpfr_f_fm() { mpfr_set_default_prec(2384); mpfr_init(r23116); mpfr_init(r23117); mpfr_init(r23118); mpfr_init(r23119); mpfr_init(r23120); mpfr_init(r23121); mpfr_init(r23122); mpfr_init_set_str(r23123, "1", 10, MPFR_RNDN); mpfr_init(r23124); mpfr_init(r23125); mpfr_init(r23126); mpfr_init(r23127); mpfr_init(r23128); mpfr_init_set_str(r23129, "-1.0000000000010343", 10, MPFR_RNDN); mpfr_init(r23130); mpfr_init(r23131); mpfr_init(r23132); mpfr_init(r23133); mpfr_init(r23134); mpfr_init(r23135); mpfr_init_set_str(r23136, "-1.0", 10, MPFR_RNDN); mpfr_init(r23137); mpfr_init(r23138); mpfr_init(r23139); } double f_fm(double x, double y, double z) { mpfr_set_d(r23116, x, MPFR_RNDN); mpfr_set_d(r23117, y, MPFR_RNDN); mpfr_set_d(r23118, z, MPFR_RNDN); mpfr_fma(r23119, r23116, r23117, r23118, MPFR_RNDN); mpfr_sub(r23120, r23119, r23118, MPFR_RNDN); mpfr_mul(r23121, r23117, r23116, MPFR_RNDN); mpfr_sub(r23122, r23120, r23121, MPFR_RNDN); ; mpfr_sub(r23124, r23122, r23123, MPFR_RNDN); mpfr_exp(r23125, r23124, MPFR_RNDN); mpfr_log(r23126, r23125, MPFR_RNDN); mpfr_exp(r23127, r23126, MPFR_RNDN); mpfr_log(r23128, r23127, MPFR_RNDN); ; mpfr_set_si(r23130, mpfr_cmp(r23128, r23129) <= 0, MPFR_RNDN); mpfr_sub(r23131, r23119, r23121, MPFR_RNDN); mpfr_sub(r23132, r23131, r23118, MPFR_RNDN); mpfr_sub(r23133, r23132, r23123, MPFR_RNDN); mpfr_exp(r23134, r23133, MPFR_RNDN); mpfr_log(r23135, r23134, MPFR_RNDN); ; mpfr_set_si(r23137, mpfr_cmp(r23128, r23136) <= 0, MPFR_RNDN); if (mpfr_get_si(r23137, MPFR_RNDN)) { mpfr_set(r23138, r23128, MPFR_RNDN); } else { mpfr_set(r23138, r23135, MPFR_RNDN); }; if (mpfr_get_si(r23130, MPFR_RNDN)) { mpfr_set(r23139, r23135, MPFR_RNDN); } else { mpfr_set(r23139, r23138, MPFR_RNDN); }; return mpfr_get_d(r23139, MPFR_RNDN); } static mpfr_t r23140, r23141, r23142, r23143, r23144, r23145, r23146, r23147, r23148, r23149, r23150, r23151, r23152, r23153, r23154, r23155, r23156, r23157, r23158, r23159, r23160, r23161, r23162, r23163; void setup_mpfr_f_dm() { mpfr_set_default_prec(2384); mpfr_init(r23140); mpfr_init(r23141); mpfr_init(r23142); mpfr_init(r23143); mpfr_init(r23144); mpfr_init(r23145); mpfr_init(r23146); mpfr_init_set_str(r23147, "1", 10, MPFR_RNDN); mpfr_init(r23148); mpfr_init(r23149); mpfr_init(r23150); mpfr_init(r23151); mpfr_init(r23152); mpfr_init_set_str(r23153, "-1.0000000000010343", 10, MPFR_RNDN); mpfr_init(r23154); mpfr_init(r23155); mpfr_init(r23156); mpfr_init(r23157); mpfr_init(r23158); mpfr_init(r23159); mpfr_init_set_str(r23160, "-1.0", 10, MPFR_RNDN); mpfr_init(r23161); mpfr_init(r23162); mpfr_init(r23163); } double f_dm(double x, double y, double z) { mpfr_set_d(r23140, x, MPFR_RNDN); mpfr_set_d(r23141, y, MPFR_RNDN); mpfr_set_d(r23142, z, MPFR_RNDN); mpfr_fma(r23143, r23140, r23141, r23142, MPFR_RNDN); mpfr_sub(r23144, r23143, r23142, MPFR_RNDN); mpfr_mul(r23145, r23141, r23140, MPFR_RNDN); mpfr_sub(r23146, r23144, r23145, MPFR_RNDN); ; mpfr_sub(r23148, r23146, r23147, MPFR_RNDN); mpfr_exp(r23149, r23148, MPFR_RNDN); mpfr_log(r23150, r23149, MPFR_RNDN); mpfr_exp(r23151, r23150, MPFR_RNDN); mpfr_log(r23152, r23151, MPFR_RNDN); ; mpfr_set_si(r23154, mpfr_cmp(r23152, r23153) <= 0, MPFR_RNDN); mpfr_sub(r23155, r23143, r23145, MPFR_RNDN); mpfr_sub(r23156, r23155, r23142, MPFR_RNDN); mpfr_sub(r23157, r23156, r23147, MPFR_RNDN); mpfr_exp(r23158, r23157, MPFR_RNDN); mpfr_log(r23159, r23158, MPFR_RNDN); ; mpfr_set_si(r23161, mpfr_cmp(r23152, r23160) <= 0, MPFR_RNDN); if (mpfr_get_si(r23161, MPFR_RNDN)) { mpfr_set(r23162, r23152, MPFR_RNDN); } else { mpfr_set(r23162, r23159, MPFR_RNDN); }; if (mpfr_get_si(r23154, MPFR_RNDN)) { mpfr_set(r23163, r23159, MPFR_RNDN); } else { mpfr_set(r23163, r23162, MPFR_RNDN); }; return mpfr_get_d(r23163, MPFR_RNDN); }