#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 r20977 = x;
        float r20978 = eps;
        float r20979 = r20977 + r20978;
        float r20980 = sin(r20979);
        float r20981 = sin(r20977);
        float r20982 = r20980 - r20981;
        return r20982;
}

double f_id(double x, double eps) {
        double r20983 = x;
        double r20984 = eps;
        double r20985 = r20983 + r20984;
        double r20986 = sin(r20985);
        double r20987 = sin(r20983);
        double r20988 = r20986 - r20987;
        return r20988;
}


double f_of(float x, float eps) {
        float r20989 = eps;
        float r20990 = -1.3331382704773037e-08f;
        bool r20991 = r20989 <= r20990;
        float r20992 = x;
        float r20993 = sin(r20992);
        float r20994 = cos(r20989);
        float r20995 = r20993 * r20994;
        float r20996 = cos(r20992);
        float r20997 = sin(r20989);
        float r20998 = r20996 * r20997;
        float r20999 = r20995 + r20998;
        float r21000 = r20999 - r20993;
        float r21001 = 8.172994376936087f;
        bool r21002 = r20989 <= r21001;
        float r21003 = 2.0f;
        float r21004 = r20989 / r21003;
        float r21005 = sin(r21004);
        float r21006 = r20992 + r20992;
        float r21007 = r20989 + r21006;
        float r21008 = r21007 / r21003;
        float r21009 = cos(r21008);
        float r21010 = r21005 * r21009;
        float r21011 = r21003 * r21010;
        float r21012 = r21002 ? r21011 : r21000;
        float r21013 = r20991 ? r21000 : r21012;
        return r21013;
}

double f_od(double x, double eps) {
        double r21014 = eps;
        double r21015 = -1.3331382704773037e-08;
        bool r21016 = r21014 <= r21015;
        double r21017 = x;
        double r21018 = sin(r21017);
        double r21019 = cos(r21014);
        double r21020 = r21018 * r21019;
        double r21021 = cos(r21017);
        double r21022 = sin(r21014);
        double r21023 = r21021 * r21022;
        double r21024 = r21020 + r21023;
        double r21025 = r21024 - r21018;
        double r21026 = 8.172994376936087;
        bool r21027 = r21014 <= r21026;
        double r21028 = 2.0;
        double r21029 = r21014 / r21028;
        double r21030 = sin(r21029);
        double r21031 = r21017 + r21017;
        double r21032 = r21014 + r21031;
        double r21033 = r21032 / r21028;
        double r21034 = cos(r21033);
        double r21035 = r21030 * r21034;
        double r21036 = r21028 * r21035;
        double r21037 = r21027 ? r21036 : r21025;
        double r21038 = r21016 ? r21025 : r21037;
        return r21038;
}

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 r21039, r21040, r21041, r21042, r21043, r21044;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2192);
        mpfr_init(r21039);
        mpfr_init(r21040);
        mpfr_init(r21041);
        mpfr_init(r21042);
        mpfr_init(r21043);
        mpfr_init(r21044);
}

double f_im(double x, double eps) {
        mpfr_set_d(r21039, x, MPFR_RNDN);
        mpfr_set_d(r21040, eps, MPFR_RNDN);
        mpfr_add(r21041, r21039, r21040, MPFR_RNDN);
        mpfr_sin(r21042, r21041, MPFR_RNDN);
        mpfr_sin(r21043, r21039, MPFR_RNDN);
        mpfr_sub(r21044, r21042, r21043, MPFR_RNDN);
        return mpfr_get_d(r21044, MPFR_RNDN);
}

static mpfr_t r21045, r21046, r21047, r21048, r21049, r21050, r21051, r21052, r21053, r21054, r21055, r21056, r21057, r21058, r21059, r21060, r21061, r21062, r21063, r21064, r21065, r21066, r21067, r21068, r21069;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2192);
        mpfr_init(r21045);
        mpfr_init_set_str(r21046, "-1.3331382704773037e-08", 10, MPFR_RNDN);
        mpfr_init(r21047);
        mpfr_init(r21048);
        mpfr_init(r21049);
        mpfr_init(r21050);
        mpfr_init(r21051);
        mpfr_init(r21052);
        mpfr_init(r21053);
        mpfr_init(r21054);
        mpfr_init(r21055);
        mpfr_init(r21056);
        mpfr_init_set_str(r21057, "8.172994376936087", 10, MPFR_RNDN);
        mpfr_init(r21058);
        mpfr_init_set_str(r21059, "2", 10, MPFR_RNDN);
        mpfr_init(r21060);
        mpfr_init(r21061);
        mpfr_init(r21062);
        mpfr_init(r21063);
        mpfr_init(r21064);
        mpfr_init(r21065);
        mpfr_init(r21066);
        mpfr_init(r21067);
        mpfr_init(r21068);
        mpfr_init(r21069);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r21045, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r21047, mpfr_cmp(r21045, r21046) <= 0, MPFR_RNDN);
        mpfr_set_d(r21048, x, MPFR_RNDN);
        mpfr_sin(r21049, r21048, MPFR_RNDN);
        mpfr_cos(r21050, r21045, MPFR_RNDN);
        mpfr_mul(r21051, r21049, r21050, MPFR_RNDN);
        mpfr_cos(r21052, r21048, MPFR_RNDN);
        mpfr_sin(r21053, r21045, MPFR_RNDN);
        mpfr_mul(r21054, r21052, r21053, MPFR_RNDN);
        mpfr_add(r21055, r21051, r21054, MPFR_RNDN);
        mpfr_sub(r21056, r21055, r21049, MPFR_RNDN);
        ;
        mpfr_set_si(r21058, mpfr_cmp(r21045, r21057) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r21060, r21045, r21059, MPFR_RNDN);
        mpfr_sin(r21061, r21060, MPFR_RNDN);
        mpfr_add(r21062, r21048, r21048, MPFR_RNDN);
        mpfr_add(r21063, r21045, r21062, MPFR_RNDN);
        mpfr_div(r21064, r21063, r21059, MPFR_RNDN);
        mpfr_cos(r21065, r21064, MPFR_RNDN);
        mpfr_mul(r21066, r21061, r21065, MPFR_RNDN);
        mpfr_mul(r21067, r21059, r21066, MPFR_RNDN);
        if (mpfr_get_si(r21058, MPFR_RNDN)) { mpfr_set(r21068, r21067, MPFR_RNDN); } else { mpfr_set(r21068, r21056, MPFR_RNDN); };
        if (mpfr_get_si(r21047, MPFR_RNDN)) { mpfr_set(r21069, r21056, MPFR_RNDN); } else { mpfr_set(r21069, r21068, MPFR_RNDN); };
        return mpfr_get_d(r21069, MPFR_RNDN);
}

static mpfr_t r21070, r21071, r21072, r21073, r21074, r21075, r21076, r21077, r21078, r21079, r21080, r21081, r21082, r21083, r21084, r21085, r21086, r21087, r21088, r21089, r21090, r21091, r21092, r21093, r21094;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2192);
        mpfr_init(r21070);
        mpfr_init_set_str(r21071, "-1.3331382704773037e-08", 10, MPFR_RNDN);
        mpfr_init(r21072);
        mpfr_init(r21073);
        mpfr_init(r21074);
        mpfr_init(r21075);
        mpfr_init(r21076);
        mpfr_init(r21077);
        mpfr_init(r21078);
        mpfr_init(r21079);
        mpfr_init(r21080);
        mpfr_init(r21081);
        mpfr_init_set_str(r21082, "8.172994376936087", 10, MPFR_RNDN);
        mpfr_init(r21083);
        mpfr_init_set_str(r21084, "2", 10, MPFR_RNDN);
        mpfr_init(r21085);
        mpfr_init(r21086);
        mpfr_init(r21087);
        mpfr_init(r21088);
        mpfr_init(r21089);
        mpfr_init(r21090);
        mpfr_init(r21091);
        mpfr_init(r21092);
        mpfr_init(r21093);
        mpfr_init(r21094);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r21070, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r21072, mpfr_cmp(r21070, r21071) <= 0, MPFR_RNDN);
        mpfr_set_d(r21073, x, MPFR_RNDN);
        mpfr_sin(r21074, r21073, MPFR_RNDN);
        mpfr_cos(r21075, r21070, MPFR_RNDN);
        mpfr_mul(r21076, r21074, r21075, MPFR_RNDN);
        mpfr_cos(r21077, r21073, MPFR_RNDN);
        mpfr_sin(r21078, r21070, MPFR_RNDN);
        mpfr_mul(r21079, r21077, r21078, MPFR_RNDN);
        mpfr_add(r21080, r21076, r21079, MPFR_RNDN);
        mpfr_sub(r21081, r21080, r21074, MPFR_RNDN);
        ;
        mpfr_set_si(r21083, mpfr_cmp(r21070, r21082) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r21085, r21070, r21084, MPFR_RNDN);
        mpfr_sin(r21086, r21085, MPFR_RNDN);
        mpfr_add(r21087, r21073, r21073, MPFR_RNDN);
        mpfr_add(r21088, r21070, r21087, MPFR_RNDN);
        mpfr_div(r21089, r21088, r21084, MPFR_RNDN);
        mpfr_cos(r21090, r21089, MPFR_RNDN);
        mpfr_mul(r21091, r21086, r21090, MPFR_RNDN);
        mpfr_mul(r21092, r21084, r21091, MPFR_RNDN);
        if (mpfr_get_si(r21083, MPFR_RNDN)) { mpfr_set(r21093, r21092, MPFR_RNDN); } else { mpfr_set(r21093, r21081, MPFR_RNDN); };
        if (mpfr_get_si(r21072, MPFR_RNDN)) { mpfr_set(r21094, r21081, MPFR_RNDN); } else { mpfr_set(r21094, r21093, MPFR_RNDN); };
        return mpfr_get_d(r21094, MPFR_RNDN);
}