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