#include #include #include #include #include char *name = "(+ x (- (tan (+ y z)) (tan a)))"; double f_if(float x, float y, float z, float a) { float r24892 = x; float r24893 = y; float r24894 = z; float r24895 = r24893 + r24894; float r24896 = tan(r24895); float r24897 = a; float r24898 = tan(r24897); float r24899 = r24896 - r24898; float r24900 = r24892 + r24899; return r24900; } double f_id(double x, double y, double z, double a) { double r24901 = x; double r24902 = y; double r24903 = z; double r24904 = r24902 + r24903; double r24905 = tan(r24904); double r24906 = a; double r24907 = tan(r24906); double r24908 = r24905 - r24907; double r24909 = r24901 + r24908; return r24909; } double f_of(float x, float y, float z, float a) { float r24910 = x; float r24911 = z; float r24912 = tan(r24911); float r24913 = y; float r24914 = tan(r24913); float r24915 = r24912 + r24914; float r24916 = a; float r24917 = cos(r24916); float r24918 = r24915 * r24917; float r24919 = sin(r24916); float r24920 = 1; float r24921 = r24912 * r24914; float r24922 = r24920 - r24921; float r24923 = r24919 * r24922; float r24924 = r24918 - r24923; float r24925 = 3; float r24926 = pow(r24921, r24925); float r24927 = r24920 - r24926; float r24928 = r24917 * r24927; float r24929 = r24924 / r24928; float r24930 = r24920 * r24920; float r24931 = r24914 * r24912; float r24932 = r24931 * r24931; float r24933 = r24920 * r24931; float r24934 = r24932 + r24933; float r24935 = r24930 + r24934; float r24936 = r24929 * r24935; float r24937 = r24910 + r24936; return r24937; } double f_od(double x, double y, double z, double a) { double r24938 = x; double r24939 = z; double r24940 = tan(r24939); double r24941 = y; double r24942 = tan(r24941); double r24943 = r24940 + r24942; double r24944 = a; double r24945 = cos(r24944); double r24946 = r24943 * r24945; double r24947 = sin(r24944); double r24948 = 1; double r24949 = r24940 * r24942; double r24950 = r24948 - r24949; double r24951 = r24947 * r24950; double r24952 = r24946 - r24951; double r24953 = 3; double r24954 = pow(r24949, r24953); double r24955 = r24948 - r24954; double r24956 = r24945 * r24955; double r24957 = r24952 / r24956; double r24958 = r24948 * r24948; double r24959 = r24942 * r24940; double r24960 = r24959 * r24959; double r24961 = r24948 * r24959; double r24962 = r24960 + r24961; double r24963 = r24958 + r24962; double r24964 = r24957 * r24963; double r24965 = r24938 + r24964; return r24965; } 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 r24966, r24967, r24968, r24969, r24970, r24971, r24972, r24973, r24974; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); 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); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r24966, x, MPFR_RNDN); mpfr_set_d(r24967, y, MPFR_RNDN); mpfr_set_d(r24968, z, MPFR_RNDN); mpfr_add(r24969, r24967, r24968, MPFR_RNDN); mpfr_tan(r24970, r24969, MPFR_RNDN); mpfr_set_d(r24971, a, MPFR_RNDN); mpfr_tan(r24972, r24971, MPFR_RNDN); mpfr_sub(r24973, r24970, r24972, MPFR_RNDN); mpfr_add(r24974, r24966, r24973, MPFR_RNDN); return mpfr_get_d(r24974, MPFR_RNDN); } static mpfr_t r24975, r24976, r24977, r24978, r24979, 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_fm() { mpfr_set_default_prec(1424); mpfr_init(r24975); mpfr_init(r24976); mpfr_init(r24977); mpfr_init(r24978); mpfr_init(r24979); mpfr_init(r24980); mpfr_init(r24981); mpfr_init(r24982); mpfr_init(r24983); mpfr_init(r24984); mpfr_init_set_str(r24985, "1", 10, MPFR_RNDN); mpfr_init(r24986); mpfr_init(r24987); mpfr_init(r24988); mpfr_init(r24989); mpfr_init_set_str(r24990, "3", 10, MPFR_RNDN); 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_fm(double x, double y, double z, double a) { mpfr_set_d(r24975, x, MPFR_RNDN); mpfr_set_d(r24976, z, MPFR_RNDN); mpfr_tan(r24977, r24976, MPFR_RNDN); mpfr_set_d(r24978, y, MPFR_RNDN); mpfr_tan(r24979, r24978, MPFR_RNDN); mpfr_add(r24980, r24977, r24979, MPFR_RNDN); mpfr_set_d(r24981, a, MPFR_RNDN); mpfr_cos(r24982, r24981, MPFR_RNDN); mpfr_mul(r24983, r24980, r24982, MPFR_RNDN); mpfr_sin(r24984, r24981, MPFR_RNDN); ; mpfr_mul(r24986, r24977, r24979, MPFR_RNDN); mpfr_sub(r24987, r24985, r24986, MPFR_RNDN); mpfr_mul(r24988, r24984, r24987, MPFR_RNDN); mpfr_sub(r24989, r24983, r24988, MPFR_RNDN); ; mpfr_pow(r24991, r24986, r24990, MPFR_RNDN); mpfr_sub(r24992, r24985, r24991, MPFR_RNDN); mpfr_mul(r24993, r24982, r24992, MPFR_RNDN); mpfr_div(r24994, r24989, r24993, MPFR_RNDN); mpfr_mul(r24995, r24985, r24985, MPFR_RNDN); mpfr_mul(r24996, r24979, r24977, MPFR_RNDN); mpfr_mul(r24997, r24996, r24996, MPFR_RNDN); mpfr_mul(r24998, r24985, r24996, MPFR_RNDN); mpfr_add(r24999, r24997, r24998, MPFR_RNDN); mpfr_add(r25000, r24995, r24999, MPFR_RNDN); mpfr_mul(r25001, r24994, r25000, MPFR_RNDN); mpfr_add(r25002, r24975, r25001, MPFR_RNDN); return mpfr_get_d(r25002, MPFR_RNDN); } static mpfr_t r25003, r25004, r25005, r25006, r25007, r25008, r25009, r25010, r25011, r25012, r25013, r25014, r25015, r25016, r25017, r25018, r25019, r25020, r25021, r25022, r25023, r25024, r25025, r25026, r25027, r25028, r25029, r25030; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r25003); mpfr_init(r25004); mpfr_init(r25005); mpfr_init(r25006); mpfr_init(r25007); mpfr_init(r25008); mpfr_init(r25009); mpfr_init(r25010); mpfr_init(r25011); mpfr_init(r25012); mpfr_init_set_str(r25013, "1", 10, MPFR_RNDN); mpfr_init(r25014); mpfr_init(r25015); mpfr_init(r25016); mpfr_init(r25017); mpfr_init_set_str(r25018, "3", 10, MPFR_RNDN); mpfr_init(r25019); mpfr_init(r25020); mpfr_init(r25021); mpfr_init(r25022); mpfr_init(r25023); mpfr_init(r25024); mpfr_init(r25025); mpfr_init(r25026); mpfr_init(r25027); mpfr_init(r25028); mpfr_init(r25029); mpfr_init(r25030); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r25003, x, MPFR_RNDN); mpfr_set_d(r25004, z, MPFR_RNDN); mpfr_tan(r25005, r25004, MPFR_RNDN); mpfr_set_d(r25006, y, MPFR_RNDN); mpfr_tan(r25007, r25006, MPFR_RNDN); mpfr_add(r25008, r25005, r25007, MPFR_RNDN); mpfr_set_d(r25009, a, MPFR_RNDN); mpfr_cos(r25010, r25009, MPFR_RNDN); mpfr_mul(r25011, r25008, r25010, MPFR_RNDN); mpfr_sin(r25012, r25009, MPFR_RNDN); ; mpfr_mul(r25014, r25005, r25007, MPFR_RNDN); mpfr_sub(r25015, r25013, r25014, MPFR_RNDN); mpfr_mul(r25016, r25012, r25015, MPFR_RNDN); mpfr_sub(r25017, r25011, r25016, MPFR_RNDN); ; mpfr_pow(r25019, r25014, r25018, MPFR_RNDN); mpfr_sub(r25020, r25013, r25019, MPFR_RNDN); mpfr_mul(r25021, r25010, r25020, MPFR_RNDN); mpfr_div(r25022, r25017, r25021, MPFR_RNDN); mpfr_mul(r25023, r25013, r25013, MPFR_RNDN); mpfr_mul(r25024, r25007, r25005, MPFR_RNDN); mpfr_mul(r25025, r25024, r25024, MPFR_RNDN); mpfr_mul(r25026, r25013, r25024, MPFR_RNDN); mpfr_add(r25027, r25025, r25026, MPFR_RNDN); mpfr_add(r25028, r25023, r25027, MPFR_RNDN); mpfr_mul(r25029, r25022, r25028, MPFR_RNDN); mpfr_add(r25030, r25003, r25029, MPFR_RNDN); return mpfr_get_d(r25030, MPFR_RNDN); }