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

char *name = "Radioactive exchange between two surfaces";

double f_if(float x, float y) {
        float r8063 = x;
        float r8064 = 4.0f;
        float r8065 = pow(r8063, r8064);
        float r8066 = y;
        float r8067 = pow(r8066, r8064);
        float r8068 = r8065 - r8067;
        return r8068;
}

double f_id(double x, double y) {
        double r8069 = x;
        double r8070 = 4.0;
        double r8071 = pow(r8069, r8070);
        double r8072 = y;
        double r8073 = pow(r8072, r8070);
        double r8074 = r8071 - r8073;
        return r8074;
}


double f_of(float x, float y) {
        float r8075 = x;
        float r8076 = 4.0f;
        float r8077 = pow(r8075, r8076);
        float r8078 = y;
        float r8079 = pow(r8078, r8076);
        float r8080 = r8077 - r8079;
        return r8080;
}

double f_od(double x, double y) {
        double r8081 = x;
        double r8082 = 4.0;
        double r8083 = pow(r8081, r8082);
        double r8084 = y;
        double r8085 = pow(r8084, r8082);
        double r8086 = r8083 - r8085;
        return r8086;
}

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 r8087, r8088, r8089, r8090, r8091, r8092;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r8087);
        mpfr_init_set_str(r8088, "4", 10, MPFR_RNDN);
        mpfr_init(r8089);
        mpfr_init(r8090);
        mpfr_init(r8091);
        mpfr_init(r8092);
}

double f_im(double x, double y) {
        mpfr_set_d(r8087, x, MPFR_RNDN);
        ;
        mpfr_pow(r8089, r8087, r8088, MPFR_RNDN);
        mpfr_set_d(r8090, y, MPFR_RNDN);
        mpfr_pow(r8091, r8090, r8088, MPFR_RNDN);
        mpfr_sub(r8092, r8089, r8091, MPFR_RNDN);
        return mpfr_get_d(r8092, MPFR_RNDN);
}

static mpfr_t r8093, r8094, r8095, r8096, r8097, r8098;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r8093);
        mpfr_init_set_str(r8094, "4", 10, MPFR_RNDN);
        mpfr_init(r8095);
        mpfr_init(r8096);
        mpfr_init(r8097);
        mpfr_init(r8098);
}

double f_fm(double x, double y) {
        mpfr_set_d(r8093, x, MPFR_RNDN);
        ;
        mpfr_pow(r8095, r8093, r8094, MPFR_RNDN);
        mpfr_set_d(r8096, y, MPFR_RNDN);
        mpfr_pow(r8097, r8096, r8094, MPFR_RNDN);
        mpfr_sub(r8098, r8095, r8097, MPFR_RNDN);
        return mpfr_get_d(r8098, MPFR_RNDN);
}

static mpfr_t r8099, r8100, r8101, r8102, r8103, r8104;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r8099);
        mpfr_init_set_str(r8100, "4", 10, MPFR_RNDN);
        mpfr_init(r8101);
        mpfr_init(r8102);
        mpfr_init(r8103);
        mpfr_init(r8104);
}

double f_dm(double x, double y) {
        mpfr_set_d(r8099, x, MPFR_RNDN);
        ;
        mpfr_pow(r8101, r8099, r8100, MPFR_RNDN);
        mpfr_set_d(r8102, y, MPFR_RNDN);
        mpfr_pow(r8103, r8102, r8100, MPFR_RNDN);
        mpfr_sub(r8104, r8101, r8103, MPFR_RNDN);
        return mpfr_get_d(r8104, MPFR_RNDN);
}