#include #include #include #include #include char *name = "Octave 3.8, jcobi/2"; double f_if(float alpha, float beta, float i) { float r9264 = alpha; float r9265 = beta; float r9266 = r9264 + r9265; float r9267 = r9265 - r9264; float r9268 = r9266 * r9267; float r9269 = 2.0f; float r9270 = i; float r9271 = r9269 * r9270; float r9272 = r9266 + r9271; float r9273 = r9268 / r9272; float r9274 = 2.0f; float r9275 = r9272 + r9274; float r9276 = r9273 / r9275; float r9277 = 1.0f; float r9278 = r9276 + r9277; float r9279 = r9278 / r9274; return r9279; } double f_id(double alpha, double beta, double i) { double r9280 = alpha; double r9281 = beta; double r9282 = r9280 + r9281; double r9283 = r9281 - r9280; double r9284 = r9282 * r9283; double r9285 = 2.0; double r9286 = i; double r9287 = r9285 * r9286; double r9288 = r9282 + r9287; double r9289 = r9284 / r9288; double r9290 = 2.0; double r9291 = r9288 + r9290; double r9292 = r9289 / r9291; double r9293 = 1.0; double r9294 = r9292 + r9293; double r9295 = r9294 / r9290; return r9295; } double f_of(float alpha, float beta, float i) { float r9296 = alpha; float r9297 = beta; float r9298 = r9296 + r9297; float r9299 = 2.0f; float r9300 = r9297 + r9299; float r9301 = 2.0f; float r9302 = i; float r9303 = fma(r9301, r9302, r9296); float r9304 = r9300 + r9303; float r9305 = sqrt(r9304); float r9306 = r9298 / r9305; float r9307 = 1.0f; float r9308 = r9307 / r9305; float r9309 = r9306 * r9308; float r9310 = r9309 * r9309; float r9311 = r9310 * r9309; float r9312 = cbrt(r9311); float r9313 = r9297 - r9296; float r9314 = r9297 + r9303; float r9315 = r9313 / r9314; float r9316 = 1.0f; float r9317 = fma(r9312, r9315, r9316); float r9318 = log(r9317); float r9319 = r9318 * r9318; float r9320 = r9319 * r9318; float r9321 = cbrt(r9320); float r9322 = exp(r9321); float r9323 = r9322 / r9299; return r9323; } double f_od(double alpha, double beta, double i) { double r9324 = alpha; double r9325 = beta; double r9326 = r9324 + r9325; double r9327 = 2.0; double r9328 = r9325 + r9327; double r9329 = 2.0; double r9330 = i; double r9331 = fma(r9329, r9330, r9324); double r9332 = r9328 + r9331; double r9333 = sqrt(r9332); double r9334 = r9326 / r9333; double r9335 = 1.0; double r9336 = r9335 / r9333; double r9337 = r9334 * r9336; double r9338 = r9337 * r9337; double r9339 = r9338 * r9337; double r9340 = cbrt(r9339); double r9341 = r9325 - r9324; double r9342 = r9325 + r9331; double r9343 = r9341 / r9342; double r9344 = 1.0; double r9345 = fma(r9340, r9343, r9344); double r9346 = log(r9345); double r9347 = r9346 * r9346; double r9348 = r9347 * r9346; double r9349 = cbrt(r9348); double r9350 = exp(r9349); double r9351 = r9350 / r9327; return r9351; } 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 r9352, r9353, r9354, r9355, r9356, r9357, r9358, r9359, r9360, r9361, r9362, r9363, r9364, r9365, r9366, r9367; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r9352); mpfr_init(r9353); mpfr_init(r9354); mpfr_init(r9355); mpfr_init(r9356); mpfr_init_set_str(r9357, "2", 10, MPFR_RNDN); mpfr_init(r9358); mpfr_init(r9359); mpfr_init(r9360); mpfr_init(r9361); mpfr_init_set_str(r9362, "2.0", 10, MPFR_RNDN); mpfr_init(r9363); mpfr_init(r9364); mpfr_init_set_str(r9365, "1.0", 10, MPFR_RNDN); mpfr_init(r9366); mpfr_init(r9367); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r9352, alpha, MPFR_RNDN); mpfr_set_d(r9353, beta, MPFR_RNDN); mpfr_add(r9354, r9352, r9353, MPFR_RNDN); mpfr_sub(r9355, r9353, r9352, MPFR_RNDN); mpfr_mul(r9356, r9354, r9355, MPFR_RNDN); ; mpfr_set_d(r9358, i, MPFR_RNDN); mpfr_mul(r9359, r9357, r9358, MPFR_RNDN); mpfr_add(r9360, r9354, r9359, MPFR_RNDN); mpfr_div(r9361, r9356, r9360, MPFR_RNDN); ; mpfr_add(r9363, r9360, r9362, MPFR_RNDN); mpfr_div(r9364, r9361, r9363, MPFR_RNDN); ; mpfr_add(r9366, r9364, r9365, MPFR_RNDN); mpfr_div(r9367, r9366, r9362, MPFR_RNDN); return mpfr_get_d(r9367, MPFR_RNDN); } static mpfr_t r9368, r9369, r9370, r9371, r9372, r9373, r9374, r9375, r9376, r9377, r9378, r9379, r9380, r9381, r9382, r9383, r9384, r9385, r9386, r9387, r9388, r9389, r9390, r9391, r9392, r9393, r9394, r9395; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r9368); mpfr_init(r9369); mpfr_init(r9370); mpfr_init_set_str(r9371, "2.0", 10, MPFR_RNDN); mpfr_init(r9372); mpfr_init_set_str(r9373, "2", 10, MPFR_RNDN); mpfr_init(r9374); mpfr_init(r9375); mpfr_init(r9376); mpfr_init(r9377); mpfr_init(r9378); mpfr_init_set_str(r9379, "1", 10, MPFR_RNDN); mpfr_init(r9380); mpfr_init(r9381); mpfr_init(r9382); mpfr_init(r9383); mpfr_init(r9384); mpfr_init(r9385); mpfr_init(r9386); mpfr_init(r9387); mpfr_init_set_str(r9388, "1.0", 10, MPFR_RNDN); mpfr_init(r9389); mpfr_init(r9390); mpfr_init(r9391); mpfr_init(r9392); mpfr_init(r9393); mpfr_init(r9394); mpfr_init(r9395); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r9368, alpha, MPFR_RNDN); mpfr_set_d(r9369, beta, MPFR_RNDN); mpfr_add(r9370, r9368, r9369, MPFR_RNDN); ; mpfr_add(r9372, r9369, r9371, MPFR_RNDN); ; mpfr_set_d(r9374, i, MPFR_RNDN); mpfr_fma(r9375, r9373, r9374, r9368, MPFR_RNDN); mpfr_add(r9376, r9372, r9375, MPFR_RNDN); mpfr_sqrt(r9377, r9376, MPFR_RNDN); mpfr_div(r9378, r9370, r9377, MPFR_RNDN); ; mpfr_div(r9380, r9379, r9377, MPFR_RNDN); mpfr_mul(r9381, r9378, r9380, MPFR_RNDN); mpfr_mul(r9382, r9381, r9381, MPFR_RNDN); mpfr_mul(r9383, r9382, r9381, MPFR_RNDN); mpfr_cbrt(r9384, r9383, MPFR_RNDN); mpfr_sub(r9385, r9369, r9368, MPFR_RNDN); mpfr_add(r9386, r9369, r9375, MPFR_RNDN); mpfr_div(r9387, r9385, r9386, MPFR_RNDN); ; mpfr_fma(r9389, r9384, r9387, r9388, MPFR_RNDN); mpfr_log(r9390, r9389, MPFR_RNDN); mpfr_mul(r9391, r9390, r9390, MPFR_RNDN); mpfr_mul(r9392, r9391, r9390, MPFR_RNDN); mpfr_cbrt(r9393, r9392, MPFR_RNDN); mpfr_exp(r9394, r9393, MPFR_RNDN); mpfr_div(r9395, r9394, r9371, MPFR_RNDN); return mpfr_get_d(r9395, MPFR_RNDN); } static mpfr_t r9396, r9397, r9398, r9399, r9400, r9401, r9402, r9403, r9404, r9405, r9406, r9407, r9408, r9409, r9410, r9411, r9412, r9413, r9414, r9415, r9416, r9417, r9418, r9419, r9420, r9421, r9422, r9423; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r9396); mpfr_init(r9397); mpfr_init(r9398); mpfr_init_set_str(r9399, "2.0", 10, MPFR_RNDN); mpfr_init(r9400); mpfr_init_set_str(r9401, "2", 10, MPFR_RNDN); mpfr_init(r9402); mpfr_init(r9403); mpfr_init(r9404); mpfr_init(r9405); mpfr_init(r9406); mpfr_init_set_str(r9407, "1", 10, MPFR_RNDN); mpfr_init(r9408); mpfr_init(r9409); mpfr_init(r9410); mpfr_init(r9411); mpfr_init(r9412); mpfr_init(r9413); mpfr_init(r9414); mpfr_init(r9415); mpfr_init_set_str(r9416, "1.0", 10, MPFR_RNDN); mpfr_init(r9417); mpfr_init(r9418); mpfr_init(r9419); mpfr_init(r9420); mpfr_init(r9421); mpfr_init(r9422); mpfr_init(r9423); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r9396, alpha, MPFR_RNDN); mpfr_set_d(r9397, beta, MPFR_RNDN); mpfr_add(r9398, r9396, r9397, MPFR_RNDN); ; mpfr_add(r9400, r9397, r9399, MPFR_RNDN); ; mpfr_set_d(r9402, i, MPFR_RNDN); mpfr_fma(r9403, r9401, r9402, r9396, MPFR_RNDN); mpfr_add(r9404, r9400, r9403, MPFR_RNDN); mpfr_sqrt(r9405, r9404, MPFR_RNDN); mpfr_div(r9406, r9398, r9405, MPFR_RNDN); ; mpfr_div(r9408, r9407, r9405, MPFR_RNDN); mpfr_mul(r9409, r9406, r9408, MPFR_RNDN); mpfr_mul(r9410, r9409, r9409, MPFR_RNDN); mpfr_mul(r9411, r9410, r9409, MPFR_RNDN); mpfr_cbrt(r9412, r9411, MPFR_RNDN); mpfr_sub(r9413, r9397, r9396, MPFR_RNDN); mpfr_add(r9414, r9397, r9403, MPFR_RNDN); mpfr_div(r9415, r9413, r9414, MPFR_RNDN); ; mpfr_fma(r9417, r9412, r9415, r9416, MPFR_RNDN); mpfr_log(r9418, r9417, MPFR_RNDN); mpfr_mul(r9419, r9418, r9418, MPFR_RNDN); mpfr_mul(r9420, r9419, r9418, MPFR_RNDN); mpfr_cbrt(r9421, r9420, MPFR_RNDN); mpfr_exp(r9422, r9421, MPFR_RNDN); mpfr_div(r9423, r9422, r9399, MPFR_RNDN); return mpfr_get_d(r9423, MPFR_RNDN); }