#include #include #include #include #include char *name = "Diagrams.Trail:splitAtParam from diagrams-lib-1.3.0.3, D"; double f_if(float x, float y) { float r52611 = 1.0; float r52612 = x; float r52613 = r52611 - r52612; float r52614 = y; float r52615 = r52613 * r52614; float r52616 = r52614 + r52611; float r52617 = r52615 / r52616; float r52618 = r52611 - r52617; return r52618; } double f_id(double x, double y) { double r52619 = 1.0; double r52620 = x; double r52621 = r52619 - r52620; double r52622 = y; double r52623 = r52621 * r52622; double r52624 = r52622 + r52619; double r52625 = r52623 / r52624; double r52626 = r52619 - r52625; return r52626; } double f_of(float x, float y) { float r52627 = y; float r52628 = -164028623.78812802; bool r52629 = r52627 <= r52628; float r52630 = 1.0; float r52631 = r52630 / r52627; float r52632 = x; float r52633 = r52631 + r52632; float r52634 = r52632 * r52630; float r52635 = r52634 / r52627; float r52636 = r52633 - r52635; float r52637 = 232001311.40530783; bool r52638 = r52627 <= r52637; float r52639 = r52630 - r52632; float r52640 = r52639 * r52627; float r52641 = r52627 + r52630; float r52642 = r52640 / r52641; float r52643 = r52630 - r52642; float r52644 = r52638 ? r52643 : r52636; float r52645 = r52629 ? r52636 : r52644; return r52645; } double f_od(double x, double y) { double r52646 = y; double r52647 = -164028623.78812802; bool r52648 = r52646 <= r52647; double r52649 = 1.0; double r52650 = r52649 / r52646; double r52651 = x; double r52652 = r52650 + r52651; double r52653 = r52651 * r52649; double r52654 = r52653 / r52646; double r52655 = r52652 - r52654; double r52656 = 232001311.40530783; bool r52657 = r52646 <= r52656; double r52658 = r52649 - r52651; double r52659 = r52658 * r52646; double r52660 = r52646 + r52649; double r52661 = r52659 / r52660; double r52662 = r52649 - r52661; double r52663 = r52657 ? r52662 : r52655; double r52664 = r52648 ? r52655 : r52663; return r52664; } 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 r52665, r52666, r52667, r52668, r52669, r52670, r52671, r52672; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r52665, "1.0", 10, MPFR_RNDN); mpfr_init(r52666); mpfr_init(r52667); mpfr_init(r52668); mpfr_init(r52669); mpfr_init(r52670); mpfr_init(r52671); mpfr_init(r52672); } double f_im(double x, double y) { ; mpfr_set_d(r52666, x, MPFR_RNDN); mpfr_sub(r52667, r52665, r52666, MPFR_RNDN); mpfr_set_d(r52668, y, MPFR_RNDN); mpfr_mul(r52669, r52667, r52668, MPFR_RNDN); mpfr_add(r52670, r52668, r52665, MPFR_RNDN); mpfr_div(r52671, r52669, r52670, MPFR_RNDN); mpfr_sub(r52672, r52665, r52671, MPFR_RNDN); return mpfr_get_d(r52672, MPFR_RNDN); } static mpfr_t r52673, r52674, r52675, r52676, r52677, r52678, r52679, r52680, r52681, r52682, r52683, r52684, r52685, r52686, r52687, r52688, r52689, r52690, r52691; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r52673); mpfr_init_set_str(r52674, "-164028623.78812802", 10, MPFR_RNDN); mpfr_init(r52675); mpfr_init_set_str(r52676, "1.0", 10, MPFR_RNDN); mpfr_init(r52677); mpfr_init(r52678); mpfr_init(r52679); mpfr_init(r52680); mpfr_init(r52681); mpfr_init(r52682); mpfr_init_set_str(r52683, "232001311.40530783", 10, MPFR_RNDN); mpfr_init(r52684); mpfr_init(r52685); mpfr_init(r52686); mpfr_init(r52687); mpfr_init(r52688); mpfr_init(r52689); mpfr_init(r52690); mpfr_init(r52691); } double f_fm(double x, double y) { mpfr_set_d(r52673, y, MPFR_RNDN); ; mpfr_set_si(r52675, mpfr_cmp(r52673, r52674) <= 0, MPFR_RNDN); ; mpfr_div(r52677, r52676, r52673, MPFR_RNDN); mpfr_set_d(r52678, x, MPFR_RNDN); mpfr_add(r52679, r52677, r52678, MPFR_RNDN); mpfr_mul(r52680, r52678, r52676, MPFR_RNDN); mpfr_div(r52681, r52680, r52673, MPFR_RNDN); mpfr_sub(r52682, r52679, r52681, MPFR_RNDN); ; mpfr_set_si(r52684, mpfr_cmp(r52673, r52683) <= 0, MPFR_RNDN); mpfr_sub(r52685, r52676, r52678, MPFR_RNDN); mpfr_mul(r52686, r52685, r52673, MPFR_RNDN); mpfr_add(r52687, r52673, r52676, MPFR_RNDN); mpfr_div(r52688, r52686, r52687, MPFR_RNDN); mpfr_sub(r52689, r52676, r52688, MPFR_RNDN); if (mpfr_get_si(r52684, MPFR_RNDN)) { mpfr_set(r52690, r52689, MPFR_RNDN); } else { mpfr_set(r52690, r52682, MPFR_RNDN); }; if (mpfr_get_si(r52675, MPFR_RNDN)) { mpfr_set(r52691, r52682, MPFR_RNDN); } else { mpfr_set(r52691, r52690, MPFR_RNDN); }; return mpfr_get_d(r52691, MPFR_RNDN); } static mpfr_t r52692, r52693, r52694, r52695, r52696, r52697, r52698, r52699, r52700, r52701, r52702, r52703, r52704, r52705, r52706, r52707, r52708, r52709, r52710; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r52692); mpfr_init_set_str(r52693, "-164028623.78812802", 10, MPFR_RNDN); mpfr_init(r52694); mpfr_init_set_str(r52695, "1.0", 10, MPFR_RNDN); mpfr_init(r52696); mpfr_init(r52697); mpfr_init(r52698); mpfr_init(r52699); mpfr_init(r52700); mpfr_init(r52701); mpfr_init_set_str(r52702, "232001311.40530783", 10, MPFR_RNDN); mpfr_init(r52703); mpfr_init(r52704); mpfr_init(r52705); mpfr_init(r52706); mpfr_init(r52707); mpfr_init(r52708); mpfr_init(r52709); mpfr_init(r52710); } double f_dm(double x, double y) { mpfr_set_d(r52692, y, MPFR_RNDN); ; mpfr_set_si(r52694, mpfr_cmp(r52692, r52693) <= 0, MPFR_RNDN); ; mpfr_div(r52696, r52695, r52692, MPFR_RNDN); mpfr_set_d(r52697, x, MPFR_RNDN); mpfr_add(r52698, r52696, r52697, MPFR_RNDN); mpfr_mul(r52699, r52697, r52695, MPFR_RNDN); mpfr_div(r52700, r52699, r52692, MPFR_RNDN); mpfr_sub(r52701, r52698, r52700, MPFR_RNDN); ; mpfr_set_si(r52703, mpfr_cmp(r52692, r52702) <= 0, MPFR_RNDN); mpfr_sub(r52704, r52695, r52697, MPFR_RNDN); mpfr_mul(r52705, r52704, r52692, MPFR_RNDN); mpfr_add(r52706, r52692, r52695, MPFR_RNDN); mpfr_div(r52707, r52705, r52706, MPFR_RNDN); mpfr_sub(r52708, r52695, r52707, MPFR_RNDN); if (mpfr_get_si(r52703, MPFR_RNDN)) { mpfr_set(r52709, r52708, MPFR_RNDN); } else { mpfr_set(r52709, r52701, MPFR_RNDN); }; if (mpfr_get_si(r52694, MPFR_RNDN)) { mpfr_set(r52710, r52701, MPFR_RNDN); } else { mpfr_set(r52710, r52709, MPFR_RNDN); }; return mpfr_get_d(r52710, MPFR_RNDN); }