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

char *name = "2atan (example 3.5)";

double f_if(float N) {
        float r24880 = N;
        float r24881 = 1;
        float r24882 = r24880 + r24881;
        float r24883 = atan(r24882);
        float r24884 = atan(r24880);
        float r24885 = r24883 - r24884;
        return r24885;
}

double f_id(double N) {
        double r24886 = N;
        double r24887 = 1;
        double r24888 = r24886 + r24887;
        double r24889 = atan(r24888);
        double r24890 = atan(r24886);
        double r24891 = r24889 - r24890;
        return r24891;
}


double f_of(float N) {
        float r24892 = 1;
        float r24893 = N;
        float r24894 = r24893 + r24892;
        float r24895 = r24894 * r24893;
        float r24896 = r24892 + r24895;
        float r24897 = atan2(r24892, r24896);
        return r24897;
}

double f_od(double N) {
        double r24898 = 1;
        double r24899 = N;
        double r24900 = r24899 + r24898;
        double r24901 = r24900 * r24899;
        double r24902 = r24898 + r24901;
        double r24903 = atan2(r24898, r24902);
        return r24903;
}

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 r24904, r24905, r24906, r24907, r24908, r24909;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r24904);
        mpfr_init_set_str(r24905, "1", 10, MPFR_RNDN);
        mpfr_init(r24906);
        mpfr_init(r24907);
        mpfr_init(r24908);
        mpfr_init(r24909);
}

double f_im(double N) {
        mpfr_set_d(r24904, N, MPFR_RNDN);
        ;
        mpfr_add(r24906, r24904, r24905, MPFR_RNDN);
        mpfr_atan(r24907, r24906, MPFR_RNDN);
        mpfr_atan(r24908, r24904, MPFR_RNDN);
        mpfr_sub(r24909, r24907, r24908, MPFR_RNDN);
        return mpfr_get_d(r24909, MPFR_RNDN);
}

static mpfr_t r24910, r24911, r24912, r24913, r24914, r24915;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r24910, "1", 10, MPFR_RNDN);
        mpfr_init(r24911);
        mpfr_init(r24912);
        mpfr_init(r24913);
        mpfr_init(r24914);
        mpfr_init(r24915);
}

double f_fm(double N) {
        ;
        mpfr_set_d(r24911, N, MPFR_RNDN);
        mpfr_add(r24912, r24911, r24910, MPFR_RNDN);
        mpfr_mul(r24913, r24912, r24911, MPFR_RNDN);
        mpfr_add(r24914, r24910, r24913, MPFR_RNDN);
        mpfr_atan2(r24915, r24910, r24914, MPFR_RNDN);
        return mpfr_get_d(r24915, MPFR_RNDN);
}

static mpfr_t r24916, r24917, r24918, r24919, r24920, r24921;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r24916, "1", 10, MPFR_RNDN);
        mpfr_init(r24917);
        mpfr_init(r24918);
        mpfr_init(r24919);
        mpfr_init(r24920);
        mpfr_init(r24921);
}

double f_dm(double N) {
        ;
        mpfr_set_d(r24917, N, MPFR_RNDN);
        mpfr_add(r24918, r24917, r24916, MPFR_RNDN);
        mpfr_mul(r24919, r24918, r24917, MPFR_RNDN);
        mpfr_add(r24920, r24916, r24919, MPFR_RNDN);
        mpfr_atan2(r24921, r24916, r24920, MPFR_RNDN);
        return mpfr_get_d(r24921, MPFR_RNDN);
}