#include #include #include #include #include char *name = "Asymptote C"; double f_if(float x) { float r26641 = x; float r26642 = 1; float r26643 = r26641 + r26642; float r26644 = r26641 / r26643; float r26645 = r26641 - r26642; float r26646 = r26643 / r26645; float r26647 = r26644 - r26646; return r26647; } double f_id(double x) { double r26648 = x; double r26649 = 1; double r26650 = r26648 + r26649; double r26651 = r26648 / r26650; double r26652 = r26648 - r26649; double r26653 = r26650 / r26652; double r26654 = r26651 - r26653; return r26654; } double f_of(float x) { float r26655 = x; float r26656 = 1; float r26657 = r26656 + r26655; float r26658 = r26655 / r26657; float r26659 = r26655 - r26656; float r26660 = r26657 / r26659; float r26661 = r26658 - r26660; float r26662 = 1.229801583672145e-07; bool r26663 = r26661 <= r26662; float r26664 = 3; float r26665 = r26664 / r26655; float r26666 = r26665 + r26656; float r26667 = -r26656; float r26668 = r26655 * r26655; float r26669 = r26667 / r26668; float r26670 = -r26664; float r26671 = r26670 / r26655; float r26672 = fma(r26666, r26669, r26671); float r26673 = -r26657; float r26674 = fma(r26668, r26655, r26667); float r26675 = r26673 / r26674; float r26676 = fma(r26657, r26655, r26656); float r26677 = fma(r26675, r26676, r26658); float r26678 = r26663 ? r26672 : r26677; return r26678; } double f_od(double x) { double r26679 = x; double r26680 = 1; double r26681 = r26680 + r26679; double r26682 = r26679 / r26681; double r26683 = r26679 - r26680; double r26684 = r26681 / r26683; double r26685 = r26682 - r26684; double r26686 = 1.229801583672145e-07; bool r26687 = r26685 <= r26686; double r26688 = 3; double r26689 = r26688 / r26679; double r26690 = r26689 + r26680; double r26691 = -r26680; double r26692 = r26679 * r26679; double r26693 = r26691 / r26692; double r26694 = -r26688; double r26695 = r26694 / r26679; double r26696 = fma(r26690, r26693, r26695); double r26697 = -r26681; double r26698 = fma(r26692, r26679, r26691); double r26699 = r26697 / r26698; double r26700 = fma(r26681, r26679, r26680); double r26701 = fma(r26699, r26700, r26682); double r26702 = r26687 ? r26696 : r26701; return r26702; } 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 r26703, r26704, r26705, r26706, r26707, r26708, r26709; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r26703); mpfr_init_set_str(r26704, "1", 10, MPFR_RNDN); mpfr_init(r26705); mpfr_init(r26706); mpfr_init(r26707); mpfr_init(r26708); mpfr_init(r26709); } double f_im(double x) { mpfr_set_d(r26703, x, MPFR_RNDN); ; mpfr_add(r26705, r26703, r26704, MPFR_RNDN); mpfr_div(r26706, r26703, r26705, MPFR_RNDN); mpfr_sub(r26707, r26703, r26704, MPFR_RNDN); mpfr_div(r26708, r26705, r26707, MPFR_RNDN); mpfr_sub(r26709, r26706, 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r26710); mpfr_init_set_str(r26711, "1", 10, MPFR_RNDN); mpfr_init(r26712); mpfr_init(r26713); mpfr_init(r26714); mpfr_init(r26715); mpfr_init(r26716); mpfr_init_set_str(r26717, "1.229801583672145e-07", 10, MPFR_RNDN); mpfr_init(r26718); mpfr_init_set_str(r26719, "3", 10, MPFR_RNDN); mpfr_init(r26720); mpfr_init(r26721); mpfr_init(r26722); mpfr_init(r26723); mpfr_init(r26724); 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); } double f_fm(double x) { mpfr_set_d(r26710, x, MPFR_RNDN); ; mpfr_add(r26712, r26711, r26710, MPFR_RNDN); mpfr_div(r26713, r26710, r26712, MPFR_RNDN); mpfr_sub(r26714, r26710, r26711, MPFR_RNDN); mpfr_div(r26715, r26712, r26714, MPFR_RNDN); mpfr_sub(r26716, r26713, r26715, MPFR_RNDN); ; mpfr_set_si(r26718, mpfr_cmp(r26716, r26717) <= 0, MPFR_RNDN); ; mpfr_div(r26720, r26719, r26710, MPFR_RNDN); mpfr_add(r26721, r26720, r26711, MPFR_RNDN); mpfr_neg(r26722, r26711, MPFR_RNDN); mpfr_mul(r26723, r26710, r26710, MPFR_RNDN); mpfr_div(r26724, r26722, r26723, MPFR_RNDN); mpfr_neg(r26725, r26719, MPFR_RNDN); mpfr_div(r26726, r26725, r26710, MPFR_RNDN); mpfr_fma(r26727, r26721, r26724, r26726, MPFR_RNDN); mpfr_neg(r26728, r26712, MPFR_RNDN); mpfr_fma(r26729, r26723, r26710, r26722, MPFR_RNDN); mpfr_div(r26730, r26728, r26729, MPFR_RNDN); mpfr_fma(r26731, r26712, r26710, r26711, MPFR_RNDN); mpfr_fma(r26732, r26730, r26731, r26713, MPFR_RNDN); if (mpfr_get_si(r26718, MPFR_RNDN)) { mpfr_set(r26733, r26727, MPFR_RNDN); } else { mpfr_set(r26733, r26732, MPFR_RNDN); }; return mpfr_get_d(r26733, MPFR_RNDN); } static mpfr_t r26734, r26735, r26736, r26737, r26738, r26739, r26740, r26741, r26742, r26743, r26744, r26745, r26746, r26747, r26748, r26749, r26750, r26751, r26752, r26753, r26754, r26755, r26756, r26757; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r26734); mpfr_init_set_str(r26735, "1", 10, MPFR_RNDN); mpfr_init(r26736); mpfr_init(r26737); mpfr_init(r26738); mpfr_init(r26739); mpfr_init(r26740); mpfr_init_set_str(r26741, "1.229801583672145e-07", 10, MPFR_RNDN); mpfr_init(r26742); mpfr_init_set_str(r26743, "3", 10, MPFR_RNDN); mpfr_init(r26744); mpfr_init(r26745); mpfr_init(r26746); mpfr_init(r26747); mpfr_init(r26748); mpfr_init(r26749); mpfr_init(r26750); mpfr_init(r26751); mpfr_init(r26752); mpfr_init(r26753); mpfr_init(r26754); mpfr_init(r26755); mpfr_init(r26756); mpfr_init(r26757); } double f_dm(double x) { mpfr_set_d(r26734, x, MPFR_RNDN); ; mpfr_add(r26736, r26735, r26734, MPFR_RNDN); mpfr_div(r26737, r26734, r26736, MPFR_RNDN); mpfr_sub(r26738, r26734, r26735, MPFR_RNDN); mpfr_div(r26739, r26736, r26738, MPFR_RNDN); mpfr_sub(r26740, r26737, r26739, MPFR_RNDN); ; mpfr_set_si(r26742, mpfr_cmp(r26740, r26741) <= 0, MPFR_RNDN); ; mpfr_div(r26744, r26743, r26734, MPFR_RNDN); mpfr_add(r26745, r26744, r26735, MPFR_RNDN); mpfr_neg(r26746, r26735, MPFR_RNDN); mpfr_mul(r26747, r26734, r26734, MPFR_RNDN); mpfr_div(r26748, r26746, r26747, MPFR_RNDN); mpfr_neg(r26749, r26743, MPFR_RNDN); mpfr_div(r26750, r26749, r26734, MPFR_RNDN); mpfr_fma(r26751, r26745, r26748, r26750, MPFR_RNDN); mpfr_neg(r26752, r26736, MPFR_RNDN); mpfr_fma(r26753, r26747, r26734, r26746, MPFR_RNDN); mpfr_div(r26754, r26752, r26753, MPFR_RNDN); mpfr_fma(r26755, r26736, r26734, r26735, MPFR_RNDN); mpfr_fma(r26756, r26754, r26755, r26737, MPFR_RNDN); if (mpfr_get_si(r26742, MPFR_RNDN)) { mpfr_set(r26757, r26751, MPFR_RNDN); } else { mpfr_set(r26757, r26756, MPFR_RNDN); }; return mpfr_get_d(r26757, MPFR_RNDN); }