#include #include #include #include #include char *name = "NMSE Section 6.1 mentioned, B"; double f_if(float a, float b) { float r20623 = atan2(1.0, 0.0); float r20624 = 2; float r20625 = r20623 / r20624; float r20626 = 1; float r20627 = b; float r20628 = r20627 * r20627; float r20629 = a; float r20630 = r20629 * r20629; float r20631 = r20628 - r20630; float r20632 = r20626 / r20631; float r20633 = r20625 * r20632; float r20634 = r20626 / r20629; float r20635 = r20626 / r20627; float r20636 = r20634 - r20635; float r20637 = r20633 * r20636; return r20637; } double f_id(double a, double b) { double r20638 = atan2(1.0, 0.0); double r20639 = 2; double r20640 = r20638 / r20639; double r20641 = 1; double r20642 = b; double r20643 = r20642 * r20642; double r20644 = a; double r20645 = r20644 * r20644; double r20646 = r20643 - r20645; double r20647 = r20641 / r20646; double r20648 = r20640 * r20647; double r20649 = r20641 / r20644; double r20650 = r20641 / r20642; double r20651 = r20649 - r20650; double r20652 = r20648 * r20651; return r20652; } double f_of(float a, float b) { float r20653 = 1; float r20654 = a; float r20655 = r20653 / r20654; float r20656 = b; float r20657 = r20653 / r20656; float r20658 = r20655 - r20657; float r20659 = r20656 - r20654; float r20660 = atan2(1.0, 0.0); float r20661 = 2; float r20662 = r20660 / r20661; float r20663 = r20659 / r20662; float r20664 = r20658 / r20663; float r20665 = r20656 + r20654; float r20666 = r20664 / r20665; return r20666; } double f_od(double a, double b) { double r20667 = 1; double r20668 = a; double r20669 = r20667 / r20668; double r20670 = b; double r20671 = r20667 / r20670; double r20672 = r20669 - r20671; double r20673 = r20670 - r20668; double r20674 = atan2(1.0, 0.0); double r20675 = 2; double r20676 = r20674 / r20675; double r20677 = r20673 / r20676; double r20678 = r20672 / r20677; double r20679 = r20670 + r20668; double r20680 = r20678 / r20679; return r20680; } 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 r20681, r20682, r20683, r20684, r20685, r20686, r20687, r20688, r20689, r20690, r20691, r20692, r20693, r20694, r20695; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r20681); mpfr_init_set_str(r20682, "2", 10, MPFR_RNDN); mpfr_init(r20683); mpfr_init_set_str(r20684, "1", 10, MPFR_RNDN); mpfr_init(r20685); mpfr_init(r20686); mpfr_init(r20687); mpfr_init(r20688); mpfr_init(r20689); mpfr_init(r20690); mpfr_init(r20691); mpfr_init(r20692); mpfr_init(r20693); mpfr_init(r20694); mpfr_init(r20695); } double f_im(double a, double b) { mpfr_const_pi(r20681, MPFR_RNDN); ; mpfr_div(r20683, r20681, r20682, MPFR_RNDN); ; mpfr_set_d(r20685, b, MPFR_RNDN); mpfr_mul(r20686, r20685, r20685, MPFR_RNDN); mpfr_set_d(r20687, a, MPFR_RNDN); mpfr_mul(r20688, r20687, r20687, MPFR_RNDN); mpfr_sub(r20689, r20686, r20688, MPFR_RNDN); mpfr_div(r20690, r20684, r20689, MPFR_RNDN); mpfr_mul(r20691, r20683, r20690, MPFR_RNDN); mpfr_div(r20692, r20684, r20687, MPFR_RNDN); mpfr_div(r20693, r20684, r20685, MPFR_RNDN); mpfr_sub(r20694, r20692, r20693, MPFR_RNDN); mpfr_mul(r20695, r20691, r20694, MPFR_RNDN); return mpfr_get_d(r20695, MPFR_RNDN); } static mpfr_t r20696, r20697, r20698, r20699, r20700, r20701, r20702, r20703, r20704, r20705, r20706, r20707, r20708, r20709; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20696, "1", 10, MPFR_RNDN); mpfr_init(r20697); mpfr_init(r20698); mpfr_init(r20699); mpfr_init(r20700); mpfr_init(r20701); mpfr_init(r20702); mpfr_init(r20703); mpfr_init_set_str(r20704, "2", 10, MPFR_RNDN); mpfr_init(r20705); mpfr_init(r20706); mpfr_init(r20707); mpfr_init(r20708); mpfr_init(r20709); } double f_fm(double a, double b) { ; mpfr_set_d(r20697, a, MPFR_RNDN); mpfr_div(r20698, r20696, r20697, MPFR_RNDN); mpfr_set_d(r20699, b, MPFR_RNDN); mpfr_div(r20700, r20696, r20699, MPFR_RNDN); mpfr_sub(r20701, r20698, r20700, MPFR_RNDN); mpfr_sub(r20702, r20699, r20697, MPFR_RNDN); mpfr_const_pi(r20703, MPFR_RNDN); ; mpfr_div(r20705, r20703, r20704, MPFR_RNDN); mpfr_div(r20706, r20702, r20705, MPFR_RNDN); mpfr_div(r20707, r20701, r20706, MPFR_RNDN); mpfr_add(r20708, r20699, r20697, MPFR_RNDN); mpfr_div(r20709, r20707, r20708, MPFR_RNDN); return mpfr_get_d(r20709, MPFR_RNDN); } static mpfr_t r20710, r20711, r20712, r20713, r20714, r20715, r20716, r20717, r20718, r20719, r20720, r20721, r20722, r20723; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20710, "1", 10, MPFR_RNDN); mpfr_init(r20711); mpfr_init(r20712); mpfr_init(r20713); mpfr_init(r20714); mpfr_init(r20715); mpfr_init(r20716); mpfr_init(r20717); mpfr_init_set_str(r20718, "2", 10, MPFR_RNDN); mpfr_init(r20719); mpfr_init(r20720); mpfr_init(r20721); mpfr_init(r20722); mpfr_init(r20723); } double f_dm(double a, double b) { ; mpfr_set_d(r20711, a, MPFR_RNDN); mpfr_div(r20712, r20710, r20711, MPFR_RNDN); mpfr_set_d(r20713, b, MPFR_RNDN); mpfr_div(r20714, r20710, r20713, MPFR_RNDN); mpfr_sub(r20715, r20712, r20714, MPFR_RNDN); mpfr_sub(r20716, r20713, r20711, MPFR_RNDN); mpfr_const_pi(r20717, MPFR_RNDN); ; mpfr_div(r20719, r20717, r20718, MPFR_RNDN); mpfr_div(r20720, r20716, r20719, MPFR_RNDN); mpfr_div(r20721, r20715, r20720, MPFR_RNDN); mpfr_add(r20722, r20713, r20711, MPFR_RNDN); mpfr_div(r20723, r20721, r20722, MPFR_RNDN); return mpfr_get_d(r20723, MPFR_RNDN); }