#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "2sin (example 3.3)";

double f_if(float x, float eps) {
        float r20851 = x;
        float r20852 = eps;
        float r20853 = r20851 + r20852;
        float r20854 = sin(r20853);
        float r20855 = sin(r20851);
        float r20856 = r20854 - r20855;
        return r20856;
}

double f_id(double x, double eps) {
        double r20857 = x;
        double r20858 = eps;
        double r20859 = r20857 + r20858;
        double r20860 = sin(r20859);
        double r20861 = sin(r20857);
        double r20862 = r20860 - r20861;
        return r20862;
}


double f_of(float x, float eps) {
        float r20863 = eps;
        float r20864 = -3.6204066776528013e-10f;
        bool r20865 = r20863 <= r20864;
        float r20866 = x;
        float r20867 = sin(r20866);
        float r20868 = cos(r20863);
        float r20869 = r20867 * r20868;
        float r20870 = cos(r20866);
        float r20871 = sin(r20863);
        float r20872 = r20870 * r20871;
        float r20873 = r20869 + r20872;
        float r20874 = r20873 - r20867;
        float r20875 = 6.909575705248122e-21f;
        bool r20876 = r20863 <= r20875;
        float r20877 = 2.0f;
        float r20878 = r20863 / r20877;
        float r20879 = sin(r20878);
        float r20880 = r20866 + r20866;
        float r20881 = r20863 + r20880;
        float r20882 = r20881 / r20877;
        float r20883 = cos(r20882);
        float r20884 = r20879 * r20883;
        float r20885 = r20877 * r20884;
        float r20886 = r20872 - r20867;
        float r20887 = r20869 + r20886;
        float r20888 = r20876 ? r20885 : r20887;
        float r20889 = r20865 ? r20874 : r20888;
        return r20889;
}

double f_od(double x, double eps) {
        double r20890 = eps;
        double r20891 = -3.6204066776528013e-10;
        bool r20892 = r20890 <= r20891;
        double r20893 = x;
        double r20894 = sin(r20893);
        double r20895 = cos(r20890);
        double r20896 = r20894 * r20895;
        double r20897 = cos(r20893);
        double r20898 = sin(r20890);
        double r20899 = r20897 * r20898;
        double r20900 = r20896 + r20899;
        double r20901 = r20900 - r20894;
        double r20902 = 6.909575705248122e-21;
        bool r20903 = r20890 <= r20902;
        double r20904 = 2.0;
        double r20905 = r20890 / r20904;
        double r20906 = sin(r20905);
        double r20907 = r20893 + r20893;
        double r20908 = r20890 + r20907;
        double r20909 = r20908 / r20904;
        double r20910 = cos(r20909);
        double r20911 = r20906 * r20910;
        double r20912 = r20904 * r20911;
        double r20913 = r20899 - r20894;
        double r20914 = r20896 + r20913;
        double r20915 = r20903 ? r20912 : r20914;
        double r20916 = r20892 ? r20901 : r20915;
        return r20916;
}

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 r20917, r20918, r20919, r20920, r20921, r20922;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init(r20917);
        mpfr_init(r20918);
        mpfr_init(r20919);
        mpfr_init(r20920);
        mpfr_init(r20921);
        mpfr_init(r20922);
}

double f_im(double x, double eps) {
        mpfr_set_d(r20917, x, MPFR_RNDN);
        mpfr_set_d(r20918, eps, MPFR_RNDN);
        mpfr_add(r20919, r20917, r20918, MPFR_RNDN);
        mpfr_sin(r20920, r20919, MPFR_RNDN);
        mpfr_sin(r20921, r20917, MPFR_RNDN);
        mpfr_sub(r20922, r20920, r20921, MPFR_RNDN);
        return mpfr_get_d(r20922, MPFR_RNDN);
}

static mpfr_t r20923, r20924, r20925, r20926, r20927, r20928, r20929, r20930, r20931, r20932, r20933, r20934, r20935, r20936, r20937, r20938, r20939, r20940, r20941, r20942, r20943, r20944, r20945, r20946, r20947, r20948, r20949;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r20923);
        mpfr_init_set_str(r20924, "-3.6204066776528013e-10", 10, MPFR_RNDN);
        mpfr_init(r20925);
        mpfr_init(r20926);
        mpfr_init(r20927);
        mpfr_init(r20928);
        mpfr_init(r20929);
        mpfr_init(r20930);
        mpfr_init(r20931);
        mpfr_init(r20932);
        mpfr_init(r20933);
        mpfr_init(r20934);
        mpfr_init_set_str(r20935, "6.909575705248122e-21", 10, MPFR_RNDN);
        mpfr_init(r20936);
        mpfr_init_set_str(r20937, "2", 10, MPFR_RNDN);
        mpfr_init(r20938);
        mpfr_init(r20939);
        mpfr_init(r20940);
        mpfr_init(r20941);
        mpfr_init(r20942);
        mpfr_init(r20943);
        mpfr_init(r20944);
        mpfr_init(r20945);
        mpfr_init(r20946);
        mpfr_init(r20947);
        mpfr_init(r20948);
        mpfr_init(r20949);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r20923, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r20925, mpfr_cmp(r20923, r20924) <= 0, MPFR_RNDN);
        mpfr_set_d(r20926, x, MPFR_RNDN);
        mpfr_sin(r20927, r20926, MPFR_RNDN);
        mpfr_cos(r20928, r20923, MPFR_RNDN);
        mpfr_mul(r20929, r20927, r20928, MPFR_RNDN);
        mpfr_cos(r20930, r20926, MPFR_RNDN);
        mpfr_sin(r20931, r20923, MPFR_RNDN);
        mpfr_mul(r20932, r20930, r20931, MPFR_RNDN);
        mpfr_add(r20933, r20929, r20932, MPFR_RNDN);
        mpfr_sub(r20934, r20933, r20927, MPFR_RNDN);
        ;
        mpfr_set_si(r20936, mpfr_cmp(r20923, r20935) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r20938, r20923, r20937, MPFR_RNDN);
        mpfr_sin(r20939, r20938, MPFR_RNDN);
        mpfr_add(r20940, r20926, r20926, MPFR_RNDN);
        mpfr_add(r20941, r20923, r20940, MPFR_RNDN);
        mpfr_div(r20942, r20941, r20937, MPFR_RNDN);
        mpfr_cos(r20943, r20942, MPFR_RNDN);
        mpfr_mul(r20944, r20939, r20943, MPFR_RNDN);
        mpfr_mul(r20945, r20937, r20944, MPFR_RNDN);
        mpfr_sub(r20946, r20932, r20927, MPFR_RNDN);
        mpfr_add(r20947, r20929, r20946, MPFR_RNDN);
        if (mpfr_get_si(r20936, MPFR_RNDN)) { mpfr_set(r20948, r20945, MPFR_RNDN); } else { mpfr_set(r20948, r20947, MPFR_RNDN); };
        if (mpfr_get_si(r20925, MPFR_RNDN)) { mpfr_set(r20949, r20934, MPFR_RNDN); } else { mpfr_set(r20949, r20948, MPFR_RNDN); };
        return mpfr_get_d(r20949, MPFR_RNDN);
}

static mpfr_t r20950, r20951, r20952, r20953, r20954, r20955, r20956, r20957, r20958, r20959, r20960, r20961, r20962, r20963, r20964, r20965, r20966, r20967, r20968, r20969, r20970, r20971, r20972, r20973, r20974, r20975, r20976;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r20950);
        mpfr_init_set_str(r20951, "-3.6204066776528013e-10", 10, MPFR_RNDN);
        mpfr_init(r20952);
        mpfr_init(r20953);
        mpfr_init(r20954);
        mpfr_init(r20955);
        mpfr_init(r20956);
        mpfr_init(r20957);
        mpfr_init(r20958);
        mpfr_init(r20959);
        mpfr_init(r20960);
        mpfr_init(r20961);
        mpfr_init_set_str(r20962, "6.909575705248122e-21", 10, MPFR_RNDN);
        mpfr_init(r20963);
        mpfr_init_set_str(r20964, "2", 10, MPFR_RNDN);
        mpfr_init(r20965);
        mpfr_init(r20966);
        mpfr_init(r20967);
        mpfr_init(r20968);
        mpfr_init(r20969);
        mpfr_init(r20970);
        mpfr_init(r20971);
        mpfr_init(r20972);
        mpfr_init(r20973);
        mpfr_init(r20974);
        mpfr_init(r20975);
        mpfr_init(r20976);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r20950, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r20952, mpfr_cmp(r20950, r20951) <= 0, MPFR_RNDN);
        mpfr_set_d(r20953, x, MPFR_RNDN);
        mpfr_sin(r20954, r20953, MPFR_RNDN);
        mpfr_cos(r20955, r20950, MPFR_RNDN);
        mpfr_mul(r20956, r20954, r20955, MPFR_RNDN);
        mpfr_cos(r20957, r20953, MPFR_RNDN);
        mpfr_sin(r20958, r20950, MPFR_RNDN);
        mpfr_mul(r20959, r20957, r20958, MPFR_RNDN);
        mpfr_add(r20960, r20956, r20959, MPFR_RNDN);
        mpfr_sub(r20961, r20960, r20954, MPFR_RNDN);
        ;
        mpfr_set_si(r20963, mpfr_cmp(r20950, r20962) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r20965, r20950, r20964, MPFR_RNDN);
        mpfr_sin(r20966, r20965, MPFR_RNDN);
        mpfr_add(r20967, r20953, r20953, MPFR_RNDN);
        mpfr_add(r20968, r20950, r20967, MPFR_RNDN);
        mpfr_div(r20969, r20968, r20964, MPFR_RNDN);
        mpfr_cos(r20970, r20969, MPFR_RNDN);
        mpfr_mul(r20971, r20966, r20970, MPFR_RNDN);
        mpfr_mul(r20972, r20964, r20971, MPFR_RNDN);
        mpfr_sub(r20973, r20959, r20954, MPFR_RNDN);
        mpfr_add(r20974, r20956, r20973, MPFR_RNDN);
        if (mpfr_get_si(r20963, MPFR_RNDN)) { mpfr_set(r20975, r20972, MPFR_RNDN); } else { mpfr_set(r20975, r20974, MPFR_RNDN); };
        if (mpfr_get_si(r20952, MPFR_RNDN)) { mpfr_set(r20976, r20961, MPFR_RNDN); } else { mpfr_set(r20976, r20975, MPFR_RNDN); };
        return mpfr_get_d(r20976, MPFR_RNDN);
}