#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 r8640 = x;
        float r8641 = 4;
        float r8642 = pow(r8640, r8641);
        float r8643 = y;
        float r8644 = pow(r8643, r8641);
        float r8645 = r8642 - r8644;
        return r8645;
}

double f_id(double x, double y) {
        double r8646 = x;
        double r8647 = 4;
        double r8648 = pow(r8646, r8647);
        double r8649 = y;
        double r8650 = pow(r8649, r8647);
        double r8651 = r8648 - r8650;
        return r8651;
}


double f_of(float x, float y) {
        float r8652 = x;
        float r8653 = 4;
        float r8654 = pow(r8652, r8653);
        float r8655 = y;
        float r8656 = pow(r8655, r8653);
        float r8657 = r8654 - r8656;
        return r8657;
}

double f_od(double x, double y) {
        double r8658 = x;
        double r8659 = 4;
        double r8660 = pow(r8658, r8659);
        double r8661 = y;
        double r8662 = pow(r8661, r8659);
        double r8663 = r8660 - r8662;
        return r8663;
}

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 r8664, r8665, r8666, r8667, r8668, r8669;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r8664);
        mpfr_init_set_str(r8665, "4", 10, MPFR_RNDN);
        mpfr_init(r8666);
        mpfr_init(r8667);
        mpfr_init(r8668);
        mpfr_init(r8669);
}

double f_im(double x, double y) {
        mpfr_set_d(r8664, x, MPFR_RNDN);
        ;
        mpfr_pow(r8666, r8664, r8665, MPFR_RNDN);
        mpfr_set_d(r8667, y, MPFR_RNDN);
        mpfr_pow(r8668, r8667, r8665, MPFR_RNDN);
        mpfr_sub(r8669, r8666, r8668, MPFR_RNDN);
        return mpfr_get_d(r8669, MPFR_RNDN);
}

static mpfr_t r8670, r8671, r8672, r8673, r8674, r8675;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r8670);
        mpfr_init_set_str(r8671, "4", 10, MPFR_RNDN);
        mpfr_init(r8672);
        mpfr_init(r8673);
        mpfr_init(r8674);
        mpfr_init(r8675);
}

double f_fm(double x, double y) {
        mpfr_set_d(r8670, x, MPFR_RNDN);
        ;
        mpfr_pow(r8672, r8670, r8671, MPFR_RNDN);
        mpfr_set_d(r8673, y, MPFR_RNDN);
        mpfr_pow(r8674, r8673, r8671, MPFR_RNDN);
        mpfr_sub(r8675, r8672, r8674, MPFR_RNDN);
        return mpfr_get_d(r8675, MPFR_RNDN);
}

static mpfr_t r8676, r8677, r8678, r8679, r8680, r8681;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r8676);
        mpfr_init_set_str(r8677, "4", 10, MPFR_RNDN);
        mpfr_init(r8678);
        mpfr_init(r8679);
        mpfr_init(r8680);
        mpfr_init(r8681);
}

double f_dm(double x, double y) {
        mpfr_set_d(r8676, x, MPFR_RNDN);
        ;
        mpfr_pow(r8678, r8676, r8677, MPFR_RNDN);
        mpfr_set_d(r8679, y, MPFR_RNDN);
        mpfr_pow(r8680, r8679, r8677, MPFR_RNDN);
        mpfr_sub(r8681, r8678, r8680, MPFR_RNDN);
        return mpfr_get_d(r8681, MPFR_RNDN);
}