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

char *name = "Data.Colour.RGBSpace.HSV:hsv from colour-2.3.3, J";

double f_if(float x, float y, float z) {
        float r56613 = x;
        float r56614 = 1.0;
        float r56615 = y;
        float r56616 = r56614 - r56615;
        float r56617 = z;
        float r56618 = r56616 * r56617;
        float r56619 = r56614 - r56618;
        float r56620 = r56613 * r56619;
        return r56620;
}

double f_id(double x, double y, double z) {
        double r56621 = x;
        double r56622 = 1.0;
        double r56623 = y;
        double r56624 = r56622 - r56623;
        double r56625 = z;
        double r56626 = r56624 * r56625;
        double r56627 = r56622 - r56626;
        double r56628 = r56621 * r56627;
        return r56628;
}


double f_of(float x, float y, float z) {
        float r56629 = x;
        float r56630 = 1.0;
        float r56631 = r56629 * r56630;
        float r56632 = z;
        float r56633 = r56632 * r56629;
        float r56634 = y;
        float r56635 = r56630 - r56634;
        float r56636 = r56633 * r56635;
        float r56637 = r56631 - r56636;
        return r56637;
}

double f_od(double x, double y, double z) {
        double r56638 = x;
        double r56639 = 1.0;
        double r56640 = r56638 * r56639;
        double r56641 = z;
        double r56642 = r56641 * r56638;
        double r56643 = y;
        double r56644 = r56639 - r56643;
        double r56645 = r56642 * r56644;
        double r56646 = r56640 - r56645;
        return r56646;
}

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 r56647, r56648, r56649, r56650, r56651, r56652, r56653, r56654;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r56647);
        mpfr_init_set_str(r56648, "1.0", 10, MPFR_RNDN);
        mpfr_init(r56649);
        mpfr_init(r56650);
        mpfr_init(r56651);
        mpfr_init(r56652);
        mpfr_init(r56653);
        mpfr_init(r56654);
}

double f_im(double x, double y, double z) {
        mpfr_set_d(r56647, x, MPFR_RNDN);
        ;
        mpfr_set_d(r56649, y, MPFR_RNDN);
        mpfr_sub(r56650, r56648, r56649, MPFR_RNDN);
        mpfr_set_d(r56651, z, MPFR_RNDN);
        mpfr_mul(r56652, r56650, r56651, MPFR_RNDN);
        mpfr_sub(r56653, r56648, r56652, MPFR_RNDN);
        mpfr_mul(r56654, r56647, r56653, MPFR_RNDN);
        return mpfr_get_d(r56654, MPFR_RNDN);
}

static mpfr_t r56655, r56656, r56657, r56658, r56659, r56660, r56661, r56662, r56663;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r56655);
        mpfr_init_set_str(r56656, "1.0", 10, MPFR_RNDN);
        mpfr_init(r56657);
        mpfr_init(r56658);
        mpfr_init(r56659);
        mpfr_init(r56660);
        mpfr_init(r56661);
        mpfr_init(r56662);
        mpfr_init(r56663);
}

double f_fm(double x, double y, double z) {
        mpfr_set_d(r56655, x, MPFR_RNDN);
        ;
        mpfr_mul(r56657, r56655, r56656, MPFR_RNDN);
        mpfr_set_d(r56658, z, MPFR_RNDN);
        mpfr_mul(r56659, r56658, r56655, MPFR_RNDN);
        mpfr_set_d(r56660, y, MPFR_RNDN);
        mpfr_sub(r56661, r56656, r56660, MPFR_RNDN);
        mpfr_mul(r56662, r56659, r56661, MPFR_RNDN);
        mpfr_sub(r56663, r56657, r56662, MPFR_RNDN);
        return mpfr_get_d(r56663, MPFR_RNDN);
}

static mpfr_t r56664, r56665, r56666, r56667, r56668, r56669, r56670, r56671, r56672;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r56664);
        mpfr_init_set_str(r56665, "1.0", 10, MPFR_RNDN);
        mpfr_init(r56666);
        mpfr_init(r56667);
        mpfr_init(r56668);
        mpfr_init(r56669);
        mpfr_init(r56670);
        mpfr_init(r56671);
        mpfr_init(r56672);
}

double f_dm(double x, double y, double z) {
        mpfr_set_d(r56664, x, MPFR_RNDN);
        ;
        mpfr_mul(r56666, r56664, r56665, MPFR_RNDN);
        mpfr_set_d(r56667, z, MPFR_RNDN);
        mpfr_mul(r56668, r56667, r56664, MPFR_RNDN);
        mpfr_set_d(r56669, y, MPFR_RNDN);
        mpfr_sub(r56670, r56665, r56669, MPFR_RNDN);
        mpfr_mul(r56671, r56668, r56670, MPFR_RNDN);
        mpfr_sub(r56672, r56666, r56671, MPFR_RNDN);
        return mpfr_get_d(r56672, MPFR_RNDN);
}