#include #include #include #include #include char *name = "Jmat.Real.gamma, branch z less than 0.5"; double f_if(float z) { float r27102 = atan2(1.0, 0.0); float r27103 = z; float r27104 = r27102 * r27103; float r27105 = sin(r27104); float r27106 = r27102 / r27105; float r27107 = 2; float r27108 = r27102 * r27107; float r27109 = sqrt(r27108); float r27110 = 1; float r27111 = r27110 - r27103; float r27112 = r27111 - r27110; float r27113 = 7; float r27114 = r27112 + r27113; float r27115 = 0.5; float r27116 = r27114 + r27115; float r27117 = r27112 + r27115; float r27118 = pow(r27116, r27117); float r27119 = r27109 * r27118; float r27120 = -r27116; float r27121 = exp(r27120); float r27122 = r27119 * r27121; float r27123 = 0.9999999999998099; float r27124 = 676.5203681218851; float r27125 = r27112 + r27110; float r27126 = r27124 / r27125; float r27127 = r27123 + r27126; float r27128 = -1259.1392167224028; float r27129 = r27112 + r27107; float r27130 = r27128 / r27129; float r27131 = r27127 + r27130; float r27132 = 771.3234287776531; float r27133 = 3; float r27134 = r27112 + r27133; float r27135 = r27132 / r27134; float r27136 = r27131 + r27135; float r27137 = -176.6150291621406; float r27138 = 4; float r27139 = r27112 + r27138; float r27140 = r27137 / r27139; float r27141 = r27136 + r27140; float r27142 = 12.507343278686905; float r27143 = 5; float r27144 = r27112 + r27143; float r27145 = r27142 / r27144; float r27146 = r27141 + r27145; float r27147 = -0.13857109526572012; float r27148 = 6; float r27149 = r27112 + r27148; float r27150 = r27147 / r27149; float r27151 = r27146 + r27150; float r27152 = 9.984369578019572e-06; float r27153 = r27152 / r27114; float r27154 = r27151 + r27153; float r27155 = 1.5056327351493116e-07; float r27156 = 8; float r27157 = r27112 + r27156; float r27158 = r27155 / r27157; float r27159 = r27154 + r27158; float r27160 = r27122 * r27159; float r27161 = r27106 * r27160; return r27161; } double f_id(double z) { double r27162 = atan2(1.0, 0.0); double r27163 = z; double r27164 = r27162 * r27163; double r27165 = sin(r27164); double r27166 = r27162 / r27165; double r27167 = 2; double r27168 = r27162 * r27167; double r27169 = sqrt(r27168); double r27170 = 1; double r27171 = r27170 - r27163; double r27172 = r27171 - r27170; double r27173 = 7; double r27174 = r27172 + r27173; double r27175 = 0.5; double r27176 = r27174 + r27175; double r27177 = r27172 + r27175; double r27178 = pow(r27176, r27177); double r27179 = r27169 * r27178; double r27180 = -r27176; double r27181 = exp(r27180); double r27182 = r27179 * r27181; double r27183 = 0.9999999999998099; double r27184 = 676.5203681218851; double r27185 = r27172 + r27170; double r27186 = r27184 / r27185; double r27187 = r27183 + r27186; double r27188 = -1259.1392167224028; double r27189 = r27172 + r27167; double r27190 = r27188 / r27189; double r27191 = r27187 + r27190; double r27192 = 771.3234287776531; double r27193 = 3; double r27194 = r27172 + r27193; double r27195 = r27192 / r27194; double r27196 = r27191 + r27195; double r27197 = -176.6150291621406; double r27198 = 4; double r27199 = r27172 + r27198; double r27200 = r27197 / r27199; double r27201 = r27196 + r27200; double r27202 = 12.507343278686905; double r27203 = 5; double r27204 = r27172 + r27203; double r27205 = r27202 / r27204; double r27206 = r27201 + r27205; double r27207 = -0.13857109526572012; double r27208 = 6; double r27209 = r27172 + r27208; double r27210 = r27207 / r27209; double r27211 = r27206 + r27210; double r27212 = 9.984369578019572e-06; double r27213 = r27212 / r27174; double r27214 = r27211 + r27213; double r27215 = 1.5056327351493116e-07; double r27216 = 8; double r27217 = r27172 + r27216; double r27218 = r27215 / r27217; double r27219 = r27214 + r27218; double r27220 = r27182 * r27219; double r27221 = r27166 * r27220; return r27221; } double f_of(float z) { float r27222 = 1.5056327351493116e-07; float r27223 = 8; float r27224 = 1; float r27225 = r27223 + r27224; float r27226 = z; float r27227 = r27224 + r27226; float r27228 = r27225 - r27227; float r27229 = r27222 / r27228; float r27230 = 771.3234287776531; float r27231 = r27224 - r27226; float r27232 = 3; float r27233 = r27224 - r27232; float r27234 = r27231 - r27233; float r27235 = r27230 / r27234; float r27236 = 0.9999999999998099; float r27237 = 676.5203681218851; float r27238 = r27237 / r27231; float r27239 = r27236 + r27238; float r27240 = -1259.1392167224028; float r27241 = 2; float r27242 = r27224 + r27241; float r27243 = r27242 - r27227; float r27244 = r27240 / r27243; float r27245 = r27239 + r27244; float r27246 = r27235 + r27245; float r27247 = 9.984369578019572e-06; float r27248 = 7; float r27249 = r27248 + r27224; float r27250 = r27249 - r27227; float r27251 = r27247 / r27250; float r27252 = -0.13857109526572012; float r27253 = 6; float r27254 = r27253 + r27224; float r27255 = r27254 - r27227; float r27256 = r27252 / r27255; float r27257 = r27251 + r27256; float r27258 = -176.6150291621406; float r27259 = 4; float r27260 = r27224 - r27259; float r27261 = r27231 - r27260; float r27262 = r27258 / r27261; float r27263 = 12.507343278686905; float r27264 = 5; float r27265 = r27224 - r27264; float r27266 = r27231 - r27265; float r27267 = r27263 / r27266; float r27268 = r27262 + r27267; float r27269 = r27257 + r27268; float r27270 = r27246 + r27269; float r27271 = r27229 + r27270; float r27272 = atan2(1.0, 0.0); float r27273 = r27226 * r27272; float r27274 = sin(r27273); float r27275 = r27272 / r27274; float r27276 = 7.5; float r27277 = r27276 - r27226; float r27278 = 0.5; float r27279 = r27278 - r27226; float r27280 = pow(r27277, r27279); float r27281 = sqrt(r27241); float r27282 = r27280 * r27281; float r27283 = r27226 - r27276; float r27284 = exp(r27283); float r27285 = sqrt(r27272); float r27286 = r27284 * r27285; float r27287 = r27282 * r27286; float r27288 = r27275 * r27287; float r27289 = r27271 * r27288; return r27289; } double f_od(double z) { double r27290 = 1.5056327351493116e-07; double r27291 = 8; double r27292 = 1; double r27293 = r27291 + r27292; double r27294 = z; double r27295 = r27292 + r27294; double r27296 = r27293 - r27295; double r27297 = r27290 / r27296; double r27298 = 771.3234287776531; double r27299 = r27292 - r27294; double r27300 = 3; double r27301 = r27292 - r27300; double r27302 = r27299 - r27301; double r27303 = r27298 / r27302; double r27304 = 0.9999999999998099; double r27305 = 676.5203681218851; double r27306 = r27305 / r27299; double r27307 = r27304 + r27306; double r27308 = -1259.1392167224028; double r27309 = 2; double r27310 = r27292 + r27309; double r27311 = r27310 - r27295; double r27312 = r27308 / r27311; double r27313 = r27307 + r27312; double r27314 = r27303 + r27313; double r27315 = 9.984369578019572e-06; double r27316 = 7; double r27317 = r27316 + r27292; double r27318 = r27317 - r27295; double r27319 = r27315 / r27318; double r27320 = -0.13857109526572012; double r27321 = 6; double r27322 = r27321 + r27292; double r27323 = r27322 - r27295; double r27324 = r27320 / r27323; double r27325 = r27319 + r27324; double r27326 = -176.6150291621406; double r27327 = 4; double r27328 = r27292 - r27327; double r27329 = r27299 - r27328; double r27330 = r27326 / r27329; double r27331 = 12.507343278686905; double r27332 = 5; double r27333 = r27292 - r27332; double r27334 = r27299 - r27333; double r27335 = r27331 / r27334; double r27336 = r27330 + r27335; double r27337 = r27325 + r27336; double r27338 = r27314 + r27337; double r27339 = r27297 + r27338; double r27340 = atan2(1.0, 0.0); double r27341 = r27294 * r27340; double r27342 = sin(r27341); double r27343 = r27340 / r27342; double r27344 = 7.5; double r27345 = r27344 - r27294; double r27346 = 0.5; double r27347 = r27346 - r27294; double r27348 = pow(r27345, r27347); double r27349 = sqrt(r27309); double r27350 = r27348 * r27349; double r27351 = r27294 - r27344; double r27352 = exp(r27351); double r27353 = sqrt(r27340); double r27354 = r27352 * r27353; double r27355 = r27350 * r27354; double r27356 = r27343 * r27355; double r27357 = r27339 * r27356; return r27357; } 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 r27358, r27359, r27360, r27361, r27362, r27363, r27364, r27365, r27366, r27367, r27368, r27369, r27370, r27371, r27372, r27373, r27374, r27375, r27376, r27377, r27378, r27379, r27380, r27381, r27382, r27383, r27384, r27385, r27386, r27387, r27388, r27389, r27390, r27391, r27392, r27393, r27394, r27395, r27396, r27397, r27398, r27399, r27400, r27401, r27402, r27403, r27404, r27405, r27406, r27407, r27408, r27409, r27410, r27411, r27412, r27413, r27414, r27415, r27416, r27417; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r27358); mpfr_init(r27359); mpfr_init(r27360); mpfr_init(r27361); mpfr_init(r27362); mpfr_init_set_str(r27363, "2", 10, MPFR_RNDN); mpfr_init(r27364); mpfr_init(r27365); mpfr_init_set_str(r27366, "1", 10, MPFR_RNDN); mpfr_init(r27367); mpfr_init(r27368); mpfr_init_set_str(r27369, "7", 10, MPFR_RNDN); mpfr_init(r27370); mpfr_init_set_str(r27371, "0.5", 10, MPFR_RNDN); mpfr_init(r27372); mpfr_init(r27373); mpfr_init(r27374); mpfr_init(r27375); mpfr_init(r27376); mpfr_init(r27377); mpfr_init(r27378); mpfr_init_set_str(r27379, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r27380, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r27381); mpfr_init(r27382); mpfr_init(r27383); mpfr_init_set_str(r27384, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init(r27385); mpfr_init(r27386); mpfr_init(r27387); mpfr_init_set_str(r27388, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r27389, "3", 10, MPFR_RNDN); mpfr_init(r27390); mpfr_init(r27391); mpfr_init(r27392); mpfr_init_set_str(r27393, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r27394, "4", 10, MPFR_RNDN); mpfr_init(r27395); mpfr_init(r27396); mpfr_init(r27397); mpfr_init_set_str(r27398, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r27399, "5", 10, MPFR_RNDN); mpfr_init(r27400); mpfr_init(r27401); mpfr_init(r27402); mpfr_init_set_str(r27403, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r27404, "6", 10, MPFR_RNDN); mpfr_init(r27405); mpfr_init(r27406); mpfr_init(r27407); mpfr_init_set_str(r27408, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init(r27409); mpfr_init(r27410); mpfr_init_set_str(r27411, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r27412, "8", 10, MPFR_RNDN); mpfr_init(r27413); mpfr_init(r27414); mpfr_init(r27415); mpfr_init(r27416); mpfr_init(r27417); } double f_im(double z) { mpfr_const_pi(r27358, MPFR_RNDN); mpfr_set_d(r27359, z, MPFR_RNDN); mpfr_mul(r27360, r27358, r27359, MPFR_RNDN); mpfr_sin(r27361, r27360, MPFR_RNDN); mpfr_div(r27362, r27358, r27361, MPFR_RNDN); ; mpfr_mul(r27364, r27358, r27363, MPFR_RNDN); mpfr_sqrt(r27365, r27364, MPFR_RNDN); ; mpfr_sub(r27367, r27366, r27359, MPFR_RNDN); mpfr_sub(r27368, r27367, r27366, MPFR_RNDN); ; mpfr_add(r27370, r27368, r27369, MPFR_RNDN); ; mpfr_add(r27372, r27370, r27371, MPFR_RNDN); mpfr_add(r27373, r27368, r27371, MPFR_RNDN); mpfr_pow(r27374, r27372, r27373, MPFR_RNDN); mpfr_mul(r27375, r27365, r27374, MPFR_RNDN); mpfr_neg(r27376, r27372, MPFR_RNDN); mpfr_exp(r27377, r27376, MPFR_RNDN); mpfr_mul(r27378, r27375, r27377, MPFR_RNDN); ; ; mpfr_add(r27381, r27368, r27366, MPFR_RNDN); mpfr_div(r27382, r27380, r27381, MPFR_RNDN); mpfr_add(r27383, r27379, r27382, MPFR_RNDN); ; mpfr_add(r27385, r27368, r27363, MPFR_RNDN); mpfr_div(r27386, r27384, r27385, MPFR_RNDN); mpfr_add(r27387, r27383, r27386, MPFR_RNDN); ; ; mpfr_add(r27390, r27368, r27389, MPFR_RNDN); mpfr_div(r27391, r27388, r27390, MPFR_RNDN); mpfr_add(r27392, r27387, r27391, MPFR_RNDN); ; ; mpfr_add(r27395, r27368, r27394, MPFR_RNDN); mpfr_div(r27396, r27393, r27395, MPFR_RNDN); mpfr_add(r27397, r27392, r27396, MPFR_RNDN); ; ; mpfr_add(r27400, r27368, r27399, MPFR_RNDN); mpfr_div(r27401, r27398, r27400, MPFR_RNDN); mpfr_add(r27402, r27397, r27401, MPFR_RNDN); ; ; mpfr_add(r27405, r27368, r27404, MPFR_RNDN); mpfr_div(r27406, r27403, r27405, MPFR_RNDN); mpfr_add(r27407, r27402, r27406, MPFR_RNDN); ; mpfr_div(r27409, r27408, r27370, MPFR_RNDN); mpfr_add(r27410, r27407, r27409, MPFR_RNDN); ; ; mpfr_add(r27413, r27368, r27412, MPFR_RNDN); mpfr_div(r27414, r27411, r27413, MPFR_RNDN); mpfr_add(r27415, r27410, r27414, MPFR_RNDN); mpfr_mul(r27416, r27378, r27415, MPFR_RNDN); mpfr_mul(r27417, r27362, r27416, MPFR_RNDN); return mpfr_get_d(r27417, MPFR_RNDN); } static mpfr_t r27418, r27419, r27420, r27421, r27422, r27423, r27424, r27425, r27426, r27427, r27428, r27429, r27430, r27431, r27432, r27433, r27434, r27435, r27436, r27437, r27438, r27439, r27440, r27441, r27442, r27443, r27444, r27445, r27446, r27447, r27448, r27449, r27450, r27451, r27452, r27453, r27454, r27455, r27456, r27457, r27458, r27459, r27460, r27461, r27462, r27463, r27464, r27465, r27466, r27467, r27468, r27469, r27470, r27471, r27472, r27473, r27474, r27475, r27476, r27477, r27478, r27479, r27480, r27481, r27482, r27483, r27484, r27485; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27418, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r27419, "8", 10, MPFR_RNDN); mpfr_init_set_str(r27420, "1", 10, MPFR_RNDN); mpfr_init(r27421); mpfr_init(r27422); mpfr_init(r27423); mpfr_init(r27424); mpfr_init(r27425); mpfr_init_set_str(r27426, "771.3234287776531", 10, MPFR_RNDN); mpfr_init(r27427); mpfr_init_set_str(r27428, "3", 10, MPFR_RNDN); mpfr_init(r27429); mpfr_init(r27430); mpfr_init(r27431); mpfr_init_set_str(r27432, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r27433, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r27434); mpfr_init(r27435); mpfr_init_set_str(r27436, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r27437, "2", 10, MPFR_RNDN); mpfr_init(r27438); mpfr_init(r27439); mpfr_init(r27440); mpfr_init(r27441); mpfr_init(r27442); mpfr_init_set_str(r27443, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init_set_str(r27444, "7", 10, MPFR_RNDN); mpfr_init(r27445); mpfr_init(r27446); mpfr_init(r27447); mpfr_init_set_str(r27448, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r27449, "6", 10, MPFR_RNDN); mpfr_init(r27450); mpfr_init(r27451); mpfr_init(r27452); mpfr_init(r27453); mpfr_init_set_str(r27454, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r27455, "4", 10, MPFR_RNDN); mpfr_init(r27456); mpfr_init(r27457); mpfr_init(r27458); mpfr_init_set_str(r27459, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r27460, "5", 10, MPFR_RNDN); mpfr_init(r27461); mpfr_init(r27462); mpfr_init(r27463); mpfr_init(r27464); mpfr_init(r27465); mpfr_init(r27466); mpfr_init(r27467); mpfr_init(r27468); mpfr_init(r27469); mpfr_init(r27470); mpfr_init(r27471); mpfr_init_set_str(r27472, "7.5", 10, MPFR_RNDN); mpfr_init(r27473); mpfr_init_set_str(r27474, "0.5", 10, MPFR_RNDN); mpfr_init(r27475); mpfr_init(r27476); mpfr_init(r27477); mpfr_init(r27478); mpfr_init(r27479); mpfr_init(r27480); mpfr_init(r27481); mpfr_init(r27482); mpfr_init(r27483); mpfr_init(r27484); mpfr_init(r27485); } double f_fm(double z) { ; ; ; mpfr_add(r27421, r27419, r27420, MPFR_RNDN); mpfr_set_d(r27422, z, MPFR_RNDN); mpfr_add(r27423, r27420, r27422, MPFR_RNDN); mpfr_sub(r27424, r27421, r27423, MPFR_RNDN); mpfr_div(r27425, r27418, r27424, MPFR_RNDN); ; mpfr_sub(r27427, r27420, r27422, MPFR_RNDN); ; mpfr_sub(r27429, r27420, r27428, MPFR_RNDN); mpfr_sub(r27430, r27427, r27429, MPFR_RNDN); mpfr_div(r27431, r27426, r27430, MPFR_RNDN); ; ; mpfr_div(r27434, r27433, r27427, MPFR_RNDN); mpfr_add(r27435, r27432, r27434, MPFR_RNDN); ; ; mpfr_add(r27438, r27420, r27437, MPFR_RNDN); mpfr_sub(r27439, r27438, r27423, MPFR_RNDN); mpfr_div(r27440, r27436, r27439, MPFR_RNDN); mpfr_add(r27441, r27435, r27440, MPFR_RNDN); mpfr_add(r27442, r27431, r27441, MPFR_RNDN); ; ; mpfr_add(r27445, r27444, r27420, MPFR_RNDN); mpfr_sub(r27446, r27445, r27423, MPFR_RNDN); mpfr_div(r27447, r27443, r27446, MPFR_RNDN); ; ; mpfr_add(r27450, r27449, r27420, MPFR_RNDN); mpfr_sub(r27451, r27450, r27423, MPFR_RNDN); mpfr_div(r27452, r27448, r27451, MPFR_RNDN); mpfr_add(r27453, r27447, r27452, MPFR_RNDN); ; ; mpfr_sub(r27456, r27420, r27455, MPFR_RNDN); mpfr_sub(r27457, r27427, r27456, MPFR_RNDN); mpfr_div(r27458, r27454, r27457, MPFR_RNDN); ; ; mpfr_sub(r27461, r27420, r27460, MPFR_RNDN); mpfr_sub(r27462, r27427, r27461, MPFR_RNDN); mpfr_div(r27463, r27459, r27462, MPFR_RNDN); mpfr_add(r27464, r27458, r27463, MPFR_RNDN); mpfr_add(r27465, r27453, r27464, MPFR_RNDN); mpfr_add(r27466, r27442, r27465, MPFR_RNDN); mpfr_add(r27467, r27425, r27466, MPFR_RNDN); mpfr_const_pi(r27468, MPFR_RNDN); mpfr_mul(r27469, r27422, r27468, MPFR_RNDN); mpfr_sin(r27470, r27469, MPFR_RNDN); mpfr_div(r27471, r27468, r27470, MPFR_RNDN); ; mpfr_sub(r27473, r27472, r27422, MPFR_RNDN); ; mpfr_sub(r27475, r27474, r27422, MPFR_RNDN); mpfr_pow(r27476, r27473, r27475, MPFR_RNDN); mpfr_sqrt(r27477, r27437, MPFR_RNDN); mpfr_mul(r27478, r27476, r27477, MPFR_RNDN); mpfr_sub(r27479, r27422, r27472, MPFR_RNDN); mpfr_exp(r27480, r27479, MPFR_RNDN); mpfr_sqrt(r27481, r27468, MPFR_RNDN); mpfr_mul(r27482, r27480, r27481, MPFR_RNDN); mpfr_mul(r27483, r27478, r27482, MPFR_RNDN); mpfr_mul(r27484, r27471, r27483, MPFR_RNDN); mpfr_mul(r27485, r27467, r27484, MPFR_RNDN); return mpfr_get_d(r27485, MPFR_RNDN); } static mpfr_t r27486, r27487, r27488, r27489, r27490, r27491, r27492, r27493, r27494, r27495, r27496, r27497, r27498, r27499, r27500, r27501, r27502, r27503, r27504, r27505, r27506, r27507, r27508, r27509, r27510, r27511, r27512, r27513, r27514, r27515, r27516, r27517, r27518, r27519, r27520, r27521, r27522, r27523, r27524, r27525, r27526, r27527, r27528, r27529, r27530, r27531, r27532, r27533, r27534, r27535, r27536, r27537, r27538, r27539, r27540, r27541, r27542, r27543, r27544, r27545, r27546, r27547, r27548, r27549, r27550, r27551, r27552, r27553; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27486, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r27487, "8", 10, MPFR_RNDN); mpfr_init_set_str(r27488, "1", 10, MPFR_RNDN); mpfr_init(r27489); mpfr_init(r27490); mpfr_init(r27491); mpfr_init(r27492); mpfr_init(r27493); mpfr_init_set_str(r27494, "771.3234287776531", 10, MPFR_RNDN); mpfr_init(r27495); mpfr_init_set_str(r27496, "3", 10, MPFR_RNDN); mpfr_init(r27497); mpfr_init(r27498); mpfr_init(r27499); mpfr_init_set_str(r27500, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r27501, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r27502); mpfr_init(r27503); mpfr_init_set_str(r27504, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r27505, "2", 10, MPFR_RNDN); mpfr_init(r27506); mpfr_init(r27507); mpfr_init(r27508); mpfr_init(r27509); mpfr_init(r27510); mpfr_init_set_str(r27511, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init_set_str(r27512, "7", 10, MPFR_RNDN); mpfr_init(r27513); mpfr_init(r27514); mpfr_init(r27515); mpfr_init_set_str(r27516, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r27517, "6", 10, MPFR_RNDN); mpfr_init(r27518); mpfr_init(r27519); mpfr_init(r27520); mpfr_init(r27521); mpfr_init_set_str(r27522, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r27523, "4", 10, MPFR_RNDN); mpfr_init(r27524); mpfr_init(r27525); mpfr_init(r27526); mpfr_init_set_str(r27527, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r27528, "5", 10, MPFR_RNDN); mpfr_init(r27529); mpfr_init(r27530); mpfr_init(r27531); mpfr_init(r27532); mpfr_init(r27533); mpfr_init(r27534); mpfr_init(r27535); mpfr_init(r27536); mpfr_init(r27537); mpfr_init(r27538); mpfr_init(r27539); mpfr_init_set_str(r27540, "7.5", 10, MPFR_RNDN); mpfr_init(r27541); mpfr_init_set_str(r27542, "0.5", 10, MPFR_RNDN); mpfr_init(r27543); mpfr_init(r27544); mpfr_init(r27545); mpfr_init(r27546); mpfr_init(r27547); mpfr_init(r27548); mpfr_init(r27549); mpfr_init(r27550); mpfr_init(r27551); mpfr_init(r27552); mpfr_init(r27553); } double f_dm(double z) { ; ; ; mpfr_add(r27489, r27487, r27488, MPFR_RNDN); mpfr_set_d(r27490, z, MPFR_RNDN); mpfr_add(r27491, r27488, r27490, MPFR_RNDN); mpfr_sub(r27492, r27489, r27491, MPFR_RNDN); mpfr_div(r27493, r27486, r27492, MPFR_RNDN); ; mpfr_sub(r27495, r27488, r27490, MPFR_RNDN); ; mpfr_sub(r27497, r27488, r27496, MPFR_RNDN); mpfr_sub(r27498, r27495, r27497, MPFR_RNDN); mpfr_div(r27499, r27494, r27498, MPFR_RNDN); ; ; mpfr_div(r27502, r27501, r27495, MPFR_RNDN); mpfr_add(r27503, r27500, r27502, MPFR_RNDN); ; ; mpfr_add(r27506, r27488, r27505, MPFR_RNDN); mpfr_sub(r27507, r27506, r27491, MPFR_RNDN); mpfr_div(r27508, r27504, r27507, MPFR_RNDN); mpfr_add(r27509, r27503, r27508, MPFR_RNDN); mpfr_add(r27510, r27499, r27509, MPFR_RNDN); ; ; mpfr_add(r27513, r27512, r27488, MPFR_RNDN); mpfr_sub(r27514, r27513, r27491, MPFR_RNDN); mpfr_div(r27515, r27511, r27514, MPFR_RNDN); ; ; mpfr_add(r27518, r27517, r27488, MPFR_RNDN); mpfr_sub(r27519, r27518, r27491, MPFR_RNDN); mpfr_div(r27520, r27516, r27519, MPFR_RNDN); mpfr_add(r27521, r27515, r27520, MPFR_RNDN); ; ; mpfr_sub(r27524, r27488, r27523, MPFR_RNDN); mpfr_sub(r27525, r27495, r27524, MPFR_RNDN); mpfr_div(r27526, r27522, r27525, MPFR_RNDN); ; ; mpfr_sub(r27529, r27488, r27528, MPFR_RNDN); mpfr_sub(r27530, r27495, r27529, MPFR_RNDN); mpfr_div(r27531, r27527, r27530, MPFR_RNDN); mpfr_add(r27532, r27526, r27531, MPFR_RNDN); mpfr_add(r27533, r27521, r27532, MPFR_RNDN); mpfr_add(r27534, r27510, r27533, MPFR_RNDN); mpfr_add(r27535, r27493, r27534, MPFR_RNDN); mpfr_const_pi(r27536, MPFR_RNDN); mpfr_mul(r27537, r27490, r27536, MPFR_RNDN); mpfr_sin(r27538, r27537, MPFR_RNDN); mpfr_div(r27539, r27536, r27538, MPFR_RNDN); ; mpfr_sub(r27541, r27540, r27490, MPFR_RNDN); ; mpfr_sub(r27543, r27542, r27490, MPFR_RNDN); mpfr_pow(r27544, r27541, r27543, MPFR_RNDN); mpfr_sqrt(r27545, r27505, MPFR_RNDN); mpfr_mul(r27546, r27544, r27545, MPFR_RNDN); mpfr_sub(r27547, r27490, r27540, MPFR_RNDN); mpfr_exp(r27548, r27547, MPFR_RNDN); mpfr_sqrt(r27549, r27536, MPFR_RNDN); mpfr_mul(r27550, r27548, r27549, MPFR_RNDN); mpfr_mul(r27551, r27546, r27550, MPFR_RNDN); mpfr_mul(r27552, r27539, r27551, MPFR_RNDN); mpfr_mul(r27553, r27535, r27552, MPFR_RNDN); return mpfr_get_d(r27553, MPFR_RNDN); }