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

char *name = "Difference of squares";

double f_if(float a, float b) {
        float r23459 = a;
        float r23460 = r23459 * r23459;
        float r23461 = b;
        float r23462 = r23461 * r23461;
        float r23463 = r23460 - r23462;
        return r23463;
}

double f_id(double a, double b) {
        double r23464 = a;
        double r23465 = r23464 * r23464;
        double r23466 = b;
        double r23467 = r23466 * r23466;
        double r23468 = r23465 - r23467;
        return r23468;
}


double f_of(float a, float b) {
        float r23469 = a;
        float r23470 = b;
        float r23471 = r23469 + r23470;
        float r23472 = r23469 - r23470;
        float r23473 = r23471 * r23472;
        return r23473;
}

double f_od(double a, double b) {
        double r23474 = a;
        double r23475 = b;
        double r23476 = r23474 + r23475;
        double r23477 = r23474 - r23475;
        double r23478 = r23476 * r23477;
        return r23478;
}

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 r23479, r23480, r23481, r23482, r23483;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r23479);
        mpfr_init(r23480);
        mpfr_init(r23481);
        mpfr_init(r23482);
        mpfr_init(r23483);
}

double f_im(double a, double b) {
        mpfr_set_d(r23479, a, MPFR_RNDN);
        mpfr_mul(r23480, r23479, r23479, MPFR_RNDN);
        mpfr_set_d(r23481, b, MPFR_RNDN);
        mpfr_mul(r23482, r23481, r23481, MPFR_RNDN);
        mpfr_sub(r23483, r23480, r23482, MPFR_RNDN);
        return mpfr_get_d(r23483, MPFR_RNDN);
}

static mpfr_t r23484, r23485, r23486, r23487, r23488;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r23484);
        mpfr_init(r23485);
        mpfr_init(r23486);
        mpfr_init(r23487);
        mpfr_init(r23488);
}

double f_fm(double a, double b) {
        mpfr_set_d(r23484, a, MPFR_RNDN);
        mpfr_set_d(r23485, b, MPFR_RNDN);
        mpfr_add(r23486, r23484, r23485, MPFR_RNDN);
        mpfr_sub(r23487, r23484, r23485, MPFR_RNDN);
        mpfr_mul(r23488, r23486, r23487, MPFR_RNDN);
        return mpfr_get_d(r23488, MPFR_RNDN);
}

static mpfr_t r23489, r23490, r23491, r23492, r23493;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r23489);
        mpfr_init(r23490);
        mpfr_init(r23491);
        mpfr_init(r23492);
        mpfr_init(r23493);
}

double f_dm(double a, double b) {
        mpfr_set_d(r23489, a, MPFR_RNDN);
        mpfr_set_d(r23490, b, MPFR_RNDN);
        mpfr_add(r23491, r23489, r23490, MPFR_RNDN);
        mpfr_sub(r23492, r23489, r23490, MPFR_RNDN);
        mpfr_mul(r23493, r23491, r23492, MPFR_RNDN);
        return mpfr_get_d(r23493, MPFR_RNDN);
}