#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 r25002 = x;
        float r25003 = eps;
        float r25004 = r25002 + r25003;
        float r25005 = sin(r25004);
        float r25006 = sin(r25002);
        float r25007 = r25005 - r25006;
        return r25007;
}

double f_id(double x, double eps) {
        double r25008 = x;
        double r25009 = eps;
        double r25010 = r25008 + r25009;
        double r25011 = sin(r25010);
        double r25012 = sin(r25008);
        double r25013 = r25011 - r25012;
        return r25013;
}


double f_of(float x, float eps) {
        float r25014 = x;
        float r25015 = sin(r25014);
        float r25016 = eps;
        float r25017 = cos(r25016);
        float r25018 = r25015 * r25017;
        float r25019 = cos(r25014);
        float r25020 = sin(r25016);
        float r25021 = r25019 * r25020;
        float r25022 = r25021 - r25015;
        float r25023 = r25018 + r25022;
        float r25024 = -0.009223036685637842;
        bool r25025 = r25023 <= r25024;
        float r25026 = 0.0024644411836328532;
        bool r25027 = r25023 <= r25026;
        float r25028 = 2;
        float r25029 = r25016 / r25028;
        float r25030 = sin(r25029);
        float r25031 = r25016 + r25014;
        float r25032 = r25014 + r25031;
        float r25033 = r25032 / r25028;
        float r25034 = cos(r25033);
        float r25035 = r25030 * r25034;
        float r25036 = r25028 * r25035;
        float r25037 = r25027 ? r25036 : r25023;
        float r25038 = r25025 ? r25023 : r25037;
        return r25038;
}

double f_od(double x, double eps) {
        double r25039 = x;
        double r25040 = sin(r25039);
        double r25041 = eps;
        double r25042 = cos(r25041);
        double r25043 = r25040 * r25042;
        double r25044 = cos(r25039);
        double r25045 = sin(r25041);
        double r25046 = r25044 * r25045;
        double r25047 = r25046 - r25040;
        double r25048 = r25043 + r25047;
        double r25049 = -0.009223036685637842;
        bool r25050 = r25048 <= r25049;
        double r25051 = 0.0024644411836328532;
        bool r25052 = r25048 <= r25051;
        double r25053 = 2;
        double r25054 = r25041 / r25053;
        double r25055 = sin(r25054);
        double r25056 = r25041 + r25039;
        double r25057 = r25039 + r25056;
        double r25058 = r25057 / r25053;
        double r25059 = cos(r25058);
        double r25060 = r25055 * r25059;
        double r25061 = r25053 * r25060;
        double r25062 = r25052 ? r25061 : r25048;
        double r25063 = r25050 ? r25048 : r25062;
        return r25063;
}

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 r25064, r25065, r25066, r25067, r25068, r25069;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r25064);
        mpfr_init(r25065);
        mpfr_init(r25066);
        mpfr_init(r25067);
        mpfr_init(r25068);
        mpfr_init(r25069);
}

double f_im(double x, double eps) {
        mpfr_set_d(r25064, x, MPFR_RNDN);
        mpfr_set_d(r25065, eps, MPFR_RNDN);
        mpfr_add(r25066, r25064, r25065, MPFR_RNDN);
        mpfr_sin(r25067, r25066, MPFR_RNDN);
        mpfr_sin(r25068, r25064, MPFR_RNDN);
        mpfr_sub(r25069, r25067, r25068, MPFR_RNDN);
        return mpfr_get_d(r25069, MPFR_RNDN);
}

static mpfr_t r25070, r25071, r25072, r25073, r25074, r25075, r25076, r25077, r25078, r25079, r25080, r25081, r25082, r25083, r25084, r25085, r25086, r25087, r25088, r25089, r25090, r25091, r25092, r25093, r25094;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r25070);
        mpfr_init(r25071);
        mpfr_init(r25072);
        mpfr_init(r25073);
        mpfr_init(r25074);
        mpfr_init(r25075);
        mpfr_init(r25076);
        mpfr_init(r25077);
        mpfr_init(r25078);
        mpfr_init(r25079);
        mpfr_init_set_str(r25080, "-0.009223036685637842", 10, MPFR_RNDN);
        mpfr_init(r25081);
        mpfr_init_set_str(r25082, "0.0024644411836328532", 10, MPFR_RNDN);
        mpfr_init(r25083);
        mpfr_init_set_str(r25084, "2", 10, MPFR_RNDN);
        mpfr_init(r25085);
        mpfr_init(r25086);
        mpfr_init(r25087);
        mpfr_init(r25088);
        mpfr_init(r25089);
        mpfr_init(r25090);
        mpfr_init(r25091);
        mpfr_init(r25092);
        mpfr_init(r25093);
        mpfr_init(r25094);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r25070, x, MPFR_RNDN);
        mpfr_sin(r25071, r25070, MPFR_RNDN);
        mpfr_set_d(r25072, eps, MPFR_RNDN);
        mpfr_cos(r25073, r25072, MPFR_RNDN);
        mpfr_mul(r25074, r25071, r25073, MPFR_RNDN);
        mpfr_cos(r25075, r25070, MPFR_RNDN);
        mpfr_sin(r25076, r25072, MPFR_RNDN);
        mpfr_mul(r25077, r25075, r25076, MPFR_RNDN);
        mpfr_sub(r25078, r25077, r25071, MPFR_RNDN);
        mpfr_add(r25079, r25074, r25078, MPFR_RNDN);
        ;
        mpfr_set_si(r25081, mpfr_cmp(r25079, r25080) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r25083, mpfr_cmp(r25079, r25082) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r25085, r25072, r25084, MPFR_RNDN);
        mpfr_sin(r25086, r25085, MPFR_RNDN);
        mpfr_add(r25087, r25072, r25070, MPFR_RNDN);
        mpfr_add(r25088, r25070, r25087, MPFR_RNDN);
        mpfr_div(r25089, r25088, r25084, MPFR_RNDN);
        mpfr_cos(r25090, r25089, MPFR_RNDN);
        mpfr_mul(r25091, r25086, r25090, MPFR_RNDN);
        mpfr_mul(r25092, r25084, r25091, MPFR_RNDN);
        if (mpfr_get_si(r25083, MPFR_RNDN)) { mpfr_set(r25093, r25092, MPFR_RNDN); } else { mpfr_set(r25093, r25079, MPFR_RNDN); };
        if (mpfr_get_si(r25081, MPFR_RNDN)) { mpfr_set(r25094, r25079, MPFR_RNDN); } else { mpfr_set(r25094, r25093, MPFR_RNDN); };
        return mpfr_get_d(r25094, MPFR_RNDN);
}

static mpfr_t r25095, r25096, r25097, r25098, r25099, r25100, r25101, r25102, r25103, r25104, r25105, r25106, r25107, r25108, r25109, r25110, r25111, r25112, r25113, r25114, r25115, r25116, r25117, r25118, r25119;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r25095);
        mpfr_init(r25096);
        mpfr_init(r25097);
        mpfr_init(r25098);
        mpfr_init(r25099);
        mpfr_init(r25100);
        mpfr_init(r25101);
        mpfr_init(r25102);
        mpfr_init(r25103);
        mpfr_init(r25104);
        mpfr_init_set_str(r25105, "-0.009223036685637842", 10, MPFR_RNDN);
        mpfr_init(r25106);
        mpfr_init_set_str(r25107, "0.0024644411836328532", 10, MPFR_RNDN);
        mpfr_init(r25108);
        mpfr_init_set_str(r25109, "2", 10, MPFR_RNDN);
        mpfr_init(r25110);
        mpfr_init(r25111);
        mpfr_init(r25112);
        mpfr_init(r25113);
        mpfr_init(r25114);
        mpfr_init(r25115);
        mpfr_init(r25116);
        mpfr_init(r25117);
        mpfr_init(r25118);
        mpfr_init(r25119);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r25095, x, MPFR_RNDN);
        mpfr_sin(r25096, r25095, MPFR_RNDN);
        mpfr_set_d(r25097, eps, MPFR_RNDN);
        mpfr_cos(r25098, r25097, MPFR_RNDN);
        mpfr_mul(r25099, r25096, r25098, MPFR_RNDN);
        mpfr_cos(r25100, r25095, MPFR_RNDN);
        mpfr_sin(r25101, r25097, MPFR_RNDN);
        mpfr_mul(r25102, r25100, r25101, MPFR_RNDN);
        mpfr_sub(r25103, r25102, r25096, MPFR_RNDN);
        mpfr_add(r25104, r25099, r25103, MPFR_RNDN);
        ;
        mpfr_set_si(r25106, mpfr_cmp(r25104, r25105) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r25108, mpfr_cmp(r25104, r25107) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r25110, r25097, r25109, MPFR_RNDN);
        mpfr_sin(r25111, r25110, MPFR_RNDN);
        mpfr_add(r25112, r25097, r25095, MPFR_RNDN);
        mpfr_add(r25113, r25095, r25112, MPFR_RNDN);
        mpfr_div(r25114, r25113, r25109, MPFR_RNDN);
        mpfr_cos(r25115, r25114, MPFR_RNDN);
        mpfr_mul(r25116, r25111, r25115, MPFR_RNDN);
        mpfr_mul(r25117, r25109, r25116, MPFR_RNDN);
        if (mpfr_get_si(r25108, MPFR_RNDN)) { mpfr_set(r25118, r25117, MPFR_RNDN); } else { mpfr_set(r25118, r25104, MPFR_RNDN); };
        if (mpfr_get_si(r25106, MPFR_RNDN)) { mpfr_set(r25119, r25104, MPFR_RNDN); } else { mpfr_set(r25119, r25118, MPFR_RNDN); };
        return mpfr_get_d(r25119, MPFR_RNDN);
}