#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 r24789 = N;
        float r24790 = 1;
        float r24791 = r24789 + r24790;
        float r24792 = atan(r24791);
        float r24793 = atan(r24789);
        float r24794 = r24792 - r24793;
        return r24794;
}

double f_id(double N) {
        double r24795 = N;
        double r24796 = 1;
        double r24797 = r24795 + r24796;
        double r24798 = atan(r24797);
        double r24799 = atan(r24795);
        double r24800 = r24798 - r24799;
        return r24800;
}


double f_of(float N) {
        float r24801 = 1;
        float r24802 = N;
        float r24803 = r24802 + r24801;
        float r24804 = fma(r24803, r24802, r24801);
        float r24805 = atan2(r24801, r24804);
        return r24805;
}

double f_od(double N) {
        double r24806 = 1;
        double r24807 = N;
        double r24808 = r24807 + r24806;
        double r24809 = fma(r24808, r24807, r24806);
        double r24810 = atan2(r24806, r24809);
        return r24810;
}

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 r24811, r24812, r24813, r24814, r24815, r24816;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r24811);
        mpfr_init_set_str(r24812, "1", 10, MPFR_RNDN);
        mpfr_init(r24813);
        mpfr_init(r24814);
        mpfr_init(r24815);
        mpfr_init(r24816);
}

double f_im(double N) {
        mpfr_set_d(r24811, N, MPFR_RNDN);
        ;
        mpfr_add(r24813, r24811, r24812, MPFR_RNDN);
        mpfr_atan(r24814, r24813, MPFR_RNDN);
        mpfr_atan(r24815, r24811, MPFR_RNDN);
        mpfr_sub(r24816, r24814, r24815, MPFR_RNDN);
        return mpfr_get_d(r24816, MPFR_RNDN);
}

static mpfr_t r24817, r24818, r24819, r24820, r24821;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r24817, "1", 10, MPFR_RNDN);
        mpfr_init(r24818);
        mpfr_init(r24819);
        mpfr_init(r24820);
        mpfr_init(r24821);
}

double f_fm(double N) {
        ;
        mpfr_set_d(r24818, N, MPFR_RNDN);
        mpfr_add(r24819, r24818, r24817, MPFR_RNDN);
        mpfr_fma(r24820, r24819, r24818, r24817, MPFR_RNDN);
        mpfr_atan2(r24821, r24817, r24820, MPFR_RNDN);
        return mpfr_get_d(r24821, MPFR_RNDN);
}

static mpfr_t r24822, r24823, r24824, r24825, r24826;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r24822, "1", 10, MPFR_RNDN);
        mpfr_init(r24823);
        mpfr_init(r24824);
        mpfr_init(r24825);
        mpfr_init(r24826);
}

double f_dm(double N) {
        ;
        mpfr_set_d(r24823, N, MPFR_RNDN);
        mpfr_add(r24824, r24823, r24822, MPFR_RNDN);
        mpfr_fma(r24825, r24824, r24823, r24822, MPFR_RNDN);
        mpfr_atan2(r24826, r24822, r24825, MPFR_RNDN);
        return mpfr_get_d(r24826, MPFR_RNDN);
}