#include #include #include #include #include char *name = "r*sin(b)/cos(a+b), A"; double f_if(float r, float a, float b) { float r14863 = r; float r14864 = b; float r14865 = sin(r14864); float r14866 = r14863 * r14865; float r14867 = a; float r14868 = r14867 + r14864; float r14869 = cos(r14868); float r14870 = r14866 / r14869; return r14870; } double f_id(double r, double a, double b) { double r14871 = r; double r14872 = b; double r14873 = sin(r14872); double r14874 = r14871 * r14873; double r14875 = a; double r14876 = r14875 + r14872; double r14877 = cos(r14876); double r14878 = r14874 / r14877; return r14878; } double f_of(float r, float a, float b) { float r14879 = r; float r14880 = 1.0f; float r14881 = r14879 / r14880; float r14882 = b; float r14883 = sin(r14882); float r14884 = a; float r14885 = cos(r14884); float r14886 = cos(r14882); float r14887 = r14885 * r14886; float r14888 = sin(r14884); float r14889 = r14888 * (r14888 * r14888); float r14890 = r14883 * (r14883 * r14883); float r14891 = r14889 * r14890; float r14892 = cbrt(r14891); float r14893 = r14887 - r14892; float r14894 = r14883 / r14893; float r14895 = r14881 * r14894; return r14895; } double f_od(double r, double a, double b) { double r14896 = r; double r14897 = 1.0; double r14898 = r14896 / r14897; double r14899 = b; double r14900 = sin(r14899); double r14901 = a; double r14902 = cos(r14901); double r14903 = cos(r14899); double r14904 = r14902 * r14903; double r14905 = sin(r14901); double r14906 = r14905 * (r14905 * r14905); double r14907 = r14900 * (r14900 * r14900); double r14908 = r14906 * r14907; double r14909 = cbrt(r14908); double r14910 = r14904 - r14909; double r14911 = r14900 / r14910; double r14912 = r14898 * r14911; return r14912; } 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 r14913, r14914, r14915, r14916, r14917, r14918, r14919, r14920; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r14913); mpfr_init(r14914); mpfr_init(r14915); mpfr_init(r14916); mpfr_init(r14917); mpfr_init(r14918); mpfr_init(r14919); mpfr_init(r14920); } double f_im(double r, double a, double b) { mpfr_set_d(r14913, r, MPFR_RNDN); mpfr_set_d(r14914, b, MPFR_RNDN); mpfr_sin(r14915, r14914, MPFR_RNDN); mpfr_mul(r14916, r14913, r14915, MPFR_RNDN); mpfr_set_d(r14917, a, MPFR_RNDN); mpfr_add(r14918, r14917, r14914, MPFR_RNDN); mpfr_cos(r14919, r14918, MPFR_RNDN); mpfr_div(r14920, r14916, r14919, MPFR_RNDN); return mpfr_get_d(r14920, MPFR_RNDN); } static mpfr_t r14921, r14922, r14923, r14924, r14925, r14926, r14927, r14928, r14929, r14930, r14931, r14932, r14933, r14934, r14935, r14936, r14937; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r14921); mpfr_init_set_str(r14922, "1", 10, MPFR_RNDN); mpfr_init(r14923); mpfr_init(r14924); mpfr_init(r14925); mpfr_init(r14926); mpfr_init(r14927); mpfr_init(r14928); mpfr_init(r14929); mpfr_init(r14930); mpfr_init(r14931); mpfr_init(r14932); mpfr_init(r14933); mpfr_init(r14934); mpfr_init(r14935); mpfr_init(r14936); mpfr_init(r14937); } double f_fm(double r, double a, double b) { mpfr_set_d(r14921, r, MPFR_RNDN); ; mpfr_div(r14923, r14921, r14922, MPFR_RNDN); mpfr_set_d(r14924, b, MPFR_RNDN); mpfr_sin(r14925, r14924, MPFR_RNDN); mpfr_set_d(r14926, a, MPFR_RNDN); mpfr_cos(r14927, r14926, MPFR_RNDN); mpfr_cos(r14928, r14924, MPFR_RNDN); mpfr_mul(r14929, r14927, r14928, MPFR_RNDN); mpfr_sin(r14930, r14926, MPFR_RNDN); mpfr_mul(r14931, r14930, r14930, MPFR_RNDN); mpfr_mul(r14931, r14931, r14930, MPFR_RNDN); mpfr_mul(r14932, r14925, r14925, MPFR_RNDN); mpfr_mul(r14932, r14932, r14925, MPFR_RNDN); mpfr_mul(r14933, r14931, r14932, MPFR_RNDN); mpfr_cbrt(r14934, r14933, MPFR_RNDN); mpfr_sub(r14935, r14929, r14934, MPFR_RNDN); mpfr_div(r14936, r14925, r14935, MPFR_RNDN); mpfr_mul(r14937, r14923, r14936, MPFR_RNDN); return mpfr_get_d(r14937, MPFR_RNDN); } static mpfr_t r14938, r14939, r14940, r14941, r14942, r14943, r14944, r14945, r14946, r14947, r14948, r14949, r14950, r14951, r14952, r14953, r14954; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r14938); mpfr_init_set_str(r14939, "1", 10, MPFR_RNDN); mpfr_init(r14940); mpfr_init(r14941); mpfr_init(r14942); mpfr_init(r14943); mpfr_init(r14944); mpfr_init(r14945); mpfr_init(r14946); mpfr_init(r14947); mpfr_init(r14948); mpfr_init(r14949); mpfr_init(r14950); mpfr_init(r14951); mpfr_init(r14952); mpfr_init(r14953); mpfr_init(r14954); } double f_dm(double r, double a, double b) { mpfr_set_d(r14938, r, MPFR_RNDN); ; mpfr_div(r14940, r14938, r14939, MPFR_RNDN); mpfr_set_d(r14941, b, MPFR_RNDN); mpfr_sin(r14942, r14941, MPFR_RNDN); mpfr_set_d(r14943, a, MPFR_RNDN); mpfr_cos(r14944, r14943, MPFR_RNDN); mpfr_cos(r14945, r14941, MPFR_RNDN); mpfr_mul(r14946, r14944, r14945, MPFR_RNDN); mpfr_sin(r14947, r14943, MPFR_RNDN); mpfr_mul(r14948, r14947, r14947, MPFR_RNDN); mpfr_mul(r14948, r14948, r14947, MPFR_RNDN); mpfr_mul(r14949, r14942, r14942, MPFR_RNDN); mpfr_mul(r14949, r14949, r14942, MPFR_RNDN); mpfr_mul(r14950, r14948, r14949, MPFR_RNDN); mpfr_cbrt(r14951, r14950, MPFR_RNDN); mpfr_sub(r14952, r14946, r14951, MPFR_RNDN); mpfr_div(r14953, r14942, r14952, MPFR_RNDN); mpfr_mul(r14954, r14940, r14953, MPFR_RNDN); return mpfr_get_d(r14954, MPFR_RNDN); }