#include #include #include #include #include char *name = "Octave 3.8, jcobi/4, as called"; double f_if(float i) { float r25868 = i; float r25869 = r25868 * r25868; float r25870 = r25869 * r25869; float r25871 = 2; float r25872 = r25871 * r25868; float r25873 = r25872 * r25872; float r25874 = r25870 / r25873; float r25875 = 1.0; float r25876 = r25873 - r25875; float r25877 = r25874 / r25876; return r25877; } double f_id(double i) { double r25878 = i; double r25879 = r25878 * r25878; double r25880 = r25879 * r25879; double r25881 = 2; double r25882 = r25881 * r25878; double r25883 = r25882 * r25882; double r25884 = r25880 / r25883; double r25885 = 1.0; double r25886 = r25883 - r25885; double r25887 = r25884 / r25886; return r25887; } double f_of(float i) { float r25888 = i; float r25889 = 225.58464822734692; bool r25890 = r25888 <= r25889; float r25891 = 1; float r25892 = 2; float r25893 = r25891 / r25892; float r25894 = r25893 * r25888; float r25895 = r25894 * r25894; float r25896 = r25892 * r25888; float r25897 = r25896 * r25896; float r25898 = 1.0; float r25899 = r25897 - r25898; float r25900 = r25895 / r25899; float r25901 = 0.015625; float r25902 = r25901 / r25888; float r25903 = r25902 / r25888; float r25904 = 0.00390625; float r25905 = 4; float r25906 = pow(r25888, r25905); float r25907 = r25904 / r25906; float r25908 = 1/16; float r25909 = r25907 + r25908; float r25910 = r25903 + r25909; float r25911 = r25890 ? r25900 : r25910; return r25911; } double f_od(double i) { double r25912 = i; double r25913 = 225.58464822734692; bool r25914 = r25912 <= r25913; double r25915 = 1; double r25916 = 2; double r25917 = r25915 / r25916; double r25918 = r25917 * r25912; double r25919 = r25918 * r25918; double r25920 = r25916 * r25912; double r25921 = r25920 * r25920; double r25922 = 1.0; double r25923 = r25921 - r25922; double r25924 = r25919 / r25923; double r25925 = 0.015625; double r25926 = r25925 / r25912; double r25927 = r25926 / r25912; double r25928 = 0.00390625; double r25929 = 4; double r25930 = pow(r25912, r25929); double r25931 = r25928 / r25930; double r25932 = 1/16; double r25933 = r25931 + r25932; double r25934 = r25927 + r25933; double r25935 = r25914 ? r25924 : r25934; return r25935; } 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 r25936, r25937, r25938, r25939, r25940, r25941, r25942, r25943, r25944, r25945; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r25936); mpfr_init(r25937); mpfr_init(r25938); mpfr_init_set_str(r25939, "2", 10, MPFR_RNDN); mpfr_init(r25940); mpfr_init(r25941); mpfr_init(r25942); mpfr_init_set_str(r25943, "1.0", 10, MPFR_RNDN); mpfr_init(r25944); mpfr_init(r25945); } double f_im(double i) { mpfr_set_d(r25936, i, MPFR_RNDN); mpfr_mul(r25937, r25936, r25936, MPFR_RNDN); mpfr_mul(r25938, r25937, r25937, MPFR_RNDN); ; mpfr_mul(r25940, r25939, r25936, MPFR_RNDN); mpfr_mul(r25941, r25940, r25940, MPFR_RNDN); mpfr_div(r25942, r25938, r25941, MPFR_RNDN); ; mpfr_sub(r25944, r25941, r25943, MPFR_RNDN); mpfr_div(r25945, r25942, r25944, MPFR_RNDN); return mpfr_get_d(r25945, MPFR_RNDN); } static mpfr_t r25946, r25947, r25948, r25949, r25950, r25951, r25952, r25953, r25954, r25955, r25956, r25957, r25958, r25959, r25960, r25961, r25962, r25963, r25964, r25965, r25966, r25967, r25968, r25969; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r25946); mpfr_init_set_str(r25947, "225.58464822734692", 10, MPFR_RNDN); mpfr_init(r25948); mpfr_init_set_str(r25949, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25950, "2", 10, MPFR_RNDN); mpfr_init(r25951); mpfr_init(r25952); mpfr_init(r25953); mpfr_init(r25954); mpfr_init(r25955); mpfr_init_set_str(r25956, "1.0", 10, MPFR_RNDN); mpfr_init(r25957); mpfr_init(r25958); mpfr_init_set_str(r25959, "0.015625", 10, MPFR_RNDN); mpfr_init(r25960); mpfr_init(r25961); mpfr_init_set_str(r25962, "0.00390625", 10, MPFR_RNDN); mpfr_init_set_str(r25963, "4", 10, MPFR_RNDN); mpfr_init(r25964); mpfr_init(r25965); mpfr_init_set_str(r25966, "1/16", 10, MPFR_RNDN); mpfr_init(r25967); mpfr_init(r25968); mpfr_init(r25969); } double f_fm(double i) { mpfr_set_d(r25946, i, MPFR_RNDN); ; mpfr_set_si(r25948, mpfr_cmp(r25946, r25947) <= 0, MPFR_RNDN); ; ; mpfr_div(r25951, r25949, r25950, MPFR_RNDN); mpfr_mul(r25952, r25951, r25946, MPFR_RNDN); mpfr_mul(r25953, r25952, r25952, MPFR_RNDN); mpfr_mul(r25954, r25950, r25946, MPFR_RNDN); mpfr_mul(r25955, r25954, r25954, MPFR_RNDN); ; mpfr_sub(r25957, r25955, r25956, MPFR_RNDN); mpfr_div(r25958, r25953, r25957, MPFR_RNDN); ; mpfr_div(r25960, r25959, r25946, MPFR_RNDN); mpfr_div(r25961, r25960, r25946, MPFR_RNDN); ; ; mpfr_pow(r25964, r25946, r25963, MPFR_RNDN); mpfr_div(r25965, r25962, r25964, MPFR_RNDN); ; mpfr_add(r25967, r25965, r25966, MPFR_RNDN); mpfr_add(r25968, r25961, r25967, MPFR_RNDN); if (mpfr_get_si(r25948, MPFR_RNDN)) { mpfr_set(r25969, r25958, MPFR_RNDN); } else { mpfr_set(r25969, r25968, MPFR_RNDN); }; return mpfr_get_d(r25969, MPFR_RNDN); } static mpfr_t r25970, r25971, r25972, r25973, r25974, r25975, r25976, r25977, r25978, r25979, r25980, r25981, r25982, r25983, r25984, r25985, r25986, r25987, r25988, r25989, r25990, r25991, r25992, r25993; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r25970); mpfr_init_set_str(r25971, "225.58464822734692", 10, MPFR_RNDN); mpfr_init(r25972); mpfr_init_set_str(r25973, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25974, "2", 10, MPFR_RNDN); mpfr_init(r25975); mpfr_init(r25976); mpfr_init(r25977); mpfr_init(r25978); mpfr_init(r25979); mpfr_init_set_str(r25980, "1.0", 10, MPFR_RNDN); mpfr_init(r25981); mpfr_init(r25982); mpfr_init_set_str(r25983, "0.015625", 10, MPFR_RNDN); mpfr_init(r25984); mpfr_init(r25985); mpfr_init_set_str(r25986, "0.00390625", 10, MPFR_RNDN); mpfr_init_set_str(r25987, "4", 10, MPFR_RNDN); mpfr_init(r25988); mpfr_init(r25989); mpfr_init_set_str(r25990, "1/16", 10, MPFR_RNDN); mpfr_init(r25991); mpfr_init(r25992); mpfr_init(r25993); } double f_dm(double i) { mpfr_set_d(r25970, i, MPFR_RNDN); ; mpfr_set_si(r25972, mpfr_cmp(r25970, r25971) <= 0, MPFR_RNDN); ; ; mpfr_div(r25975, r25973, r25974, MPFR_RNDN); mpfr_mul(r25976, r25975, r25970, MPFR_RNDN); mpfr_mul(r25977, r25976, r25976, MPFR_RNDN); mpfr_mul(r25978, r25974, r25970, MPFR_RNDN); mpfr_mul(r25979, r25978, r25978, MPFR_RNDN); ; mpfr_sub(r25981, r25979, r25980, MPFR_RNDN); mpfr_div(r25982, r25977, r25981, MPFR_RNDN); ; mpfr_div(r25984, r25983, r25970, MPFR_RNDN); mpfr_div(r25985, r25984, r25970, MPFR_RNDN); ; ; mpfr_pow(r25988, r25970, r25987, MPFR_RNDN); mpfr_div(r25989, r25986, r25988, MPFR_RNDN); ; mpfr_add(r25991, r25989, r25990, MPFR_RNDN); mpfr_add(r25992, r25985, r25991, MPFR_RNDN); if (mpfr_get_si(r25972, MPFR_RNDN)) { mpfr_set(r25993, r25982, MPFR_RNDN); } else { mpfr_set(r25993, r25992, MPFR_RNDN); }; return mpfr_get_d(r25993, MPFR_RNDN); }