#include #include #include #include #include char *name = "Octave 3.8, jcobi/2"; double f_if(float alpha, float beta, float i) { float r25923 = alpha; float r25924 = beta; float r25925 = r25923 + r25924; float r25926 = r25924 - r25923; float r25927 = r25925 * r25926; float r25928 = 2; float r25929 = i; float r25930 = r25928 * r25929; float r25931 = r25925 + r25930; float r25932 = r25927 / r25931; float r25933 = 2.0; float r25934 = r25931 + r25933; float r25935 = r25932 / r25934; float r25936 = 1.0; float r25937 = r25935 + r25936; float r25938 = r25937 / r25933; return r25938; } double f_id(double alpha, double beta, double i) { double r25939 = alpha; double r25940 = beta; double r25941 = r25939 + r25940; double r25942 = r25940 - r25939; double r25943 = r25941 * r25942; double r25944 = 2; double r25945 = i; double r25946 = r25944 * r25945; double r25947 = r25941 + r25946; double r25948 = r25943 / r25947; double r25949 = 2.0; double r25950 = r25947 + r25949; double r25951 = r25948 / r25950; double r25952 = 1.0; double r25953 = r25951 + r25952; double r25954 = r25953 / r25949; return r25954; } double f_of(float alpha, float beta, float i) { float r25955 = 1; float r25956 = i; float r25957 = 2; float r25958 = beta; float r25959 = fma(r25956, r25957, r25958); float r25960 = 2.0; float r25961 = alpha; float r25962 = r25960 + r25961; float r25963 = r25959 + r25962; float r25964 = sqrt(r25963); float r25965 = r25955 / r25964; float r25966 = r25958 + r25961; float r25967 = r25966 / r25964; float r25968 = r25965 * r25967; float r25969 = r25958 - r25961; float r25970 = fma(r25957, r25956, r25966); float r25971 = r25969 / r25970; float r25972 = 1.0; float r25973 = fma(r25968, r25971, r25972); float r25974 = exp(r25973); float r25975 = log(r25974); float r25976 = r25975 / r25960; return r25976; } double f_od(double alpha, double beta, double i) { double r25977 = 1; double r25978 = i; double r25979 = 2; double r25980 = beta; double r25981 = fma(r25978, r25979, r25980); double r25982 = 2.0; double r25983 = alpha; double r25984 = r25982 + r25983; double r25985 = r25981 + r25984; double r25986 = sqrt(r25985); double r25987 = r25977 / r25986; double r25988 = r25980 + r25983; double r25989 = r25988 / r25986; double r25990 = r25987 * r25989; double r25991 = r25980 - r25983; double r25992 = fma(r25979, r25978, r25988); double r25993 = r25991 / r25992; double r25994 = 1.0; double r25995 = fma(r25990, r25993, r25994); double r25996 = exp(r25995); double r25997 = log(r25996); double r25998 = r25997 / r25982; return r25998; } 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 r25999, r26000, r26001, r26002, r26003, r26004, r26005, r26006, r26007, r26008, r26009, r26010, r26011, r26012, r26013, r26014; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r25999); mpfr_init(r26000); mpfr_init(r26001); mpfr_init(r26002); mpfr_init(r26003); mpfr_init_set_str(r26004, "2", 10, MPFR_RNDN); mpfr_init(r26005); mpfr_init(r26006); mpfr_init(r26007); mpfr_init(r26008); mpfr_init_set_str(r26009, "2.0", 10, MPFR_RNDN); mpfr_init(r26010); mpfr_init(r26011); mpfr_init_set_str(r26012, "1.0", 10, MPFR_RNDN); mpfr_init(r26013); mpfr_init(r26014); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r25999, alpha, MPFR_RNDN); mpfr_set_d(r26000, beta, MPFR_RNDN); mpfr_add(r26001, r25999, r26000, MPFR_RNDN); mpfr_sub(r26002, r26000, r25999, MPFR_RNDN); mpfr_mul(r26003, r26001, r26002, MPFR_RNDN); ; mpfr_set_d(r26005, i, MPFR_RNDN); mpfr_mul(r26006, r26004, r26005, MPFR_RNDN); mpfr_add(r26007, r26001, r26006, MPFR_RNDN); mpfr_div(r26008, r26003, r26007, MPFR_RNDN); ; mpfr_add(r26010, r26007, r26009, MPFR_RNDN); mpfr_div(r26011, r26008, r26010, MPFR_RNDN); ; mpfr_add(r26013, r26011, r26012, MPFR_RNDN); mpfr_div(r26014, r26013, r26009, MPFR_RNDN); return mpfr_get_d(r26014, MPFR_RNDN); } static mpfr_t r26015, r26016, r26017, r26018, r26019, r26020, r26021, r26022, r26023, r26024, r26025, r26026, r26027, r26028, r26029, r26030, r26031, r26032, r26033, r26034, r26035, r26036; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r26015, "1", 10, MPFR_RNDN); mpfr_init(r26016); mpfr_init_set_str(r26017, "2", 10, MPFR_RNDN); mpfr_init(r26018); mpfr_init(r26019); mpfr_init_set_str(r26020, "2.0", 10, MPFR_RNDN); mpfr_init(r26021); mpfr_init(r26022); mpfr_init(r26023); mpfr_init(r26024); mpfr_init(r26025); mpfr_init(r26026); mpfr_init(r26027); mpfr_init(r26028); mpfr_init(r26029); mpfr_init(r26030); mpfr_init(r26031); mpfr_init_set_str(r26032, "1.0", 10, MPFR_RNDN); mpfr_init(r26033); mpfr_init(r26034); mpfr_init(r26035); mpfr_init(r26036); } double f_fm(double alpha, double beta, double i) { ; mpfr_set_d(r26016, i, MPFR_RNDN); ; mpfr_set_d(r26018, beta, MPFR_RNDN); mpfr_fma(r26019, r26016, r26017, r26018, MPFR_RNDN); ; mpfr_set_d(r26021, alpha, MPFR_RNDN); mpfr_add(r26022, r26020, r26021, MPFR_RNDN); mpfr_add(r26023, r26019, r26022, MPFR_RNDN); mpfr_sqrt(r26024, r26023, MPFR_RNDN); mpfr_div(r26025, r26015, r26024, MPFR_RNDN); mpfr_add(r26026, r26018, r26021, MPFR_RNDN); mpfr_div(r26027, r26026, r26024, MPFR_RNDN); mpfr_mul(r26028, r26025, r26027, MPFR_RNDN); mpfr_sub(r26029, r26018, r26021, MPFR_RNDN); mpfr_fma(r26030, r26017, r26016, r26026, MPFR_RNDN); mpfr_div(r26031, r26029, r26030, MPFR_RNDN); ; mpfr_fma(r26033, r26028, r26031, r26032, MPFR_RNDN); mpfr_exp(r26034, r26033, MPFR_RNDN); mpfr_log(r26035, r26034, MPFR_RNDN); mpfr_div(r26036, r26035, r26020, MPFR_RNDN); return mpfr_get_d(r26036, MPFR_RNDN); } static mpfr_t r26037, r26038, r26039, r26040, r26041, r26042, r26043, r26044, r26045, r26046, r26047, r26048, r26049, r26050, r26051, r26052, r26053, r26054, r26055, r26056, r26057, r26058; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r26037, "1", 10, MPFR_RNDN); mpfr_init(r26038); mpfr_init_set_str(r26039, "2", 10, MPFR_RNDN); mpfr_init(r26040); mpfr_init(r26041); mpfr_init_set_str(r26042, "2.0", 10, MPFR_RNDN); mpfr_init(r26043); mpfr_init(r26044); mpfr_init(r26045); mpfr_init(r26046); mpfr_init(r26047); mpfr_init(r26048); mpfr_init(r26049); mpfr_init(r26050); mpfr_init(r26051); mpfr_init(r26052); mpfr_init(r26053); mpfr_init_set_str(r26054, "1.0", 10, MPFR_RNDN); mpfr_init(r26055); mpfr_init(r26056); mpfr_init(r26057); mpfr_init(r26058); } double f_dm(double alpha, double beta, double i) { ; mpfr_set_d(r26038, i, MPFR_RNDN); ; mpfr_set_d(r26040, beta, MPFR_RNDN); mpfr_fma(r26041, r26038, r26039, r26040, MPFR_RNDN); ; mpfr_set_d(r26043, alpha, MPFR_RNDN); mpfr_add(r26044, r26042, r26043, MPFR_RNDN); mpfr_add(r26045, r26041, r26044, MPFR_RNDN); mpfr_sqrt(r26046, r26045, MPFR_RNDN); mpfr_div(r26047, r26037, r26046, MPFR_RNDN); mpfr_add(r26048, r26040, r26043, MPFR_RNDN); mpfr_div(r26049, r26048, r26046, MPFR_RNDN); mpfr_mul(r26050, r26047, r26049, MPFR_RNDN); mpfr_sub(r26051, r26040, r26043, MPFR_RNDN); mpfr_fma(r26052, r26039, r26038, r26048, MPFR_RNDN); mpfr_div(r26053, r26051, r26052, MPFR_RNDN); ; mpfr_fma(r26055, r26050, r26053, r26054, MPFR_RNDN); mpfr_exp(r26056, r26055, MPFR_RNDN); mpfr_log(r26057, r26056, MPFR_RNDN); mpfr_div(r26058, r26057, r26042, MPFR_RNDN); return mpfr_get_d(r26058, MPFR_RNDN); }