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

char *name = "Diagrams.Color.HSV:lerp  from diagrams-contrib-1.3.0.5";

double f_if(float x, float y, float z) {
        float r54981 = 1.0;
        float r54982 = x;
        float r54983 = r54981 - r54982;
        float r54984 = y;
        float r54985 = r54983 * r54984;
        float r54986 = z;
        float r54987 = r54982 * r54986;
        float r54988 = r54985 + r54987;
        return r54988;
}

double f_id(double x, double y, double z) {
        double r54989 = 1.0;
        double r54990 = x;
        double r54991 = r54989 - r54990;
        double r54992 = y;
        double r54993 = r54991 * r54992;
        double r54994 = z;
        double r54995 = r54990 * r54994;
        double r54996 = r54993 + r54995;
        return r54996;
}


double f_of(float x, float y, float z) {
        float r54997 = 1.0;
        float r54998 = x;
        float r54999 = r54997 - r54998;
        float r55000 = y;
        float r55001 = r54999 * r55000;
        float r55002 = z;
        float r55003 = r54998 * r55002;
        float r55004 = r55001 + r55003;
        return r55004;
}

double f_od(double x, double y, double z) {
        double r55005 = 1.0;
        double r55006 = x;
        double r55007 = r55005 - r55006;
        double r55008 = y;
        double r55009 = r55007 * r55008;
        double r55010 = z;
        double r55011 = r55006 * r55010;
        double r55012 = r55009 + r55011;
        return r55012;
}

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 r55013, r55014, r55015, r55016, r55017, r55018, r55019, r55020;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r55013, "1.0", 10, MPFR_RNDN);
        mpfr_init(r55014);
        mpfr_init(r55015);
        mpfr_init(r55016);
        mpfr_init(r55017);
        mpfr_init(r55018);
        mpfr_init(r55019);
        mpfr_init(r55020);
}

double f_im(double x, double y, double z) {
        ;
        mpfr_set_d(r55014, x, MPFR_RNDN);
        mpfr_sub(r55015, r55013, r55014, MPFR_RNDN);
        mpfr_set_d(r55016, y, MPFR_RNDN);
        mpfr_mul(r55017, r55015, r55016, MPFR_RNDN);
        mpfr_set_d(r55018, z, MPFR_RNDN);
        mpfr_mul(r55019, r55014, r55018, MPFR_RNDN);
        mpfr_add(r55020, r55017, r55019, MPFR_RNDN);
        return mpfr_get_d(r55020, MPFR_RNDN);
}

static mpfr_t r55021, r55022, r55023, r55024, r55025, r55026, r55027, r55028;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r55021, "1.0", 10, MPFR_RNDN);
        mpfr_init(r55022);
        mpfr_init(r55023);
        mpfr_init(r55024);
        mpfr_init(r55025);
        mpfr_init(r55026);
        mpfr_init(r55027);
        mpfr_init(r55028);
}

double f_fm(double x, double y, double z) {
        ;
        mpfr_set_d(r55022, x, MPFR_RNDN);
        mpfr_sub(r55023, r55021, r55022, MPFR_RNDN);
        mpfr_set_d(r55024, y, MPFR_RNDN);
        mpfr_mul(r55025, r55023, r55024, MPFR_RNDN);
        mpfr_set_d(r55026, z, MPFR_RNDN);
        mpfr_mul(r55027, r55022, r55026, MPFR_RNDN);
        mpfr_add(r55028, r55025, r55027, MPFR_RNDN);
        return mpfr_get_d(r55028, MPFR_RNDN);
}

static mpfr_t r55029, r55030, r55031, r55032, r55033, r55034, r55035, r55036;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r55029, "1.0", 10, MPFR_RNDN);
        mpfr_init(r55030);
        mpfr_init(r55031);
        mpfr_init(r55032);
        mpfr_init(r55033);
        mpfr_init(r55034);
        mpfr_init(r55035);
        mpfr_init(r55036);
}

double f_dm(double x, double y, double z) {
        ;
        mpfr_set_d(r55030, x, MPFR_RNDN);
        mpfr_sub(r55031, r55029, r55030, MPFR_RNDN);
        mpfr_set_d(r55032, y, MPFR_RNDN);
        mpfr_mul(r55033, r55031, r55032, MPFR_RNDN);
        mpfr_set_d(r55034, z, MPFR_RNDN);
        mpfr_mul(r55035, r55030, r55034, MPFR_RNDN);
        mpfr_add(r55036, r55033, r55035, MPFR_RNDN);
        return mpfr_get_d(r55036, MPFR_RNDN);
}