#include #include #include #include #include char *name = "Cubic critical, wide range"; double f_if(float a, float b, float c, float __attribute__((unused)) d) { float r9888 = b; float r9889 = -r9888; float r9890 = r9888 * r9888; float r9891 = 3; float r9892 = a; float r9893 = r9891 * r9892; float r9894 = c; float r9895 = r9893 * r9894; float r9896 = r9890 - r9895; float r9897 = sqrt(r9896); float r9898 = r9889 + r9897; float r9899 = r9898 / r9893; return r9899; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r9900 = b; double r9901 = -r9900; double r9902 = r9900 * r9900; double r9903 = 3; double r9904 = a; double r9905 = r9903 * r9904; double r9906 = c; double r9907 = r9905 * r9906; double r9908 = r9902 - r9907; double r9909 = sqrt(r9908); double r9910 = r9901 + r9909; double r9911 = r9910 / r9905; return r9911; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r9912 = 3; float r9913 = c; float r9914 = r9912 * r9913; float r9915 = 1; float r9916 = r9914 / r9915; float r9917 = -r9912; float r9918 = a; float r9919 = r9918 * r9912; float r9920 = -r9913; float r9921 = b; float r9922 = r9921 * r9921; float r9923 = fma(r9919, r9920, r9922); float r9924 = sqrt(r9923); float r9925 = r9917 * r9921; float r9926 = fma(r9917, r9924, r9925); float r9927 = r9916 / r9926; return r9927; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r9928 = 3; double r9929 = c; double r9930 = r9928 * r9929; double r9931 = 1; double r9932 = r9930 / r9931; double r9933 = -r9928; double r9934 = a; double r9935 = r9934 * r9928; double r9936 = -r9929; double r9937 = b; double r9938 = r9937 * r9937; double r9939 = fma(r9935, r9936, r9938); double r9940 = sqrt(r9939); double r9941 = r9933 * r9937; double r9942 = fma(r9933, r9940, r9941); double r9943 = r9932 / r9942; return r9943; } 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 r9944, r9945, r9946, r9947, r9948, r9949, r9950, r9951, r9952, r9953, r9954, r9955; void setup_mpfr_f_im() { mpfr_set_default_prec(848); mpfr_init(r9944); mpfr_init(r9945); mpfr_init(r9946); mpfr_init_set_str(r9947, "3", 10, MPFR_RNDN); mpfr_init(r9948); mpfr_init(r9949); mpfr_init(r9950); mpfr_init(r9951); mpfr_init(r9952); mpfr_init(r9953); mpfr_init(r9954); mpfr_init(r9955); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r9944, b, MPFR_RNDN); mpfr_neg(r9945, r9944, MPFR_RNDN); mpfr_mul(r9946, r9944, r9944, MPFR_RNDN); ; mpfr_set_d(r9948, a, MPFR_RNDN); mpfr_mul(r9949, r9947, r9948, MPFR_RNDN); mpfr_set_d(r9950, c, MPFR_RNDN); mpfr_mul(r9951, r9949, r9950, MPFR_RNDN); mpfr_sub(r9952, r9946, r9951, MPFR_RNDN); mpfr_sqrt(r9953, r9952, MPFR_RNDN); mpfr_add(r9954, r9945, r9953, MPFR_RNDN); mpfr_div(r9955, r9954, r9949, MPFR_RNDN); return mpfr_get_d(r9955, MPFR_RNDN); } static mpfr_t r9956, r9957, r9958, r9959, r9960, r9961, r9962, r9963, r9964, r9965, r9966, r9967, r9968, r9969, r9970, r9971; void setup_mpfr_f_fm() { mpfr_set_default_prec(848); mpfr_init_set_str(r9956, "3", 10, MPFR_RNDN); mpfr_init(r9957); mpfr_init(r9958); mpfr_init_set_str(r9959, "1", 10, MPFR_RNDN); mpfr_init(r9960); mpfr_init(r9961); mpfr_init(r9962); mpfr_init(r9963); mpfr_init(r9964); mpfr_init(r9965); mpfr_init(r9966); mpfr_init(r9967); mpfr_init(r9968); mpfr_init(r9969); mpfr_init(r9970); mpfr_init(r9971); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { ; mpfr_set_d(r9957, c, MPFR_RNDN); mpfr_mul(r9958, r9956, r9957, MPFR_RNDN); ; mpfr_div(r9960, r9958, r9959, MPFR_RNDN); mpfr_neg(r9961, r9956, MPFR_RNDN); mpfr_set_d(r9962, a, MPFR_RNDN); mpfr_mul(r9963, r9962, r9956, MPFR_RNDN); mpfr_neg(r9964, r9957, MPFR_RNDN); mpfr_set_d(r9965, b, MPFR_RNDN); mpfr_mul(r9966, r9965, r9965, MPFR_RNDN); mpfr_fma(r9967, r9963, r9964, r9966, MPFR_RNDN); mpfr_sqrt(r9968, r9967, MPFR_RNDN); mpfr_mul(r9969, r9961, r9965, MPFR_RNDN); mpfr_fma(r9970, r9961, r9968, r9969, MPFR_RNDN); mpfr_div(r9971, r9960, r9970, MPFR_RNDN); return mpfr_get_d(r9971, MPFR_RNDN); } static mpfr_t r9972, r9973, r9974, r9975, r9976, r9977, r9978, r9979, r9980, r9981, r9982, r9983, r9984, r9985, r9986, r9987; void setup_mpfr_f_dm() { mpfr_set_default_prec(848); mpfr_init_set_str(r9972, "3", 10, MPFR_RNDN); mpfr_init(r9973); mpfr_init(r9974); mpfr_init_set_str(r9975, "1", 10, MPFR_RNDN); mpfr_init(r9976); mpfr_init(r9977); mpfr_init(r9978); mpfr_init(r9979); mpfr_init(r9980); mpfr_init(r9981); mpfr_init(r9982); mpfr_init(r9983); mpfr_init(r9984); mpfr_init(r9985); mpfr_init(r9986); mpfr_init(r9987); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { ; mpfr_set_d(r9973, c, MPFR_RNDN); mpfr_mul(r9974, r9972, r9973, MPFR_RNDN); ; mpfr_div(r9976, r9974, r9975, MPFR_RNDN); mpfr_neg(r9977, r9972, MPFR_RNDN); mpfr_set_d(r9978, a, MPFR_RNDN); mpfr_mul(r9979, r9978, r9972, MPFR_RNDN); mpfr_neg(r9980, r9973, MPFR_RNDN); mpfr_set_d(r9981, b, MPFR_RNDN); mpfr_mul(r9982, r9981, r9981, MPFR_RNDN); mpfr_fma(r9983, r9979, r9980, r9982, MPFR_RNDN); mpfr_sqrt(r9984, r9983, MPFR_RNDN); mpfr_mul(r9985, r9977, r9981, MPFR_RNDN); mpfr_fma(r9986, r9977, r9984, r9985, MPFR_RNDN); mpfr_div(r9987, r9976, r9986, MPFR_RNDN); return mpfr_get_d(r9987, MPFR_RNDN); }