#include #include #include #include #include char *name = "Octave 3.8, oct_fill_randg"; double f_if(float a, float rand) { float r16293 = a; float r16294 = 1.0f; float r16295 = 3.0f; float r16296 = r16294 / r16295; float r16297 = r16293 - r16296; float r16298 = 1.0f; float r16299 = 9.0f; float r16300 = r16299 * r16297; float r16301 = sqrt(r16300); float r16302 = r16298 / r16301; float r16303 = rand; float r16304 = r16302 * r16303; float r16305 = r16298 + r16304; float r16306 = r16297 * r16305; return r16306; } double f_id(double a, double rand) { double r16307 = a; double r16308 = 1.0; double r16309 = 3.0; double r16310 = r16308 / r16309; double r16311 = r16307 - r16310; double r16312 = 1.0; double r16313 = 9.0; double r16314 = r16313 * r16311; double r16315 = sqrt(r16314); double r16316 = r16312 / r16315; double r16317 = rand; double r16318 = r16316 * r16317; double r16319 = r16312 + r16318; double r16320 = r16311 * r16319; return r16320; } double f_of(float a, float rand) { float r16321 = a; float r16322 = 1.0f; float r16323 = 3.0f; float r16324 = r16322 / r16323; float r16325 = r16321 - r16324; float r16326 = 1.0f; float r16327 = 9.0f; float r16328 = sqrt(r16327); float r16329 = sqrt(r16325); float r16330 = r16328 * r16329; float r16331 = r16326 / r16330; float r16332 = rand; float r16333 = r16331 * r16332; float r16334 = r16326 + r16333; float r16335 = r16325 * r16334; return r16335; } double f_od(double a, double rand) { double r16336 = a; double r16337 = 1.0; double r16338 = 3.0; double r16339 = r16337 / r16338; double r16340 = r16336 - r16339; double r16341 = 1.0; double r16342 = 9.0; double r16343 = sqrt(r16342); double r16344 = sqrt(r16340); double r16345 = r16343 * r16344; double r16346 = r16341 / r16345; double r16347 = rand; double r16348 = r16346 * r16347; double r16349 = r16341 + r16348; double r16350 = r16340 * r16349; return r16350; } 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 r16351, r16352, r16353, r16354, r16355, r16356, r16357, r16358, r16359, r16360, r16361, r16362, r16363, r16364; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16351); mpfr_init_set_str(r16352, "1.0", 10, MPFR_RNDN); mpfr_init_set_str(r16353, "3.0", 10, MPFR_RNDN); mpfr_init(r16354); mpfr_init(r16355); mpfr_init_set_str(r16356, "1", 10, MPFR_RNDN); mpfr_init_set_str(r16357, "9", 10, MPFR_RNDN); mpfr_init(r16358); mpfr_init(r16359); mpfr_init(r16360); mpfr_init(r16361); mpfr_init(r16362); mpfr_init(r16363); mpfr_init(r16364); } double f_im(double a, double rand) { mpfr_set_d(r16351, a, MPFR_RNDN); ; ; mpfr_div(r16354, r16352, r16353, MPFR_RNDN); mpfr_sub(r16355, r16351, r16354, MPFR_RNDN); ; ; mpfr_mul(r16358, r16357, r16355, MPFR_RNDN); mpfr_sqrt(r16359, r16358, MPFR_RNDN); mpfr_div(r16360, r16356, r16359, MPFR_RNDN); mpfr_set_d(r16361, rand, MPFR_RNDN); mpfr_mul(r16362, r16360, r16361, MPFR_RNDN); mpfr_add(r16363, r16356, r16362, MPFR_RNDN); mpfr_mul(r16364, r16355, r16363, MPFR_RNDN); return mpfr_get_d(r16364, MPFR_RNDN); } static mpfr_t r16365, r16366, r16367, r16368, r16369, r16370, r16371, r16372, r16373, r16374, r16375, r16376, r16377, r16378, r16379; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r16365); mpfr_init_set_str(r16366, "1.0", 10, MPFR_RNDN); mpfr_init_set_str(r16367, "3.0", 10, MPFR_RNDN); mpfr_init(r16368); mpfr_init(r16369); mpfr_init_set_str(r16370, "1", 10, MPFR_RNDN); mpfr_init_set_str(r16371, "9", 10, MPFR_RNDN); mpfr_init(r16372); mpfr_init(r16373); mpfr_init(r16374); mpfr_init(r16375); mpfr_init(r16376); mpfr_init(r16377); mpfr_init(r16378); mpfr_init(r16379); } double f_fm(double a, double rand) { mpfr_set_d(r16365, a, MPFR_RNDN); ; ; mpfr_div(r16368, r16366, r16367, MPFR_RNDN); mpfr_sub(r16369, r16365, r16368, MPFR_RNDN); ; ; mpfr_sqrt(r16372, r16371, MPFR_RNDN); mpfr_sqrt(r16373, r16369, MPFR_RNDN); mpfr_mul(r16374, r16372, r16373, MPFR_RNDN); mpfr_div(r16375, r16370, r16374, MPFR_RNDN); mpfr_set_d(r16376, rand, MPFR_RNDN); mpfr_mul(r16377, r16375, r16376, MPFR_RNDN); mpfr_add(r16378, r16370, r16377, MPFR_RNDN); mpfr_mul(r16379, r16369, r16378, MPFR_RNDN); return mpfr_get_d(r16379, MPFR_RNDN); } static mpfr_t r16380, r16381, r16382, r16383, r16384, r16385, r16386, r16387, r16388, r16389, r16390, r16391, r16392, r16393, r16394; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r16380); mpfr_init_set_str(r16381, "1.0", 10, MPFR_RNDN); mpfr_init_set_str(r16382, "3.0", 10, MPFR_RNDN); mpfr_init(r16383); mpfr_init(r16384); mpfr_init_set_str(r16385, "1", 10, MPFR_RNDN); mpfr_init_set_str(r16386, "9", 10, MPFR_RNDN); mpfr_init(r16387); mpfr_init(r16388); mpfr_init(r16389); mpfr_init(r16390); mpfr_init(r16391); mpfr_init(r16392); mpfr_init(r16393); mpfr_init(r16394); } double f_dm(double a, double rand) { mpfr_set_d(r16380, a, MPFR_RNDN); ; ; mpfr_div(r16383, r16381, r16382, MPFR_RNDN); mpfr_sub(r16384, r16380, r16383, MPFR_RNDN); ; ; mpfr_sqrt(r16387, r16386, MPFR_RNDN); mpfr_sqrt(r16388, r16384, MPFR_RNDN); mpfr_mul(r16389, r16387, r16388, MPFR_RNDN); mpfr_div(r16390, r16385, r16389, MPFR_RNDN); mpfr_set_d(r16391, rand, MPFR_RNDN); mpfr_mul(r16392, r16390, r16391, MPFR_RNDN); mpfr_add(r16393, r16385, r16392, MPFR_RNDN); mpfr_mul(r16394, r16384, r16393, MPFR_RNDN); return mpfr_get_d(r16394, MPFR_RNDN); }