#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r21865 = i; float r21866 = alpha; float r21867 = beta; float r21868 = r21866 + r21867; float r21869 = r21868 + r21865; float r21870 = r21865 * r21869; float r21871 = r21867 * r21866; float r21872 = r21871 + r21870; float r21873 = r21870 * r21872; float r21874 = 2; float r21875 = r21874 * r21865; float r21876 = r21868 + r21875; float r21877 = r21876 * r21876; float r21878 = r21873 / r21877; float r21879 = 1.0; float r21880 = r21877 - r21879; float r21881 = r21878 / r21880; return r21881; } double f_id(double alpha, double beta, double i) { double r21882 = i; double r21883 = alpha; double r21884 = beta; double r21885 = r21883 + r21884; double r21886 = r21885 + r21882; double r21887 = r21882 * r21886; double r21888 = r21884 * r21883; double r21889 = r21888 + r21887; double r21890 = r21887 * r21889; double r21891 = 2; double r21892 = r21891 * r21882; double r21893 = r21885 + r21892; double r21894 = r21893 * r21893; double r21895 = r21890 / r21894; double r21896 = 1.0; double r21897 = r21894 - r21896; double r21898 = r21895 / r21897; return r21898; } double f_of(float alpha, float beta, float i) { float r21899 = alpha; float r21900 = 1.5352145518006176e+162; bool r21901 = r21899 <= r21900; float r21902 = i; float r21903 = beta; float r21904 = r21899 + r21903; float r21905 = r21904 + r21902; float r21906 = r21902 * r21905; float r21907 = 2; float r21908 = r21907 * r21902; float r21909 = r21904 + r21908; float r21910 = r21909 * r21909; float r21911 = r21903 * r21899; float r21912 = r21911 + r21906; float r21913 = r21910 / r21912; float r21914 = r21906 / r21913; float r21915 = 1.0; float r21916 = r21910 - r21915; float r21917 = r21914 / r21916; float r21918 = 0; float r21919 = r21901 ? r21917 : r21918; return r21919; } double f_od(double alpha, double beta, double i) { double r21920 = alpha; double r21921 = 1.5352145518006176e+162; bool r21922 = r21920 <= r21921; double r21923 = i; double r21924 = beta; double r21925 = r21920 + r21924; double r21926 = r21925 + r21923; double r21927 = r21923 * r21926; double r21928 = 2; double r21929 = r21928 * r21923; double r21930 = r21925 + r21929; double r21931 = r21930 * r21930; double r21932 = r21924 * r21920; double r21933 = r21932 + r21927; double r21934 = r21931 / r21933; double r21935 = r21927 / r21934; double r21936 = 1.0; double r21937 = r21931 - r21936; double r21938 = r21935 / r21937; double r21939 = 0; double r21940 = r21922 ? r21938 : r21939; return r21940; } 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 r21941, r21942, r21943, r21944, r21945, r21946, r21947, r21948, r21949, r21950, r21951, r21952, r21953, r21954, r21955, r21956, r21957; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r21941); mpfr_init(r21942); mpfr_init(r21943); mpfr_init(r21944); mpfr_init(r21945); mpfr_init(r21946); mpfr_init(r21947); mpfr_init(r21948); mpfr_init(r21949); mpfr_init_set_str(r21950, "2", 10, MPFR_RNDN); mpfr_init(r21951); mpfr_init(r21952); mpfr_init(r21953); mpfr_init(r21954); mpfr_init_set_str(r21955, "1.0", 10, MPFR_RNDN); mpfr_init(r21956); mpfr_init(r21957); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r21941, i, MPFR_RNDN); mpfr_set_d(r21942, alpha, MPFR_RNDN); mpfr_set_d(r21943, beta, MPFR_RNDN); mpfr_add(r21944, r21942, r21943, MPFR_RNDN); mpfr_add(r21945, r21944, r21941, MPFR_RNDN); mpfr_mul(r21946, r21941, r21945, MPFR_RNDN); mpfr_mul(r21947, r21943, r21942, MPFR_RNDN); mpfr_add(r21948, r21947, r21946, MPFR_RNDN); mpfr_mul(r21949, r21946, r21948, MPFR_RNDN); ; mpfr_mul(r21951, r21950, r21941, MPFR_RNDN); mpfr_add(r21952, r21944, r21951, MPFR_RNDN); mpfr_mul(r21953, r21952, r21952, MPFR_RNDN); mpfr_div(r21954, r21949, r21953, MPFR_RNDN); ; mpfr_sub(r21956, r21953, r21955, MPFR_RNDN); mpfr_div(r21957, r21954, r21956, MPFR_RNDN); return mpfr_get_d(r21957, MPFR_RNDN); } static mpfr_t r21958, r21959, r21960, r21961, r21962, r21963, r21964, r21965, r21966, r21967, r21968, r21969, r21970, r21971, r21972, r21973, r21974, r21975, r21976, r21977, r21978; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r21958); mpfr_init_set_str(r21959, "1.5352145518006176e+162", 10, MPFR_RNDN); mpfr_init(r21960); mpfr_init(r21961); mpfr_init(r21962); mpfr_init(r21963); mpfr_init(r21964); mpfr_init(r21965); mpfr_init_set_str(r21966, "2", 10, MPFR_RNDN); mpfr_init(r21967); mpfr_init(r21968); mpfr_init(r21969); mpfr_init(r21970); mpfr_init(r21971); mpfr_init(r21972); mpfr_init(r21973); mpfr_init_set_str(r21974, "1.0", 10, MPFR_RNDN); mpfr_init(r21975); mpfr_init(r21976); mpfr_init_set_str(r21977, "0", 10, MPFR_RNDN); mpfr_init(r21978); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r21958, alpha, MPFR_RNDN); ; mpfr_set_si(r21960, mpfr_cmp(r21958, r21959) <= 0, MPFR_RNDN); mpfr_set_d(r21961, i, MPFR_RNDN); mpfr_set_d(r21962, beta, MPFR_RNDN); mpfr_add(r21963, r21958, r21962, MPFR_RNDN); mpfr_add(r21964, r21963, r21961, MPFR_RNDN); mpfr_mul(r21965, r21961, r21964, MPFR_RNDN); ; mpfr_mul(r21967, r21966, r21961, MPFR_RNDN); mpfr_add(r21968, r21963, r21967, MPFR_RNDN); mpfr_mul(r21969, r21968, r21968, MPFR_RNDN); mpfr_mul(r21970, r21962, r21958, MPFR_RNDN); mpfr_add(r21971, r21970, r21965, MPFR_RNDN); mpfr_div(r21972, r21969, r21971, MPFR_RNDN); mpfr_div(r21973, r21965, r21972, MPFR_RNDN); ; mpfr_sub(r21975, r21969, r21974, MPFR_RNDN); mpfr_div(r21976, r21973, r21975, MPFR_RNDN); ; if (mpfr_get_si(r21960, MPFR_RNDN)) { mpfr_set(r21978, r21976, MPFR_RNDN); } else { mpfr_set(r21978, r21977, MPFR_RNDN); }; return mpfr_get_d(r21978, MPFR_RNDN); } static mpfr_t r21979, r21980, r21981, r21982, r21983, r21984, r21985, r21986, r21987, r21988, r21989, r21990, r21991, r21992, r21993, r21994, r21995, r21996, r21997, r21998, r21999; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r21979); mpfr_init_set_str(r21980, "1.5352145518006176e+162", 10, MPFR_RNDN); mpfr_init(r21981); mpfr_init(r21982); mpfr_init(r21983); mpfr_init(r21984); mpfr_init(r21985); mpfr_init(r21986); mpfr_init_set_str(r21987, "2", 10, MPFR_RNDN); mpfr_init(r21988); mpfr_init(r21989); mpfr_init(r21990); mpfr_init(r21991); mpfr_init(r21992); mpfr_init(r21993); mpfr_init(r21994); mpfr_init_set_str(r21995, "1.0", 10, MPFR_RNDN); mpfr_init(r21996); mpfr_init(r21997); mpfr_init_set_str(r21998, "0", 10, MPFR_RNDN); mpfr_init(r21999); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r21979, alpha, MPFR_RNDN); ; mpfr_set_si(r21981, mpfr_cmp(r21979, r21980) <= 0, MPFR_RNDN); mpfr_set_d(r21982, i, MPFR_RNDN); mpfr_set_d(r21983, beta, MPFR_RNDN); mpfr_add(r21984, r21979, r21983, MPFR_RNDN); mpfr_add(r21985, r21984, r21982, MPFR_RNDN); mpfr_mul(r21986, r21982, r21985, MPFR_RNDN); ; mpfr_mul(r21988, r21987, r21982, MPFR_RNDN); mpfr_add(r21989, r21984, r21988, MPFR_RNDN); mpfr_mul(r21990, r21989, r21989, MPFR_RNDN); mpfr_mul(r21991, r21983, r21979, MPFR_RNDN); mpfr_add(r21992, r21991, r21986, MPFR_RNDN); mpfr_div(r21993, r21990, r21992, MPFR_RNDN); mpfr_div(r21994, r21986, r21993, MPFR_RNDN); ; mpfr_sub(r21996, r21990, r21995, MPFR_RNDN); mpfr_div(r21997, r21994, r21996, MPFR_RNDN); ; if (mpfr_get_si(r21981, MPFR_RNDN)) { mpfr_set(r21999, r21997, MPFR_RNDN); } else { mpfr_set(r21999, r21998, MPFR_RNDN); }; return mpfr_get_d(r21999, MPFR_RNDN); }