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

char *name = "Diagrams.Backend.Rasterific:rasterificRadialGradient from diagrams-rasterific-1.3.1.3";

double f_if(float x, float y, float z) {
        float r55270 = x;
        float r55271 = y;
        float r55272 = z;
        float r55273 = r55272 - r55270;
        float r55274 = r55271 * r55273;
        float r55275 = r55270 + r55274;
        float r55276 = r55275 / r55272;
        return r55276;
}

double f_id(double x, double y, double z) {
        double r55277 = x;
        double r55278 = y;
        double r55279 = z;
        double r55280 = r55279 - r55277;
        double r55281 = r55278 * r55280;
        double r55282 = r55277 + r55281;
        double r55283 = r55282 / r55279;
        return r55283;
}


double f_of(float x, float y, float z) {
        float r55284 = y;
        float r55285 = x;
        float r55286 = z;
        float r55287 = r55285 / r55286;
        float r55288 = r55284 + r55287;
        float r55289 = r55286 / r55285;
        float r55290 = r55284 / r55289;
        float r55291 = r55288 - r55290;
        return r55291;
}

double f_od(double x, double y, double z) {
        double r55292 = y;
        double r55293 = x;
        double r55294 = z;
        double r55295 = r55293 / r55294;
        double r55296 = r55292 + r55295;
        double r55297 = r55294 / r55293;
        double r55298 = r55292 / r55297;
        double r55299 = r55296 - r55298;
        return r55299;
}

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 r55300, r55301, r55302, r55303, r55304, r55305, r55306;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r55300);
        mpfr_init(r55301);
        mpfr_init(r55302);
        mpfr_init(r55303);
        mpfr_init(r55304);
        mpfr_init(r55305);
        mpfr_init(r55306);
}

double f_im(double x, double y, double z) {
        mpfr_set_d(r55300, x, MPFR_RNDN);
        mpfr_set_d(r55301, y, MPFR_RNDN);
        mpfr_set_d(r55302, z, MPFR_RNDN);
        mpfr_sub(r55303, r55302, r55300, MPFR_RNDN);
        mpfr_mul(r55304, r55301, r55303, MPFR_RNDN);
        mpfr_add(r55305, r55300, r55304, MPFR_RNDN);
        mpfr_div(r55306, r55305, r55302, MPFR_RNDN);
        return mpfr_get_d(r55306, MPFR_RNDN);
}

static mpfr_t r55307, r55308, r55309, r55310, r55311, r55312, r55313, r55314;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r55307);
        mpfr_init(r55308);
        mpfr_init(r55309);
        mpfr_init(r55310);
        mpfr_init(r55311);
        mpfr_init(r55312);
        mpfr_init(r55313);
        mpfr_init(r55314);
}

double f_fm(double x, double y, double z) {
        mpfr_set_d(r55307, y, MPFR_RNDN);
        mpfr_set_d(r55308, x, MPFR_RNDN);
        mpfr_set_d(r55309, z, MPFR_RNDN);
        mpfr_div(r55310, r55308, r55309, MPFR_RNDN);
        mpfr_add(r55311, r55307, r55310, MPFR_RNDN);
        mpfr_div(r55312, r55309, r55308, MPFR_RNDN);
        mpfr_div(r55313, r55307, r55312, MPFR_RNDN);
        mpfr_sub(r55314, r55311, r55313, MPFR_RNDN);
        return mpfr_get_d(r55314, MPFR_RNDN);
}

static mpfr_t r55315, r55316, r55317, r55318, r55319, r55320, r55321, r55322;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r55315);
        mpfr_init(r55316);
        mpfr_init(r55317);
        mpfr_init(r55318);
        mpfr_init(r55319);
        mpfr_init(r55320);
        mpfr_init(r55321);
        mpfr_init(r55322);
}

double f_dm(double x, double y, double z) {
        mpfr_set_d(r55315, y, MPFR_RNDN);
        mpfr_set_d(r55316, x, MPFR_RNDN);
        mpfr_set_d(r55317, z, MPFR_RNDN);
        mpfr_div(r55318, r55316, r55317, MPFR_RNDN);
        mpfr_add(r55319, r55315, r55318, MPFR_RNDN);
        mpfr_div(r55320, r55317, r55316, MPFR_RNDN);
        mpfr_div(r55321, r55315, r55320, MPFR_RNDN);
        mpfr_sub(r55322, r55319, r55321, MPFR_RNDN);
        return mpfr_get_d(r55322, MPFR_RNDN);
}