#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 r21125 = x;
        float r21126 = 4;
        float r21127 = pow(r21125, r21126);
        float r21128 = y;
        float r21129 = pow(r21128, r21126);
        float r21130 = r21127 - r21129;
        return r21130;
}

double f_id(double x, double y) {
        double r21131 = x;
        double r21132 = 4;
        double r21133 = pow(r21131, r21132);
        double r21134 = y;
        double r21135 = pow(r21134, r21132);
        double r21136 = r21133 - r21135;
        return r21136;
}


double f_of(float x, float y) {
        float r21137 = x;
        float r21138 = 4;
        float r21139 = pow(r21137, r21138);
        float r21140 = y;
        float r21141 = pow(r21140, r21138);
        float r21142 = r21139 - r21141;
        return r21142;
}

double f_od(double x, double y) {
        double r21143 = x;
        double r21144 = 4;
        double r21145 = pow(r21143, r21144);
        double r21146 = y;
        double r21147 = pow(r21146, r21144);
        double r21148 = r21145 - r21147;
        return r21148;
}

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 r21149, r21150, r21151, r21152, r21153, r21154;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r21149);
        mpfr_init_set_str(r21150, "4", 10, MPFR_RNDN);
        mpfr_init(r21151);
        mpfr_init(r21152);
        mpfr_init(r21153);
        mpfr_init(r21154);
}

double f_im(double x, double y) {
        mpfr_set_d(r21149, x, MPFR_RNDN);
        ;
        mpfr_pow(r21151, r21149, r21150, MPFR_RNDN);
        mpfr_set_d(r21152, y, MPFR_RNDN);
        mpfr_pow(r21153, r21152, r21150, MPFR_RNDN);
        mpfr_sub(r21154, r21151, r21153, MPFR_RNDN);
        return mpfr_get_d(r21154, MPFR_RNDN);
}

static mpfr_t r21155, r21156, r21157, r21158, r21159, r21160;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r21155);
        mpfr_init_set_str(r21156, "4", 10, MPFR_RNDN);
        mpfr_init(r21157);
        mpfr_init(r21158);
        mpfr_init(r21159);
        mpfr_init(r21160);
}

double f_fm(double x, double y) {
        mpfr_set_d(r21155, x, MPFR_RNDN);
        ;
        mpfr_pow(r21157, r21155, r21156, MPFR_RNDN);
        mpfr_set_d(r21158, y, MPFR_RNDN);
        mpfr_pow(r21159, r21158, r21156, MPFR_RNDN);
        mpfr_sub(r21160, r21157, r21159, MPFR_RNDN);
        return mpfr_get_d(r21160, MPFR_RNDN);
}

static mpfr_t r21161, r21162, r21163, r21164, r21165, r21166;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r21161);
        mpfr_init_set_str(r21162, "4", 10, MPFR_RNDN);
        mpfr_init(r21163);
        mpfr_init(r21164);
        mpfr_init(r21165);
        mpfr_init(r21166);
}

double f_dm(double x, double y) {
        mpfr_set_d(r21161, x, MPFR_RNDN);
        ;
        mpfr_pow(r21163, r21161, r21162, MPFR_RNDN);
        mpfr_set_d(r21164, y, MPFR_RNDN);
        mpfr_pow(r21165, r21164, r21162, MPFR_RNDN);
        mpfr_sub(r21166, r21163, r21165, MPFR_RNDN);
        return mpfr_get_d(r21166, MPFR_RNDN);
}