#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r22645 = i; float r22646 = alpha; float r22647 = beta; float r22648 = r22646 + r22647; float r22649 = r22648 + r22645; float r22650 = r22645 * r22649; float r22651 = r22647 * r22646; float r22652 = r22651 + r22650; float r22653 = r22650 * r22652; float r22654 = 2; float r22655 = r22654 * r22645; float r22656 = r22648 + r22655; float r22657 = r22656 * r22656; float r22658 = r22653 / r22657; float r22659 = 1.0; float r22660 = r22657 - r22659; float r22661 = r22658 / r22660; return r22661; } double f_id(double alpha, double beta, double i) { double r22662 = i; double r22663 = alpha; double r22664 = beta; double r22665 = r22663 + r22664; double r22666 = r22665 + r22662; double r22667 = r22662 * r22666; double r22668 = r22664 * r22663; double r22669 = r22668 + r22667; double r22670 = r22667 * r22669; double r22671 = 2; double r22672 = r22671 * r22662; double r22673 = r22665 + r22672; double r22674 = r22673 * r22673; double r22675 = r22670 / r22674; double r22676 = 1.0; double r22677 = r22674 - r22676; double r22678 = r22675 / r22677; return r22678; } double f_of(float alpha, float beta, float i) { float r22679 = beta; float r22680 = 2.2950245454369153e+151; bool r22681 = r22679 <= r22680; float r22682 = alpha; float r22683 = i; float r22684 = r22682 + r22683; float r22685 = r22679 + r22684; float r22686 = 2; float r22687 = fma(r22683, r22686, r22682); float r22688 = r22687 + r22679; float r22689 = r22685 / r22688; float r22690 = r22679 * r22682; float r22691 = fma(r22685, r22683, r22690); float r22692 = r22689 * r22691; float r22693 = r22683 / r22688; float r22694 = r22692 * r22693; float r22695 = 1.0; float r22696 = -r22695; float r22697 = fma(r22688, r22688, r22696); float r22698 = r22694 / r22697; float r22699 = sqrt(r22698); float r22700 = r22699 * r22699; float r22701 = +inf.0; bool r22702 = r22679 <= r22701; float r22703 = 0; float r22704 = r22702 ? r22703 : r22703; float r22705 = r22681 ? r22700 : r22704; return r22705; } double f_od(double alpha, double beta, double i) { double r22706 = beta; double r22707 = 2.2950245454369153e+151; bool r22708 = r22706 <= r22707; double r22709 = alpha; double r22710 = i; double r22711 = r22709 + r22710; double r22712 = r22706 + r22711; double r22713 = 2; double r22714 = fma(r22710, r22713, r22709); double r22715 = r22714 + r22706; double r22716 = r22712 / r22715; double r22717 = r22706 * r22709; double r22718 = fma(r22712, r22710, r22717); double r22719 = r22716 * r22718; double r22720 = r22710 / r22715; double r22721 = r22719 * r22720; double r22722 = 1.0; double r22723 = -r22722; double r22724 = fma(r22715, r22715, r22723); double r22725 = r22721 / r22724; double r22726 = sqrt(r22725); double r22727 = r22726 * r22726; double r22728 = +inf.0; bool r22729 = r22706 <= r22728; double r22730 = 0; double r22731 = r22729 ? r22730 : r22730; double r22732 = r22708 ? r22727 : r22731; return r22732; } 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 r22733, r22734, r22735, r22736, r22737, r22738, r22739, r22740, r22741, r22742, r22743, r22744, r22745, r22746, r22747, r22748, r22749; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22733); mpfr_init(r22734); mpfr_init(r22735); mpfr_init(r22736); mpfr_init(r22737); mpfr_init(r22738); mpfr_init(r22739); mpfr_init(r22740); mpfr_init(r22741); mpfr_init_set_str(r22742, "2", 10, MPFR_RNDN); mpfr_init(r22743); mpfr_init(r22744); mpfr_init(r22745); mpfr_init(r22746); mpfr_init_set_str(r22747, "1.0", 10, MPFR_RNDN); mpfr_init(r22748); mpfr_init(r22749); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r22733, i, MPFR_RNDN); mpfr_set_d(r22734, alpha, MPFR_RNDN); mpfr_set_d(r22735, beta, MPFR_RNDN); mpfr_add(r22736, r22734, r22735, MPFR_RNDN); mpfr_add(r22737, r22736, r22733, MPFR_RNDN); mpfr_mul(r22738, r22733, r22737, MPFR_RNDN); mpfr_mul(r22739, r22735, r22734, MPFR_RNDN); mpfr_add(r22740, r22739, r22738, MPFR_RNDN); mpfr_mul(r22741, r22738, r22740, MPFR_RNDN); ; mpfr_mul(r22743, r22742, r22733, MPFR_RNDN); mpfr_add(r22744, r22736, r22743, MPFR_RNDN); mpfr_mul(r22745, r22744, r22744, MPFR_RNDN); mpfr_div(r22746, r22741, r22745, MPFR_RNDN); ; mpfr_sub(r22748, r22745, r22747, MPFR_RNDN); mpfr_div(r22749, r22746, r22748, MPFR_RNDN); return mpfr_get_d(r22749, MPFR_RNDN); } static mpfr_t r22750, r22751, r22752, r22753, r22754, r22755, r22756, r22757, r22758, r22759, r22760, r22761, r22762, r22763, r22764, r22765, r22766, r22767, r22768, r22769, r22770, r22771, r22772, r22773, r22774, r22775, r22776; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22750); mpfr_init_set_str(r22751, "2.2950245454369153e+151", 10, MPFR_RNDN); mpfr_init(r22752); mpfr_init(r22753); mpfr_init(r22754); mpfr_init(r22755); mpfr_init(r22756); mpfr_init_set_str(r22757, "2", 10, MPFR_RNDN); mpfr_init(r22758); mpfr_init(r22759); mpfr_init(r22760); mpfr_init(r22761); mpfr_init(r22762); mpfr_init(r22763); mpfr_init(r22764); mpfr_init(r22765); mpfr_init_set_str(r22766, "1.0", 10, MPFR_RNDN); mpfr_init(r22767); mpfr_init(r22768); mpfr_init(r22769); mpfr_init(r22770); mpfr_init(r22771); mpfr_init_set_str(r22772, "+inf.0", 10, MPFR_RNDN); mpfr_init(r22773); mpfr_init_set_str(r22774, "0", 10, MPFR_RNDN); mpfr_init(r22775); mpfr_init(r22776); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r22750, beta, MPFR_RNDN); ; mpfr_set_si(r22752, mpfr_cmp(r22750, r22751) <= 0, MPFR_RNDN); mpfr_set_d(r22753, alpha, MPFR_RNDN); mpfr_set_d(r22754, i, MPFR_RNDN); mpfr_add(r22755, r22753, r22754, MPFR_RNDN); mpfr_add(r22756, r22750, r22755, MPFR_RNDN); ; mpfr_fma(r22758, r22754, r22757, r22753, MPFR_RNDN); mpfr_add(r22759, r22758, r22750, MPFR_RNDN); mpfr_div(r22760, r22756, r22759, MPFR_RNDN); mpfr_mul(r22761, r22750, r22753, MPFR_RNDN); mpfr_fma(r22762, r22756, r22754, r22761, MPFR_RNDN); mpfr_mul(r22763, r22760, r22762, MPFR_RNDN); mpfr_div(r22764, r22754, r22759, MPFR_RNDN); mpfr_mul(r22765, r22763, r22764, MPFR_RNDN); ; mpfr_neg(r22767, r22766, MPFR_RNDN); mpfr_fma(r22768, r22759, r22759, r22767, MPFR_RNDN); mpfr_div(r22769, r22765, r22768, MPFR_RNDN); mpfr_sqrt(r22770, r22769, MPFR_RNDN); mpfr_mul(r22771, r22770, r22770, MPFR_RNDN); ; mpfr_set_si(r22773, mpfr_cmp(r22750, r22772) <= 0, MPFR_RNDN); ; if (mpfr_get_si(r22773, MPFR_RNDN)) { mpfr_set(r22775, r22774, MPFR_RNDN); } else { mpfr_set(r22775, r22774, MPFR_RNDN); }; if (mpfr_get_si(r22752, MPFR_RNDN)) { mpfr_set(r22776, r22771, MPFR_RNDN); } else { mpfr_set(r22776, r22775, MPFR_RNDN); }; return mpfr_get_d(r22776, MPFR_RNDN); } static mpfr_t r22777, r22778, r22779, r22780, r22781, r22782, r22783, r22784, r22785, r22786, r22787, r22788, r22789, r22790, r22791, r22792, r22793, r22794, r22795, r22796, r22797, r22798, r22799, r22800, r22801, r22802, r22803; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22777); mpfr_init_set_str(r22778, "2.2950245454369153e+151", 10, MPFR_RNDN); mpfr_init(r22779); mpfr_init(r22780); mpfr_init(r22781); mpfr_init(r22782); mpfr_init(r22783); mpfr_init_set_str(r22784, "2", 10, MPFR_RNDN); mpfr_init(r22785); mpfr_init(r22786); mpfr_init(r22787); mpfr_init(r22788); mpfr_init(r22789); mpfr_init(r22790); mpfr_init(r22791); mpfr_init(r22792); mpfr_init_set_str(r22793, "1.0", 10, MPFR_RNDN); mpfr_init(r22794); mpfr_init(r22795); mpfr_init(r22796); mpfr_init(r22797); mpfr_init(r22798); mpfr_init_set_str(r22799, "+inf.0", 10, MPFR_RNDN); mpfr_init(r22800); mpfr_init_set_str(r22801, "0", 10, MPFR_RNDN); mpfr_init(r22802); mpfr_init(r22803); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r22777, beta, MPFR_RNDN); ; mpfr_set_si(r22779, mpfr_cmp(r22777, r22778) <= 0, MPFR_RNDN); mpfr_set_d(r22780, alpha, MPFR_RNDN); mpfr_set_d(r22781, i, MPFR_RNDN); mpfr_add(r22782, r22780, r22781, MPFR_RNDN); mpfr_add(r22783, r22777, r22782, MPFR_RNDN); ; mpfr_fma(r22785, r22781, r22784, r22780, MPFR_RNDN); mpfr_add(r22786, r22785, r22777, MPFR_RNDN); mpfr_div(r22787, r22783, r22786, MPFR_RNDN); mpfr_mul(r22788, r22777, r22780, MPFR_RNDN); mpfr_fma(r22789, r22783, r22781, r22788, MPFR_RNDN); mpfr_mul(r22790, r22787, r22789, MPFR_RNDN); mpfr_div(r22791, r22781, r22786, MPFR_RNDN); mpfr_mul(r22792, r22790, r22791, MPFR_RNDN); ; mpfr_neg(r22794, r22793, MPFR_RNDN); mpfr_fma(r22795, r22786, r22786, r22794, MPFR_RNDN); mpfr_div(r22796, r22792, r22795, MPFR_RNDN); mpfr_sqrt(r22797, r22796, MPFR_RNDN); mpfr_mul(r22798, r22797, r22797, MPFR_RNDN); ; mpfr_set_si(r22800, mpfr_cmp(r22777, r22799) <= 0, MPFR_RNDN); ; if (mpfr_get_si(r22800, MPFR_RNDN)) { mpfr_set(r22802, r22801, MPFR_RNDN); } else { mpfr_set(r22802, r22801, MPFR_RNDN); }; if (mpfr_get_si(r22779, MPFR_RNDN)) { mpfr_set(r22803, r22798, MPFR_RNDN); } else { mpfr_set(r22803, r22802, MPFR_RNDN); }; return mpfr_get_d(r22803, MPFR_RNDN); }