#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 r24884 = x; float r24885 = y; float r24886 = z; float r24887 = r24885 + r24886; float r24888 = tan(r24887); float r24889 = a; float r24890 = tan(r24889); float r24891 = r24888 - r24890; float r24892 = r24884 + r24891; return r24892; } double f_id(double x, double y, double z, double a) { double r24893 = x; double r24894 = y; double r24895 = z; double r24896 = r24894 + r24895; double r24897 = tan(r24896); double r24898 = a; double r24899 = tan(r24898); double r24900 = r24897 - r24899; double r24901 = r24893 + r24900; return r24901; } double f_of(float x, float y, float z, float a) { float r24902 = x; float r24903 = y; float r24904 = tan(r24903); float r24905 = z; float r24906 = tan(r24905); float r24907 = r24904 + r24906; float r24908 = a; float r24909 = cos(r24908); float r24910 = r24907 * r24909; float r24911 = 1; float r24912 = r24904 * r24906; float r24913 = r24911 - r24912; float r24914 = sin(r24908); float r24915 = r24913 * r24914; float r24916 = r24910 - r24915; float r24917 = r24913 * r24909; float r24918 = r24916 / r24917; float r24919 = r24902 + r24918; return r24919; } double f_od(double x, double y, double z, double a) { double r24920 = x; double r24921 = y; double r24922 = tan(r24921); double r24923 = z; double r24924 = tan(r24923); double r24925 = r24922 + r24924; double r24926 = a; double r24927 = cos(r24926); double r24928 = r24925 * r24927; double r24929 = 1; double r24930 = r24922 * r24924; double r24931 = r24929 - r24930; double r24932 = sin(r24926); double r24933 = r24931 * r24932; double r24934 = r24928 - r24933; double r24935 = r24931 * r24927; double r24936 = r24934 / r24935; double r24937 = r24920 + r24936; return r24937; } 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 r24938, r24939, r24940, r24941, r24942, r24943, r24944, r24945, r24946; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); 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); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r24938, x, MPFR_RNDN); mpfr_set_d(r24939, y, MPFR_RNDN); mpfr_set_d(r24940, z, MPFR_RNDN); mpfr_add(r24941, r24939, r24940, MPFR_RNDN); mpfr_tan(r24942, r24941, MPFR_RNDN); mpfr_set_d(r24943, a, MPFR_RNDN); mpfr_tan(r24944, r24943, MPFR_RNDN); mpfr_sub(r24945, r24942, r24944, MPFR_RNDN); mpfr_add(r24946, r24938, r24945, MPFR_RNDN); return mpfr_get_d(r24946, MPFR_RNDN); } static mpfr_t r24947, r24948, r24949, r24950, r24951, r24952, r24953, r24954, r24955, r24956, r24957, r24958, r24959, r24960, r24961, r24962, r24963, r24964; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); 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_set_str(r24956, "1", 10, MPFR_RNDN); mpfr_init(r24957); mpfr_init(r24958); mpfr_init(r24959); mpfr_init(r24960); mpfr_init(r24961); mpfr_init(r24962); mpfr_init(r24963); mpfr_init(r24964); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r24947, x, MPFR_RNDN); mpfr_set_d(r24948, y, MPFR_RNDN); mpfr_tan(r24949, r24948, MPFR_RNDN); mpfr_set_d(r24950, z, MPFR_RNDN); mpfr_tan(r24951, r24950, MPFR_RNDN); mpfr_add(r24952, r24949, r24951, MPFR_RNDN); mpfr_set_d(r24953, a, MPFR_RNDN); mpfr_cos(r24954, r24953, MPFR_RNDN); mpfr_mul(r24955, r24952, r24954, MPFR_RNDN); ; mpfr_mul(r24957, r24949, r24951, MPFR_RNDN); mpfr_sub(r24958, r24956, r24957, MPFR_RNDN); mpfr_sin(r24959, r24953, MPFR_RNDN); mpfr_mul(r24960, r24958, r24959, MPFR_RNDN); mpfr_sub(r24961, r24955, r24960, MPFR_RNDN); mpfr_mul(r24962, r24958, r24954, MPFR_RNDN); mpfr_div(r24963, r24961, r24962, MPFR_RNDN); mpfr_add(r24964, r24947, r24963, MPFR_RNDN); return mpfr_get_d(r24964, MPFR_RNDN); } static mpfr_t r24965, r24966, r24967, r24968, r24969, r24970, r24971, r24972, r24973, r24974, r24975, r24976, r24977, r24978, r24979, r24980, r24981, r24982; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); 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_set_str(r24974, "1", 10, MPFR_RNDN); mpfr_init(r24975); mpfr_init(r24976); mpfr_init(r24977); mpfr_init(r24978); mpfr_init(r24979); mpfr_init(r24980); mpfr_init(r24981); mpfr_init(r24982); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r24965, x, MPFR_RNDN); mpfr_set_d(r24966, y, MPFR_RNDN); mpfr_tan(r24967, r24966, MPFR_RNDN); mpfr_set_d(r24968, z, MPFR_RNDN); mpfr_tan(r24969, r24968, MPFR_RNDN); mpfr_add(r24970, r24967, r24969, MPFR_RNDN); mpfr_set_d(r24971, a, MPFR_RNDN); mpfr_cos(r24972, r24971, MPFR_RNDN); mpfr_mul(r24973, r24970, r24972, MPFR_RNDN); ; mpfr_mul(r24975, r24967, r24969, MPFR_RNDN); mpfr_sub(r24976, r24974, r24975, MPFR_RNDN); mpfr_sin(r24977, r24971, MPFR_RNDN); mpfr_mul(r24978, r24976, r24977, MPFR_RNDN); mpfr_sub(r24979, r24973, r24978, MPFR_RNDN); mpfr_mul(r24980, r24976, r24972, MPFR_RNDN); mpfr_div(r24981, r24979, r24980, MPFR_RNDN); mpfr_add(r24982, r24965, r24981, MPFR_RNDN); return mpfr_get_d(r24982, MPFR_RNDN); }