#include #include #include #include #include char *name = "Hyperbolic arcsine"; double f_if(float x) { float r19722 = x; float r19723 = r19722 * r19722; float r19724 = 1.0f; float r19725 = r19723 + r19724; float r19726 = sqrt(r19725); float r19727 = r19722 + r19726; float r19728 = log(r19727); return r19728; } double f_id(double x) { double r19729 = x; double r19730 = r19729 * r19729; double r19731 = 1.0; double r19732 = r19730 + r19731; double r19733 = sqrt(r19732); double r19734 = r19729 + r19733; double r19735 = log(r19734); return r19735; } double f_of(float x) { float r19736 = x; float r19737 = -0.11215751184115397f; bool r19738 = r19736 <= r19737; float r19739 = 0.125f; float r19740 = r19736 * (r19736 * r19736); float r19741 = r19739 / r19740; float r19742 = 0.5f; float r19743 = r19742 / r19736; float r19744 = r19741 - r19743; float r19745 = log(r19744); float r19746 = 7.988163388129748f; bool r19747 = r19736 <= r19746; float r19748 = 0.075f; float r19749 = 5.0f; float r19750 = pow(r19736, r19749); float r19751 = r19748 * r19750; float r19752 = r19751 + r19736; float r19753 = 0.16666666666666666f; float r19754 = 3.0f; float r19755 = pow(r19736, r19754); float r19756 = r19753 * r19755; float r19757 = r19752 - r19756; float r19758 = r19736 + r19736; float r19759 = r19758 + r19743; float r19760 = r19759 - r19741; float r19761 = log(r19760); float r19762 = r19747 ? r19757 : r19761; float r19763 = r19738 ? r19745 : r19762; return r19763; } double f_od(double x) { double r19764 = x; double r19765 = -0.11215751184115397; bool r19766 = r19764 <= r19765; double r19767 = 0.125; double r19768 = r19764 * (r19764 * r19764); double r19769 = r19767 / r19768; double r19770 = 0.5; double r19771 = r19770 / r19764; double r19772 = r19769 - r19771; double r19773 = log(r19772); double r19774 = 7.988163388129748; bool r19775 = r19764 <= r19774; double r19776 = 0.075; double r19777 = 5.0; double r19778 = pow(r19764, r19777); double r19779 = r19776 * r19778; double r19780 = r19779 + r19764; double r19781 = 0.16666666666666666; double r19782 = 3.0; double r19783 = pow(r19764, r19782); double r19784 = r19781 * r19783; double r19785 = r19780 - r19784; double r19786 = r19764 + r19764; double r19787 = r19786 + r19771; double r19788 = r19787 - r19769; double r19789 = log(r19788); double r19790 = r19775 ? r19785 : r19789; double r19791 = r19766 ? r19773 : r19790; return r19791; } 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 r19792, r19793, r19794, r19795, r19796, r19797, r19798; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init(r19792); mpfr_init(r19793); mpfr_init_set_str(r19794, "1", 10, MPFR_RNDN); mpfr_init(r19795); mpfr_init(r19796); mpfr_init(r19797); mpfr_init(r19798); } double f_im(double x) { mpfr_set_d(r19792, x, MPFR_RNDN); mpfr_sqr(r19793, r19792, MPFR_RNDN); ; mpfr_add(r19795, r19793, r19794, MPFR_RNDN); mpfr_sqrt(r19796, r19795, MPFR_RNDN); mpfr_add(r19797, r19792, r19796, MPFR_RNDN); mpfr_log(r19798, r19797, MPFR_RNDN); return mpfr_get_d(r19798, MPFR_RNDN); } static mpfr_t r19799, r19800, r19801, r19802, r19803, r19804, r19805, r19806, r19807, r19808, r19809, r19810, r19811, r19812, r19813, r19814, r19815, r19816, r19817, r19818, r19819, r19820, r19821, r19822, r19823, r19824, r19825, r19826; void setup_mpfr_f_fm() { mpfr_set_default_prec(2448); mpfr_init(r19799); mpfr_init_set_str(r19800, "-0.11215751184115397", 10, MPFR_RNDN); mpfr_init(r19801); mpfr_init_set_str(r19802, "1/8", 10, MPFR_RNDN); mpfr_init(r19803); mpfr_init(r19804); mpfr_init_set_str(r19805, "1/2", 10, MPFR_RNDN); mpfr_init(r19806); mpfr_init(r19807); mpfr_init(r19808); mpfr_init_set_str(r19809, "7.988163388129748", 10, MPFR_RNDN); mpfr_init(r19810); mpfr_init_set_str(r19811, "3/40", 10, MPFR_RNDN); mpfr_init_set_str(r19812, "5", 10, MPFR_RNDN); mpfr_init(r19813); mpfr_init(r19814); mpfr_init(r19815); mpfr_init_set_str(r19816, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r19817, "3", 10, MPFR_RNDN); mpfr_init(r19818); mpfr_init(r19819); mpfr_init(r19820); mpfr_init(r19821); mpfr_init(r19822); mpfr_init(r19823); mpfr_init(r19824); mpfr_init(r19825); mpfr_init(r19826); } double f_fm(double x) { mpfr_set_d(r19799, x, MPFR_RNDN); ; mpfr_set_si(r19801, mpfr_cmp(r19799, r19800) <= 0, MPFR_RNDN); ; mpfr_mul(r19803, r19799, r19799, MPFR_RNDN); mpfr_mul(r19803, r19803, r19799, MPFR_RNDN); mpfr_div(r19804, r19802, r19803, MPFR_RNDN); ; mpfr_div(r19806, r19805, r19799, MPFR_RNDN); mpfr_sub(r19807, r19804, r19806, MPFR_RNDN); mpfr_log(r19808, r19807, MPFR_RNDN); ; mpfr_set_si(r19810, mpfr_cmp(r19799, r19809) <= 0, MPFR_RNDN); ; ; mpfr_pow(r19813, r19799, r19812, MPFR_RNDN); mpfr_mul(r19814, r19811, r19813, MPFR_RNDN); mpfr_add(r19815, r19814, r19799, MPFR_RNDN); ; ; mpfr_pow(r19818, r19799, r19817, MPFR_RNDN); mpfr_mul(r19819, r19816, r19818, MPFR_RNDN); mpfr_sub(r19820, r19815, r19819, MPFR_RNDN); mpfr_add(r19821, r19799, r19799, MPFR_RNDN); mpfr_add(r19822, r19821, r19806, MPFR_RNDN); mpfr_sub(r19823, r19822, r19804, MPFR_RNDN); mpfr_log(r19824, r19823, MPFR_RNDN); if (mpfr_get_si(r19810, MPFR_RNDN)) { mpfr_set(r19825, r19820, MPFR_RNDN); } else { mpfr_set(r19825, r19824, MPFR_RNDN); }; if (mpfr_get_si(r19801, MPFR_RNDN)) { mpfr_set(r19826, r19808, MPFR_RNDN); } else { mpfr_set(r19826, r19825, MPFR_RNDN); }; return mpfr_get_d(r19826, MPFR_RNDN); } static mpfr_t r19827, r19828, r19829, r19830, r19831, r19832, r19833, r19834, r19835, r19836, r19837, r19838, r19839, r19840, r19841, r19842, r19843, r19844, r19845, r19846, r19847, r19848, r19849, r19850, r19851, r19852, r19853, r19854; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init(r19827); mpfr_init_set_str(r19828, "-0.11215751184115397", 10, MPFR_RNDN); mpfr_init(r19829); mpfr_init_set_str(r19830, "1/8", 10, MPFR_RNDN); mpfr_init(r19831); mpfr_init(r19832); mpfr_init_set_str(r19833, "1/2", 10, MPFR_RNDN); mpfr_init(r19834); mpfr_init(r19835); mpfr_init(r19836); mpfr_init_set_str(r19837, "7.988163388129748", 10, MPFR_RNDN); mpfr_init(r19838); mpfr_init_set_str(r19839, "3/40", 10, MPFR_RNDN); mpfr_init_set_str(r19840, "5", 10, MPFR_RNDN); mpfr_init(r19841); mpfr_init(r19842); mpfr_init(r19843); mpfr_init_set_str(r19844, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r19845, "3", 10, MPFR_RNDN); mpfr_init(r19846); mpfr_init(r19847); mpfr_init(r19848); mpfr_init(r19849); mpfr_init(r19850); mpfr_init(r19851); mpfr_init(r19852); mpfr_init(r19853); mpfr_init(r19854); } double f_dm(double x) { mpfr_set_d(r19827, x, MPFR_RNDN); ; mpfr_set_si(r19829, mpfr_cmp(r19827, r19828) <= 0, MPFR_RNDN); ; mpfr_mul(r19831, r19827, r19827, MPFR_RNDN); mpfr_mul(r19831, r19831, r19827, MPFR_RNDN); mpfr_div(r19832, r19830, r19831, MPFR_RNDN); ; mpfr_div(r19834, r19833, r19827, MPFR_RNDN); mpfr_sub(r19835, r19832, r19834, MPFR_RNDN); mpfr_log(r19836, r19835, MPFR_RNDN); ; mpfr_set_si(r19838, mpfr_cmp(r19827, r19837) <= 0, MPFR_RNDN); ; ; mpfr_pow(r19841, r19827, r19840, MPFR_RNDN); mpfr_mul(r19842, r19839, r19841, MPFR_RNDN); mpfr_add(r19843, r19842, r19827, MPFR_RNDN); ; ; mpfr_pow(r19846, r19827, r19845, MPFR_RNDN); mpfr_mul(r19847, r19844, r19846, MPFR_RNDN); mpfr_sub(r19848, r19843, r19847, MPFR_RNDN); mpfr_add(r19849, r19827, r19827, MPFR_RNDN); mpfr_add(r19850, r19849, r19834, MPFR_RNDN); mpfr_sub(r19851, r19850, r19832, MPFR_RNDN); mpfr_log(r19852, r19851, MPFR_RNDN); if (mpfr_get_si(r19838, MPFR_RNDN)) { mpfr_set(r19853, r19848, MPFR_RNDN); } else { mpfr_set(r19853, r19852, MPFR_RNDN); }; if (mpfr_get_si(r19829, MPFR_RNDN)) { mpfr_set(r19854, r19836, MPFR_RNDN); } else { mpfr_set(r19854, r19853, MPFR_RNDN); }; return mpfr_get_d(r19854, MPFR_RNDN); }