#include #include #include #include #include char *name = "Asymptote C"; double f_if(float x) { float r26608 = x; float r26609 = 1; float r26610 = r26608 + r26609; float r26611 = r26608 / r26610; float r26612 = r26608 - r26609; float r26613 = r26610 / r26612; float r26614 = r26611 - r26613; return r26614; } double f_id(double x) { double r26615 = x; double r26616 = 1; double r26617 = r26615 + r26616; double r26618 = r26615 / r26617; double r26619 = r26615 - r26616; double r26620 = r26617 / r26619; double r26621 = r26618 - r26620; return r26621; } double f_of(float x) { float r26622 = 1; float r26623 = 3; float r26624 = x; float r26625 = r26623 / r26624; float r26626 = r26622 + r26625; float r26627 = -r26622; float r26628 = r26624 * r26624; float r26629 = r26627 / r26628; float r26630 = -r26625; float r26631 = fma(r26626, r26629, r26630); float r26632 = -1.6839866920797528e-11; bool r26633 = r26631 <= r26632; float r26634 = r26623 + r26624; float r26635 = fma(r26624, r26634, r26622); float r26636 = 1.4542143759088224e-05; bool r26637 = r26631 <= r26636; float r26638 = r26624 + r26622; float r26639 = r26624 / r26638; float r26640 = r26639 * r26639; float r26641 = r26624 - r26622; float r26642 = r26638 / r26641; float r26643 = r26642 * r26642; float r26644 = r26640 - r26643; float r26645 = r26639 + r26642; float r26646 = r26644 / r26645; float r26647 = r26637 ? r26631 : r26646; float r26648 = r26633 ? r26635 : r26647; return r26648; } double f_od(double x) { double r26649 = 1; double r26650 = 3; double r26651 = x; double r26652 = r26650 / r26651; double r26653 = r26649 + r26652; double r26654 = -r26649; double r26655 = r26651 * r26651; double r26656 = r26654 / r26655; double r26657 = -r26652; double r26658 = fma(r26653, r26656, r26657); double r26659 = -1.6839866920797528e-11; bool r26660 = r26658 <= r26659; double r26661 = r26650 + r26651; double r26662 = fma(r26651, r26661, r26649); double r26663 = 1.4542143759088224e-05; bool r26664 = r26658 <= r26663; double r26665 = r26651 + r26649; double r26666 = r26651 / r26665; double r26667 = r26666 * r26666; double r26668 = r26651 - r26649; double r26669 = r26665 / r26668; double r26670 = r26669 * r26669; double r26671 = r26667 - r26670; double r26672 = r26666 + r26669; double r26673 = r26671 / r26672; double r26674 = r26664 ? r26658 : r26673; double r26675 = r26660 ? r26662 : r26674; return r26675; } 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 r26676, r26677, r26678, r26679, r26680, r26681, r26682; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r26676); mpfr_init_set_str(r26677, "1", 10, MPFR_RNDN); mpfr_init(r26678); mpfr_init(r26679); mpfr_init(r26680); mpfr_init(r26681); mpfr_init(r26682); } double f_im(double x) { mpfr_set_d(r26676, x, MPFR_RNDN); ; mpfr_add(r26678, r26676, r26677, MPFR_RNDN); mpfr_div(r26679, r26676, r26678, MPFR_RNDN); mpfr_sub(r26680, r26676, r26677, MPFR_RNDN); mpfr_div(r26681, r26678, r26680, MPFR_RNDN); mpfr_sub(r26682, r26679, r26681, MPFR_RNDN); return mpfr_get_d(r26682, MPFR_RNDN); } static mpfr_t r26683, r26684, r26685, r26686, r26687, r26688, r26689, r26690, r26691, r26692, r26693, r26694, r26695, r26696, r26697, r26698, r26699, r26700, r26701, r26702, r26703, r26704, r26705, r26706, r26707, r26708, r26709; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r26683, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26684, "3", 10, MPFR_RNDN); mpfr_init(r26685); mpfr_init(r26686); mpfr_init(r26687); mpfr_init(r26688); mpfr_init(r26689); mpfr_init(r26690); mpfr_init(r26691); mpfr_init(r26692); mpfr_init_set_str(r26693, "-1.6839866920797528e-11", 10, MPFR_RNDN); mpfr_init(r26694); mpfr_init(r26695); mpfr_init(r26696); mpfr_init_set_str(r26697, "1.4542143759088224e-05", 10, MPFR_RNDN); mpfr_init(r26698); mpfr_init(r26699); mpfr_init(r26700); mpfr_init(r26701); mpfr_init(r26702); mpfr_init(r26703); mpfr_init(r26704); mpfr_init(r26705); mpfr_init(r26706); mpfr_init(r26707); mpfr_init(r26708); mpfr_init(r26709); } double f_fm(double x) { ; ; mpfr_set_d(r26685, x, MPFR_RNDN); mpfr_div(r26686, r26684, r26685, MPFR_RNDN); mpfr_add(r26687, r26683, r26686, MPFR_RNDN); mpfr_neg(r26688, r26683, MPFR_RNDN); mpfr_mul(r26689, r26685, r26685, MPFR_RNDN); mpfr_div(r26690, r26688, r26689, MPFR_RNDN); mpfr_neg(r26691, r26686, MPFR_RNDN); mpfr_fma(r26692, r26687, r26690, r26691, MPFR_RNDN); ; mpfr_set_si(r26694, mpfr_cmp(r26692, r26693) <= 0, MPFR_RNDN); mpfr_add(r26695, r26684, r26685, MPFR_RNDN); mpfr_fma(r26696, r26685, r26695, r26683, MPFR_RNDN); ; mpfr_set_si(r26698, mpfr_cmp(r26692, r26697) <= 0, MPFR_RNDN); mpfr_add(r26699, r26685, r26683, MPFR_RNDN); mpfr_div(r26700, r26685, r26699, MPFR_RNDN); mpfr_mul(r26701, r26700, r26700, MPFR_RNDN); mpfr_sub(r26702, r26685, r26683, MPFR_RNDN); mpfr_div(r26703, r26699, r26702, MPFR_RNDN); mpfr_mul(r26704, r26703, r26703, MPFR_RNDN); mpfr_sub(r26705, r26701, r26704, MPFR_RNDN); mpfr_add(r26706, r26700, r26703, MPFR_RNDN); mpfr_div(r26707, r26705, r26706, MPFR_RNDN); if (mpfr_get_si(r26698, MPFR_RNDN)) { mpfr_set(r26708, r26692, MPFR_RNDN); } else { mpfr_set(r26708, r26707, MPFR_RNDN); }; if (mpfr_get_si(r26694, MPFR_RNDN)) { mpfr_set(r26709, r26696, MPFR_RNDN); } else { mpfr_set(r26709, r26708, MPFR_RNDN); }; return mpfr_get_d(r26709, MPFR_RNDN); } static mpfr_t r26710, r26711, r26712, r26713, r26714, r26715, r26716, r26717, r26718, r26719, r26720, r26721, r26722, r26723, r26724, r26725, r26726, r26727, r26728, r26729, r26730, r26731, r26732, r26733, r26734, r26735, r26736; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r26710, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26711, "3", 10, MPFR_RNDN); mpfr_init(r26712); mpfr_init(r26713); mpfr_init(r26714); mpfr_init(r26715); mpfr_init(r26716); mpfr_init(r26717); mpfr_init(r26718); mpfr_init(r26719); mpfr_init_set_str(r26720, "-1.6839866920797528e-11", 10, MPFR_RNDN); mpfr_init(r26721); mpfr_init(r26722); mpfr_init(r26723); mpfr_init_set_str(r26724, "1.4542143759088224e-05", 10, MPFR_RNDN); mpfr_init(r26725); mpfr_init(r26726); mpfr_init(r26727); mpfr_init(r26728); mpfr_init(r26729); mpfr_init(r26730); mpfr_init(r26731); mpfr_init(r26732); mpfr_init(r26733); mpfr_init(r26734); mpfr_init(r26735); mpfr_init(r26736); } double f_dm(double x) { ; ; mpfr_set_d(r26712, x, MPFR_RNDN); mpfr_div(r26713, r26711, r26712, MPFR_RNDN); mpfr_add(r26714, r26710, r26713, MPFR_RNDN); mpfr_neg(r26715, r26710, MPFR_RNDN); mpfr_mul(r26716, r26712, r26712, MPFR_RNDN); mpfr_div(r26717, r26715, r26716, MPFR_RNDN); mpfr_neg(r26718, r26713, MPFR_RNDN); mpfr_fma(r26719, r26714, r26717, r26718, MPFR_RNDN); ; mpfr_set_si(r26721, mpfr_cmp(r26719, r26720) <= 0, MPFR_RNDN); mpfr_add(r26722, r26711, r26712, MPFR_RNDN); mpfr_fma(r26723, r26712, r26722, r26710, MPFR_RNDN); ; mpfr_set_si(r26725, mpfr_cmp(r26719, r26724) <= 0, MPFR_RNDN); mpfr_add(r26726, r26712, r26710, MPFR_RNDN); mpfr_div(r26727, r26712, r26726, MPFR_RNDN); mpfr_mul(r26728, r26727, r26727, MPFR_RNDN); mpfr_sub(r26729, r26712, r26710, MPFR_RNDN); mpfr_div(r26730, r26726, r26729, MPFR_RNDN); mpfr_mul(r26731, r26730, r26730, MPFR_RNDN); mpfr_sub(r26732, r26728, r26731, MPFR_RNDN); mpfr_add(r26733, r26727, r26730, MPFR_RNDN); mpfr_div(r26734, r26732, r26733, MPFR_RNDN); if (mpfr_get_si(r26725, MPFR_RNDN)) { mpfr_set(r26735, r26719, MPFR_RNDN); } else { mpfr_set(r26735, r26734, MPFR_RNDN); }; if (mpfr_get_si(r26721, MPFR_RNDN)) { mpfr_set(r26736, r26723, MPFR_RNDN); } else { mpfr_set(r26736, r26735, MPFR_RNDN); }; return mpfr_get_d(r26736, MPFR_RNDN); }