#include #include #include #include #include char *name = "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2"; double f_if(float x, float y, float z, float t, float a, float b, float c) { float r23037 = x; float r23038 = y; float r23039 = 2.0; float r23040 = z; float r23041 = t; float r23042 = a; float r23043 = r23041 + r23042; float r23044 = sqrt(r23043); float r23045 = r23040 * r23044; float r23046 = r23045 / r23041; float r23047 = b; float r23048 = c; float r23049 = r23047 - r23048; float r23050 = 5.0; float r23051 = 6.0; float r23052 = r23050 / r23051; float r23053 = r23042 + r23052; float r23054 = 3.0; float r23055 = r23041 * r23054; float r23056 = r23039 / r23055; float r23057 = r23053 - r23056; float r23058 = r23049 * r23057; float r23059 = r23046 - r23058; float r23060 = r23039 * r23059; float r23061 = exp(r23060); float r23062 = r23038 * r23061; float r23063 = r23037 + r23062; float r23064 = r23037 / r23063; return r23064; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r23065 = x; double r23066 = y; double r23067 = 2.0; double r23068 = z; double r23069 = t; double r23070 = a; double r23071 = r23069 + r23070; double r23072 = sqrt(r23071); double r23073 = r23068 * r23072; double r23074 = r23073 / r23069; double r23075 = b; double r23076 = c; double r23077 = r23075 - r23076; double r23078 = 5.0; double r23079 = 6.0; double r23080 = r23078 / r23079; double r23081 = r23070 + r23080; double r23082 = 3.0; double r23083 = r23069 * r23082; double r23084 = r23067 / r23083; double r23085 = r23081 - r23084; double r23086 = r23077 * r23085; double r23087 = r23074 - r23086; double r23088 = r23067 * r23087; double r23089 = exp(r23088); double r23090 = r23066 * r23089; double r23091 = r23065 + r23090; double r23092 = r23065 / r23091; return r23092; } double f_of(float x, float y, float z, float t, float a, float b, float c) { float r23093 = x; float r23094 = y; float r23095 = 2.0; float r23096 = a; float r23097 = 5.0; float r23098 = 6.0; float r23099 = r23097 / r23098; float r23100 = r23096 + r23099; float r23101 = 3.0; float r23102 = t; float r23103 = r23101 * r23102; float r23104 = r23095 / r23103; float r23105 = r23100 - r23104; float r23106 = c; float r23107 = b; float r23108 = r23106 - r23107; float r23109 = r23102 + r23096; float r23110 = sqrt(r23109); float r23111 = z; float r23112 = r23111 / r23102; float r23113 = r23110 * r23112; float r23114 = fma(r23105, r23108, r23113); float r23115 = 0; float r23116 = r23114 + r23115; float r23117 = r23095 * r23116; float r23118 = exp(r23117); float r23119 = r23094 * r23118; float r23120 = r23093 + r23119; float r23121 = r23093 / r23120; return r23121; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r23122 = x; double r23123 = y; double r23124 = 2.0; double r23125 = a; double r23126 = 5.0; double r23127 = 6.0; double r23128 = r23126 / r23127; double r23129 = r23125 + r23128; double r23130 = 3.0; double r23131 = t; double r23132 = r23130 * r23131; double r23133 = r23124 / r23132; double r23134 = r23129 - r23133; double r23135 = c; double r23136 = b; double r23137 = r23135 - r23136; double r23138 = r23131 + r23125; double r23139 = sqrt(r23138); double r23140 = z; double r23141 = r23140 / r23131; double r23142 = r23139 * r23141; double r23143 = fma(r23134, r23137, r23142); double r23144 = 0; double r23145 = r23143 + r23144; double r23146 = r23124 * r23145; double r23147 = exp(r23146); double r23148 = r23123 * r23147; double r23149 = r23122 + r23148; double r23150 = r23122 / r23149; return r23150; } 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 r23151, r23152, r23153, r23154, r23155, r23156, r23157, r23158, r23159, r23160, r23161, r23162, r23163, r23164, r23165, r23166, r23167, r23168, r23169, r23170, r23171, r23172, r23173, r23174, r23175, r23176, r23177, r23178; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r23151); mpfr_init(r23152); mpfr_init_set_str(r23153, "2.0", 10, MPFR_RNDN); mpfr_init(r23154); mpfr_init(r23155); mpfr_init(r23156); mpfr_init(r23157); mpfr_init(r23158); mpfr_init(r23159); mpfr_init(r23160); mpfr_init(r23161); mpfr_init(r23162); mpfr_init(r23163); mpfr_init_set_str(r23164, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23165, "6.0", 10, MPFR_RNDN); mpfr_init(r23166); mpfr_init(r23167); mpfr_init_set_str(r23168, "3.0", 10, MPFR_RNDN); mpfr_init(r23169); mpfr_init(r23170); mpfr_init(r23171); mpfr_init(r23172); mpfr_init(r23173); mpfr_init(r23174); mpfr_init(r23175); mpfr_init(r23176); mpfr_init(r23177); mpfr_init(r23178); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23151, x, MPFR_RNDN); mpfr_set_d(r23152, y, MPFR_RNDN); ; mpfr_set_d(r23154, z, MPFR_RNDN); mpfr_set_d(r23155, t, MPFR_RNDN); mpfr_set_d(r23156, a, MPFR_RNDN); mpfr_add(r23157, r23155, r23156, MPFR_RNDN); mpfr_sqrt(r23158, r23157, MPFR_RNDN); mpfr_mul(r23159, r23154, r23158, MPFR_RNDN); mpfr_div(r23160, r23159, r23155, MPFR_RNDN); mpfr_set_d(r23161, b, MPFR_RNDN); mpfr_set_d(r23162, c, MPFR_RNDN); mpfr_sub(r23163, r23161, r23162, MPFR_RNDN); ; ; mpfr_div(r23166, r23164, r23165, MPFR_RNDN); mpfr_add(r23167, r23156, r23166, MPFR_RNDN); ; mpfr_mul(r23169, r23155, r23168, MPFR_RNDN); mpfr_div(r23170, r23153, r23169, MPFR_RNDN); mpfr_sub(r23171, r23167, r23170, MPFR_RNDN); mpfr_mul(r23172, r23163, r23171, MPFR_RNDN); mpfr_sub(r23173, r23160, r23172, MPFR_RNDN); mpfr_mul(r23174, r23153, r23173, MPFR_RNDN); mpfr_exp(r23175, r23174, MPFR_RNDN); mpfr_mul(r23176, r23152, r23175, MPFR_RNDN); mpfr_add(r23177, r23151, r23176, MPFR_RNDN); mpfr_div(r23178, r23151, r23177, MPFR_RNDN); return mpfr_get_d(r23178, MPFR_RNDN); } static mpfr_t r23179, r23180, r23181, r23182, r23183, r23184, r23185, r23186, r23187, r23188, r23189, r23190, r23191, r23192, r23193, r23194, r23195, r23196, r23197, r23198, r23199, r23200, r23201, r23202, r23203, r23204, r23205, r23206, r23207; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r23179); mpfr_init(r23180); mpfr_init_set_str(r23181, "2.0", 10, MPFR_RNDN); mpfr_init(r23182); mpfr_init_set_str(r23183, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23184, "6.0", 10, MPFR_RNDN); mpfr_init(r23185); mpfr_init(r23186); mpfr_init_set_str(r23187, "3.0", 10, MPFR_RNDN); mpfr_init(r23188); mpfr_init(r23189); mpfr_init(r23190); mpfr_init(r23191); mpfr_init(r23192); mpfr_init(r23193); mpfr_init(r23194); mpfr_init(r23195); mpfr_init(r23196); mpfr_init(r23197); mpfr_init(r23198); mpfr_init(r23199); mpfr_init(r23200); mpfr_init_set_str(r23201, "0", 10, MPFR_RNDN); mpfr_init(r23202); mpfr_init(r23203); mpfr_init(r23204); mpfr_init(r23205); mpfr_init(r23206); mpfr_init(r23207); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23179, x, MPFR_RNDN); mpfr_set_d(r23180, y, MPFR_RNDN); ; mpfr_set_d(r23182, a, MPFR_RNDN); ; ; mpfr_div(r23185, r23183, r23184, MPFR_RNDN); mpfr_add(r23186, r23182, r23185, MPFR_RNDN); ; mpfr_set_d(r23188, t, MPFR_RNDN); mpfr_mul(r23189, r23187, r23188, MPFR_RNDN); mpfr_div(r23190, r23181, r23189, MPFR_RNDN); mpfr_sub(r23191, r23186, r23190, MPFR_RNDN); mpfr_set_d(r23192, c, MPFR_RNDN); mpfr_set_d(r23193, b, MPFR_RNDN); mpfr_sub(r23194, r23192, r23193, MPFR_RNDN); mpfr_add(r23195, r23188, r23182, MPFR_RNDN); mpfr_sqrt(r23196, r23195, MPFR_RNDN); mpfr_set_d(r23197, z, MPFR_RNDN); mpfr_div(r23198, r23197, r23188, MPFR_RNDN); mpfr_mul(r23199, r23196, r23198, MPFR_RNDN); mpfr_fma(r23200, r23191, r23194, r23199, MPFR_RNDN); ; mpfr_add(r23202, r23200, r23201, MPFR_RNDN); mpfr_mul(r23203, r23181, r23202, MPFR_RNDN); mpfr_exp(r23204, r23203, MPFR_RNDN); mpfr_mul(r23205, r23180, r23204, MPFR_RNDN); mpfr_add(r23206, r23179, r23205, MPFR_RNDN); mpfr_div(r23207, r23179, r23206, MPFR_RNDN); return mpfr_get_d(r23207, MPFR_RNDN); } static mpfr_t r23208, r23209, r23210, r23211, r23212, r23213, r23214, r23215, r23216, r23217, r23218, r23219, r23220, r23221, r23222, r23223, r23224, r23225, r23226, r23227, r23228, r23229, r23230, r23231, r23232, r23233, r23234, r23235, r23236; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r23208); mpfr_init(r23209); mpfr_init_set_str(r23210, "2.0", 10, MPFR_RNDN); mpfr_init(r23211); mpfr_init_set_str(r23212, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23213, "6.0", 10, MPFR_RNDN); mpfr_init(r23214); mpfr_init(r23215); mpfr_init_set_str(r23216, "3.0", 10, MPFR_RNDN); mpfr_init(r23217); mpfr_init(r23218); mpfr_init(r23219); mpfr_init(r23220); mpfr_init(r23221); mpfr_init(r23222); mpfr_init(r23223); mpfr_init(r23224); mpfr_init(r23225); mpfr_init(r23226); mpfr_init(r23227); mpfr_init(r23228); mpfr_init(r23229); mpfr_init_set_str(r23230, "0", 10, MPFR_RNDN); mpfr_init(r23231); mpfr_init(r23232); mpfr_init(r23233); mpfr_init(r23234); mpfr_init(r23235); mpfr_init(r23236); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23208, x, MPFR_RNDN); mpfr_set_d(r23209, y, MPFR_RNDN); ; mpfr_set_d(r23211, a, MPFR_RNDN); ; ; mpfr_div(r23214, r23212, r23213, MPFR_RNDN); mpfr_add(r23215, r23211, r23214, MPFR_RNDN); ; mpfr_set_d(r23217, t, MPFR_RNDN); mpfr_mul(r23218, r23216, r23217, MPFR_RNDN); mpfr_div(r23219, r23210, r23218, MPFR_RNDN); mpfr_sub(r23220, r23215, r23219, MPFR_RNDN); mpfr_set_d(r23221, c, MPFR_RNDN); mpfr_set_d(r23222, b, MPFR_RNDN); mpfr_sub(r23223, r23221, r23222, MPFR_RNDN); mpfr_add(r23224, r23217, r23211, MPFR_RNDN); mpfr_sqrt(r23225, r23224, MPFR_RNDN); mpfr_set_d(r23226, z, MPFR_RNDN); mpfr_div(r23227, r23226, r23217, MPFR_RNDN); mpfr_mul(r23228, r23225, r23227, MPFR_RNDN); mpfr_fma(r23229, r23220, r23223, r23228, MPFR_RNDN); ; mpfr_add(r23231, r23229, r23230, MPFR_RNDN); mpfr_mul(r23232, r23210, r23231, MPFR_RNDN); mpfr_exp(r23233, r23232, MPFR_RNDN); mpfr_mul(r23234, r23209, r23233, MPFR_RNDN); mpfr_add(r23235, r23208, r23234, MPFR_RNDN); mpfr_div(r23236, r23208, r23235, MPFR_RNDN); return mpfr_get_d(r23236, MPFR_RNDN); }