#include #include #include #include #include char *name = "Harley's example"; double f_if(float c_p, float c_n, float t, float s) { float r24845 = 1; float r24846 = s; float r24847 = -r24846; float r24848 = exp(r24847); float r24849 = r24845 + r24848; float r24850 = r24845 / r24849; float r24851 = c_p; float r24852 = pow(r24850, r24851); float r24853 = r24845 - r24850; float r24854 = c_n; float r24855 = pow(r24853, r24854); float r24856 = r24852 * r24855; float r24857 = t; float r24858 = -r24857; float r24859 = exp(r24858); float r24860 = r24845 + r24859; float r24861 = r24845 / r24860; float r24862 = pow(r24861, r24851); float r24863 = r24845 - r24861; float r24864 = pow(r24863, r24854); float r24865 = r24862 * r24864; float r24866 = r24856 / r24865; return r24866; } double f_id(double c_p, double c_n, double t, double s) { double r24867 = 1; double r24868 = s; double r24869 = -r24868; double r24870 = exp(r24869); double r24871 = r24867 + r24870; double r24872 = r24867 / r24871; double r24873 = c_p; double r24874 = pow(r24872, r24873); double r24875 = r24867 - r24872; double r24876 = c_n; double r24877 = pow(r24875, r24876); double r24878 = r24874 * r24877; double r24879 = t; double r24880 = -r24879; double r24881 = exp(r24880); double r24882 = r24867 + r24881; double r24883 = r24867 / r24882; double r24884 = pow(r24883, r24873); double r24885 = r24867 - r24883; double r24886 = pow(r24885, r24876); double r24887 = r24884 * r24886; double r24888 = r24878 / r24887; return r24888; } double f_of(float c_p, float c_n, float t, float s) { float r24889 = 1; float r24890 = -r24889; float r24891 = s; float r24892 = -r24891; float r24893 = exp(r24892); float r24894 = r24889 + r24893; float r24895 = r24890 / r24894; float r24896 = log1p(r24895); float r24897 = t; float r24898 = -r24897; float r24899 = exp(r24898); float r24900 = r24899 + r24889; float r24901 = r24890 / r24900; float r24902 = log1p(r24901); float r24903 = r24896 - r24902; float r24904 = c_n; float r24905 = c_p; float r24906 = log1p(r24899); float r24907 = log1p(r24893); float r24908 = r24906 - r24907; float r24909 = r24905 * r24908; float r24910 = fma(r24903, r24904, r24909); float r24911 = exp(r24910); return r24911; } double f_od(double c_p, double c_n, double t, double s) { double r24912 = 1; double r24913 = -r24912; double r24914 = s; double r24915 = -r24914; double r24916 = exp(r24915); double r24917 = r24912 + r24916; double r24918 = r24913 / r24917; double r24919 = log1p(r24918); double r24920 = t; double r24921 = -r24920; double r24922 = exp(r24921); double r24923 = r24922 + r24912; double r24924 = r24913 / r24923; double r24925 = log1p(r24924); double r24926 = r24919 - r24925; double r24927 = c_n; double r24928 = c_p; double r24929 = log1p(r24922); double r24930 = log1p(r24916); double r24931 = r24929 - r24930; double r24932 = r24928 * r24931; double r24933 = fma(r24926, r24927, r24932); double r24934 = exp(r24933); return r24934; } 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 r24935, r24936, r24937, r24938, r24939, r24940, r24941, r24942, r24943, r24944, r24945, r24946, r24947, r24948, r24949, r24950, r24951, r24952, r24953, r24954, r24955, r24956; void setup_mpfr_f_im() { mpfr_set_default_prec(1616); mpfr_init_set_str(r24935, "1", 10, MPFR_RNDN); mpfr_init(r24936); mpfr_init(r24937); mpfr_init(r24938); mpfr_init(r24939); mpfr_init(r24940); mpfr_init(r24941); mpfr_init(r24942); mpfr_init(r24943); mpfr_init(r24944); mpfr_init(r24945); mpfr_init(r24946); mpfr_init(r24947); mpfr_init(r24948); mpfr_init(r24949); mpfr_init(r24950); mpfr_init(r24951); mpfr_init(r24952); mpfr_init(r24953); mpfr_init(r24954); mpfr_init(r24955); mpfr_init(r24956); } double f_im(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r24936, s, MPFR_RNDN); mpfr_neg(r24937, r24936, MPFR_RNDN); mpfr_exp(r24938, r24937, MPFR_RNDN); mpfr_add(r24939, r24935, r24938, MPFR_RNDN); mpfr_div(r24940, r24935, r24939, MPFR_RNDN); mpfr_set_d(r24941, c_p, MPFR_RNDN); mpfr_pow(r24942, r24940, r24941, MPFR_RNDN); mpfr_sub(r24943, r24935, r24940, MPFR_RNDN); mpfr_set_d(r24944, c_n, MPFR_RNDN); mpfr_pow(r24945, r24943, r24944, MPFR_RNDN); mpfr_mul(r24946, r24942, r24945, MPFR_RNDN); mpfr_set_d(r24947, t, MPFR_RNDN); mpfr_neg(r24948, r24947, MPFR_RNDN); mpfr_exp(r24949, r24948, MPFR_RNDN); mpfr_add(r24950, r24935, r24949, MPFR_RNDN); mpfr_div(r24951, r24935, r24950, MPFR_RNDN); mpfr_pow(r24952, r24951, r24941, MPFR_RNDN); mpfr_sub(r24953, r24935, r24951, MPFR_RNDN); mpfr_pow(r24954, r24953, r24944, MPFR_RNDN); mpfr_mul(r24955, r24952, r24954, MPFR_RNDN); mpfr_div(r24956, r24946, r24955, MPFR_RNDN); return mpfr_get_d(r24956, MPFR_RNDN); } static mpfr_t r24957, r24958, r24959, r24960, r24961, r24962, r24963, r24964, r24965, r24966, r24967, r24968, r24969, r24970, r24971, r24972, r24973, r24974, r24975, r24976, r24977, r24978, r24979; void setup_mpfr_f_fm() { mpfr_set_default_prec(1616); mpfr_init_set_str(r24957, "1", 10, MPFR_RNDN); mpfr_init(r24958); mpfr_init(r24959); mpfr_init(r24960); mpfr_init(r24961); mpfr_init(r24962); mpfr_init(r24963); mpfr_init(r24964); mpfr_init(r24965); mpfr_init(r24966); mpfr_init(r24967); mpfr_init(r24968); mpfr_init(r24969); mpfr_init(r24970); mpfr_init(r24971); mpfr_init(r24972); mpfr_init(r24973); mpfr_init(r24974); mpfr_init(r24975); mpfr_init(r24976); mpfr_init(r24977); mpfr_init(r24978); mpfr_init(r24979); } double f_fm(double c_p, double c_n, double t, double s) { ; mpfr_neg(r24958, r24957, MPFR_RNDN); mpfr_set_d(r24959, s, MPFR_RNDN); mpfr_neg(r24960, r24959, MPFR_RNDN); mpfr_exp(r24961, r24960, MPFR_RNDN); mpfr_add(r24962, r24957, r24961, MPFR_RNDN); mpfr_div(r24963, r24958, r24962, MPFR_RNDN); mpfr_log1p(r24964, r24963, MPFR_RNDN); mpfr_set_d(r24965, t, MPFR_RNDN); mpfr_neg(r24966, r24965, MPFR_RNDN); mpfr_exp(r24967, r24966, MPFR_RNDN); mpfr_add(r24968, r24967, r24957, MPFR_RNDN); mpfr_div(r24969, r24958, r24968, MPFR_RNDN); mpfr_log1p(r24970, r24969, MPFR_RNDN); mpfr_sub(r24971, r24964, r24970, MPFR_RNDN); mpfr_set_d(r24972, c_n, MPFR_RNDN); mpfr_set_d(r24973, c_p, MPFR_RNDN); mpfr_log1p(r24974, r24967, MPFR_RNDN); mpfr_log1p(r24975, r24961, MPFR_RNDN); mpfr_sub(r24976, r24974, r24975, MPFR_RNDN); mpfr_mul(r24977, r24973, r24976, MPFR_RNDN); mpfr_fma(r24978, r24971, r24972, r24977, MPFR_RNDN); mpfr_exp(r24979, r24978, MPFR_RNDN); return mpfr_get_d(r24979, MPFR_RNDN); } static mpfr_t r24980, r24981, r24982, r24983, r24984, r24985, r24986, r24987, r24988, r24989, r24990, r24991, r24992, r24993, r24994, r24995, r24996, r24997, r24998, r24999, r25000, r25001, r25002; void setup_mpfr_f_dm() { mpfr_set_default_prec(1616); mpfr_init_set_str(r24980, "1", 10, MPFR_RNDN); mpfr_init(r24981); mpfr_init(r24982); mpfr_init(r24983); mpfr_init(r24984); mpfr_init(r24985); mpfr_init(r24986); mpfr_init(r24987); mpfr_init(r24988); mpfr_init(r24989); mpfr_init(r24990); mpfr_init(r24991); mpfr_init(r24992); mpfr_init(r24993); mpfr_init(r24994); mpfr_init(r24995); mpfr_init(r24996); mpfr_init(r24997); mpfr_init(r24998); mpfr_init(r24999); mpfr_init(r25000); mpfr_init(r25001); mpfr_init(r25002); } double f_dm(double c_p, double c_n, double t, double s) { ; mpfr_neg(r24981, r24980, MPFR_RNDN); mpfr_set_d(r24982, s, MPFR_RNDN); mpfr_neg(r24983, r24982, MPFR_RNDN); mpfr_exp(r24984, r24983, MPFR_RNDN); mpfr_add(r24985, r24980, r24984, MPFR_RNDN); mpfr_div(r24986, r24981, r24985, MPFR_RNDN); mpfr_log1p(r24987, r24986, MPFR_RNDN); mpfr_set_d(r24988, t, MPFR_RNDN); mpfr_neg(r24989, r24988, MPFR_RNDN); mpfr_exp(r24990, r24989, MPFR_RNDN); mpfr_add(r24991, r24990, r24980, MPFR_RNDN); mpfr_div(r24992, r24981, r24991, MPFR_RNDN); mpfr_log1p(r24993, r24992, MPFR_RNDN); mpfr_sub(r24994, r24987, r24993, MPFR_RNDN); mpfr_set_d(r24995, c_n, MPFR_RNDN); mpfr_set_d(r24996, c_p, MPFR_RNDN); mpfr_log1p(r24997, r24990, MPFR_RNDN); mpfr_log1p(r24998, r24984, MPFR_RNDN); mpfr_sub(r24999, r24997, r24998, MPFR_RNDN); mpfr_mul(r25000, r24996, r24999, MPFR_RNDN); mpfr_fma(r25001, r24994, r24995, r25000, MPFR_RNDN); mpfr_exp(r25002, r25001, MPFR_RNDN); return mpfr_get_d(r25002, MPFR_RNDN); }