#include #include #include #include #include char *name = "Numeric.SpecFunctions:incompleteBetaApprox from math-functions-0.1.5.2, B"; double f_if(float x, float y, float z, float t, float a, float b) { float r44974 = x; float r44975 = y; float r44976 = z; float r44977 = log(r44976); float r44978 = t; float r44979 = r44977 - r44978; float r44980 = r44975 * r44979; float r44981 = a; float r44982 = 1.0; float r44983 = r44982 - r44976; float r44984 = log(r44983); float r44985 = b; float r44986 = r44984 - r44985; float r44987 = r44981 * r44986; float r44988 = r44980 + r44987; float r44989 = exp(r44988); float r44990 = r44974 * r44989; return r44990; } double f_id(double x, double y, double z, double t, double a, double b) { double r44991 = x; double r44992 = y; double r44993 = z; double r44994 = log(r44993); double r44995 = t; double r44996 = r44994 - r44995; double r44997 = r44992 * r44996; double r44998 = a; double r44999 = 1.0; double r45000 = r44999 - r44993; double r45001 = log(r45000); double r45002 = b; double r45003 = r45001 - r45002; double r45004 = r44998 * r45003; double r45005 = r44997 + r45004; double r45006 = exp(r45005); double r45007 = r44991 * r45006; return r45007; } double f_of(float x, float y, float z, float t, float a, float b) { float r45008 = y; float r45009 = -2.1800019238241644e+116; bool r45010 = r45008 <= r45009; float r45011 = x; float r45012 = 1.0; float r45013 = z; float r45014 = r45012 - r45013; float r45015 = b; float r45016 = exp(r45015); float r45017 = r45014 / r45016; float r45018 = a; float r45019 = pow(r45017, r45018); float r45020 = r45011 * r45019; float r45021 = log(r45013); float r45022 = r45021 / r45008; float r45023 = 1; float r45024 = t; float r45025 = r45008 * r45024; float r45026 = r45023 / r45025; float r45027 = r45022 + r45026; float r45028 = exp(r45027); float r45029 = r45020 / r45028; float r45030 = 7.870676857854245e-219; bool r45031 = r45008 <= r45030; float r45032 = log(r45012); float r45033 = r45012 * r45013; float r45034 = 1/2; float r45035 = 2; float r45036 = pow(r45013, r45035); float r45037 = pow(r45012, r45035); float r45038 = r45036 / r45037; float r45039 = r45034 * r45038; float r45040 = r45033 + r45039; float r45041 = r45032 - r45040; float r45042 = r45041 - r45015; float r45043 = r45042 * r45018; float r45044 = exp(r45043); float r45045 = r45011 * r45044; float r45046 = r45021 - r45024; float r45047 = r45008 * r45046; float r45048 = exp(r45047); float r45049 = r45045 * r45048; float r45050 = 4.4393079520886555e+231; bool r45051 = r45008 <= r45050; float r45052 = log(r45014); float r45053 = r45052 - r45015; float r45054 = r45053 * r45018; float r45055 = exp(r45054); float r45056 = r45011 * r45055; float r45057 = r45056 * r45048; float r45058 = r45023 / r45015; float r45059 = exp(r45058); float r45060 = r45012 / r45059; float r45061 = r45013 / r45059; float r45062 = r45060 - r45061; float r45063 = pow(r45062, r45018); float r45064 = exp(r45024); float r45065 = r45013 / r45064; float r45066 = pow(r45065, r45008); float r45067 = r45011 * r45066; float r45068 = r45063 * r45067; float r45069 = r45051 ? r45057 : r45068; float r45070 = r45031 ? r45049 : r45069; float r45071 = r45010 ? r45029 : r45070; return r45071; } double f_od(double x, double y, double z, double t, double a, double b) { double r45072 = y; double r45073 = -2.1800019238241644e+116; bool r45074 = r45072 <= r45073; double r45075 = x; double r45076 = 1.0; double r45077 = z; double r45078 = r45076 - r45077; double r45079 = b; double r45080 = exp(r45079); double r45081 = r45078 / r45080; double r45082 = a; double r45083 = pow(r45081, r45082); double r45084 = r45075 * r45083; double r45085 = log(r45077); double r45086 = r45085 / r45072; double r45087 = 1; double r45088 = t; double r45089 = r45072 * r45088; double r45090 = r45087 / r45089; double r45091 = r45086 + r45090; double r45092 = exp(r45091); double r45093 = r45084 / r45092; double r45094 = 7.870676857854245e-219; bool r45095 = r45072 <= r45094; double r45096 = log(r45076); double r45097 = r45076 * r45077; double r45098 = 1/2; double r45099 = 2; double r45100 = pow(r45077, r45099); double r45101 = pow(r45076, r45099); double r45102 = r45100 / r45101; double r45103 = r45098 * r45102; double r45104 = r45097 + r45103; double r45105 = r45096 - r45104; double r45106 = r45105 - r45079; double r45107 = r45106 * r45082; double r45108 = exp(r45107); double r45109 = r45075 * r45108; double r45110 = r45085 - r45088; double r45111 = r45072 * r45110; double r45112 = exp(r45111); double r45113 = r45109 * r45112; double r45114 = 4.4393079520886555e+231; bool r45115 = r45072 <= r45114; double r45116 = log(r45078); double r45117 = r45116 - r45079; double r45118 = r45117 * r45082; double r45119 = exp(r45118); double r45120 = r45075 * r45119; double r45121 = r45120 * r45112; double r45122 = r45087 / r45079; double r45123 = exp(r45122); double r45124 = r45076 / r45123; double r45125 = r45077 / r45123; double r45126 = r45124 - r45125; double r45127 = pow(r45126, r45082); double r45128 = exp(r45088); double r45129 = r45077 / r45128; double r45130 = pow(r45129, r45072); double r45131 = r45075 * r45130; double r45132 = r45127 * r45131; double r45133 = r45115 ? r45121 : r45132; double r45134 = r45095 ? r45113 : r45133; double r45135 = r45074 ? r45093 : r45134; return r45135; } 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 r45136, r45137, r45138, r45139, r45140, r45141, r45142, r45143, r45144, r45145, r45146, r45147, r45148, r45149, r45150, r45151, r45152; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r45136); mpfr_init(r45137); mpfr_init(r45138); mpfr_init(r45139); mpfr_init(r45140); mpfr_init(r45141); mpfr_init(r45142); mpfr_init(r45143); mpfr_init_set_str(r45144, "1.0", 10, MPFR_RNDN); mpfr_init(r45145); mpfr_init(r45146); mpfr_init(r45147); mpfr_init(r45148); mpfr_init(r45149); mpfr_init(r45150); mpfr_init(r45151); mpfr_init(r45152); } double f_im(double x, double y, double z, double t, double a, double b) { mpfr_set_d(r45136, x, MPFR_RNDN); mpfr_set_d(r45137, y, MPFR_RNDN); mpfr_set_d(r45138, z, MPFR_RNDN); mpfr_log(r45139, r45138, MPFR_RNDN); mpfr_set_d(r45140, t, MPFR_RNDN); mpfr_sub(r45141, r45139, r45140, MPFR_RNDN); mpfr_mul(r45142, r45137, r45141, MPFR_RNDN); mpfr_set_d(r45143, a, MPFR_RNDN); ; mpfr_sub(r45145, r45144, r45138, MPFR_RNDN); mpfr_log(r45146, r45145, MPFR_RNDN); mpfr_set_d(r45147, b, MPFR_RNDN); mpfr_sub(r45148, r45146, r45147, MPFR_RNDN); mpfr_mul(r45149, r45143, r45148, MPFR_RNDN); mpfr_add(r45150, r45142, r45149, MPFR_RNDN); mpfr_exp(r45151, r45150, MPFR_RNDN); mpfr_mul(r45152, r45136, r45151, MPFR_RNDN); return mpfr_get_d(r45152, MPFR_RNDN); } static mpfr_t r45153, r45154, r45155, r45156, r45157, r45158, r45159, r45160, r45161, r45162, r45163, r45164, r45165, r45166, r45167, r45168, r45169, r45170, r45171, r45172, r45173, r45174, r45175, r45176, r45177, r45178, r45179, r45180, r45181, r45182, r45183, r45184, r45185, r45186, r45187, r45188, r45189, r45190, r45191, r45192, r45193, r45194, r45195, r45196, r45197, r45198, r45199, r45200, r45201, r45202, r45203, r45204, r45205, r45206, r45207, r45208, r45209, r45210, r45211, r45212, r45213, r45214, r45215, r45216; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r45153); mpfr_init_set_str(r45154, "-2.1800019238241644e+116", 10, MPFR_RNDN); mpfr_init(r45155); mpfr_init(r45156); mpfr_init_set_str(r45157, "1.0", 10, MPFR_RNDN); mpfr_init(r45158); mpfr_init(r45159); mpfr_init(r45160); mpfr_init(r45161); mpfr_init(r45162); mpfr_init(r45163); mpfr_init(r45164); mpfr_init(r45165); mpfr_init(r45166); mpfr_init(r45167); mpfr_init_set_str(r45168, "1", 10, MPFR_RNDN); mpfr_init(r45169); mpfr_init(r45170); mpfr_init(r45171); mpfr_init(r45172); mpfr_init(r45173); mpfr_init(r45174); mpfr_init_set_str(r45175, "7.870676857854245e-219", 10, MPFR_RNDN); mpfr_init(r45176); mpfr_init(r45177); mpfr_init(r45178); mpfr_init_set_str(r45179, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r45180, "2", 10, MPFR_RNDN); mpfr_init(r45181); mpfr_init(r45182); mpfr_init(r45183); mpfr_init(r45184); mpfr_init(r45185); mpfr_init(r45186); mpfr_init(r45187); mpfr_init(r45188); mpfr_init(r45189); mpfr_init(r45190); mpfr_init(r45191); mpfr_init(r45192); mpfr_init(r45193); mpfr_init(r45194); mpfr_init_set_str(r45195, "4.4393079520886555e+231", 10, MPFR_RNDN); mpfr_init(r45196); mpfr_init(r45197); mpfr_init(r45198); mpfr_init(r45199); mpfr_init(r45200); mpfr_init(r45201); mpfr_init(r45202); mpfr_init(r45203); mpfr_init(r45204); mpfr_init(r45205); mpfr_init(r45206); mpfr_init(r45207); mpfr_init(r45208); mpfr_init(r45209); mpfr_init(r45210); mpfr_init(r45211); mpfr_init(r45212); mpfr_init(r45213); mpfr_init(r45214); mpfr_init(r45215); mpfr_init(r45216); } double f_fm(double x, double y, double z, double t, double a, double b) { mpfr_set_d(r45153, y, MPFR_RNDN); ; mpfr_set_si(r45155, mpfr_cmp(r45153, r45154) <= 0, MPFR_RNDN); mpfr_set_d(r45156, x, MPFR_RNDN); ; mpfr_set_d(r45158, z, MPFR_RNDN); mpfr_sub(r45159, r45157, r45158, MPFR_RNDN); mpfr_set_d(r45160, b, MPFR_RNDN); mpfr_exp(r45161, r45160, MPFR_RNDN); mpfr_div(r45162, r45159, r45161, MPFR_RNDN); mpfr_set_d(r45163, a, MPFR_RNDN); mpfr_pow(r45164, r45162, r45163, MPFR_RNDN); mpfr_mul(r45165, r45156, r45164, MPFR_RNDN); mpfr_log(r45166, r45158, MPFR_RNDN); mpfr_div(r45167, r45166, r45153, MPFR_RNDN); ; mpfr_set_d(r45169, t, MPFR_RNDN); mpfr_mul(r45170, r45153, r45169, MPFR_RNDN); mpfr_div(r45171, r45168, r45170, MPFR_RNDN); mpfr_add(r45172, r45167, r45171, MPFR_RNDN); mpfr_exp(r45173, r45172, MPFR_RNDN); mpfr_div(r45174, r45165, r45173, MPFR_RNDN); ; mpfr_set_si(r45176, mpfr_cmp(r45153, r45175) <= 0, MPFR_RNDN); mpfr_log(r45177, r45157, MPFR_RNDN); mpfr_mul(r45178, r45157, r45158, MPFR_RNDN); ; ; mpfr_pow(r45181, r45158, r45180, MPFR_RNDN); mpfr_pow(r45182, r45157, r45180, MPFR_RNDN); mpfr_div(r45183, r45181, r45182, MPFR_RNDN); mpfr_mul(r45184, r45179, r45183, MPFR_RNDN); mpfr_add(r45185, r45178, r45184, MPFR_RNDN); mpfr_sub(r45186, r45177, r45185, MPFR_RNDN); mpfr_sub(r45187, r45186, r45160, MPFR_RNDN); mpfr_mul(r45188, r45187, r45163, MPFR_RNDN); mpfr_exp(r45189, r45188, MPFR_RNDN); mpfr_mul(r45190, r45156, r45189, MPFR_RNDN); mpfr_sub(r45191, r45166, r45169, MPFR_RNDN); mpfr_mul(r45192, r45153, r45191, MPFR_RNDN); mpfr_exp(r45193, r45192, MPFR_RNDN); mpfr_mul(r45194, r45190, r45193, MPFR_RNDN); ; mpfr_set_si(r45196, mpfr_cmp(r45153, r45195) <= 0, MPFR_RNDN); mpfr_log(r45197, r45159, MPFR_RNDN); mpfr_sub(r45198, r45197, r45160, MPFR_RNDN); mpfr_mul(r45199, r45198, r45163, MPFR_RNDN); mpfr_exp(r45200, r45199, MPFR_RNDN); mpfr_mul(r45201, r45156, r45200, MPFR_RNDN); mpfr_mul(r45202, r45201, r45193, MPFR_RNDN); mpfr_div(r45203, r45168, r45160, MPFR_RNDN); mpfr_exp(r45204, r45203, MPFR_RNDN); mpfr_div(r45205, r45157, r45204, MPFR_RNDN); mpfr_div(r45206, r45158, r45204, MPFR_RNDN); mpfr_sub(r45207, r45205, r45206, MPFR_RNDN); mpfr_pow(r45208, r45207, r45163, MPFR_RNDN); mpfr_exp(r45209, r45169, MPFR_RNDN); mpfr_div(r45210, r45158, r45209, MPFR_RNDN); mpfr_pow(r45211, r45210, r45153, MPFR_RNDN); mpfr_mul(r45212, r45156, r45211, MPFR_RNDN); mpfr_mul(r45213, r45208, r45212, MPFR_RNDN); if (mpfr_get_si(r45196, MPFR_RNDN)) { mpfr_set(r45214, r45202, MPFR_RNDN); } else { mpfr_set(r45214, r45213, MPFR_RNDN); }; if (mpfr_get_si(r45176, MPFR_RNDN)) { mpfr_set(r45215, r45194, MPFR_RNDN); } else { mpfr_set(r45215, r45214, MPFR_RNDN); }; if (mpfr_get_si(r45155, MPFR_RNDN)) { mpfr_set(r45216, r45174, MPFR_RNDN); } else { mpfr_set(r45216, r45215, MPFR_RNDN); }; return mpfr_get_d(r45216, MPFR_RNDN); } static mpfr_t r45217, r45218, r45219, r45220, r45221, r45222, r45223, r45224, r45225, r45226, r45227, r45228, r45229, r45230, r45231, r45232, r45233, r45234, r45235, r45236, r45237, r45238, r45239, r45240, r45241, r45242, r45243, r45244, r45245, r45246, r45247, r45248, r45249, r45250, r45251, r45252, r45253, r45254, r45255, r45256, r45257, r45258, r45259, r45260, r45261, r45262, r45263, r45264, r45265, r45266, r45267, r45268, r45269, r45270, r45271, r45272, r45273, r45274, r45275, r45276, r45277, r45278, r45279, r45280; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r45217); mpfr_init_set_str(r45218, "-2.1800019238241644e+116", 10, MPFR_RNDN); mpfr_init(r45219); mpfr_init(r45220); mpfr_init_set_str(r45221, "1.0", 10, MPFR_RNDN); mpfr_init(r45222); mpfr_init(r45223); mpfr_init(r45224); mpfr_init(r45225); mpfr_init(r45226); mpfr_init(r45227); mpfr_init(r45228); mpfr_init(r45229); mpfr_init(r45230); mpfr_init(r45231); mpfr_init_set_str(r45232, "1", 10, MPFR_RNDN); mpfr_init(r45233); mpfr_init(r45234); mpfr_init(r45235); mpfr_init(r45236); mpfr_init(r45237); mpfr_init(r45238); mpfr_init_set_str(r45239, "7.870676857854245e-219", 10, MPFR_RNDN); mpfr_init(r45240); mpfr_init(r45241); mpfr_init(r45242); mpfr_init_set_str(r45243, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r45244, "2", 10, MPFR_RNDN); mpfr_init(r45245); mpfr_init(r45246); mpfr_init(r45247); mpfr_init(r45248); mpfr_init(r45249); mpfr_init(r45250); mpfr_init(r45251); mpfr_init(r45252); mpfr_init(r45253); mpfr_init(r45254); mpfr_init(r45255); mpfr_init(r45256); mpfr_init(r45257); mpfr_init(r45258); mpfr_init_set_str(r45259, "4.4393079520886555e+231", 10, MPFR_RNDN); mpfr_init(r45260); mpfr_init(r45261); mpfr_init(r45262); mpfr_init(r45263); mpfr_init(r45264); mpfr_init(r45265); mpfr_init(r45266); mpfr_init(r45267); mpfr_init(r45268); mpfr_init(r45269); mpfr_init(r45270); mpfr_init(r45271); mpfr_init(r45272); mpfr_init(r45273); mpfr_init(r45274); mpfr_init(r45275); mpfr_init(r45276); mpfr_init(r45277); mpfr_init(r45278); mpfr_init(r45279); mpfr_init(r45280); } double f_dm(double x, double y, double z, double t, double a, double b) { mpfr_set_d(r45217, y, MPFR_RNDN); ; mpfr_set_si(r45219, mpfr_cmp(r45217, r45218) <= 0, MPFR_RNDN); mpfr_set_d(r45220, x, MPFR_RNDN); ; mpfr_set_d(r45222, z, MPFR_RNDN); mpfr_sub(r45223, r45221, r45222, MPFR_RNDN); mpfr_set_d(r45224, b, MPFR_RNDN); mpfr_exp(r45225, r45224, MPFR_RNDN); mpfr_div(r45226, r45223, r45225, MPFR_RNDN); mpfr_set_d(r45227, a, MPFR_RNDN); mpfr_pow(r45228, r45226, r45227, MPFR_RNDN); mpfr_mul(r45229, r45220, r45228, MPFR_RNDN); mpfr_log(r45230, r45222, MPFR_RNDN); mpfr_div(r45231, r45230, r45217, MPFR_RNDN); ; mpfr_set_d(r45233, t, MPFR_RNDN); mpfr_mul(r45234, r45217, r45233, MPFR_RNDN); mpfr_div(r45235, r45232, r45234, MPFR_RNDN); mpfr_add(r45236, r45231, r45235, MPFR_RNDN); mpfr_exp(r45237, r45236, MPFR_RNDN); mpfr_div(r45238, r45229, r45237, MPFR_RNDN); ; mpfr_set_si(r45240, mpfr_cmp(r45217, r45239) <= 0, MPFR_RNDN); mpfr_log(r45241, r45221, MPFR_RNDN); mpfr_mul(r45242, r45221, r45222, MPFR_RNDN); ; ; mpfr_pow(r45245, r45222, r45244, MPFR_RNDN); mpfr_pow(r45246, r45221, r45244, MPFR_RNDN); mpfr_div(r45247, r45245, r45246, MPFR_RNDN); mpfr_mul(r45248, r45243, r45247, MPFR_RNDN); mpfr_add(r45249, r45242, r45248, MPFR_RNDN); mpfr_sub(r45250, r45241, r45249, MPFR_RNDN); mpfr_sub(r45251, r45250, r45224, MPFR_RNDN); mpfr_mul(r45252, r45251, r45227, MPFR_RNDN); mpfr_exp(r45253, r45252, MPFR_RNDN); mpfr_mul(r45254, r45220, r45253, MPFR_RNDN); mpfr_sub(r45255, r45230, r45233, MPFR_RNDN); mpfr_mul(r45256, r45217, r45255, MPFR_RNDN); mpfr_exp(r45257, r45256, MPFR_RNDN); mpfr_mul(r45258, r45254, r45257, MPFR_RNDN); ; mpfr_set_si(r45260, mpfr_cmp(r45217, r45259) <= 0, MPFR_RNDN); mpfr_log(r45261, r45223, MPFR_RNDN); mpfr_sub(r45262, r45261, r45224, MPFR_RNDN); mpfr_mul(r45263, r45262, r45227, MPFR_RNDN); mpfr_exp(r45264, r45263, MPFR_RNDN); mpfr_mul(r45265, r45220, r45264, MPFR_RNDN); mpfr_mul(r45266, r45265, r45257, MPFR_RNDN); mpfr_div(r45267, r45232, r45224, MPFR_RNDN); mpfr_exp(r45268, r45267, MPFR_RNDN); mpfr_div(r45269, r45221, r45268, MPFR_RNDN); mpfr_div(r45270, r45222, r45268, MPFR_RNDN); mpfr_sub(r45271, r45269, r45270, MPFR_RNDN); mpfr_pow(r45272, r45271, r45227, MPFR_RNDN); mpfr_exp(r45273, r45233, MPFR_RNDN); mpfr_div(r45274, r45222, r45273, MPFR_RNDN); mpfr_pow(r45275, r45274, r45217, MPFR_RNDN); mpfr_mul(r45276, r45220, r45275, MPFR_RNDN); mpfr_mul(r45277, r45272, r45276, MPFR_RNDN); if (mpfr_get_si(r45260, MPFR_RNDN)) { mpfr_set(r45278, r45266, MPFR_RNDN); } else { mpfr_set(r45278, r45277, MPFR_RNDN); }; if (mpfr_get_si(r45240, MPFR_RNDN)) { mpfr_set(r45279, r45258, MPFR_RNDN); } else { mpfr_set(r45279, r45278, MPFR_RNDN); }; if (mpfr_get_si(r45219, MPFR_RNDN)) { mpfr_set(r45280, r45238, MPFR_RNDN); } else { mpfr_set(r45280, r45279, MPFR_RNDN); }; return mpfr_get_d(r45280, MPFR_RNDN); }