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

char *name = "Graphics.Rasterific.Shading:$sgradientColorAt from Rasterific-0.6.1";

double f_if(float x, float y, float z) {
        float r50916 = x;
        float r50917 = y;
        float r50918 = r50916 - r50917;
        float r50919 = z;
        float r50920 = r50919 - r50917;
        float r50921 = r50918 / r50920;
        return r50921;
}

double f_id(double x, double y, double z) {
        double r50922 = x;
        double r50923 = y;
        double r50924 = r50922 - r50923;
        double r50925 = z;
        double r50926 = r50925 - r50923;
        double r50927 = r50924 / r50926;
        return r50927;
}


double f_of(float x, float y, float z) {
        float r50928 = x;
        float r50929 = y;
        float r50930 = r50928 - r50929;
        float r50931 = z;
        float r50932 = r50931 - r50929;
        float r50933 = r50930 / r50932;
        return r50933;
}

double f_od(double x, double y, double z) {
        double r50934 = x;
        double r50935 = y;
        double r50936 = r50934 - r50935;
        double r50937 = z;
        double r50938 = r50937 - r50935;
        double r50939 = r50936 / r50938;
        return r50939;
}

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 r50940, r50941, r50942, r50943, r50944, r50945;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r50940);
        mpfr_init(r50941);
        mpfr_init(r50942);
        mpfr_init(r50943);
        mpfr_init(r50944);
        mpfr_init(r50945);
}

double f_im(double x, double y, double z) {
        mpfr_set_d(r50940, x, MPFR_RNDN);
        mpfr_set_d(r50941, y, MPFR_RNDN);
        mpfr_sub(r50942, r50940, r50941, MPFR_RNDN);
        mpfr_set_d(r50943, z, MPFR_RNDN);
        mpfr_sub(r50944, r50943, r50941, MPFR_RNDN);
        mpfr_div(r50945, r50942, r50944, MPFR_RNDN);
        return mpfr_get_d(r50945, MPFR_RNDN);
}

static mpfr_t r50946, r50947, r50948, r50949, r50950, r50951;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r50946);
        mpfr_init(r50947);
        mpfr_init(r50948);
        mpfr_init(r50949);
        mpfr_init(r50950);
        mpfr_init(r50951);
}

double f_fm(double x, double y, double z) {
        mpfr_set_d(r50946, x, MPFR_RNDN);
        mpfr_set_d(r50947, y, MPFR_RNDN);
        mpfr_sub(r50948, r50946, r50947, MPFR_RNDN);
        mpfr_set_d(r50949, z, MPFR_RNDN);
        mpfr_sub(r50950, r50949, r50947, MPFR_RNDN);
        mpfr_div(r50951, r50948, r50950, MPFR_RNDN);
        return mpfr_get_d(r50951, MPFR_RNDN);
}

static mpfr_t r50952, r50953, r50954, r50955, r50956, r50957;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r50952);
        mpfr_init(r50953);
        mpfr_init(r50954);
        mpfr_init(r50955);
        mpfr_init(r50956);
        mpfr_init(r50957);
}

double f_dm(double x, double y, double z) {
        mpfr_set_d(r50952, x, MPFR_RNDN);
        mpfr_set_d(r50953, y, MPFR_RNDN);
        mpfr_sub(r50954, r50952, r50953, MPFR_RNDN);
        mpfr_set_d(r50955, z, MPFR_RNDN);
        mpfr_sub(r50956, r50955, r50953, MPFR_RNDN);
        mpfr_div(r50957, r50954, r50956, MPFR_RNDN);
        return mpfr_get_d(r50957, MPFR_RNDN);
}