#include #include #include #include #include char *name = "Harley's example"; double f_if(float c_p, float c_n, float t, float s) { float r24191 = 1; float r24192 = s; float r24193 = -r24192; float r24194 = exp(r24193); float r24195 = r24191 + r24194; float r24196 = r24191 / r24195; float r24197 = c_p; float r24198 = pow(r24196, r24197); float r24199 = r24191 - r24196; float r24200 = c_n; float r24201 = pow(r24199, r24200); float r24202 = r24198 * r24201; float r24203 = t; float r24204 = -r24203; float r24205 = exp(r24204); float r24206 = r24191 + r24205; float r24207 = r24191 / r24206; float r24208 = pow(r24207, r24197); float r24209 = r24191 - r24207; float r24210 = pow(r24209, r24200); float r24211 = r24208 * r24210; float r24212 = r24202 / r24211; return r24212; } double f_id(double c_p, double c_n, double t, double s) { double r24213 = 1; double r24214 = s; double r24215 = -r24214; double r24216 = exp(r24215); double r24217 = r24213 + r24216; double r24218 = r24213 / r24217; double r24219 = c_p; double r24220 = pow(r24218, r24219); double r24221 = r24213 - r24218; double r24222 = c_n; double r24223 = pow(r24221, r24222); double r24224 = r24220 * r24223; double r24225 = t; double r24226 = -r24225; double r24227 = exp(r24226); double r24228 = r24213 + r24227; double r24229 = r24213 / r24228; double r24230 = pow(r24229, r24219); double r24231 = r24213 - r24229; double r24232 = pow(r24231, r24222); double r24233 = r24230 * r24232; double r24234 = r24224 / r24233; return r24234; } double f_of(float c_p, float c_n, float t, float s) { float r24235 = 1; float r24236 = s; float r24237 = -r24236; float r24238 = exp(r24237); float r24239 = r24235 + r24238; float r24240 = r24235 / r24239; float r24241 = r24235 - r24240; float r24242 = c_n; float r24243 = pow(r24241, r24242); float r24244 = t; float r24245 = -r24244; float r24246 = exp(r24245); float r24247 = r24235 + r24246; float r24248 = r24235 / r24247; float r24249 = r24235 - r24248; float r24250 = pow(r24249, r24242); float r24251 = r24243 / r24250; float r24252 = c_p; float r24253 = pow(r24240, r24252); float r24254 = 1/2; float r24255 = log(r24254); float r24256 = r24244 * r24254; float r24257 = r24255 + r24256; float r24258 = r24257 * r24252; float r24259 = r24258 + r24235; float r24260 = r24253 / r24259; float r24261 = r24251 * r24260; return r24261; } double f_od(double c_p, double c_n, double t, double s) { double r24262 = 1; double r24263 = s; double r24264 = -r24263; double r24265 = exp(r24264); double r24266 = r24262 + r24265; double r24267 = r24262 / r24266; double r24268 = r24262 - r24267; double r24269 = c_n; double r24270 = pow(r24268, r24269); double r24271 = t; double r24272 = -r24271; double r24273 = exp(r24272); double r24274 = r24262 + r24273; double r24275 = r24262 / r24274; double r24276 = r24262 - r24275; double r24277 = pow(r24276, r24269); double r24278 = r24270 / r24277; double r24279 = c_p; double r24280 = pow(r24267, r24279); double r24281 = 1/2; double r24282 = log(r24281); double r24283 = r24271 * r24281; double r24284 = r24282 + r24283; double r24285 = r24284 * r24279; double r24286 = r24285 + r24262; double r24287 = r24280 / r24286; double r24288 = r24278 * r24287; return r24288; } 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 r24289, r24290, r24291, r24292, r24293, r24294, r24295, r24296, r24297, r24298, r24299, r24300, r24301, r24302, r24303, r24304, r24305, r24306, r24307, r24308, r24309, r24310; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r24289, "1", 10, MPFR_RNDN); mpfr_init(r24290); mpfr_init(r24291); mpfr_init(r24292); mpfr_init(r24293); mpfr_init(r24294); mpfr_init(r24295); mpfr_init(r24296); mpfr_init(r24297); mpfr_init(r24298); mpfr_init(r24299); mpfr_init(r24300); mpfr_init(r24301); mpfr_init(r24302); mpfr_init(r24303); mpfr_init(r24304); mpfr_init(r24305); mpfr_init(r24306); mpfr_init(r24307); mpfr_init(r24308); mpfr_init(r24309); mpfr_init(r24310); } double f_im(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r24290, s, MPFR_RNDN); mpfr_neg(r24291, r24290, MPFR_RNDN); mpfr_exp(r24292, r24291, MPFR_RNDN); mpfr_add(r24293, r24289, r24292, MPFR_RNDN); mpfr_div(r24294, r24289, r24293, MPFR_RNDN); mpfr_set_d(r24295, c_p, MPFR_RNDN); mpfr_pow(r24296, r24294, r24295, MPFR_RNDN); mpfr_sub(r24297, r24289, r24294, MPFR_RNDN); mpfr_set_d(r24298, c_n, MPFR_RNDN); mpfr_pow(r24299, r24297, r24298, MPFR_RNDN); mpfr_mul(r24300, r24296, r24299, MPFR_RNDN); mpfr_set_d(r24301, t, MPFR_RNDN); mpfr_neg(r24302, r24301, MPFR_RNDN); mpfr_exp(r24303, r24302, MPFR_RNDN); mpfr_add(r24304, r24289, r24303, MPFR_RNDN); mpfr_div(r24305, r24289, r24304, MPFR_RNDN); mpfr_pow(r24306, r24305, r24295, MPFR_RNDN); mpfr_sub(r24307, r24289, r24305, MPFR_RNDN); mpfr_pow(r24308, r24307, r24298, MPFR_RNDN); mpfr_mul(r24309, r24306, r24308, MPFR_RNDN); mpfr_div(r24310, r24300, r24309, MPFR_RNDN); return mpfr_get_d(r24310, MPFR_RNDN); } static mpfr_t r24311, r24312, r24313, r24314, r24315, r24316, r24317, r24318, r24319, r24320, r24321, r24322, r24323, r24324, r24325, r24326, r24327, r24328, r24329, r24330, r24331, r24332, r24333, r24334, r24335, r24336, r24337; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r24311, "1", 10, MPFR_RNDN); mpfr_init(r24312); mpfr_init(r24313); mpfr_init(r24314); mpfr_init(r24315); mpfr_init(r24316); mpfr_init(r24317); mpfr_init(r24318); mpfr_init(r24319); mpfr_init(r24320); mpfr_init(r24321); mpfr_init(r24322); mpfr_init(r24323); mpfr_init(r24324); mpfr_init(r24325); mpfr_init(r24326); mpfr_init(r24327); mpfr_init(r24328); mpfr_init(r24329); mpfr_init_set_str(r24330, "1/2", 10, MPFR_RNDN); mpfr_init(r24331); mpfr_init(r24332); mpfr_init(r24333); mpfr_init(r24334); mpfr_init(r24335); mpfr_init(r24336); mpfr_init(r24337); } double f_fm(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r24312, s, MPFR_RNDN); mpfr_neg(r24313, r24312, MPFR_RNDN); mpfr_exp(r24314, r24313, MPFR_RNDN); mpfr_add(r24315, r24311, r24314, MPFR_RNDN); mpfr_div(r24316, r24311, r24315, MPFR_RNDN); mpfr_sub(r24317, r24311, r24316, MPFR_RNDN); mpfr_set_d(r24318, c_n, MPFR_RNDN); mpfr_pow(r24319, r24317, r24318, MPFR_RNDN); mpfr_set_d(r24320, t, MPFR_RNDN); mpfr_neg(r24321, r24320, MPFR_RNDN); mpfr_exp(r24322, r24321, MPFR_RNDN); mpfr_add(r24323, r24311, r24322, MPFR_RNDN); mpfr_div(r24324, r24311, r24323, MPFR_RNDN); mpfr_sub(r24325, r24311, r24324, MPFR_RNDN); mpfr_pow(r24326, r24325, r24318, MPFR_RNDN); mpfr_div(r24327, r24319, r24326, MPFR_RNDN); mpfr_set_d(r24328, c_p, MPFR_RNDN); mpfr_pow(r24329, r24316, r24328, MPFR_RNDN); ; mpfr_log(r24331, r24330, MPFR_RNDN); mpfr_mul(r24332, r24320, r24330, MPFR_RNDN); mpfr_add(r24333, r24331, r24332, MPFR_RNDN); mpfr_mul(r24334, r24333, r24328, MPFR_RNDN); mpfr_add(r24335, r24334, r24311, MPFR_RNDN); mpfr_div(r24336, r24329, r24335, MPFR_RNDN); mpfr_mul(r24337, r24327, r24336, MPFR_RNDN); return mpfr_get_d(r24337, MPFR_RNDN); } static mpfr_t r24338, r24339, r24340, r24341, r24342, r24343, r24344, r24345, r24346, r24347, r24348, r24349, r24350, r24351, r24352, r24353, r24354, r24355, r24356, r24357, r24358, r24359, r24360, r24361, r24362, r24363, r24364; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r24338, "1", 10, MPFR_RNDN); mpfr_init(r24339); mpfr_init(r24340); mpfr_init(r24341); mpfr_init(r24342); mpfr_init(r24343); mpfr_init(r24344); mpfr_init(r24345); mpfr_init(r24346); mpfr_init(r24347); mpfr_init(r24348); mpfr_init(r24349); mpfr_init(r24350); mpfr_init(r24351); mpfr_init(r24352); mpfr_init(r24353); mpfr_init(r24354); mpfr_init(r24355); mpfr_init(r24356); mpfr_init_set_str(r24357, "1/2", 10, MPFR_RNDN); mpfr_init(r24358); mpfr_init(r24359); mpfr_init(r24360); mpfr_init(r24361); mpfr_init(r24362); mpfr_init(r24363); mpfr_init(r24364); } double f_dm(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r24339, s, MPFR_RNDN); mpfr_neg(r24340, r24339, MPFR_RNDN); mpfr_exp(r24341, r24340, MPFR_RNDN); mpfr_add(r24342, r24338, r24341, MPFR_RNDN); mpfr_div(r24343, r24338, r24342, MPFR_RNDN); mpfr_sub(r24344, r24338, r24343, MPFR_RNDN); mpfr_set_d(r24345, c_n, MPFR_RNDN); mpfr_pow(r24346, r24344, r24345, MPFR_RNDN); mpfr_set_d(r24347, t, MPFR_RNDN); mpfr_neg(r24348, r24347, MPFR_RNDN); mpfr_exp(r24349, r24348, MPFR_RNDN); mpfr_add(r24350, r24338, r24349, MPFR_RNDN); mpfr_div(r24351, r24338, r24350, MPFR_RNDN); mpfr_sub(r24352, r24338, r24351, MPFR_RNDN); mpfr_pow(r24353, r24352, r24345, MPFR_RNDN); mpfr_div(r24354, r24346, r24353, MPFR_RNDN); mpfr_set_d(r24355, c_p, MPFR_RNDN); mpfr_pow(r24356, r24343, r24355, MPFR_RNDN); ; mpfr_log(r24358, r24357, MPFR_RNDN); mpfr_mul(r24359, r24347, r24357, MPFR_RNDN); mpfr_add(r24360, r24358, r24359, MPFR_RNDN); mpfr_mul(r24361, r24360, r24355, MPFR_RNDN); mpfr_add(r24362, r24361, r24338, MPFR_RNDN); mpfr_div(r24363, r24356, r24362, MPFR_RNDN); mpfr_mul(r24364, r24354, r24363, MPFR_RNDN); return mpfr_get_d(r24364, MPFR_RNDN); }