#include #include #include #include #include char *name = "Linear.V4:$cdot from linear-1.19.1.3"; double f_if(float x, float y, float z, float t, float a, float b, float c, float i) { float r22761 = x; float r22762 = y; float r22763 = r22761 * r22762; float r22764 = z; float r22765 = t; float r22766 = r22764 * r22765; float r22767 = r22763 + r22766; float r22768 = a; float r22769 = b; float r22770 = r22768 * r22769; float r22771 = r22767 + r22770; float r22772 = c; float r22773 = i; float r22774 = r22772 * r22773; float r22775 = r22771 + r22774; return r22775; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r22776 = x; double r22777 = y; double r22778 = r22776 * r22777; double r22779 = z; double r22780 = t; double r22781 = r22779 * r22780; double r22782 = r22778 + r22781; double r22783 = a; double r22784 = b; double r22785 = r22783 * r22784; double r22786 = r22782 + r22785; double r22787 = c; double r22788 = i; double r22789 = r22787 * r22788; double r22790 = r22786 + r22789; return r22790; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r22791 = c; float r22792 = i; float r22793 = x; float r22794 = y; float r22795 = r22793 * r22794; float r22796 = fma(r22791, r22792, r22795); float r22797 = b; float r22798 = a; float r22799 = t; float r22800 = z; float r22801 = r22799 * r22800; float r22802 = fma(r22797, r22798, r22801); float r22803 = r22796 + r22802; return r22803; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r22804 = c; double r22805 = i; double r22806 = x; double r22807 = y; double r22808 = r22806 * r22807; double r22809 = fma(r22804, r22805, r22808); double r22810 = b; double r22811 = a; double r22812 = t; double r22813 = z; double r22814 = r22812 * r22813; double r22815 = fma(r22810, r22811, r22814); double r22816 = r22809 + r22815; return r22816; } 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 r22817, r22818, r22819, r22820, r22821, r22822, r22823, r22824, r22825, r22826, r22827, r22828, r22829, r22830, r22831; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22817); mpfr_init(r22818); mpfr_init(r22819); mpfr_init(r22820); mpfr_init(r22821); mpfr_init(r22822); mpfr_init(r22823); mpfr_init(r22824); mpfr_init(r22825); mpfr_init(r22826); mpfr_init(r22827); mpfr_init(r22828); mpfr_init(r22829); mpfr_init(r22830); mpfr_init(r22831); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r22817, x, MPFR_RNDN); mpfr_set_d(r22818, y, MPFR_RNDN); mpfr_mul(r22819, r22817, r22818, MPFR_RNDN); mpfr_set_d(r22820, z, MPFR_RNDN); mpfr_set_d(r22821, t, MPFR_RNDN); mpfr_mul(r22822, r22820, r22821, MPFR_RNDN); mpfr_add(r22823, r22819, r22822, MPFR_RNDN); mpfr_set_d(r22824, a, MPFR_RNDN); mpfr_set_d(r22825, b, MPFR_RNDN); mpfr_mul(r22826, r22824, r22825, MPFR_RNDN); mpfr_add(r22827, r22823, r22826, MPFR_RNDN); mpfr_set_d(r22828, c, MPFR_RNDN); mpfr_set_d(r22829, i, MPFR_RNDN); mpfr_mul(r22830, r22828, r22829, MPFR_RNDN); mpfr_add(r22831, r22827, r22830, MPFR_RNDN); return mpfr_get_d(r22831, MPFR_RNDN); } static mpfr_t r22832, r22833, r22834, r22835, r22836, r22837, r22838, r22839, r22840, r22841, r22842, r22843, r22844; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22832); mpfr_init(r22833); mpfr_init(r22834); mpfr_init(r22835); mpfr_init(r22836); mpfr_init(r22837); mpfr_init(r22838); mpfr_init(r22839); mpfr_init(r22840); mpfr_init(r22841); mpfr_init(r22842); mpfr_init(r22843); mpfr_init(r22844); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r22832, c, MPFR_RNDN); mpfr_set_d(r22833, i, MPFR_RNDN); mpfr_set_d(r22834, x, MPFR_RNDN); mpfr_set_d(r22835, y, MPFR_RNDN); mpfr_mul(r22836, r22834, r22835, MPFR_RNDN); mpfr_fma(r22837, r22832, r22833, r22836, MPFR_RNDN); mpfr_set_d(r22838, b, MPFR_RNDN); mpfr_set_d(r22839, a, MPFR_RNDN); mpfr_set_d(r22840, t, MPFR_RNDN); mpfr_set_d(r22841, z, MPFR_RNDN); mpfr_mul(r22842, r22840, r22841, MPFR_RNDN); mpfr_fma(r22843, r22838, r22839, r22842, MPFR_RNDN); mpfr_add(r22844, r22837, r22843, MPFR_RNDN); return mpfr_get_d(r22844, MPFR_RNDN); } static mpfr_t r22845, r22846, r22847, r22848, r22849, r22850, r22851, r22852, r22853, r22854, r22855, r22856, r22857; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22845); mpfr_init(r22846); mpfr_init(r22847); mpfr_init(r22848); mpfr_init(r22849); mpfr_init(r22850); mpfr_init(r22851); mpfr_init(r22852); mpfr_init(r22853); mpfr_init(r22854); mpfr_init(r22855); mpfr_init(r22856); mpfr_init(r22857); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r22845, c, MPFR_RNDN); mpfr_set_d(r22846, i, MPFR_RNDN); mpfr_set_d(r22847, x, MPFR_RNDN); mpfr_set_d(r22848, y, MPFR_RNDN); mpfr_mul(r22849, r22847, r22848, MPFR_RNDN); mpfr_fma(r22850, r22845, r22846, r22849, MPFR_RNDN); mpfr_set_d(r22851, b, MPFR_RNDN); mpfr_set_d(r22852, a, MPFR_RNDN); mpfr_set_d(r22853, t, MPFR_RNDN); mpfr_set_d(r22854, z, MPFR_RNDN); mpfr_mul(r22855, r22853, r22854, MPFR_RNDN); mpfr_fma(r22856, r22851, r22852, r22855, MPFR_RNDN); mpfr_add(r22857, r22850, r22856, MPFR_RNDN); return mpfr_get_d(r22857, MPFR_RNDN); }