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

char *name = "NMSE Section 6.1 mentioned, B";

double f_if(float a, float b) {
        float r8474 = atan2(1.0, 0.0);
        float r8475 = 2;
        float r8476 = r8474 / r8475;
        float r8477 = 1;
        float r8478 = b;
        float r8479 = r8478 * r8478;
        float r8480 = a;
        float r8481 = r8480 * r8480;
        float r8482 = r8479 - r8481;
        float r8483 = r8477 / r8482;
        float r8484 = r8476 * r8483;
        float r8485 = r8477 / r8480;
        float r8486 = r8477 / r8478;
        float r8487 = r8485 - r8486;
        float r8488 = r8484 * r8487;
        return r8488;
}

double f_id(double a, double b) {
        double r8489 = atan2(1.0, 0.0);
        double r8490 = 2;
        double r8491 = r8489 / r8490;
        double r8492 = 1;
        double r8493 = b;
        double r8494 = r8493 * r8493;
        double r8495 = a;
        double r8496 = r8495 * r8495;
        double r8497 = r8494 - r8496;
        double r8498 = r8492 / r8497;
        double r8499 = r8491 * r8498;
        double r8500 = r8492 / r8495;
        double r8501 = r8492 / r8493;
        double r8502 = r8500 - r8501;
        double r8503 = r8499 * r8502;
        return r8503;
}


double f_of(float a, float b) {
        float r8504 = 1;
        float r8505 = a;
        float r8506 = r8504 / r8505;
        float r8507 = b;
        float r8508 = r8504 / r8507;
        float r8509 = r8506 - r8508;
        float r8510 = r8505 + r8507;
        float r8511 = 2;
        float r8512 = r8510 * r8511;
        float r8513 = r8509 / r8512;
        float r8514 = r8507 - r8505;
        float r8515 = atan2(1.0, 0.0);
        float r8516 = r8514 / r8515;
        float r8517 = r8513 / r8516;
        return r8517;
}

double f_od(double a, double b) {
        double r8518 = 1;
        double r8519 = a;
        double r8520 = r8518 / r8519;
        double r8521 = b;
        double r8522 = r8518 / r8521;
        double r8523 = r8520 - r8522;
        double r8524 = r8519 + r8521;
        double r8525 = 2;
        double r8526 = r8524 * r8525;
        double r8527 = r8523 / r8526;
        double r8528 = r8521 - r8519;
        double r8529 = atan2(1.0, 0.0);
        double r8530 = r8528 / r8529;
        double r8531 = r8527 / r8530;
        return r8531;
}

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 r8532, r8533, r8534, r8535, r8536, r8537, r8538, r8539, r8540, r8541, r8542, r8543, r8544, r8545, r8546;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r8532);
        mpfr_init_set_str(r8533, "2", 10, MPFR_RNDN);
        mpfr_init(r8534);
        mpfr_init_set_str(r8535, "1", 10, MPFR_RNDN);
        mpfr_init(r8536);
        mpfr_init(r8537);
        mpfr_init(r8538);
        mpfr_init(r8539);
        mpfr_init(r8540);
        mpfr_init(r8541);
        mpfr_init(r8542);
        mpfr_init(r8543);
        mpfr_init(r8544);
        mpfr_init(r8545);
        mpfr_init(r8546);
}

double f_im(double a, double b) {
        mpfr_const_pi(r8532, MPFR_RNDN);
        ;
        mpfr_div(r8534, r8532, r8533, MPFR_RNDN);
        ;
        mpfr_set_d(r8536, b, MPFR_RNDN);
        mpfr_mul(r8537, r8536, r8536, MPFR_RNDN);
        mpfr_set_d(r8538, a, MPFR_RNDN);
        mpfr_mul(r8539, r8538, r8538, MPFR_RNDN);
        mpfr_sub(r8540, r8537, r8539, MPFR_RNDN);
        mpfr_div(r8541, r8535, r8540, MPFR_RNDN);
        mpfr_mul(r8542, r8534, r8541, MPFR_RNDN);
        mpfr_div(r8543, r8535, r8538, MPFR_RNDN);
        mpfr_div(r8544, r8535, r8536, MPFR_RNDN);
        mpfr_sub(r8545, r8543, r8544, MPFR_RNDN);
        mpfr_mul(r8546, r8542, r8545, MPFR_RNDN);
        return mpfr_get_d(r8546, MPFR_RNDN);
}

static mpfr_t r8547, r8548, r8549, r8550, r8551, r8552, r8553, r8554, r8555, r8556, r8557, r8558, r8559, r8560;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r8547, "1", 10, MPFR_RNDN);
        mpfr_init(r8548);
        mpfr_init(r8549);
        mpfr_init(r8550);
        mpfr_init(r8551);
        mpfr_init(r8552);
        mpfr_init(r8553);
        mpfr_init_set_str(r8554, "2", 10, MPFR_RNDN);
        mpfr_init(r8555);
        mpfr_init(r8556);
        mpfr_init(r8557);
        mpfr_init(r8558);
        mpfr_init(r8559);
        mpfr_init(r8560);
}

double f_fm(double a, double b) {
        ;
        mpfr_set_d(r8548, a, MPFR_RNDN);
        mpfr_div(r8549, r8547, r8548, MPFR_RNDN);
        mpfr_set_d(r8550, b, MPFR_RNDN);
        mpfr_div(r8551, r8547, r8550, MPFR_RNDN);
        mpfr_sub(r8552, r8549, r8551, MPFR_RNDN);
        mpfr_add(r8553, r8548, r8550, MPFR_RNDN);
        ;
        mpfr_mul(r8555, r8553, r8554, MPFR_RNDN);
        mpfr_div(r8556, r8552, r8555, MPFR_RNDN);
        mpfr_sub(r8557, r8550, r8548, MPFR_RNDN);
        mpfr_const_pi(r8558, MPFR_RNDN);
        mpfr_div(r8559, r8557, r8558, MPFR_RNDN);
        mpfr_div(r8560, r8556, r8559, MPFR_RNDN);
        return mpfr_get_d(r8560, MPFR_RNDN);
}

static mpfr_t r8561, r8562, r8563, r8564, r8565, r8566, r8567, r8568, r8569, r8570, r8571, r8572, r8573, r8574;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r8561, "1", 10, MPFR_RNDN);
        mpfr_init(r8562);
        mpfr_init(r8563);
        mpfr_init(r8564);
        mpfr_init(r8565);
        mpfr_init(r8566);
        mpfr_init(r8567);
        mpfr_init_set_str(r8568, "2", 10, MPFR_RNDN);
        mpfr_init(r8569);
        mpfr_init(r8570);
        mpfr_init(r8571);
        mpfr_init(r8572);
        mpfr_init(r8573);
        mpfr_init(r8574);
}

double f_dm(double a, double b) {
        ;
        mpfr_set_d(r8562, a, MPFR_RNDN);
        mpfr_div(r8563, r8561, r8562, MPFR_RNDN);
        mpfr_set_d(r8564, b, MPFR_RNDN);
        mpfr_div(r8565, r8561, r8564, MPFR_RNDN);
        mpfr_sub(r8566, r8563, r8565, MPFR_RNDN);
        mpfr_add(r8567, r8562, r8564, MPFR_RNDN);
        ;
        mpfr_mul(r8569, r8567, r8568, MPFR_RNDN);
        mpfr_div(r8570, r8566, r8569, MPFR_RNDN);
        mpfr_sub(r8571, r8564, r8562, MPFR_RNDN);
        mpfr_const_pi(r8572, MPFR_RNDN);
        mpfr_div(r8573, r8571, r8572, MPFR_RNDN);
        mpfr_div(r8574, r8570, r8573, MPFR_RNDN);
        return mpfr_get_d(r8574, MPFR_RNDN);
}