#include #include #include #include #include char *name = "Toniolo and Linder, Equation (7)"; double f_if(float x, float l, float t) { float r19959 = 2; float r19960 = sqrt(r19959); float r19961 = t; float r19962 = r19960 * r19961; float r19963 = x; float r19964 = 1; float r19965 = r19963 + r19964; float r19966 = r19963 - r19964; float r19967 = r19965 / r19966; float r19968 = l; float r19969 = r19968 * r19968; float r19970 = r19961 * r19961; float r19971 = r19959 * r19970; float r19972 = r19969 + r19971; float r19973 = r19967 * r19972; float r19974 = r19973 - r19969; float r19975 = sqrt(r19974); float r19976 = r19962 / r19975; return r19976; } double f_id(double x, double l, double t) { double r19977 = 2; double r19978 = sqrt(r19977); double r19979 = t; double r19980 = r19978 * r19979; double r19981 = x; double r19982 = 1; double r19983 = r19981 + r19982; double r19984 = r19981 - r19982; double r19985 = r19983 / r19984; double r19986 = l; double r19987 = r19986 * r19986; double r19988 = r19979 * r19979; double r19989 = r19977 * r19988; double r19990 = r19987 + r19989; double r19991 = r19985 * r19990; double r19992 = r19991 - r19987; double r19993 = sqrt(r19992); double r19994 = r19980 / r19993; return r19994; } double f_of(float x, float l, float t) { float r19995 = t; float r19996 = -4.655030284999508e+62; bool r19997 = r19995 <= r19996; float r19998 = 2; float r19999 = sqrt(r19998); float r20000 = r19995 * r19999; float r20001 = x; float r20002 = r19995 / r20001; float r20003 = r20002 / r20001; float r20004 = r20003 / r19999; float r20005 = r20004 - r20000; float r20006 = r19998 / r20001; float r20007 = r20006 / r19999; float r20008 = r20002 + r19995; float r20009 = r20007 * r20008; float r20010 = r20005 - r20009; float r20011 = r20000 / r20010; float r20012 = -2.6020978138459956e-176; bool r20013 = r19995 <= r20012; float r20014 = cbrt(r19999); float r20015 = r20014 * r20014; float r20016 = r20014 * r19995; float r20017 = r20015 * r20016; float r20018 = pow(r19995, r19998); float r20019 = r19998 * r20018; float r20020 = l; float r20021 = r20001 / r20020; float r20022 = r20020 / r20021; float r20023 = r19998 * r20022; float r20024 = 4; float r20025 = r20018 / r20001; float r20026 = r20024 * r20025; float r20027 = r20023 + r20026; float r20028 = r20019 + r20027; float r20029 = sqrt(r20028); float r20030 = r20017 / r20029; float r20031 = -3.816162453277027e-251; bool r20032 = r19995 <= r20031; float r20033 = 6.267275199579288e-307; bool r20034 = r19995 <= r20033; float r20035 = sqrt(r19999); float r20036 = r20035 * r19995; float r20037 = r20035 * r20036; float r20038 = r20037 / r20029; float r20039 = 2.2638065892910295e-177; bool r20040 = r19995 <= r20039; float r20041 = r20009 + r20000; float r20042 = r20041 - r20004; float r20043 = r20000 / r20042; float r20044 = 3.9224840747176794e+133; bool r20045 = r19995 <= r20044; float r20046 = r19999 * r19995; float r20047 = r20046 / r20029; float r20048 = r20045 ? r20047 : r20043; float r20049 = r20040 ? r20043 : r20048; float r20050 = r20034 ? r20038 : r20049; float r20051 = r20032 ? r20011 : r20050; float r20052 = r20013 ? r20030 : r20051; float r20053 = r19997 ? r20011 : r20052; return r20053; } double f_od(double x, double l, double t) { double r20054 = t; double r20055 = -4.655030284999508e+62; bool r20056 = r20054 <= r20055; double r20057 = 2; double r20058 = sqrt(r20057); double r20059 = r20054 * r20058; double r20060 = x; double r20061 = r20054 / r20060; double r20062 = r20061 / r20060; double r20063 = r20062 / r20058; double r20064 = r20063 - r20059; double r20065 = r20057 / r20060; double r20066 = r20065 / r20058; double r20067 = r20061 + r20054; double r20068 = r20066 * r20067; double r20069 = r20064 - r20068; double r20070 = r20059 / r20069; double r20071 = -2.6020978138459956e-176; bool r20072 = r20054 <= r20071; double r20073 = cbrt(r20058); double r20074 = r20073 * r20073; double r20075 = r20073 * r20054; double r20076 = r20074 * r20075; double r20077 = pow(r20054, r20057); double r20078 = r20057 * r20077; double r20079 = l; double r20080 = r20060 / r20079; double r20081 = r20079 / r20080; double r20082 = r20057 * r20081; double r20083 = 4; double r20084 = r20077 / r20060; double r20085 = r20083 * r20084; double r20086 = r20082 + r20085; double r20087 = r20078 + r20086; double r20088 = sqrt(r20087); double r20089 = r20076 / r20088; double r20090 = -3.816162453277027e-251; bool r20091 = r20054 <= r20090; double r20092 = 6.267275199579288e-307; bool r20093 = r20054 <= r20092; double r20094 = sqrt(r20058); double r20095 = r20094 * r20054; double r20096 = r20094 * r20095; double r20097 = r20096 / r20088; double r20098 = 2.2638065892910295e-177; bool r20099 = r20054 <= r20098; double r20100 = r20068 + r20059; double r20101 = r20100 - r20063; double r20102 = r20059 / r20101; double r20103 = 3.9224840747176794e+133; bool r20104 = r20054 <= r20103; double r20105 = r20058 * r20054; double r20106 = r20105 / r20088; double r20107 = r20104 ? r20106 : r20102; double r20108 = r20099 ? r20102 : r20107; double r20109 = r20093 ? r20097 : r20108; double r20110 = r20091 ? r20070 : r20109; double r20111 = r20072 ? r20089 : r20110; double r20112 = r20056 ? r20070 : r20111; return r20112; } 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 r20113, r20114, r20115, r20116, r20117, r20118, r20119, r20120, r20121, r20122, r20123, r20124, r20125, r20126, r20127, r20128, r20129, r20130; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r20113, "2", 10, MPFR_RNDN); mpfr_init(r20114); mpfr_init(r20115); mpfr_init(r20116); mpfr_init(r20117); mpfr_init_set_str(r20118, "1", 10, MPFR_RNDN); mpfr_init(r20119); mpfr_init(r20120); mpfr_init(r20121); mpfr_init(r20122); mpfr_init(r20123); mpfr_init(r20124); mpfr_init(r20125); mpfr_init(r20126); mpfr_init(r20127); mpfr_init(r20128); mpfr_init(r20129); mpfr_init(r20130); } double f_im(double x, double l, double t) { ; mpfr_sqrt(r20114, r20113, MPFR_RNDN); mpfr_set_d(r20115, t, MPFR_RNDN); mpfr_mul(r20116, r20114, r20115, MPFR_RNDN); mpfr_set_d(r20117, x, MPFR_RNDN); ; mpfr_add(r20119, r20117, r20118, MPFR_RNDN); mpfr_sub(r20120, r20117, r20118, MPFR_RNDN); mpfr_div(r20121, r20119, r20120, MPFR_RNDN); mpfr_set_d(r20122, l, MPFR_RNDN); mpfr_mul(r20123, r20122, r20122, MPFR_RNDN); mpfr_mul(r20124, r20115, r20115, MPFR_RNDN); mpfr_mul(r20125, r20113, r20124, MPFR_RNDN); mpfr_add(r20126, r20123, r20125, MPFR_RNDN); mpfr_mul(r20127, r20121, r20126, MPFR_RNDN); mpfr_sub(r20128, r20127, r20123, MPFR_RNDN); mpfr_sqrt(r20129, r20128, MPFR_RNDN); mpfr_div(r20130, r20116, r20129, MPFR_RNDN); return mpfr_get_d(r20130, MPFR_RNDN); } static mpfr_t r20131, r20132, r20133, r20134, r20135, r20136, r20137, r20138, r20139, r20140, r20141, r20142, r20143, r20144, r20145, r20146, r20147, r20148, r20149, r20150, r20151, r20152, r20153, r20154, r20155, r20156, r20157, r20158, r20159, r20160, r20161, r20162, r20163, r20164, r20165, r20166, r20167, r20168, r20169, r20170, r20171, r20172, r20173, r20174, r20175, r20176, r20177, r20178, r20179, r20180, r20181, r20182, r20183, r20184, r20185, r20186, r20187, r20188, r20189; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r20131); mpfr_init_set_str(r20132, "-4.655030284999508e+62", 10, MPFR_RNDN); mpfr_init(r20133); mpfr_init_set_str(r20134, "2", 10, MPFR_RNDN); mpfr_init(r20135); mpfr_init(r20136); mpfr_init(r20137); mpfr_init(r20138); mpfr_init(r20139); mpfr_init(r20140); mpfr_init(r20141); mpfr_init(r20142); mpfr_init(r20143); mpfr_init(r20144); mpfr_init(r20145); mpfr_init(r20146); mpfr_init(r20147); mpfr_init_set_str(r20148, "-2.6020978138459956e-176", 10, MPFR_RNDN); mpfr_init(r20149); mpfr_init(r20150); mpfr_init(r20151); mpfr_init(r20152); mpfr_init(r20153); mpfr_init(r20154); mpfr_init(r20155); mpfr_init(r20156); mpfr_init(r20157); mpfr_init(r20158); mpfr_init(r20159); mpfr_init_set_str(r20160, "4", 10, MPFR_RNDN); mpfr_init(r20161); mpfr_init(r20162); mpfr_init(r20163); mpfr_init(r20164); mpfr_init(r20165); mpfr_init(r20166); mpfr_init_set_str(r20167, "-3.816162453277027e-251", 10, MPFR_RNDN); mpfr_init(r20168); mpfr_init_set_str(r20169, "6.267275199579288e-307", 10, MPFR_RNDN); mpfr_init(r20170); mpfr_init(r20171); mpfr_init(r20172); mpfr_init(r20173); mpfr_init(r20174); mpfr_init_set_str(r20175, "2.2638065892910295e-177", 10, MPFR_RNDN); mpfr_init(r20176); mpfr_init(r20177); mpfr_init(r20178); mpfr_init(r20179); mpfr_init_set_str(r20180, "3.9224840747176794e+133", 10, MPFR_RNDN); mpfr_init(r20181); mpfr_init(r20182); mpfr_init(r20183); mpfr_init(r20184); mpfr_init(r20185); mpfr_init(r20186); mpfr_init(r20187); mpfr_init(r20188); mpfr_init(r20189); } double f_fm(double x, double l, double t) { mpfr_set_d(r20131, t, MPFR_RNDN); ; mpfr_set_si(r20133, mpfr_cmp(r20131, r20132) <= 0, MPFR_RNDN); ; mpfr_sqrt(r20135, r20134, MPFR_RNDN); mpfr_mul(r20136, r20131, r20135, MPFR_RNDN); mpfr_set_d(r20137, x, MPFR_RNDN); mpfr_div(r20138, r20131, r20137, MPFR_RNDN); mpfr_div(r20139, r20138, r20137, MPFR_RNDN); mpfr_div(r20140, r20139, r20135, MPFR_RNDN); mpfr_sub(r20141, r20140, r20136, MPFR_RNDN); mpfr_div(r20142, r20134, r20137, MPFR_RNDN); mpfr_div(r20143, r20142, r20135, MPFR_RNDN); mpfr_add(r20144, r20138, r20131, MPFR_RNDN); mpfr_mul(r20145, r20143, r20144, MPFR_RNDN); mpfr_sub(r20146, r20141, r20145, MPFR_RNDN); mpfr_div(r20147, r20136, r20146, MPFR_RNDN); ; mpfr_set_si(r20149, mpfr_cmp(r20131, r20148) <= 0, MPFR_RNDN); mpfr_cbrt(r20150, r20135, MPFR_RNDN); mpfr_mul(r20151, r20150, r20150, MPFR_RNDN); mpfr_mul(r20152, r20150, r20131, MPFR_RNDN); mpfr_mul(r20153, r20151, r20152, MPFR_RNDN); mpfr_pow(r20154, r20131, r20134, MPFR_RNDN); mpfr_mul(r20155, r20134, r20154, MPFR_RNDN); mpfr_set_d(r20156, l, MPFR_RNDN); mpfr_div(r20157, r20137, r20156, MPFR_RNDN); mpfr_div(r20158, r20156, r20157, MPFR_RNDN); mpfr_mul(r20159, r20134, r20158, MPFR_RNDN); ; mpfr_div(r20161, r20154, r20137, MPFR_RNDN); mpfr_mul(r20162, r20160, r20161, MPFR_RNDN); mpfr_add(r20163, r20159, r20162, MPFR_RNDN); mpfr_add(r20164, r20155, r20163, MPFR_RNDN); mpfr_sqrt(r20165, r20164, MPFR_RNDN); mpfr_div(r20166, r20153, r20165, MPFR_RNDN); ; mpfr_set_si(r20168, mpfr_cmp(r20131, r20167) <= 0, MPFR_RNDN); ; mpfr_set_si(r20170, mpfr_cmp(r20131, r20169) <= 0, MPFR_RNDN); mpfr_sqrt(r20171, r20135, MPFR_RNDN); mpfr_mul(r20172, r20171, r20131, MPFR_RNDN); mpfr_mul(r20173, r20171, r20172, MPFR_RNDN); mpfr_div(r20174, r20173, r20165, MPFR_RNDN); ; mpfr_set_si(r20176, mpfr_cmp(r20131, r20175) <= 0, MPFR_RNDN); mpfr_add(r20177, r20145, r20136, MPFR_RNDN); mpfr_sub(r20178, r20177, r20140, MPFR_RNDN); mpfr_div(r20179, r20136, r20178, MPFR_RNDN); ; mpfr_set_si(r20181, mpfr_cmp(r20131, r20180) <= 0, MPFR_RNDN); mpfr_mul(r20182, r20135, r20131, MPFR_RNDN); mpfr_div(r20183, r20182, r20165, MPFR_RNDN); if (mpfr_get_si(r20181, MPFR_RNDN)) { mpfr_set(r20184, r20183, MPFR_RNDN); } else { mpfr_set(r20184, r20179, MPFR_RNDN); }; if (mpfr_get_si(r20176, MPFR_RNDN)) { mpfr_set(r20185, r20179, MPFR_RNDN); } else { mpfr_set(r20185, r20184, MPFR_RNDN); }; if (mpfr_get_si(r20170, MPFR_RNDN)) { mpfr_set(r20186, r20174, MPFR_RNDN); } else { mpfr_set(r20186, r20185, MPFR_RNDN); }; if (mpfr_get_si(r20168, MPFR_RNDN)) { mpfr_set(r20187, r20147, MPFR_RNDN); } else { mpfr_set(r20187, r20186, MPFR_RNDN); }; if (mpfr_get_si(r20149, MPFR_RNDN)) { mpfr_set(r20188, r20166, MPFR_RNDN); } else { mpfr_set(r20188, r20187, MPFR_RNDN); }; if (mpfr_get_si(r20133, MPFR_RNDN)) { mpfr_set(r20189, r20147, MPFR_RNDN); } else { mpfr_set(r20189, r20188, MPFR_RNDN); }; return mpfr_get_d(r20189, MPFR_RNDN); } static mpfr_t r20190, r20191, r20192, r20193, r20194, r20195, r20196, r20197, r20198, r20199, r20200, r20201, r20202, r20203, r20204, r20205, r20206, r20207, r20208, r20209, r20210, r20211, r20212, r20213, r20214, r20215, r20216, r20217, r20218, r20219, r20220, r20221, r20222, r20223, r20224, r20225, r20226, r20227, r20228, r20229, r20230, r20231, r20232, r20233, r20234, r20235, r20236, r20237, r20238, r20239, r20240, r20241, r20242, r20243, r20244, r20245, r20246, r20247, r20248; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r20190); mpfr_init_set_str(r20191, "-4.655030284999508e+62", 10, MPFR_RNDN); mpfr_init(r20192); mpfr_init_set_str(r20193, "2", 10, MPFR_RNDN); mpfr_init(r20194); mpfr_init(r20195); mpfr_init(r20196); mpfr_init(r20197); mpfr_init(r20198); mpfr_init(r20199); mpfr_init(r20200); mpfr_init(r20201); mpfr_init(r20202); mpfr_init(r20203); mpfr_init(r20204); mpfr_init(r20205); mpfr_init(r20206); mpfr_init_set_str(r20207, "-2.6020978138459956e-176", 10, MPFR_RNDN); mpfr_init(r20208); mpfr_init(r20209); mpfr_init(r20210); mpfr_init(r20211); mpfr_init(r20212); mpfr_init(r20213); mpfr_init(r20214); mpfr_init(r20215); mpfr_init(r20216); mpfr_init(r20217); mpfr_init(r20218); mpfr_init_set_str(r20219, "4", 10, MPFR_RNDN); mpfr_init(r20220); mpfr_init(r20221); mpfr_init(r20222); mpfr_init(r20223); mpfr_init(r20224); mpfr_init(r20225); mpfr_init_set_str(r20226, "-3.816162453277027e-251", 10, MPFR_RNDN); mpfr_init(r20227); mpfr_init_set_str(r20228, "6.267275199579288e-307", 10, MPFR_RNDN); mpfr_init(r20229); mpfr_init(r20230); mpfr_init(r20231); mpfr_init(r20232); mpfr_init(r20233); mpfr_init_set_str(r20234, "2.2638065892910295e-177", 10, MPFR_RNDN); mpfr_init(r20235); mpfr_init(r20236); mpfr_init(r20237); mpfr_init(r20238); mpfr_init_set_str(r20239, "3.9224840747176794e+133", 10, MPFR_RNDN); mpfr_init(r20240); mpfr_init(r20241); mpfr_init(r20242); mpfr_init(r20243); mpfr_init(r20244); mpfr_init(r20245); mpfr_init(r20246); mpfr_init(r20247); mpfr_init(r20248); } double f_dm(double x, double l, double t) { mpfr_set_d(r20190, t, MPFR_RNDN); ; mpfr_set_si(r20192, mpfr_cmp(r20190, r20191) <= 0, MPFR_RNDN); ; mpfr_sqrt(r20194, r20193, MPFR_RNDN); mpfr_mul(r20195, r20190, r20194, MPFR_RNDN); mpfr_set_d(r20196, x, MPFR_RNDN); mpfr_div(r20197, r20190, r20196, MPFR_RNDN); mpfr_div(r20198, r20197, r20196, MPFR_RNDN); mpfr_div(r20199, r20198, r20194, MPFR_RNDN); mpfr_sub(r20200, r20199, r20195, MPFR_RNDN); mpfr_div(r20201, r20193, r20196, MPFR_RNDN); mpfr_div(r20202, r20201, r20194, MPFR_RNDN); mpfr_add(r20203, r20197, r20190, MPFR_RNDN); mpfr_mul(r20204, r20202, r20203, MPFR_RNDN); mpfr_sub(r20205, r20200, r20204, MPFR_RNDN); mpfr_div(r20206, r20195, r20205, MPFR_RNDN); ; mpfr_set_si(r20208, mpfr_cmp(r20190, r20207) <= 0, MPFR_RNDN); mpfr_cbrt(r20209, r20194, MPFR_RNDN); mpfr_mul(r20210, r20209, r20209, MPFR_RNDN); mpfr_mul(r20211, r20209, r20190, MPFR_RNDN); mpfr_mul(r20212, r20210, r20211, MPFR_RNDN); mpfr_pow(r20213, r20190, r20193, MPFR_RNDN); mpfr_mul(r20214, r20193, r20213, MPFR_RNDN); mpfr_set_d(r20215, l, MPFR_RNDN); mpfr_div(r20216, r20196, r20215, MPFR_RNDN); mpfr_div(r20217, r20215, r20216, MPFR_RNDN); mpfr_mul(r20218, r20193, r20217, MPFR_RNDN); ; mpfr_div(r20220, r20213, r20196, MPFR_RNDN); mpfr_mul(r20221, r20219, r20220, MPFR_RNDN); mpfr_add(r20222, r20218, r20221, MPFR_RNDN); mpfr_add(r20223, r20214, r20222, MPFR_RNDN); mpfr_sqrt(r20224, r20223, MPFR_RNDN); mpfr_div(r20225, r20212, r20224, MPFR_RNDN); ; mpfr_set_si(r20227, mpfr_cmp(r20190, r20226) <= 0, MPFR_RNDN); ; mpfr_set_si(r20229, mpfr_cmp(r20190, r20228) <= 0, MPFR_RNDN); mpfr_sqrt(r20230, r20194, MPFR_RNDN); mpfr_mul(r20231, r20230, r20190, MPFR_RNDN); mpfr_mul(r20232, r20230, r20231, MPFR_RNDN); mpfr_div(r20233, r20232, r20224, MPFR_RNDN); ; mpfr_set_si(r20235, mpfr_cmp(r20190, r20234) <= 0, MPFR_RNDN); mpfr_add(r20236, r20204, r20195, MPFR_RNDN); mpfr_sub(r20237, r20236, r20199, MPFR_RNDN); mpfr_div(r20238, r20195, r20237, MPFR_RNDN); ; mpfr_set_si(r20240, mpfr_cmp(r20190, r20239) <= 0, MPFR_RNDN); mpfr_mul(r20241, r20194, r20190, MPFR_RNDN); mpfr_div(r20242, r20241, r20224, MPFR_RNDN); if (mpfr_get_si(r20240, MPFR_RNDN)) { mpfr_set(r20243, r20242, MPFR_RNDN); } else { mpfr_set(r20243, r20238, MPFR_RNDN); }; if (mpfr_get_si(r20235, MPFR_RNDN)) { mpfr_set(r20244, r20238, MPFR_RNDN); } else { mpfr_set(r20244, r20243, MPFR_RNDN); }; if (mpfr_get_si(r20229, MPFR_RNDN)) { mpfr_set(r20245, r20233, MPFR_RNDN); } else { mpfr_set(r20245, r20244, MPFR_RNDN); }; if (mpfr_get_si(r20227, MPFR_RNDN)) { mpfr_set(r20246, r20206, MPFR_RNDN); } else { mpfr_set(r20246, r20245, MPFR_RNDN); }; if (mpfr_get_si(r20208, MPFR_RNDN)) { mpfr_set(r20247, r20225, MPFR_RNDN); } else { mpfr_set(r20247, r20246, MPFR_RNDN); }; if (mpfr_get_si(r20192, MPFR_RNDN)) { mpfr_set(r20248, r20206, MPFR_RNDN); } else { mpfr_set(r20248, r20247, MPFR_RNDN); }; return mpfr_get_d(r20248, MPFR_RNDN); }