#include #include #include #include #include char *name = "Octave 3.8, jcobi/3"; double f_if(float alpha, float beta) { float r9843 = alpha; float r9844 = beta; float r9845 = r9843 + r9844; float r9846 = r9844 * r9843; float r9847 = r9845 + r9846; float r9848 = 1.0; float r9849 = r9847 + r9848; float r9850 = 2; float r9851 = 1; float r9852 = r9850 * r9851; float r9853 = r9845 + r9852; float r9854 = r9849 / r9853; float r9855 = r9854 / r9853; float r9856 = r9853 + r9848; float r9857 = r9855 / r9856; return r9857; } double f_id(double alpha, double beta) { double r9858 = alpha; double r9859 = beta; double r9860 = r9858 + r9859; double r9861 = r9859 * r9858; double r9862 = r9860 + r9861; double r9863 = 1.0; double r9864 = r9862 + r9863; double r9865 = 2; double r9866 = 1; double r9867 = r9865 * r9866; double r9868 = r9860 + r9867; double r9869 = r9864 / r9868; double r9870 = r9869 / r9868; double r9871 = r9868 + r9863; double r9872 = r9870 / r9871; return r9872; } double f_of(float alpha, float beta) { float r9873 = 1.0; float r9874 = alpha; float r9875 = beta; float r9876 = r9874 * r9875; float r9877 = r9875 + r9874; float r9878 = r9876 + r9877; float r9879 = r9873 + r9878; float r9880 = 2; float r9881 = r9880 + r9877; float r9882 = r9879 / r9881; float r9883 = r9882 / r9881; float r9884 = r9881 + r9873; float r9885 = r9883 / r9884; float r9886 = 7.774571895904472e+307; bool r9887 = r9885 <= r9886; float r9888 = r9875 + r9873; float r9889 = fma(r9874, r9875, r9874); float r9890 = r9888 + r9889; float r9891 = r9874 + r9880; float r9892 = r9891 + r9875; float r9893 = r9890 / r9892; float r9894 = sqrt(r9893); float r9895 = r9873 + r9874; float r9896 = r9880 + r9875; float r9897 = r9895 + r9896; float r9898 = sqrt(r9892); float r9899 = r9897 * r9898; float r9900 = r9894 / r9899; float r9901 = fma(r9874, r9875, r9875); float r9902 = r9901 + r9895; float r9903 = r9902 / r9892; float r9904 = sqrt(r9903); float r9905 = r9904 / r9898; float r9906 = r9900 * r9905; float r9907 = 1; float r9908 = r9907 / r9874; float r9909 = 2.0; float r9910 = r9909 / r9874; float r9911 = r9910 - r9873; float r9912 = fma(r9908, r9911, r9907); float r9913 = r9897 * r9892; float r9914 = r9912 / r9913; float r9915 = r9887 ? r9906 : r9914; return r9915; } double f_od(double alpha, double beta) { double r9916 = 1.0; double r9917 = alpha; double r9918 = beta; double r9919 = r9917 * r9918; double r9920 = r9918 + r9917; double r9921 = r9919 + r9920; double r9922 = r9916 + r9921; double r9923 = 2; double r9924 = r9923 + r9920; double r9925 = r9922 / r9924; double r9926 = r9925 / r9924; double r9927 = r9924 + r9916; double r9928 = r9926 / r9927; double r9929 = 7.774571895904472e+307; bool r9930 = r9928 <= r9929; double r9931 = r9918 + r9916; double r9932 = fma(r9917, r9918, r9917); double r9933 = r9931 + r9932; double r9934 = r9917 + r9923; double r9935 = r9934 + r9918; double r9936 = r9933 / r9935; double r9937 = sqrt(r9936); double r9938 = r9916 + r9917; double r9939 = r9923 + r9918; double r9940 = r9938 + r9939; double r9941 = sqrt(r9935); double r9942 = r9940 * r9941; double r9943 = r9937 / r9942; double r9944 = fma(r9917, r9918, r9918); double r9945 = r9944 + r9938; double r9946 = r9945 / r9935; double r9947 = sqrt(r9946); double r9948 = r9947 / r9941; double r9949 = r9943 * r9948; double r9950 = 1; double r9951 = r9950 / r9917; double r9952 = 2.0; double r9953 = r9952 / r9917; double r9954 = r9953 - r9916; double r9955 = fma(r9951, r9954, r9950); double r9956 = r9940 * r9935; double r9957 = r9955 / r9956; double r9958 = r9930 ? r9949 : r9957; return r9958; } 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 r9959, r9960, r9961, r9962, r9963, r9964, r9965, r9966, r9967, r9968, r9969, r9970, r9971, r9972, r9973; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r9959); mpfr_init(r9960); mpfr_init(r9961); mpfr_init(r9962); mpfr_init(r9963); mpfr_init_set_str(r9964, "1.0", 10, MPFR_RNDN); mpfr_init(r9965); mpfr_init_set_str(r9966, "2", 10, MPFR_RNDN); mpfr_init_set_str(r9967, "1", 10, MPFR_RNDN); mpfr_init(r9968); mpfr_init(r9969); mpfr_init(r9970); mpfr_init(r9971); mpfr_init(r9972); mpfr_init(r9973); } double f_im(double alpha, double beta) { mpfr_set_d(r9959, alpha, MPFR_RNDN); mpfr_set_d(r9960, beta, MPFR_RNDN); mpfr_add(r9961, r9959, r9960, MPFR_RNDN); mpfr_mul(r9962, r9960, r9959, MPFR_RNDN); mpfr_add(r9963, r9961, r9962, MPFR_RNDN); ; mpfr_add(r9965, r9963, r9964, MPFR_RNDN); ; ; mpfr_mul(r9968, r9966, r9967, MPFR_RNDN); mpfr_add(r9969, r9961, r9968, MPFR_RNDN); mpfr_div(r9970, r9965, r9969, MPFR_RNDN); mpfr_div(r9971, r9970, r9969, MPFR_RNDN); mpfr_add(r9972, r9969, r9964, MPFR_RNDN); mpfr_div(r9973, r9971, r9972, MPFR_RNDN); return mpfr_get_d(r9973, MPFR_RNDN); } static mpfr_t r9974, r9975, r9976, r9977, r9978, r9979, r9980, r9981, r9982, r9983, r9984, r9985, r9986, r9987, r9988, r9989, r9990, r9991, r9992, r9993, r9994, r9995, r9996, r9997, r9998, r9999, r10000, r10001, r10002, r10003, r10004, r10005, r10006, r10007, r10008, r10009, r10010, r10011, r10012, r10013, r10014, r10015, r10016; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r9974, "1.0", 10, MPFR_RNDN); mpfr_init(r9975); mpfr_init(r9976); mpfr_init(r9977); mpfr_init(r9978); mpfr_init(r9979); mpfr_init(r9980); mpfr_init_set_str(r9981, "2", 10, MPFR_RNDN); mpfr_init(r9982); mpfr_init(r9983); mpfr_init(r9984); mpfr_init(r9985); mpfr_init(r9986); mpfr_init_set_str(r9987, "7.774571895904472e+307", 10, MPFR_RNDN); mpfr_init(r9988); mpfr_init(r9989); mpfr_init(r9990); mpfr_init(r9991); mpfr_init(r9992); mpfr_init(r9993); mpfr_init(r9994); mpfr_init(r9995); mpfr_init(r9996); mpfr_init(r9997); mpfr_init(r9998); mpfr_init(r9999); mpfr_init(r10000); mpfr_init(r10001); mpfr_init(r10002); mpfr_init(r10003); mpfr_init(r10004); mpfr_init(r10005); mpfr_init(r10006); mpfr_init(r10007); mpfr_init_set_str(r10008, "1", 10, MPFR_RNDN); mpfr_init(r10009); mpfr_init_set_str(r10010, "2.0", 10, MPFR_RNDN); mpfr_init(r10011); mpfr_init(r10012); mpfr_init(r10013); mpfr_init(r10014); mpfr_init(r10015); mpfr_init(r10016); } double f_fm(double alpha, double beta) { ; mpfr_set_d(r9975, alpha, MPFR_RNDN); mpfr_set_d(r9976, beta, MPFR_RNDN); mpfr_mul(r9977, r9975, r9976, MPFR_RNDN); mpfr_add(r9978, r9976, r9975, MPFR_RNDN); mpfr_add(r9979, r9977, r9978, MPFR_RNDN); mpfr_add(r9980, r9974, r9979, MPFR_RNDN); ; mpfr_add(r9982, r9981, r9978, MPFR_RNDN); mpfr_div(r9983, r9980, r9982, MPFR_RNDN); mpfr_div(r9984, r9983, r9982, MPFR_RNDN); mpfr_add(r9985, r9982, r9974, MPFR_RNDN); mpfr_div(r9986, r9984, r9985, MPFR_RNDN); ; mpfr_set_si(r9988, mpfr_cmp(r9986, r9987) <= 0, MPFR_RNDN); mpfr_add(r9989, r9976, r9974, MPFR_RNDN); mpfr_fma(r9990, r9975, r9976, r9975, MPFR_RNDN); mpfr_add(r9991, r9989, r9990, MPFR_RNDN); mpfr_add(r9992, r9975, r9981, MPFR_RNDN); mpfr_add(r9993, r9992, r9976, MPFR_RNDN); mpfr_div(r9994, r9991, r9993, MPFR_RNDN); mpfr_sqrt(r9995, r9994, MPFR_RNDN); mpfr_add(r9996, r9974, r9975, MPFR_RNDN); mpfr_add(r9997, r9981, r9976, MPFR_RNDN); mpfr_add(r9998, r9996, r9997, MPFR_RNDN); mpfr_sqrt(r9999, r9993, MPFR_RNDN); mpfr_mul(r10000, r9998, r9999, MPFR_RNDN); mpfr_div(r10001, r9995, r10000, MPFR_RNDN); mpfr_fma(r10002, r9975, r9976, r9976, MPFR_RNDN); mpfr_add(r10003, r10002, r9996, MPFR_RNDN); mpfr_div(r10004, r10003, r9993, MPFR_RNDN); mpfr_sqrt(r10005, r10004, MPFR_RNDN); mpfr_div(r10006, r10005, r9999, MPFR_RNDN); mpfr_mul(r10007, r10001, r10006, MPFR_RNDN); ; mpfr_div(r10009, r10008, r9975, MPFR_RNDN); ; mpfr_div(r10011, r10010, r9975, MPFR_RNDN); mpfr_sub(r10012, r10011, r9974, MPFR_RNDN); mpfr_fma(r10013, r10009, r10012, r10008, MPFR_RNDN); mpfr_mul(r10014, r9998, r9993, MPFR_RNDN); mpfr_div(r10015, r10013, r10014, MPFR_RNDN); if (mpfr_get_si(r9988, MPFR_RNDN)) { mpfr_set(r10016, r10007, MPFR_RNDN); } else { mpfr_set(r10016, r10015, MPFR_RNDN); }; return mpfr_get_d(r10016, MPFR_RNDN); } static mpfr_t r10017, r10018, r10019, r10020, r10021, r10022, r10023, r10024, r10025, r10026, r10027, r10028, r10029, r10030, r10031, r10032, r10033, r10034, r10035, r10036, r10037, r10038, r10039, r10040, r10041, r10042, r10043, r10044, r10045, r10046, r10047, r10048, r10049, r10050, r10051, r10052, r10053, r10054, r10055, r10056, r10057, r10058, r10059; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r10017, "1.0", 10, MPFR_RNDN); mpfr_init(r10018); mpfr_init(r10019); mpfr_init(r10020); mpfr_init(r10021); mpfr_init(r10022); mpfr_init(r10023); mpfr_init_set_str(r10024, "2", 10, MPFR_RNDN); mpfr_init(r10025); mpfr_init(r10026); mpfr_init(r10027); mpfr_init(r10028); mpfr_init(r10029); mpfr_init_set_str(r10030, "7.774571895904472e+307", 10, MPFR_RNDN); mpfr_init(r10031); mpfr_init(r10032); mpfr_init(r10033); mpfr_init(r10034); mpfr_init(r10035); mpfr_init(r10036); mpfr_init(r10037); mpfr_init(r10038); mpfr_init(r10039); mpfr_init(r10040); mpfr_init(r10041); mpfr_init(r10042); mpfr_init(r10043); mpfr_init(r10044); mpfr_init(r10045); mpfr_init(r10046); mpfr_init(r10047); mpfr_init(r10048); mpfr_init(r10049); mpfr_init(r10050); mpfr_init_set_str(r10051, "1", 10, MPFR_RNDN); mpfr_init(r10052); mpfr_init_set_str(r10053, "2.0", 10, MPFR_RNDN); mpfr_init(r10054); mpfr_init(r10055); mpfr_init(r10056); mpfr_init(r10057); mpfr_init(r10058); mpfr_init(r10059); } double f_dm(double alpha, double beta) { ; mpfr_set_d(r10018, alpha, MPFR_RNDN); mpfr_set_d(r10019, beta, MPFR_RNDN); mpfr_mul(r10020, r10018, r10019, MPFR_RNDN); mpfr_add(r10021, r10019, r10018, MPFR_RNDN); mpfr_add(r10022, r10020, r10021, MPFR_RNDN); mpfr_add(r10023, r10017, r10022, MPFR_RNDN); ; mpfr_add(r10025, r10024, r10021, MPFR_RNDN); mpfr_div(r10026, r10023, r10025, MPFR_RNDN); mpfr_div(r10027, r10026, r10025, MPFR_RNDN); mpfr_add(r10028, r10025, r10017, MPFR_RNDN); mpfr_div(r10029, r10027, r10028, MPFR_RNDN); ; mpfr_set_si(r10031, mpfr_cmp(r10029, r10030) <= 0, MPFR_RNDN); mpfr_add(r10032, r10019, r10017, MPFR_RNDN); mpfr_fma(r10033, r10018, r10019, r10018, MPFR_RNDN); mpfr_add(r10034, r10032, r10033, MPFR_RNDN); mpfr_add(r10035, r10018, r10024, MPFR_RNDN); mpfr_add(r10036, r10035, r10019, MPFR_RNDN); mpfr_div(r10037, r10034, r10036, MPFR_RNDN); mpfr_sqrt(r10038, r10037, MPFR_RNDN); mpfr_add(r10039, r10017, r10018, MPFR_RNDN); mpfr_add(r10040, r10024, r10019, MPFR_RNDN); mpfr_add(r10041, r10039, r10040, MPFR_RNDN); mpfr_sqrt(r10042, r10036, MPFR_RNDN); mpfr_mul(r10043, r10041, r10042, MPFR_RNDN); mpfr_div(r10044, r10038, r10043, MPFR_RNDN); mpfr_fma(r10045, r10018, r10019, r10019, MPFR_RNDN); mpfr_add(r10046, r10045, r10039, MPFR_RNDN); mpfr_div(r10047, r10046, r10036, MPFR_RNDN); mpfr_sqrt(r10048, r10047, MPFR_RNDN); mpfr_div(r10049, r10048, r10042, MPFR_RNDN); mpfr_mul(r10050, r10044, r10049, MPFR_RNDN); ; mpfr_div(r10052, r10051, r10018, MPFR_RNDN); ; mpfr_div(r10054, r10053, r10018, MPFR_RNDN); mpfr_sub(r10055, r10054, r10017, MPFR_RNDN); mpfr_fma(r10056, r10052, r10055, r10051, MPFR_RNDN); mpfr_mul(r10057, r10041, r10036, MPFR_RNDN); mpfr_div(r10058, r10056, r10057, MPFR_RNDN); if (mpfr_get_si(r10031, MPFR_RNDN)) { mpfr_set(r10059, r10050, MPFR_RNDN); } else { mpfr_set(r10059, r10058, MPFR_RNDN); }; return mpfr_get_d(r10059, MPFR_RNDN); }