#include #include #include #include #include char *name = "Octave 3.8, jcobi/4, as called"; double f_if(float i) { float r26217 = i; float r26218 = r26217 * r26217; float r26219 = r26218 * r26218; float r26220 = 2; float r26221 = r26220 * r26217; float r26222 = r26221 * r26221; float r26223 = r26219 / r26222; float r26224 = 1.0; float r26225 = r26222 - r26224; float r26226 = r26223 / r26225; return r26226; } double f_id(double i) { double r26227 = i; double r26228 = r26227 * r26227; double r26229 = r26228 * r26228; double r26230 = 2; double r26231 = r26230 * r26227; double r26232 = r26231 * r26231; double r26233 = r26229 / r26232; double r26234 = 1.0; double r26235 = r26232 - r26234; double r26236 = r26233 / r26235; return r26236; } double f_of(float i) { float r26237 = i; float r26238 = 222.38347807310183; bool r26239 = r26237 <= r26238; float r26240 = 2; float r26241 = r26237 / r26240; float r26242 = r26241 * r26241; float r26243 = r26240 * r26237; float r26244 = 1.0; float r26245 = -r26244; float r26246 = fma(r26243, r26243, r26245); float r26247 = r26242 / r26246; float r26248 = expm1(r26247); float r26249 = log1p(r26248); float r26250 = 0.00390625; float r26251 = 4; float r26252 = pow(r26237, r26251); float r26253 = r26250 / r26252; float r26254 = 0.015625; float r26255 = r26254 / r26237; float r26256 = r26255 / r26237; float r26257 = r26253 + r26256; float r26258 = 1/16; float r26259 = r26257 + r26258; float r26260 = r26239 ? r26249 : r26259; return r26260; } double f_od(double i) { double r26261 = i; double r26262 = 222.38347807310183; bool r26263 = r26261 <= r26262; double r26264 = 2; double r26265 = r26261 / r26264; double r26266 = r26265 * r26265; double r26267 = r26264 * r26261; double r26268 = 1.0; double r26269 = -r26268; double r26270 = fma(r26267, r26267, r26269); double r26271 = r26266 / r26270; double r26272 = expm1(r26271); double r26273 = log1p(r26272); double r26274 = 0.00390625; double r26275 = 4; double r26276 = pow(r26261, r26275); double r26277 = r26274 / r26276; double r26278 = 0.015625; double r26279 = r26278 / r26261; double r26280 = r26279 / r26261; double r26281 = r26277 + r26280; double r26282 = 1/16; double r26283 = r26281 + r26282; double r26284 = r26263 ? r26273 : r26283; return r26284; } 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 r26285, r26286, r26287, r26288, r26289, r26290, r26291, r26292, r26293, r26294; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r26285); mpfr_init(r26286); mpfr_init(r26287); mpfr_init_set_str(r26288, "2", 10, MPFR_RNDN); mpfr_init(r26289); mpfr_init(r26290); mpfr_init(r26291); mpfr_init_set_str(r26292, "1.0", 10, MPFR_RNDN); mpfr_init(r26293); mpfr_init(r26294); } double f_im(double i) { mpfr_set_d(r26285, i, MPFR_RNDN); mpfr_mul(r26286, r26285, r26285, MPFR_RNDN); mpfr_mul(r26287, r26286, r26286, MPFR_RNDN); ; mpfr_mul(r26289, r26288, r26285, MPFR_RNDN); mpfr_mul(r26290, r26289, r26289, MPFR_RNDN); mpfr_div(r26291, r26287, r26290, MPFR_RNDN); ; mpfr_sub(r26293, r26290, r26292, MPFR_RNDN); mpfr_div(r26294, r26291, r26293, MPFR_RNDN); return mpfr_get_d(r26294, MPFR_RNDN); } static mpfr_t r26295, r26296, r26297, r26298, r26299, r26300, r26301, r26302, r26303, r26304, r26305, r26306, r26307, r26308, r26309, r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r26295); mpfr_init_set_str(r26296, "222.38347807310183", 10, MPFR_RNDN); mpfr_init(r26297); mpfr_init_set_str(r26298, "2", 10, MPFR_RNDN); mpfr_init(r26299); mpfr_init(r26300); mpfr_init(r26301); mpfr_init_set_str(r26302, "1.0", 10, MPFR_RNDN); mpfr_init(r26303); mpfr_init(r26304); mpfr_init(r26305); mpfr_init(r26306); mpfr_init(r26307); mpfr_init_set_str(r26308, "0.00390625", 10, MPFR_RNDN); mpfr_init_set_str(r26309, "4", 10, MPFR_RNDN); mpfr_init(r26310); mpfr_init(r26311); mpfr_init_set_str(r26312, "0.015625", 10, MPFR_RNDN); mpfr_init(r26313); mpfr_init(r26314); mpfr_init(r26315); mpfr_init_set_str(r26316, "1/16", 10, MPFR_RNDN); mpfr_init(r26317); mpfr_init(r26318); } double f_fm(double i) { mpfr_set_d(r26295, i, MPFR_RNDN); ; mpfr_set_si(r26297, mpfr_cmp(r26295, r26296) <= 0, MPFR_RNDN); ; mpfr_div(r26299, r26295, r26298, MPFR_RNDN); mpfr_mul(r26300, r26299, r26299, MPFR_RNDN); mpfr_mul(r26301, r26298, r26295, MPFR_RNDN); ; mpfr_neg(r26303, r26302, MPFR_RNDN); mpfr_fma(r26304, r26301, r26301, r26303, MPFR_RNDN); mpfr_div(r26305, r26300, r26304, MPFR_RNDN); mpfr_expm1(r26306, r26305, MPFR_RNDN); mpfr_log1p(r26307, r26306, MPFR_RNDN); ; ; mpfr_pow(r26310, r26295, r26309, MPFR_RNDN); mpfr_div(r26311, r26308, r26310, MPFR_RNDN); ; mpfr_div(r26313, r26312, r26295, MPFR_RNDN); mpfr_div(r26314, r26313, r26295, MPFR_RNDN); mpfr_add(r26315, r26311, r26314, MPFR_RNDN); ; mpfr_add(r26317, r26315, r26316, MPFR_RNDN); if (mpfr_get_si(r26297, MPFR_RNDN)) { mpfr_set(r26318, r26307, MPFR_RNDN); } else { mpfr_set(r26318, r26317, MPFR_RNDN); }; return mpfr_get_d(r26318, MPFR_RNDN); } static mpfr_t r26319, r26320, r26321, r26322, r26323, r26324, r26325, r26326, r26327, r26328, r26329, r26330, r26331, r26332, r26333, r26334, r26335, r26336, r26337, r26338, r26339, r26340, r26341, r26342; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r26319); mpfr_init_set_str(r26320, "222.38347807310183", 10, MPFR_RNDN); mpfr_init(r26321); mpfr_init_set_str(r26322, "2", 10, MPFR_RNDN); mpfr_init(r26323); mpfr_init(r26324); mpfr_init(r26325); mpfr_init_set_str(r26326, "1.0", 10, MPFR_RNDN); mpfr_init(r26327); mpfr_init(r26328); mpfr_init(r26329); mpfr_init(r26330); mpfr_init(r26331); mpfr_init_set_str(r26332, "0.00390625", 10, MPFR_RNDN); mpfr_init_set_str(r26333, "4", 10, MPFR_RNDN); mpfr_init(r26334); mpfr_init(r26335); mpfr_init_set_str(r26336, "0.015625", 10, MPFR_RNDN); mpfr_init(r26337); mpfr_init(r26338); mpfr_init(r26339); mpfr_init_set_str(r26340, "1/16", 10, MPFR_RNDN); mpfr_init(r26341); mpfr_init(r26342); } double f_dm(double i) { mpfr_set_d(r26319, i, MPFR_RNDN); ; mpfr_set_si(r26321, mpfr_cmp(r26319, r26320) <= 0, MPFR_RNDN); ; mpfr_div(r26323, r26319, r26322, MPFR_RNDN); mpfr_mul(r26324, r26323, r26323, MPFR_RNDN); mpfr_mul(r26325, r26322, r26319, MPFR_RNDN); ; mpfr_neg(r26327, r26326, MPFR_RNDN); mpfr_fma(r26328, r26325, r26325, r26327, MPFR_RNDN); mpfr_div(r26329, r26324, r26328, MPFR_RNDN); mpfr_expm1(r26330, r26329, MPFR_RNDN); mpfr_log1p(r26331, r26330, MPFR_RNDN); ; ; mpfr_pow(r26334, r26319, r26333, MPFR_RNDN); mpfr_div(r26335, r26332, r26334, MPFR_RNDN); ; mpfr_div(r26337, r26336, r26319, MPFR_RNDN); mpfr_div(r26338, r26337, r26319, MPFR_RNDN); mpfr_add(r26339, r26335, r26338, MPFR_RNDN); ; mpfr_add(r26341, r26339, r26340, MPFR_RNDN); if (mpfr_get_si(r26321, MPFR_RNDN)) { mpfr_set(r26342, r26331, MPFR_RNDN); } else { mpfr_set(r26342, r26341, MPFR_RNDN); }; return mpfr_get_d(r26342, MPFR_RNDN); }