#include #include #include #include #include char *name = "Toniolo and Linder, Equation (7)"; double f_if(float x, float l, float t) { float r7853 = 2.0f; float r7854 = sqrt(r7853); float r7855 = t; float r7856 = r7854 * r7855; float r7857 = x; float r7858 = 1.0f; float r7859 = r7857 + r7858; float r7860 = r7857 - r7858; float r7861 = r7859 / r7860; float r7862 = l; float r7863 = r7862 * r7862; float r7864 = r7855 * r7855; float r7865 = r7853 * r7864; float r7866 = r7863 + r7865; float r7867 = r7861 * r7866; float r7868 = r7867 - r7863; float r7869 = sqrt(r7868); float r7870 = r7856 / r7869; return r7870; } double f_id(double x, double l, double t) { double r7871 = 2.0; double r7872 = sqrt(r7871); double r7873 = t; double r7874 = r7872 * r7873; double r7875 = x; double r7876 = 1.0; double r7877 = r7875 + r7876; double r7878 = r7875 - r7876; double r7879 = r7877 / r7878; double r7880 = l; double r7881 = r7880 * r7880; double r7882 = r7873 * r7873; double r7883 = r7871 * r7882; double r7884 = r7881 + r7883; double r7885 = r7879 * r7884; double r7886 = r7885 - r7881; double r7887 = sqrt(r7886); double r7888 = r7874 / r7887; return r7888; } double f_of(float x, float l, float t) { float r7889 = t; float r7890 = -6.446840727595846e+56f; bool r7891 = r7889 <= r7890; float r7892 = 2.0f; float r7893 = sqrt(r7892); float r7894 = r7893 * r7889; float r7895 = x; float r7896 = r7895 * r7895; float r7897 = r7889 / r7896; float r7898 = 1.0f; float r7899 = r7898 / r7893; float r7900 = r7892 / r7893; float r7901 = r7899 - r7900; float r7902 = r7892 / r7895; float r7903 = -r7889; float r7904 = r7903 / r7893; float r7905 = r7903 * r7893; float r7906 = fma(r7902, r7904, r7905); float r7907 = fma(r7897, r7901, r7906); float r7908 = r7894 / r7907; float r7909 = -1.9073716957841407e-160f; bool r7910 = r7889 <= r7909; float r7911 = cbrt(r7893); float r7912 = fabs(r7911); float r7913 = sqrt(r7911); float r7914 = r7889 * r7913; float r7915 = r7912 * r7914; float r7916 = sqrt(r7893); float r7917 = r7915 * r7916; float r7918 = l; float r7919 = r7918 / r7895; float r7920 = r7889 * r7889; float r7921 = fma(r7919, r7918, r7920); float r7922 = r7889 / r7895; float r7923 = 4.0f; float r7924 = r7889 * r7923; float r7925 = r7922 * r7924; float r7926 = fma(r7892, r7921, r7925); float r7927 = sqrt(r7926); float r7928 = r7917 / r7927; float r7929 = -1.6577109360186375e-242f; bool r7930 = r7889 <= r7929; float r7931 = 1.7103591102507113e-246f; bool r7932 = r7889 <= r7931; float r7933 = 4.6275459220879563e-169f; bool r7934 = r7889 <= r7933; float r7935 = !r7934; float r7936 = 1.8364323594358619e+145f; bool r7937 = r7889 <= r7936; bool r7938 = r7935 && r7937; bool r7939 = r7932 || r7938; float r7940 = r7922 / r7895; float r7941 = r7940 + r7922; float r7942 = r7903 / r7896; float r7943 = fma(r7899, r7942, r7894); float r7944 = fma(r7900, r7941, r7943); float r7945 = r7894 / r7944; float r7946 = r7939 ? r7928 : r7945; float r7947 = r7930 ? r7908 : r7946; float r7948 = r7910 ? r7928 : r7947; float r7949 = r7891 ? r7908 : r7948; return r7949; } double f_od(double x, double l, double t) { double r7950 = t; double r7951 = -6.446840727595846e+56; bool r7952 = r7950 <= r7951; double r7953 = 2.0; double r7954 = sqrt(r7953); double r7955 = r7954 * r7950; double r7956 = x; double r7957 = r7956 * r7956; double r7958 = r7950 / r7957; double r7959 = 1.0; double r7960 = r7959 / r7954; double r7961 = r7953 / r7954; double r7962 = r7960 - r7961; double r7963 = r7953 / r7956; double r7964 = -r7950; double r7965 = r7964 / r7954; double r7966 = r7964 * r7954; double r7967 = fma(r7963, r7965, r7966); double r7968 = fma(r7958, r7962, r7967); double r7969 = r7955 / r7968; double r7970 = -1.9073716957841407e-160; bool r7971 = r7950 <= r7970; double r7972 = cbrt(r7954); double r7973 = fabs(r7972); double r7974 = sqrt(r7972); double r7975 = r7950 * r7974; double r7976 = r7973 * r7975; double r7977 = sqrt(r7954); double r7978 = r7976 * r7977; double r7979 = l; double r7980 = r7979 / r7956; double r7981 = r7950 * r7950; double r7982 = fma(r7980, r7979, r7981); double r7983 = r7950 / r7956; double r7984 = 4.0; double r7985 = r7950 * r7984; double r7986 = r7983 * r7985; double r7987 = fma(r7953, r7982, r7986); double r7988 = sqrt(r7987); double r7989 = r7978 / r7988; double r7990 = -1.6577109360186375e-242; bool r7991 = r7950 <= r7990; double r7992 = 1.7103591102507113e-246; bool r7993 = r7950 <= r7992; double r7994 = 4.6275459220879563e-169; bool r7995 = r7950 <= r7994; double r7996 = !r7995; double r7997 = 1.8364323594358619e+145; bool r7998 = r7950 <= r7997; bool r7999 = r7996 && r7998; bool r8000 = r7993 || r7999; double r8001 = r7983 / r7956; double r8002 = r8001 + r7983; double r8003 = r7964 / r7957; double r8004 = fma(r7960, r8003, r7955); double r8005 = fma(r7961, r8002, r8004); double r8006 = r7955 / r8005; double r8007 = r8000 ? r7989 : r8006; double r8008 = r7991 ? r7969 : r8007; double r8009 = r7971 ? r7989 : r8008; double r8010 = r7952 ? r7969 : r8009; return r8010; } 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 r8011, r8012, r8013, r8014, r8015, r8016, r8017, r8018, r8019, r8020, r8021, r8022, r8023, r8024, r8025, r8026, r8027, r8028; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r8011, "2", 10, MPFR_RNDN); mpfr_init(r8012); mpfr_init(r8013); mpfr_init(r8014); mpfr_init(r8015); mpfr_init_set_str(r8016, "1", 10, MPFR_RNDN); mpfr_init(r8017); mpfr_init(r8018); mpfr_init(r8019); mpfr_init(r8020); mpfr_init(r8021); mpfr_init(r8022); mpfr_init(r8023); mpfr_init(r8024); mpfr_init(r8025); mpfr_init(r8026); mpfr_init(r8027); mpfr_init(r8028); } double f_im(double x, double l, double t) { ; mpfr_sqrt(r8012, r8011, MPFR_RNDN); mpfr_set_d(r8013, t, MPFR_RNDN); mpfr_mul(r8014, r8012, r8013, MPFR_RNDN); mpfr_set_d(r8015, x, MPFR_RNDN); ; mpfr_add(r8017, r8015, r8016, MPFR_RNDN); mpfr_sub(r8018, r8015, r8016, MPFR_RNDN); mpfr_div(r8019, r8017, r8018, MPFR_RNDN); mpfr_set_d(r8020, l, MPFR_RNDN); mpfr_mul(r8021, r8020, r8020, MPFR_RNDN); mpfr_mul(r8022, r8013, r8013, MPFR_RNDN); mpfr_mul(r8023, r8011, r8022, MPFR_RNDN); mpfr_add(r8024, r8021, r8023, MPFR_RNDN); mpfr_mul(r8025, r8019, r8024, MPFR_RNDN); mpfr_sub(r8026, r8025, r8021, MPFR_RNDN); mpfr_sqrt(r8027, r8026, MPFR_RNDN); mpfr_div(r8028, r8014, r8027, MPFR_RNDN); return mpfr_get_d(r8028, MPFR_RNDN); } static mpfr_t r8029, r8030, r8031, r8032, r8033, r8034, r8035, r8036, r8037, r8038, r8039, r8040, r8041, r8042, r8043, r8044, r8045, r8046, r8047, r8048, r8049, r8050, r8051, r8052, r8053, r8054, r8055, r8056, r8057, r8058, r8059, r8060, r8061, r8062, r8063, r8064, r8065, r8066, r8067, r8068, r8069, r8070, r8071, r8072, r8073, r8074, r8075, r8076, r8077, r8078, r8079, r8080, r8081, r8082, r8083, r8084, r8085, r8086, r8087, r8088, r8089; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r8029); mpfr_init_set_str(r8030, "-6.446840727595846e+56", 10, MPFR_RNDN); mpfr_init(r8031); mpfr_init_set_str(r8032, "2", 10, MPFR_RNDN); mpfr_init(r8033); mpfr_init(r8034); mpfr_init(r8035); mpfr_init(r8036); mpfr_init(r8037); mpfr_init_set_str(r8038, "1", 10, MPFR_RNDN); mpfr_init(r8039); mpfr_init(r8040); mpfr_init(r8041); mpfr_init(r8042); mpfr_init(r8043); mpfr_init(r8044); mpfr_init(r8045); mpfr_init(r8046); mpfr_init(r8047); mpfr_init(r8048); mpfr_init_set_str(r8049, "-1.9073716957841407e-160", 10, MPFR_RNDN); mpfr_init(r8050); mpfr_init(r8051); mpfr_init(r8052); mpfr_init(r8053); mpfr_init(r8054); mpfr_init(r8055); mpfr_init(r8056); mpfr_init(r8057); mpfr_init(r8058); mpfr_init(r8059); mpfr_init(r8060); mpfr_init(r8061); mpfr_init(r8062); mpfr_init_set_str(r8063, "4", 10, MPFR_RNDN); mpfr_init(r8064); mpfr_init(r8065); mpfr_init(r8066); mpfr_init(r8067); mpfr_init(r8068); mpfr_init_set_str(r8069, "-1.6577109360186375e-242", 10, MPFR_RNDN); mpfr_init(r8070); mpfr_init_set_str(r8071, "1.7103591102507113e-246", 10, MPFR_RNDN); mpfr_init(r8072); mpfr_init_set_str(r8073, "4.6275459220879563e-169", 10, MPFR_RNDN); mpfr_init(r8074); mpfr_init(r8075); mpfr_init_set_str(r8076, "1.8364323594358619e+145", 10, MPFR_RNDN); mpfr_init(r8077); mpfr_init(r8078); mpfr_init(r8079); mpfr_init(r8080); mpfr_init(r8081); mpfr_init(r8082); mpfr_init(r8083); mpfr_init(r8084); mpfr_init(r8085); mpfr_init(r8086); mpfr_init(r8087); mpfr_init(r8088); mpfr_init(r8089); } double f_fm(double x, double l, double t) { mpfr_set_d(r8029, t, MPFR_RNDN); ; mpfr_set_si(r8031, mpfr_cmp(r8029, r8030) <= 0, MPFR_RNDN); ; mpfr_sqrt(r8033, r8032, MPFR_RNDN); mpfr_mul(r8034, r8033, r8029, MPFR_RNDN); mpfr_set_d(r8035, x, MPFR_RNDN); mpfr_mul(r8036, r8035, r8035, MPFR_RNDN); mpfr_div(r8037, r8029, r8036, MPFR_RNDN); ; mpfr_div(r8039, r8038, r8033, MPFR_RNDN); mpfr_div(r8040, r8032, r8033, MPFR_RNDN); mpfr_sub(r8041, r8039, r8040, MPFR_RNDN); mpfr_div(r8042, r8032, r8035, MPFR_RNDN); mpfr_neg(r8043, r8029, MPFR_RNDN); mpfr_div(r8044, r8043, r8033, MPFR_RNDN); mpfr_mul(r8045, r8043, r8033, MPFR_RNDN); mpfr_fma(r8046, r8042, r8044, r8045, MPFR_RNDN); mpfr_fma(r8047, r8037, r8041, r8046, MPFR_RNDN); mpfr_div(r8048, r8034, r8047, MPFR_RNDN); ; mpfr_set_si(r8050, mpfr_cmp(r8029, r8049) <= 0, MPFR_RNDN); mpfr_cbrt(r8051, r8033, MPFR_RNDN); mpfr_abs(r8052, r8051, MPFR_RNDN); mpfr_sqrt(r8053, r8051, MPFR_RNDN); mpfr_mul(r8054, r8029, r8053, MPFR_RNDN); mpfr_mul(r8055, r8052, r8054, MPFR_RNDN); mpfr_sqrt(r8056, r8033, MPFR_RNDN); mpfr_mul(r8057, r8055, r8056, MPFR_RNDN); mpfr_set_d(r8058, l, MPFR_RNDN); mpfr_div(r8059, r8058, r8035, MPFR_RNDN); mpfr_mul(r8060, r8029, r8029, MPFR_RNDN); mpfr_fma(r8061, r8059, r8058, r8060, MPFR_RNDN); mpfr_div(r8062, r8029, r8035, MPFR_RNDN); ; mpfr_mul(r8064, r8029, r8063, MPFR_RNDN); mpfr_mul(r8065, r8062, r8064, MPFR_RNDN); mpfr_fma(r8066, r8032, r8061, r8065, MPFR_RNDN); mpfr_sqrt(r8067, r8066, MPFR_RNDN); mpfr_div(r8068, r8057, r8067, MPFR_RNDN); ; mpfr_set_si(r8070, mpfr_cmp(r8029, r8069) <= 0, MPFR_RNDN); ; mpfr_set_si(r8072, mpfr_cmp(r8029, r8071) <= 0, MPFR_RNDN); ; mpfr_set_si(r8074, mpfr_cmp(r8029, r8073) <= 0, MPFR_RNDN); mpfr_set_si(r8075, !mpfr_get_si(r8074, MPFR_RNDN), MPFR_RNDN); ; mpfr_set_si(r8077, mpfr_cmp(r8029, r8076) <= 0, MPFR_RNDN); mpfr_set_si(r8078, mpfr_get_si(r8075, MPFR_RNDN) && mpfr_get_si(r8077, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r8079, mpfr_get_si(r8072, MPFR_RNDN) || mpfr_get_si(r8078, MPFR_RNDN), MPFR_RNDN); mpfr_div(r8080, r8062, r8035, MPFR_RNDN); mpfr_add(r8081, r8080, r8062, MPFR_RNDN); mpfr_div(r8082, r8043, r8036, MPFR_RNDN); mpfr_fma(r8083, r8039, r8082, r8034, MPFR_RNDN); mpfr_fma(r8084, r8040, r8081, r8083, MPFR_RNDN); mpfr_div(r8085, r8034, r8084, MPFR_RNDN); if (mpfr_get_si(r8079, MPFR_RNDN)) { mpfr_set(r8086, r8068, MPFR_RNDN); } else { mpfr_set(r8086, r8085, MPFR_RNDN); }; if (mpfr_get_si(r8070, MPFR_RNDN)) { mpfr_set(r8087, r8048, MPFR_RNDN); } else { mpfr_set(r8087, r8086, MPFR_RNDN); }; if (mpfr_get_si(r8050, MPFR_RNDN)) { mpfr_set(r8088, r8068, MPFR_RNDN); } else { mpfr_set(r8088, r8087, MPFR_RNDN); }; if (mpfr_get_si(r8031, MPFR_RNDN)) { mpfr_set(r8089, r8048, MPFR_RNDN); } else { mpfr_set(r8089, r8088, MPFR_RNDN); }; return mpfr_get_d(r8089, MPFR_RNDN); } static mpfr_t r8090, r8091, r8092, r8093, r8094, r8095, r8096, r8097, r8098, r8099, r8100, r8101, r8102, r8103, r8104, r8105, r8106, r8107, r8108, r8109, r8110, r8111, r8112, r8113, r8114, r8115, r8116, r8117, r8118, r8119, r8120, r8121, r8122, r8123, r8124, r8125, r8126, r8127, r8128, r8129, r8130, r8131, r8132, r8133, r8134, r8135, r8136, r8137, r8138, r8139, r8140, r8141, r8142, r8143, r8144, r8145, r8146, r8147, r8148, r8149, r8150; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r8090); mpfr_init_set_str(r8091, "-6.446840727595846e+56", 10, MPFR_RNDN); mpfr_init(r8092); mpfr_init_set_str(r8093, "2", 10, MPFR_RNDN); mpfr_init(r8094); mpfr_init(r8095); mpfr_init(r8096); mpfr_init(r8097); mpfr_init(r8098); mpfr_init_set_str(r8099, "1", 10, MPFR_RNDN); mpfr_init(r8100); mpfr_init(r8101); mpfr_init(r8102); mpfr_init(r8103); mpfr_init(r8104); mpfr_init(r8105); mpfr_init(r8106); mpfr_init(r8107); mpfr_init(r8108); mpfr_init(r8109); mpfr_init_set_str(r8110, "-1.9073716957841407e-160", 10, MPFR_RNDN); mpfr_init(r8111); mpfr_init(r8112); mpfr_init(r8113); mpfr_init(r8114); mpfr_init(r8115); mpfr_init(r8116); mpfr_init(r8117); mpfr_init(r8118); mpfr_init(r8119); mpfr_init(r8120); mpfr_init(r8121); mpfr_init(r8122); mpfr_init(r8123); mpfr_init_set_str(r8124, "4", 10, MPFR_RNDN); mpfr_init(r8125); mpfr_init(r8126); mpfr_init(r8127); mpfr_init(r8128); mpfr_init(r8129); mpfr_init_set_str(r8130, "-1.6577109360186375e-242", 10, MPFR_RNDN); mpfr_init(r8131); mpfr_init_set_str(r8132, "1.7103591102507113e-246", 10, MPFR_RNDN); mpfr_init(r8133); mpfr_init_set_str(r8134, "4.6275459220879563e-169", 10, MPFR_RNDN); mpfr_init(r8135); mpfr_init(r8136); mpfr_init_set_str(r8137, "1.8364323594358619e+145", 10, MPFR_RNDN); mpfr_init(r8138); mpfr_init(r8139); mpfr_init(r8140); mpfr_init(r8141); mpfr_init(r8142); mpfr_init(r8143); mpfr_init(r8144); mpfr_init(r8145); mpfr_init(r8146); mpfr_init(r8147); mpfr_init(r8148); mpfr_init(r8149); mpfr_init(r8150); } double f_dm(double x, double l, double t) { mpfr_set_d(r8090, t, MPFR_RNDN); ; mpfr_set_si(r8092, mpfr_cmp(r8090, r8091) <= 0, MPFR_RNDN); ; mpfr_sqrt(r8094, r8093, MPFR_RNDN); mpfr_mul(r8095, r8094, r8090, MPFR_RNDN); mpfr_set_d(r8096, x, MPFR_RNDN); mpfr_mul(r8097, r8096, r8096, MPFR_RNDN); mpfr_div(r8098, r8090, r8097, MPFR_RNDN); ; mpfr_div(r8100, r8099, r8094, MPFR_RNDN); mpfr_div(r8101, r8093, r8094, MPFR_RNDN); mpfr_sub(r8102, r8100, r8101, MPFR_RNDN); mpfr_div(r8103, r8093, r8096, MPFR_RNDN); mpfr_neg(r8104, r8090, MPFR_RNDN); mpfr_div(r8105, r8104, r8094, MPFR_RNDN); mpfr_mul(r8106, r8104, r8094, MPFR_RNDN); mpfr_fma(r8107, r8103, r8105, r8106, MPFR_RNDN); mpfr_fma(r8108, r8098, r8102, r8107, MPFR_RNDN); mpfr_div(r8109, r8095, r8108, MPFR_RNDN); ; mpfr_set_si(r8111, mpfr_cmp(r8090, r8110) <= 0, MPFR_RNDN); mpfr_cbrt(r8112, r8094, MPFR_RNDN); mpfr_abs(r8113, r8112, MPFR_RNDN); mpfr_sqrt(r8114, r8112, MPFR_RNDN); mpfr_mul(r8115, r8090, r8114, MPFR_RNDN); mpfr_mul(r8116, r8113, r8115, MPFR_RNDN); mpfr_sqrt(r8117, r8094, MPFR_RNDN); mpfr_mul(r8118, r8116, r8117, MPFR_RNDN); mpfr_set_d(r8119, l, MPFR_RNDN); mpfr_div(r8120, r8119, r8096, MPFR_RNDN); mpfr_mul(r8121, r8090, r8090, MPFR_RNDN); mpfr_fma(r8122, r8120, r8119, r8121, MPFR_RNDN); mpfr_div(r8123, r8090, r8096, MPFR_RNDN); ; mpfr_mul(r8125, r8090, r8124, MPFR_RNDN); mpfr_mul(r8126, r8123, r8125, MPFR_RNDN); mpfr_fma(r8127, r8093, r8122, r8126, MPFR_RNDN); mpfr_sqrt(r8128, r8127, MPFR_RNDN); mpfr_div(r8129, r8118, r8128, MPFR_RNDN); ; mpfr_set_si(r8131, mpfr_cmp(r8090, r8130) <= 0, MPFR_RNDN); ; mpfr_set_si(r8133, mpfr_cmp(r8090, r8132) <= 0, MPFR_RNDN); ; mpfr_set_si(r8135, mpfr_cmp(r8090, r8134) <= 0, MPFR_RNDN); mpfr_set_si(r8136, !mpfr_get_si(r8135, MPFR_RNDN), MPFR_RNDN); ; mpfr_set_si(r8138, mpfr_cmp(r8090, r8137) <= 0, MPFR_RNDN); mpfr_set_si(r8139, mpfr_get_si(r8136, MPFR_RNDN) && mpfr_get_si(r8138, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r8140, mpfr_get_si(r8133, MPFR_RNDN) || mpfr_get_si(r8139, MPFR_RNDN), MPFR_RNDN); mpfr_div(r8141, r8123, r8096, MPFR_RNDN); mpfr_add(r8142, r8141, r8123, MPFR_RNDN); mpfr_div(r8143, r8104, r8097, MPFR_RNDN); mpfr_fma(r8144, r8100, r8143, r8095, MPFR_RNDN); mpfr_fma(r8145, r8101, r8142, r8144, MPFR_RNDN); mpfr_div(r8146, r8095, r8145, MPFR_RNDN); if (mpfr_get_si(r8140, MPFR_RNDN)) { mpfr_set(r8147, r8129, MPFR_RNDN); } else { mpfr_set(r8147, r8146, MPFR_RNDN); }; if (mpfr_get_si(r8131, MPFR_RNDN)) { mpfr_set(r8148, r8109, MPFR_RNDN); } else { mpfr_set(r8148, r8147, MPFR_RNDN); }; if (mpfr_get_si(r8111, MPFR_RNDN)) { mpfr_set(r8149, r8129, MPFR_RNDN); } else { mpfr_set(r8149, r8148, MPFR_RNDN); }; if (mpfr_get_si(r8092, MPFR_RNDN)) { mpfr_set(r8150, r8109, MPFR_RNDN); } else { mpfr_set(r8150, r8149, MPFR_RNDN); }; return mpfr_get_d(r8150, MPFR_RNDN); }