#include #include #include #include #include char *name = "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2"; double f_if(float x, float y, float z, float t, float a, float b) { float r19980 = x; float r19981 = y; float r19982 = z; float r19983 = log(r19982); float r19984 = r19981 * r19983; float r19985 = t; float r19986 = 1.0f; float r19987 = r19985 - r19986; float r19988 = a; float r19989 = log(r19988); float r19990 = r19987 * r19989; float r19991 = r19984 + r19990; float r19992 = b; float r19993 = r19991 - r19992; float r19994 = exp(r19993); float r19995 = r19980 * r19994; float r19996 = r19995 / r19981; return r19996; } double f_id(double x, double y, double z, double t, double a, double b) { double r19997 = x; double r19998 = y; double r19999 = z; double r20000 = log(r19999); double r20001 = r19998 * r20000; double r20002 = t; double r20003 = 1.0; double r20004 = r20002 - r20003; double r20005 = a; double r20006 = log(r20005); double r20007 = r20004 * r20006; double r20008 = r20001 + r20007; double r20009 = b; double r20010 = r20008 - r20009; double r20011 = exp(r20010); double r20012 = r19997 * r20011; double r20013 = r20012 / r19998; return r20013; } double f_of(float x, float y, float z, float t, float a, float b) { float r20014 = y; float r20015 = z; float r20016 = log(r20015); float r20017 = r20014 * r20016; float r20018 = -1.1154917005965662e+243f; bool r20019 = r20017 <= r20018; float r20020 = x; float r20021 = r20020 / r20014; float r20022 = log(r20021); float r20023 = a; float r20024 = log(r20023); float r20025 = t; float r20026 = 1.0f; float r20027 = r20025 - r20026; float r20028 = r20024 * r20027; float r20029 = b; float r20030 = r20016 * r20014; float r20031 = r20029 - r20030; float r20032 = r20028 - r20031; float r20033 = r20022 + r20032; float r20034 = exp(r20033); float r20035 = -6.82937553380511e+147f; bool r20036 = r20017 <= r20035; float r20037 = pow(r20015, r20014); float r20038 = 1.0f; float r20039 = pow(r20023, r20026); float r20040 = r20038 / r20039; float r20041 = pow(r20040, r20026); float r20042 = r20037 * r20041; float r20043 = r20024 * r20024; float r20044 = r20025 * r20025; float r20045 = 0.5f; float r20046 = r20044 * r20045; float r20047 = r20043 * r20046; float r20048 = r20025 * r20024; float r20049 = r20047 + r20048; float r20050 = r20042 * r20049; float r20051 = -1.0f; float r20052 = pow(r20023, r20051); float r20053 = r20052 * r20037; float r20054 = r20050 + r20053; float r20055 = exp(r20029); float r20056 = r20055 / r20021; float r20057 = r20054 / r20056; float r20058 = -2.850576188986565e-16f; bool r20059 = r20017 <= r20058; float r20060 = pow(r20023, r20027); float r20061 = exp(r20031); float r20062 = r20060 / r20061; float r20063 = r20021 * r20062; float r20064 = 0.0922372319414478f; bool r20065 = r20017 <= r20064; float r20066 = r20037 / r20014; float r20067 = r20028 - r20029; float r20068 = exp(r20067); float r20069 = r20066 * r20068; float r20070 = r20020 * r20069; float r20071 = r20065 ? r20070 : r20034; float r20072 = r20059 ? r20063 : r20071; float r20073 = r20036 ? r20057 : r20072; float r20074 = r20019 ? r20034 : r20073; return r20074; } double f_od(double x, double y, double z, double t, double a, double b) { double r20075 = y; double r20076 = z; double r20077 = log(r20076); double r20078 = r20075 * r20077; double r20079 = -1.1154917005965662e+243; bool r20080 = r20078 <= r20079; double r20081 = x; double r20082 = r20081 / r20075; double r20083 = log(r20082); double r20084 = a; double r20085 = log(r20084); double r20086 = t; double r20087 = 1.0; double r20088 = r20086 - r20087; double r20089 = r20085 * r20088; double r20090 = b; double r20091 = r20077 * r20075; double r20092 = r20090 - r20091; double r20093 = r20089 - r20092; double r20094 = r20083 + r20093; double r20095 = exp(r20094); double r20096 = -6.82937553380511e+147; bool r20097 = r20078 <= r20096; double r20098 = pow(r20076, r20075); double r20099 = 1.0; double r20100 = pow(r20084, r20087); double r20101 = r20099 / r20100; double r20102 = pow(r20101, r20087); double r20103 = r20098 * r20102; double r20104 = r20085 * r20085; double r20105 = r20086 * r20086; double r20106 = 0.5; double r20107 = r20105 * r20106; double r20108 = r20104 * r20107; double r20109 = r20086 * r20085; double r20110 = r20108 + r20109; double r20111 = r20103 * r20110; double r20112 = -1.0; double r20113 = pow(r20084, r20112); double r20114 = r20113 * r20098; double r20115 = r20111 + r20114; double r20116 = exp(r20090); double r20117 = r20116 / r20082; double r20118 = r20115 / r20117; double r20119 = -2.850576188986565e-16; bool r20120 = r20078 <= r20119; double r20121 = pow(r20084, r20088); double r20122 = exp(r20092); double r20123 = r20121 / r20122; double r20124 = r20082 * r20123; double r20125 = 0.0922372319414478; bool r20126 = r20078 <= r20125; double r20127 = r20098 / r20075; double r20128 = r20089 - r20090; double r20129 = exp(r20128); double r20130 = r20127 * r20129; double r20131 = r20081 * r20130; double r20132 = r20126 ? r20131 : r20095; double r20133 = r20120 ? r20124 : r20132; double r20134 = r20097 ? r20118 : r20133; double r20135 = r20080 ? r20095 : r20134; return r20135; } 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 r20136, r20137, r20138, r20139, r20140, r20141, r20142, r20143, r20144, r20145, r20146, r20147, r20148, r20149, r20150, r20151, r20152; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r20136); mpfr_init(r20137); mpfr_init(r20138); mpfr_init(r20139); mpfr_init(r20140); mpfr_init(r20141); mpfr_init_set_str(r20142, "1.0", 10, MPFR_RNDN); mpfr_init(r20143); mpfr_init(r20144); mpfr_init(r20145); mpfr_init(r20146); mpfr_init(r20147); mpfr_init(r20148); mpfr_init(r20149); mpfr_init(r20150); mpfr_init(r20151); mpfr_init(r20152); } double f_im(double x, double y, double z, double t, double a, double b) { mpfr_set_d(r20136, x, MPFR_RNDN); mpfr_set_d(r20137, y, MPFR_RNDN); mpfr_set_d(r20138, z, MPFR_RNDN); mpfr_log(r20139, r20138, MPFR_RNDN); mpfr_mul(r20140, r20137, r20139, MPFR_RNDN); mpfr_set_d(r20141, t, MPFR_RNDN); ; mpfr_sub(r20143, r20141, r20142, MPFR_RNDN); mpfr_set_d(r20144, a, MPFR_RNDN); mpfr_log(r20145, r20144, MPFR_RNDN); mpfr_mul(r20146, r20143, r20145, MPFR_RNDN); mpfr_add(r20147, r20140, r20146, MPFR_RNDN); mpfr_set_d(r20148, b, MPFR_RNDN); mpfr_sub(r20149, r20147, r20148, MPFR_RNDN); mpfr_exp(r20150, r20149, MPFR_RNDN); mpfr_mul(r20151, r20136, r20150, MPFR_RNDN); mpfr_div(r20152, r20151, r20137, MPFR_RNDN); return mpfr_get_d(r20152, MPFR_RNDN); } static mpfr_t 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, r20190, r20191, r20192, r20193, r20194, r20195, r20196, r20197, r20198, r20199, r20200, r20201, r20202, r20203, r20204, r20205, r20206, r20207, r20208, r20209, r20210, r20211, r20212, r20213; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r20153); mpfr_init(r20154); mpfr_init(r20155); mpfr_init(r20156); mpfr_init_set_str(r20157, "-1.1154917005965662e+243", 10, MPFR_RNDN); mpfr_init(r20158); mpfr_init(r20159); mpfr_init(r20160); mpfr_init(r20161); mpfr_init(r20162); mpfr_init(r20163); mpfr_init(r20164); mpfr_init_set_str(r20165, "1.0", 10, MPFR_RNDN); mpfr_init(r20166); mpfr_init(r20167); mpfr_init(r20168); mpfr_init(r20169); mpfr_init(r20170); mpfr_init(r20171); mpfr_init(r20172); mpfr_init(r20173); mpfr_init_set_str(r20174, "-6.82937553380511e+147", 10, MPFR_RNDN); mpfr_init(r20175); mpfr_init(r20176); mpfr_init_set_str(r20177, "1", 10, MPFR_RNDN); mpfr_init(r20178); mpfr_init(r20179); mpfr_init(r20180); mpfr_init(r20181); mpfr_init(r20182); mpfr_init(r20183); mpfr_init_set_str(r20184, "1/2", 10, MPFR_RNDN); mpfr_init(r20185); mpfr_init(r20186); mpfr_init(r20187); mpfr_init(r20188); mpfr_init(r20189); mpfr_init_set_str(r20190, "-1.0", 10, MPFR_RNDN); mpfr_init(r20191); mpfr_init(r20192); mpfr_init(r20193); mpfr_init(r20194); mpfr_init(r20195); mpfr_init(r20196); mpfr_init_set_str(r20197, "-2.850576188986565e-16", 10, MPFR_RNDN); mpfr_init(r20198); mpfr_init(r20199); mpfr_init(r20200); mpfr_init(r20201); mpfr_init(r20202); mpfr_init_set_str(r20203, "0.0922372319414478", 10, MPFR_RNDN); mpfr_init(r20204); mpfr_init(r20205); mpfr_init(r20206); mpfr_init(r20207); mpfr_init(r20208); mpfr_init(r20209); mpfr_init(r20210); mpfr_init(r20211); mpfr_init(r20212); mpfr_init(r20213); } double f_fm(double x, double y, double z, double t, double a, double b) { mpfr_set_d(r20153, y, MPFR_RNDN); mpfr_set_d(r20154, z, MPFR_RNDN); mpfr_log(r20155, r20154, MPFR_RNDN); mpfr_mul(r20156, r20153, r20155, MPFR_RNDN); ; mpfr_set_si(r20158, mpfr_cmp(r20156, r20157) <= 0, MPFR_RNDN); mpfr_set_d(r20159, x, MPFR_RNDN); mpfr_div(r20160, r20159, r20153, MPFR_RNDN); mpfr_log(r20161, r20160, MPFR_RNDN); mpfr_set_d(r20162, a, MPFR_RNDN); mpfr_log(r20163, r20162, MPFR_RNDN); mpfr_set_d(r20164, t, MPFR_RNDN); ; mpfr_sub(r20166, r20164, r20165, MPFR_RNDN); mpfr_mul(r20167, r20163, r20166, MPFR_RNDN); mpfr_set_d(r20168, b, MPFR_RNDN); mpfr_mul(r20169, r20155, r20153, MPFR_RNDN); mpfr_sub(r20170, r20168, r20169, MPFR_RNDN); mpfr_sub(r20171, r20167, r20170, MPFR_RNDN); mpfr_add(r20172, r20161, r20171, MPFR_RNDN); mpfr_exp(r20173, r20172, MPFR_RNDN); ; mpfr_set_si(r20175, mpfr_cmp(r20156, r20174) <= 0, MPFR_RNDN); mpfr_pow(r20176, r20154, r20153, MPFR_RNDN); ; mpfr_pow(r20178, r20162, r20165, MPFR_RNDN); mpfr_div(r20179, r20177, r20178, MPFR_RNDN); mpfr_pow(r20180, r20179, r20165, MPFR_RNDN); mpfr_mul(r20181, r20176, r20180, MPFR_RNDN); mpfr_mul(r20182, r20163, r20163, MPFR_RNDN); mpfr_mul(r20183, r20164, r20164, MPFR_RNDN); ; mpfr_mul(r20185, r20183, r20184, MPFR_RNDN); mpfr_mul(r20186, r20182, r20185, MPFR_RNDN); mpfr_mul(r20187, r20164, r20163, MPFR_RNDN); mpfr_add(r20188, r20186, r20187, MPFR_RNDN); mpfr_mul(r20189, r20181, r20188, MPFR_RNDN); ; mpfr_pow(r20191, r20162, r20190, MPFR_RNDN); mpfr_mul(r20192, r20191, r20176, MPFR_RNDN); mpfr_add(r20193, r20189, r20192, MPFR_RNDN); mpfr_exp(r20194, r20168, MPFR_RNDN); mpfr_div(r20195, r20194, r20160, MPFR_RNDN); mpfr_div(r20196, r20193, r20195, MPFR_RNDN); ; mpfr_set_si(r20198, mpfr_cmp(r20156, r20197) <= 0, MPFR_RNDN); mpfr_pow(r20199, r20162, r20166, MPFR_RNDN); mpfr_exp(r20200, r20170, MPFR_RNDN); mpfr_div(r20201, r20199, r20200, MPFR_RNDN); mpfr_mul(r20202, r20160, r20201, MPFR_RNDN); ; mpfr_set_si(r20204, mpfr_cmp(r20156, r20203) <= 0, MPFR_RNDN); mpfr_div(r20205, r20176, r20153, MPFR_RNDN); mpfr_sub(r20206, r20167, r20168, MPFR_RNDN); mpfr_exp(r20207, r20206, MPFR_RNDN); mpfr_mul(r20208, r20205, r20207, MPFR_RNDN); mpfr_mul(r20209, r20159, r20208, MPFR_RNDN); if (mpfr_get_si(r20204, MPFR_RNDN)) { mpfr_set(r20210, r20209, MPFR_RNDN); } else { mpfr_set(r20210, r20173, MPFR_RNDN); }; if (mpfr_get_si(r20198, MPFR_RNDN)) { mpfr_set(r20211, r20202, MPFR_RNDN); } else { mpfr_set(r20211, r20210, MPFR_RNDN); }; if (mpfr_get_si(r20175, MPFR_RNDN)) { mpfr_set(r20212, r20196, MPFR_RNDN); } else { mpfr_set(r20212, r20211, MPFR_RNDN); }; if (mpfr_get_si(r20158, MPFR_RNDN)) { mpfr_set(r20213, r20173, MPFR_RNDN); } else { mpfr_set(r20213, r20212, MPFR_RNDN); }; return mpfr_get_d(r20213, MPFR_RNDN); } static mpfr_t 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, r20249, r20250, r20251, r20252, r20253, r20254, r20255, r20256, r20257, r20258, r20259, r20260, r20261, r20262, r20263, r20264, r20265, r20266, r20267, r20268, r20269, r20270, r20271, r20272, r20273, r20274; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r20214); mpfr_init(r20215); mpfr_init(r20216); mpfr_init(r20217); mpfr_init_set_str(r20218, "-1.1154917005965662e+243", 10, MPFR_RNDN); mpfr_init(r20219); mpfr_init(r20220); mpfr_init(r20221); mpfr_init(r20222); mpfr_init(r20223); mpfr_init(r20224); mpfr_init(r20225); mpfr_init_set_str(r20226, "1.0", 10, MPFR_RNDN); mpfr_init(r20227); mpfr_init(r20228); mpfr_init(r20229); mpfr_init(r20230); mpfr_init(r20231); mpfr_init(r20232); mpfr_init(r20233); mpfr_init(r20234); mpfr_init_set_str(r20235, "-6.82937553380511e+147", 10, MPFR_RNDN); mpfr_init(r20236); mpfr_init(r20237); mpfr_init_set_str(r20238, "1", 10, MPFR_RNDN); mpfr_init(r20239); mpfr_init(r20240); mpfr_init(r20241); mpfr_init(r20242); mpfr_init(r20243); mpfr_init(r20244); mpfr_init_set_str(r20245, "1/2", 10, MPFR_RNDN); mpfr_init(r20246); mpfr_init(r20247); mpfr_init(r20248); mpfr_init(r20249); mpfr_init(r20250); mpfr_init_set_str(r20251, "-1.0", 10, MPFR_RNDN); mpfr_init(r20252); mpfr_init(r20253); mpfr_init(r20254); mpfr_init(r20255); mpfr_init(r20256); mpfr_init(r20257); mpfr_init_set_str(r20258, "-2.850576188986565e-16", 10, MPFR_RNDN); mpfr_init(r20259); mpfr_init(r20260); mpfr_init(r20261); mpfr_init(r20262); mpfr_init(r20263); mpfr_init_set_str(r20264, "0.0922372319414478", 10, MPFR_RNDN); mpfr_init(r20265); mpfr_init(r20266); mpfr_init(r20267); mpfr_init(r20268); mpfr_init(r20269); mpfr_init(r20270); mpfr_init(r20271); mpfr_init(r20272); mpfr_init(r20273); mpfr_init(r20274); } double f_dm(double x, double y, double z, double t, double a, double b) { mpfr_set_d(r20214, y, MPFR_RNDN); mpfr_set_d(r20215, z, MPFR_RNDN); mpfr_log(r20216, r20215, MPFR_RNDN); mpfr_mul(r20217, r20214, r20216, MPFR_RNDN); ; mpfr_set_si(r20219, mpfr_cmp(r20217, r20218) <= 0, MPFR_RNDN); mpfr_set_d(r20220, x, MPFR_RNDN); mpfr_div(r20221, r20220, r20214, MPFR_RNDN); mpfr_log(r20222, r20221, MPFR_RNDN); mpfr_set_d(r20223, a, MPFR_RNDN); mpfr_log(r20224, r20223, MPFR_RNDN); mpfr_set_d(r20225, t, MPFR_RNDN); ; mpfr_sub(r20227, r20225, r20226, MPFR_RNDN); mpfr_mul(r20228, r20224, r20227, MPFR_RNDN); mpfr_set_d(r20229, b, MPFR_RNDN); mpfr_mul(r20230, r20216, r20214, MPFR_RNDN); mpfr_sub(r20231, r20229, r20230, MPFR_RNDN); mpfr_sub(r20232, r20228, r20231, MPFR_RNDN); mpfr_add(r20233, r20222, r20232, MPFR_RNDN); mpfr_exp(r20234, r20233, MPFR_RNDN); ; mpfr_set_si(r20236, mpfr_cmp(r20217, r20235) <= 0, MPFR_RNDN); mpfr_pow(r20237, r20215, r20214, MPFR_RNDN); ; mpfr_pow(r20239, r20223, r20226, MPFR_RNDN); mpfr_div(r20240, r20238, r20239, MPFR_RNDN); mpfr_pow(r20241, r20240, r20226, MPFR_RNDN); mpfr_mul(r20242, r20237, r20241, MPFR_RNDN); mpfr_mul(r20243, r20224, r20224, MPFR_RNDN); mpfr_mul(r20244, r20225, r20225, MPFR_RNDN); ; mpfr_mul(r20246, r20244, r20245, MPFR_RNDN); mpfr_mul(r20247, r20243, r20246, MPFR_RNDN); mpfr_mul(r20248, r20225, r20224, MPFR_RNDN); mpfr_add(r20249, r20247, r20248, MPFR_RNDN); mpfr_mul(r20250, r20242, r20249, MPFR_RNDN); ; mpfr_pow(r20252, r20223, r20251, MPFR_RNDN); mpfr_mul(r20253, r20252, r20237, MPFR_RNDN); mpfr_add(r20254, r20250, r20253, MPFR_RNDN); mpfr_exp(r20255, r20229, MPFR_RNDN); mpfr_div(r20256, r20255, r20221, MPFR_RNDN); mpfr_div(r20257, r20254, r20256, MPFR_RNDN); ; mpfr_set_si(r20259, mpfr_cmp(r20217, r20258) <= 0, MPFR_RNDN); mpfr_pow(r20260, r20223, r20227, MPFR_RNDN); mpfr_exp(r20261, r20231, MPFR_RNDN); mpfr_div(r20262, r20260, r20261, MPFR_RNDN); mpfr_mul(r20263, r20221, r20262, MPFR_RNDN); ; mpfr_set_si(r20265, mpfr_cmp(r20217, r20264) <= 0, MPFR_RNDN); mpfr_div(r20266, r20237, r20214, MPFR_RNDN); mpfr_sub(r20267, r20228, r20229, MPFR_RNDN); mpfr_exp(r20268, r20267, MPFR_RNDN); mpfr_mul(r20269, r20266, r20268, MPFR_RNDN); mpfr_mul(r20270, r20220, r20269, MPFR_RNDN); if (mpfr_get_si(r20265, MPFR_RNDN)) { mpfr_set(r20271, r20270, MPFR_RNDN); } else { mpfr_set(r20271, r20234, MPFR_RNDN); }; if (mpfr_get_si(r20259, MPFR_RNDN)) { mpfr_set(r20272, r20263, MPFR_RNDN); } else { mpfr_set(r20272, r20271, MPFR_RNDN); }; if (mpfr_get_si(r20236, MPFR_RNDN)) { mpfr_set(r20273, r20257, MPFR_RNDN); } else { mpfr_set(r20273, r20272, MPFR_RNDN); }; if (mpfr_get_si(r20219, MPFR_RNDN)) { mpfr_set(r20274, r20234, MPFR_RNDN); } else { mpfr_set(r20274, r20273, MPFR_RNDN); }; return mpfr_get_d(r20274, MPFR_RNDN); }