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

char *name = "Octave 3.8, oct_fill_randg";

double f_if(float a, float rand) {
        float r26004 = a;
        float r26005 = 1.0;
        float r26006 = 3.0;
        float r26007 = r26005 / r26006;
        float r26008 = r26004 - r26007;
        float r26009 = 1;
        float r26010 = 9;
        float r26011 = r26010 * r26008;
        float r26012 = sqrt(r26011);
        float r26013 = r26009 / r26012;
        float r26014 = rand;
        float r26015 = r26013 * r26014;
        float r26016 = r26009 + r26015;
        float r26017 = r26008 * r26016;
        return r26017;
}

double f_id(double a, double rand) {
        double r26018 = a;
        double r26019 = 1.0;
        double r26020 = 3.0;
        double r26021 = r26019 / r26020;
        double r26022 = r26018 - r26021;
        double r26023 = 1;
        double r26024 = 9;
        double r26025 = r26024 * r26022;
        double r26026 = sqrt(r26025);
        double r26027 = r26023 / r26026;
        double r26028 = rand;
        double r26029 = r26027 * r26028;
        double r26030 = r26023 + r26029;
        double r26031 = r26022 * r26030;
        return r26031;
}


double f_of(float a, float rand) {
        float r26032 = a;
        float r26033 = 1.0;
        float r26034 = 3.0;
        float r26035 = r26033 / r26034;
        float r26036 = r26032 - r26035;
        float r26037 = 1;
        float r26038 = 9;
        float r26039 = r26038 * r26036;
        float r26040 = 1/2;
        float r26041 = -r26040;
        float r26042 = pow(r26039, r26041);
        float r26043 = rand;
        float r26044 = r26042 * r26043;
        float r26045 = r26037 + r26044;
        float r26046 = r26036 * r26045;
        return r26046;
}

double f_od(double a, double rand) {
        double r26047 = a;
        double r26048 = 1.0;
        double r26049 = 3.0;
        double r26050 = r26048 / r26049;
        double r26051 = r26047 - r26050;
        double r26052 = 1;
        double r26053 = 9;
        double r26054 = r26053 * r26051;
        double r26055 = 1/2;
        double r26056 = -r26055;
        double r26057 = pow(r26054, r26056);
        double r26058 = rand;
        double r26059 = r26057 * r26058;
        double r26060 = r26052 + r26059;
        double r26061 = r26051 * r26060;
        return r26061;
}

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 r26062, r26063, r26064, r26065, r26066, r26067, r26068, r26069, r26070, r26071, r26072, r26073, r26074, r26075;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r26062);
        mpfr_init_set_str(r26063, "1.0", 10, MPFR_RNDN);
        mpfr_init_set_str(r26064, "3.0", 10, MPFR_RNDN);
        mpfr_init(r26065);
        mpfr_init(r26066);
        mpfr_init_set_str(r26067, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r26068, "9", 10, MPFR_RNDN);
        mpfr_init(r26069);
        mpfr_init(r26070);
        mpfr_init(r26071);
        mpfr_init(r26072);
        mpfr_init(r26073);
        mpfr_init(r26074);
        mpfr_init(r26075);
}

double f_im(double a, double rand) {
        mpfr_set_d(r26062, a, MPFR_RNDN);
        ;
        ;
        mpfr_div(r26065, r26063, r26064, MPFR_RNDN);
        mpfr_sub(r26066, r26062, r26065, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r26069, r26068, r26066, MPFR_RNDN);
        mpfr_sqrt(r26070, r26069, MPFR_RNDN);
        mpfr_div(r26071, r26067, r26070, MPFR_RNDN);
        mpfr_set_d(r26072, rand, MPFR_RNDN);
        mpfr_mul(r26073, r26071, r26072, MPFR_RNDN);
        mpfr_add(r26074, r26067, r26073, MPFR_RNDN);
        mpfr_mul(r26075, r26066, r26074, MPFR_RNDN);
        return mpfr_get_d(r26075, MPFR_RNDN);
}

static mpfr_t r26076, r26077, r26078, r26079, r26080, r26081, r26082, r26083, r26084, r26085, r26086, r26087, r26088, r26089, r26090;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r26076);
        mpfr_init_set_str(r26077, "1.0", 10, MPFR_RNDN);
        mpfr_init_set_str(r26078, "3.0", 10, MPFR_RNDN);
        mpfr_init(r26079);
        mpfr_init(r26080);
        mpfr_init_set_str(r26081, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r26082, "9", 10, MPFR_RNDN);
        mpfr_init(r26083);
        mpfr_init_set_str(r26084, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26085);
        mpfr_init(r26086);
        mpfr_init(r26087);
        mpfr_init(r26088);
        mpfr_init(r26089);
        mpfr_init(r26090);
}

double f_fm(double a, double rand) {
        mpfr_set_d(r26076, a, MPFR_RNDN);
        ;
        ;
        mpfr_div(r26079, r26077, r26078, MPFR_RNDN);
        mpfr_sub(r26080, r26076, r26079, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r26083, r26082, r26080, MPFR_RNDN);
        ;
        mpfr_neg(r26085, r26084, MPFR_RNDN);
        mpfr_pow(r26086, r26083, r26085, MPFR_RNDN);
        mpfr_set_d(r26087, rand, MPFR_RNDN);
        mpfr_mul(r26088, r26086, r26087, MPFR_RNDN);
        mpfr_add(r26089, r26081, r26088, MPFR_RNDN);
        mpfr_mul(r26090, r26080, r26089, MPFR_RNDN);
        return mpfr_get_d(r26090, MPFR_RNDN);
}

static mpfr_t r26091, r26092, r26093, r26094, r26095, r26096, r26097, r26098, r26099, r26100, r26101, r26102, r26103, r26104, r26105;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r26091);
        mpfr_init_set_str(r26092, "1.0", 10, MPFR_RNDN);
        mpfr_init_set_str(r26093, "3.0", 10, MPFR_RNDN);
        mpfr_init(r26094);
        mpfr_init(r26095);
        mpfr_init_set_str(r26096, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r26097, "9", 10, MPFR_RNDN);
        mpfr_init(r26098);
        mpfr_init_set_str(r26099, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26100);
        mpfr_init(r26101);
        mpfr_init(r26102);
        mpfr_init(r26103);
        mpfr_init(r26104);
        mpfr_init(r26105);
}

double f_dm(double a, double rand) {
        mpfr_set_d(r26091, a, MPFR_RNDN);
        ;
        ;
        mpfr_div(r26094, r26092, r26093, MPFR_RNDN);
        mpfr_sub(r26095, r26091, r26094, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r26098, r26097, r26095, MPFR_RNDN);
        ;
        mpfr_neg(r26100, r26099, MPFR_RNDN);
        mpfr_pow(r26101, r26098, r26100, MPFR_RNDN);
        mpfr_set_d(r26102, rand, MPFR_RNDN);
        mpfr_mul(r26103, r26101, r26102, MPFR_RNDN);
        mpfr_add(r26104, r26096, r26103, MPFR_RNDN);
        mpfr_mul(r26105, r26095, r26104, MPFR_RNDN);
        return mpfr_get_d(r26105, MPFR_RNDN);
}