#include #include #include #include #include char *name = "AI.Clustering.Hierarchical.Internal:ward from clustering-0.2.1"; double f_if(float x, float y, float z, float t, float a, float b) { float r58255 = x; float r58256 = y; float r58257 = r58255 + r58256; float r58258 = z; float r58259 = r58257 * r58258; float r58260 = t; float r58261 = r58260 + r58256; float r58262 = a; float r58263 = r58261 * r58262; float r58264 = r58259 + r58263; float r58265 = b; float r58266 = r58256 * r58265; float r58267 = r58264 - r58266; float r58268 = r58255 + r58260; float r58269 = r58268 + r58256; float r58270 = r58267 / r58269; return r58270; } double f_id(double x, double y, double z, double t, double a, double b) { double r58271 = x; double r58272 = y; double r58273 = r58271 + r58272; double r58274 = z; double r58275 = r58273 * r58274; double r58276 = t; double r58277 = r58276 + r58272; double r58278 = a; double r58279 = r58277 * r58278; double r58280 = r58275 + r58279; double r58281 = b; double r58282 = r58272 * r58281; double r58283 = r58280 - r58282; double r58284 = r58271 + r58276; double r58285 = r58284 + r58272; double r58286 = r58283 / r58285; return r58286; } double f_of(float x, float y, float z, float t, float a, float b) { float r58287 = y; float r58288 = -2.4207781418622835e+64; bool r58289 = r58287 <= r58288; float r58290 = z; float r58291 = a; float r58292 = r58290 + r58291; float r58293 = b; float r58294 = r58292 - r58293; float r58295 = 4.990162669943169e+125; bool r58296 = r58287 <= r58295; float r58297 = 1; float r58298 = x; float r58299 = t; float r58300 = r58298 + r58299; float r58301 = r58300 + r58287; float r58302 = r58298 + r58287; float r58303 = r58302 * r58290; float r58304 = r58299 + r58287; float r58305 = r58304 * r58291; float r58306 = r58303 + r58305; float r58307 = r58287 * r58293; float r58308 = r58306 - r58307; float r58309 = r58301 / r58308; float r58310 = r58297 / r58309; float r58311 = r58296 ? r58310 : r58294; float r58312 = r58289 ? r58294 : r58311; return r58312; } double f_od(double x, double y, double z, double t, double a, double b) { double r58313 = y; double r58314 = -2.4207781418622835e+64; bool r58315 = r58313 <= r58314; double r58316 = z; double r58317 = a; double r58318 = r58316 + r58317; double r58319 = b; double r58320 = r58318 - r58319; double r58321 = 4.990162669943169e+125; bool r58322 = r58313 <= r58321; double r58323 = 1; double r58324 = x; double r58325 = t; double r58326 = r58324 + r58325; double r58327 = r58326 + r58313; double r58328 = r58324 + r58313; double r58329 = r58328 * r58316; double r58330 = r58325 + r58313; double r58331 = r58330 * r58317; double r58332 = r58329 + r58331; double r58333 = r58313 * r58319; double r58334 = r58332 - r58333; double r58335 = r58327 / r58334; double r58336 = r58323 / r58335; double r58337 = r58322 ? r58336 : r58320; double r58338 = r58315 ? r58320 : r58337; return r58338; } 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 r58339, r58340, r58341, r58342, r58343, r58344, r58345, r58346, r58347, r58348, r58349, r58350, r58351, r58352, r58353, r58354; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r58339); mpfr_init(r58340); mpfr_init(r58341); mpfr_init(r58342); mpfr_init(r58343); mpfr_init(r58344); mpfr_init(r58345); mpfr_init(r58346); mpfr_init(r58347); mpfr_init(r58348); mpfr_init(r58349); mpfr_init(r58350); mpfr_init(r58351); mpfr_init(r58352); mpfr_init(r58353); mpfr_init(r58354); } double f_im(double x, double y, double z, double t, double a, double b) { mpfr_set_d(r58339, x, MPFR_RNDN); mpfr_set_d(r58340, y, MPFR_RNDN); mpfr_add(r58341, r58339, r58340, MPFR_RNDN); mpfr_set_d(r58342, z, MPFR_RNDN); mpfr_mul(r58343, r58341, r58342, MPFR_RNDN); mpfr_set_d(r58344, t, MPFR_RNDN); mpfr_add(r58345, r58344, r58340, MPFR_RNDN); mpfr_set_d(r58346, a, MPFR_RNDN); mpfr_mul(r58347, r58345, r58346, MPFR_RNDN); mpfr_add(r58348, r58343, r58347, MPFR_RNDN); mpfr_set_d(r58349, b, MPFR_RNDN); mpfr_mul(r58350, r58340, r58349, MPFR_RNDN); mpfr_sub(r58351, r58348, r58350, MPFR_RNDN); mpfr_add(r58352, r58339, r58344, MPFR_RNDN); mpfr_add(r58353, r58352, r58340, MPFR_RNDN); mpfr_div(r58354, r58351, r58353, MPFR_RNDN); return mpfr_get_d(r58354, MPFR_RNDN); } static mpfr_t r58355, r58356, r58357, r58358, r58359, r58360, r58361, r58362, r58363, r58364, r58365, r58366, r58367, r58368, r58369, r58370, r58371, r58372, r58373, r58374, r58375, r58376, r58377, r58378, r58379, r58380; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r58355); mpfr_init_set_str(r58356, "-2.4207781418622835e+64", 10, MPFR_RNDN); mpfr_init(r58357); mpfr_init(r58358); mpfr_init(r58359); mpfr_init(r58360); mpfr_init(r58361); mpfr_init(r58362); mpfr_init_set_str(r58363, "4.990162669943169e+125", 10, MPFR_RNDN); mpfr_init(r58364); mpfr_init_set_str(r58365, "1", 10, MPFR_RNDN); mpfr_init(r58366); mpfr_init(r58367); mpfr_init(r58368); mpfr_init(r58369); mpfr_init(r58370); mpfr_init(r58371); mpfr_init(r58372); mpfr_init(r58373); mpfr_init(r58374); mpfr_init(r58375); mpfr_init(r58376); mpfr_init(r58377); mpfr_init(r58378); mpfr_init(r58379); mpfr_init(r58380); } double f_fm(double x, double y, double z, double t, double a, double b) { mpfr_set_d(r58355, y, MPFR_RNDN); ; mpfr_set_si(r58357, mpfr_cmp(r58355, r58356) <= 0, MPFR_RNDN); mpfr_set_d(r58358, z, MPFR_RNDN); mpfr_set_d(r58359, a, MPFR_RNDN); mpfr_add(r58360, r58358, r58359, MPFR_RNDN); mpfr_set_d(r58361, b, MPFR_RNDN); mpfr_sub(r58362, r58360, r58361, MPFR_RNDN); ; mpfr_set_si(r58364, mpfr_cmp(r58355, r58363) <= 0, MPFR_RNDN); ; mpfr_set_d(r58366, x, MPFR_RNDN); mpfr_set_d(r58367, t, MPFR_RNDN); mpfr_add(r58368, r58366, r58367, MPFR_RNDN); mpfr_add(r58369, r58368, r58355, MPFR_RNDN); mpfr_add(r58370, r58366, r58355, MPFR_RNDN); mpfr_mul(r58371, r58370, r58358, MPFR_RNDN); mpfr_add(r58372, r58367, r58355, MPFR_RNDN); mpfr_mul(r58373, r58372, r58359, MPFR_RNDN); mpfr_add(r58374, r58371, r58373, MPFR_RNDN); mpfr_mul(r58375, r58355, r58361, MPFR_RNDN); mpfr_sub(r58376, r58374, r58375, MPFR_RNDN); mpfr_div(r58377, r58369, r58376, MPFR_RNDN); mpfr_div(r58378, r58365, r58377, MPFR_RNDN); if (mpfr_get_si(r58364, MPFR_RNDN)) { mpfr_set(r58379, r58378, MPFR_RNDN); } else { mpfr_set(r58379, r58362, MPFR_RNDN); }; if (mpfr_get_si(r58357, MPFR_RNDN)) { mpfr_set(r58380, r58362, MPFR_RNDN); } else { mpfr_set(r58380, r58379, MPFR_RNDN); }; return mpfr_get_d(r58380, MPFR_RNDN); } static mpfr_t r58381, r58382, r58383, r58384, r58385, r58386, r58387, r58388, r58389, r58390, r58391, r58392, r58393, r58394, r58395, r58396, r58397, r58398, r58399, r58400, r58401, r58402, r58403, r58404, r58405, r58406; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r58381); mpfr_init_set_str(r58382, "-2.4207781418622835e+64", 10, MPFR_RNDN); mpfr_init(r58383); mpfr_init(r58384); mpfr_init(r58385); mpfr_init(r58386); mpfr_init(r58387); mpfr_init(r58388); mpfr_init_set_str(r58389, "4.990162669943169e+125", 10, MPFR_RNDN); mpfr_init(r58390); mpfr_init_set_str(r58391, "1", 10, MPFR_RNDN); mpfr_init(r58392); mpfr_init(r58393); mpfr_init(r58394); mpfr_init(r58395); mpfr_init(r58396); mpfr_init(r58397); mpfr_init(r58398); mpfr_init(r58399); mpfr_init(r58400); mpfr_init(r58401); mpfr_init(r58402); mpfr_init(r58403); mpfr_init(r58404); mpfr_init(r58405); mpfr_init(r58406); } double f_dm(double x, double y, double z, double t, double a, double b) { mpfr_set_d(r58381, y, MPFR_RNDN); ; mpfr_set_si(r58383, mpfr_cmp(r58381, r58382) <= 0, MPFR_RNDN); mpfr_set_d(r58384, z, MPFR_RNDN); mpfr_set_d(r58385, a, MPFR_RNDN); mpfr_add(r58386, r58384, r58385, MPFR_RNDN); mpfr_set_d(r58387, b, MPFR_RNDN); mpfr_sub(r58388, r58386, r58387, MPFR_RNDN); ; mpfr_set_si(r58390, mpfr_cmp(r58381, r58389) <= 0, MPFR_RNDN); ; mpfr_set_d(r58392, x, MPFR_RNDN); mpfr_set_d(r58393, t, MPFR_RNDN); mpfr_add(r58394, r58392, r58393, MPFR_RNDN); mpfr_add(r58395, r58394, r58381, MPFR_RNDN); mpfr_add(r58396, r58392, r58381, MPFR_RNDN); mpfr_mul(r58397, r58396, r58384, MPFR_RNDN); mpfr_add(r58398, r58393, r58381, MPFR_RNDN); mpfr_mul(r58399, r58398, r58385, MPFR_RNDN); mpfr_add(r58400, r58397, r58399, MPFR_RNDN); mpfr_mul(r58401, r58381, r58387, MPFR_RNDN); mpfr_sub(r58402, r58400, r58401, MPFR_RNDN); mpfr_div(r58403, r58395, r58402, MPFR_RNDN); mpfr_div(r58404, r58391, r58403, MPFR_RNDN); if (mpfr_get_si(r58390, MPFR_RNDN)) { mpfr_set(r58405, r58404, MPFR_RNDN); } else { mpfr_set(r58405, r58388, MPFR_RNDN); }; if (mpfr_get_si(r58383, MPFR_RNDN)) { mpfr_set(r58406, r58388, MPFR_RNDN); } else { mpfr_set(r58406, r58405, MPFR_RNDN); }; return mpfr_get_d(r58406, MPFR_RNDN); }