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

char *name = "Graphics.Rasterific.CubicBezier:cachedBezierAt from Rasterific-0.6.1";

double f_if(float x, float y, float z, float t, float a, float b) {
        float r50902 = x;
        float r50903 = y;
        float r50904 = z;
        float r50905 = r50903 * r50904;
        float r50906 = r50902 + r50905;
        float r50907 = t;
        float r50908 = a;
        float r50909 = r50907 * r50908;
        float r50910 = r50906 + r50909;
        float r50911 = r50908 * r50904;
        float r50912 = b;
        float r50913 = r50911 * r50912;
        float r50914 = r50910 + r50913;
        return r50914;
}

double f_id(double x, double y, double z, double t, double a, double b) {
        double r50915 = x;
        double r50916 = y;
        double r50917 = z;
        double r50918 = r50916 * r50917;
        double r50919 = r50915 + r50918;
        double r50920 = t;
        double r50921 = a;
        double r50922 = r50920 * r50921;
        double r50923 = r50919 + r50922;
        double r50924 = r50921 * r50917;
        double r50925 = b;
        double r50926 = r50924 * r50925;
        double r50927 = r50923 + r50926;
        return r50927;
}


double f_of(float x, float y, float z, float t, float a, float b) {
        float r50928 = x;
        float r50929 = a;
        float r50930 = t;
        float r50931 = r50929 * r50930;
        float r50932 = r50928 + r50931;
        float r50933 = z;
        float r50934 = b;
        float r50935 = r50934 * r50929;
        float r50936 = y;
        float r50937 = r50935 + r50936;
        float r50938 = r50933 * r50937;
        float r50939 = r50932 + r50938;
        return r50939;
}

double f_od(double x, double y, double z, double t, double a, double b) {
        double r50940 = x;
        double r50941 = a;
        double r50942 = t;
        double r50943 = r50941 * r50942;
        double r50944 = r50940 + r50943;
        double r50945 = z;
        double r50946 = b;
        double r50947 = r50946 * r50941;
        double r50948 = y;
        double r50949 = r50947 + r50948;
        double r50950 = r50945 * r50949;
        double r50951 = r50944 + r50950;
        return r50951;
}

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 r50952, r50953, r50954, r50955, r50956, r50957, r50958, r50959, r50960, r50961, r50962, r50963, r50964;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r50952);
        mpfr_init(r50953);
        mpfr_init(r50954);
        mpfr_init(r50955);
        mpfr_init(r50956);
        mpfr_init(r50957);
        mpfr_init(r50958);
        mpfr_init(r50959);
        mpfr_init(r50960);
        mpfr_init(r50961);
        mpfr_init(r50962);
        mpfr_init(r50963);
        mpfr_init(r50964);
}

double f_im(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r50952, x, MPFR_RNDN);
        mpfr_set_d(r50953, y, MPFR_RNDN);
        mpfr_set_d(r50954, z, MPFR_RNDN);
        mpfr_mul(r50955, r50953, r50954, MPFR_RNDN);
        mpfr_add(r50956, r50952, r50955, MPFR_RNDN);
        mpfr_set_d(r50957, t, MPFR_RNDN);
        mpfr_set_d(r50958, a, MPFR_RNDN);
        mpfr_mul(r50959, r50957, r50958, MPFR_RNDN);
        mpfr_add(r50960, r50956, r50959, MPFR_RNDN);
        mpfr_mul(r50961, r50958, r50954, MPFR_RNDN);
        mpfr_set_d(r50962, b, MPFR_RNDN);
        mpfr_mul(r50963, r50961, r50962, MPFR_RNDN);
        mpfr_add(r50964, r50960, r50963, MPFR_RNDN);
        return mpfr_get_d(r50964, MPFR_RNDN);
}

static mpfr_t r50965, r50966, r50967, r50968, r50969, r50970, r50971, r50972, r50973, r50974, r50975, r50976;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r50965);
        mpfr_init(r50966);
        mpfr_init(r50967);
        mpfr_init(r50968);
        mpfr_init(r50969);
        mpfr_init(r50970);
        mpfr_init(r50971);
        mpfr_init(r50972);
        mpfr_init(r50973);
        mpfr_init(r50974);
        mpfr_init(r50975);
        mpfr_init(r50976);
}

double f_fm(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r50965, x, MPFR_RNDN);
        mpfr_set_d(r50966, a, MPFR_RNDN);
        mpfr_set_d(r50967, t, MPFR_RNDN);
        mpfr_mul(r50968, r50966, r50967, MPFR_RNDN);
        mpfr_add(r50969, r50965, r50968, MPFR_RNDN);
        mpfr_set_d(r50970, z, MPFR_RNDN);
        mpfr_set_d(r50971, b, MPFR_RNDN);
        mpfr_mul(r50972, r50971, r50966, MPFR_RNDN);
        mpfr_set_d(r50973, y, MPFR_RNDN);
        mpfr_add(r50974, r50972, r50973, MPFR_RNDN);
        mpfr_mul(r50975, r50970, r50974, MPFR_RNDN);
        mpfr_add(r50976, r50969, r50975, MPFR_RNDN);
        return mpfr_get_d(r50976, MPFR_RNDN);
}

static mpfr_t r50977, r50978, r50979, r50980, r50981, r50982, r50983, r50984, r50985, r50986, r50987, r50988;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r50977);
        mpfr_init(r50978);
        mpfr_init(r50979);
        mpfr_init(r50980);
        mpfr_init(r50981);
        mpfr_init(r50982);
        mpfr_init(r50983);
        mpfr_init(r50984);
        mpfr_init(r50985);
        mpfr_init(r50986);
        mpfr_init(r50987);
        mpfr_init(r50988);
}

double f_dm(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r50977, x, MPFR_RNDN);
        mpfr_set_d(r50978, a, MPFR_RNDN);
        mpfr_set_d(r50979, t, MPFR_RNDN);
        mpfr_mul(r50980, r50978, r50979, MPFR_RNDN);
        mpfr_add(r50981, r50977, r50980, MPFR_RNDN);
        mpfr_set_d(r50982, z, MPFR_RNDN);
        mpfr_set_d(r50983, b, MPFR_RNDN);
        mpfr_mul(r50984, r50983, r50978, MPFR_RNDN);
        mpfr_set_d(r50985, y, MPFR_RNDN);
        mpfr_add(r50986, r50984, r50985, MPFR_RNDN);
        mpfr_mul(r50987, r50982, r50986, MPFR_RNDN);
        mpfr_add(r50988, r50981, r50987, MPFR_RNDN);
        return mpfr_get_d(r50988, MPFR_RNDN);
}