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

char *name = "Octave 3.8, jcobi/4, as called";

double f_if(float i) {
        float r23287 = i;
        float r23288 = r23287 * r23287;
        float r23289 = r23288 * r23288;
        float r23290 = 2;
        float r23291 = r23290 * r23287;
        float r23292 = r23291 * r23291;
        float r23293 = r23289 / r23292;
        float r23294 = 1.0;
        float r23295 = r23292 - r23294;
        float r23296 = r23293 / r23295;
        return r23296;
}

double f_id(double i) {
        double r23297 = i;
        double r23298 = r23297 * r23297;
        double r23299 = r23298 * r23298;
        double r23300 = 2;
        double r23301 = r23300 * r23297;
        double r23302 = r23301 * r23301;
        double r23303 = r23299 / r23302;
        double r23304 = 1.0;
        double r23305 = r23302 - r23304;
        double r23306 = r23303 / r23305;
        return r23306;
}


double f_of(float i) {
        float r23307 = i;
        float r23308 = 2;
        float r23309 = r23307 / r23308;
        float r23310 = 8;
        float r23311 = r23310 * r23307;
        float r23312 = 2.0;
        float r23313 = r23312 / r23307;
        float r23314 = r23311 - r23313;
        float r23315 = r23309 / r23314;
        return r23315;
}

double f_od(double i) {
        double r23316 = i;
        double r23317 = 2;
        double r23318 = r23316 / r23317;
        double r23319 = 8;
        double r23320 = r23319 * r23316;
        double r23321 = 2.0;
        double r23322 = r23321 / r23316;
        double r23323 = r23320 - r23322;
        double r23324 = r23318 / r23323;
        return r23324;
}

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 r23325, r23326, r23327, r23328, r23329, r23330, r23331, r23332, r23333, r23334;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r23325);
        mpfr_init(r23326);
        mpfr_init(r23327);
        mpfr_init_set_str(r23328, "2", 10, MPFR_RNDN);
        mpfr_init(r23329);
        mpfr_init(r23330);
        mpfr_init(r23331);
        mpfr_init_set_str(r23332, "1.0", 10, MPFR_RNDN);
        mpfr_init(r23333);
        mpfr_init(r23334);
}

double f_im(double i) {
        mpfr_set_d(r23325, i, MPFR_RNDN);
        mpfr_mul(r23326, r23325, r23325, MPFR_RNDN);
        mpfr_mul(r23327, r23326, r23326, MPFR_RNDN);
        ;
        mpfr_mul(r23329, r23328, r23325, MPFR_RNDN);
        mpfr_mul(r23330, r23329, r23329, MPFR_RNDN);
        mpfr_div(r23331, r23327, r23330, MPFR_RNDN);
        ;
        mpfr_sub(r23333, r23330, r23332, MPFR_RNDN);
        mpfr_div(r23334, r23331, r23333, MPFR_RNDN);
        return mpfr_get_d(r23334, MPFR_RNDN);
}

static mpfr_t r23335, r23336, r23337, r23338, r23339, r23340, r23341, r23342, r23343;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r23335);
        mpfr_init_set_str(r23336, "2", 10, MPFR_RNDN);
        mpfr_init(r23337);
        mpfr_init_set_str(r23338, "8", 10, MPFR_RNDN);
        mpfr_init(r23339);
        mpfr_init_set_str(r23340, "2.0", 10, MPFR_RNDN);
        mpfr_init(r23341);
        mpfr_init(r23342);
        mpfr_init(r23343);
}

double f_fm(double i) {
        mpfr_set_d(r23335, i, MPFR_RNDN);
        ;
        mpfr_div(r23337, r23335, r23336, MPFR_RNDN);
        ;
        mpfr_mul(r23339, r23338, r23335, MPFR_RNDN);
        ;
        mpfr_div(r23341, r23340, r23335, MPFR_RNDN);
        mpfr_sub(r23342, r23339, r23341, MPFR_RNDN);
        mpfr_div(r23343, r23337, r23342, MPFR_RNDN);
        return mpfr_get_d(r23343, MPFR_RNDN);
}

static mpfr_t r23344, r23345, r23346, r23347, r23348, r23349, r23350, r23351, r23352;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r23344);
        mpfr_init_set_str(r23345, "2", 10, MPFR_RNDN);
        mpfr_init(r23346);
        mpfr_init_set_str(r23347, "8", 10, MPFR_RNDN);
        mpfr_init(r23348);
        mpfr_init_set_str(r23349, "2.0", 10, MPFR_RNDN);
        mpfr_init(r23350);
        mpfr_init(r23351);
        mpfr_init(r23352);
}

double f_dm(double i) {
        mpfr_set_d(r23344, i, MPFR_RNDN);
        ;
        mpfr_div(r23346, r23344, r23345, MPFR_RNDN);
        ;
        mpfr_mul(r23348, r23347, r23344, MPFR_RNDN);
        ;
        mpfr_div(r23350, r23349, r23344, MPFR_RNDN);
        mpfr_sub(r23351, r23348, r23350, MPFR_RNDN);
        mpfr_div(r23352, r23346, r23351, MPFR_RNDN);
        return mpfr_get_d(r23352, MPFR_RNDN);
}