\left(\left(\sqrt{\pi \cdot 2} \cdot {\left(\left(\left(z - 1\right) + 7\right) + 0.5\right)}^{\left(\left(z - 1\right) + 0.5\right)}\right) \cdot e^{-\left(\left(\left(z - 1\right) + 7\right) + 0.5\right)}\right) \cdot \left(\left(\left(\left(\left(\left(\left(\left(0.9999999999998099298181841732002794742584 + \frac{676.5203681218850988443591631948947906494}{\left(z - 1\right) + 1}\right) + \frac{-1259.139216722402807135949842631816864014}{\left(z - 1\right) + 2}\right) + \frac{771.3234287776531346025876700878143310547}{\left(z - 1\right) + 3}\right) + \frac{-176.6150291621405870046146446838974952698}{\left(z - 1\right) + 4}\right) + \frac{12.50734327868690520801919774385169148445}{\left(z - 1\right) + 5}\right) + \frac{-0.1385710952657201178173096423051902092993}{\left(z - 1\right) + 6}\right) + \frac{9.984369578019571583242346146658263705831 \cdot 10^{-6}}{\left(z - 1\right) + 7}\right) + \frac{1.505632735149311617592788074479481785772 \cdot 10^{-7}}{\left(z - 1\right) + 8}\right)double f(double z) {
double r11146 = atan2(1.0, 0.0);
double r11147 = 2.0;
double r11148 = r11146 * r11147;
double r11149 = sqrt(r11148);
double r11150 = z;
double r11151 = 1.0;
double r11152 = r11150 - r11151;
double r11153 = 7.0;
double r11154 = r11152 + r11153;
double r11155 = 0.5;
double r11156 = r11154 + r11155;
double r11157 = r11152 + r11155;
double r11158 = pow(r11156, r11157);
double r11159 = r11149 * r11158;
double r11160 = -r11156;
double r11161 = exp(r11160);
double r11162 = r11159 * r11161;
double r11163 = 0.9999999999998099;
double r11164 = 676.5203681218851;
double r11165 = r11152 + r11151;
double r11166 = r11164 / r11165;
double r11167 = r11163 + r11166;
double r11168 = -1259.1392167224028;
double r11169 = r11152 + r11147;
double r11170 = r11168 / r11169;
double r11171 = r11167 + r11170;
double r11172 = 771.3234287776531;
double r11173 = 3.0;
double r11174 = r11152 + r11173;
double r11175 = r11172 / r11174;
double r11176 = r11171 + r11175;
double r11177 = -176.6150291621406;
double r11178 = 4.0;
double r11179 = r11152 + r11178;
double r11180 = r11177 / r11179;
double r11181 = r11176 + r11180;
double r11182 = 12.507343278686905;
double r11183 = 5.0;
double r11184 = r11152 + r11183;
double r11185 = r11182 / r11184;
double r11186 = r11181 + r11185;
double r11187 = -0.13857109526572012;
double r11188 = 6.0;
double r11189 = r11152 + r11188;
double r11190 = r11187 / r11189;
double r11191 = r11186 + r11190;
double r11192 = 9.984369578019572e-06;
double r11193 = r11192 / r11154;
double r11194 = r11191 + r11193;
double r11195 = 1.5056327351493116e-07;
double r11196 = 8.0;
double r11197 = r11152 + r11196;
double r11198 = r11195 / r11197;
double r11199 = r11194 + r11198;
double r11200 = r11162 * r11199;
return r11200;
}
Please include this information when filing a bug report:
herbie shell --seed 2019315 +o rules:numerics
(FPCore (z)
:name "Jmat.Real.gamma, branch z greater than 0.5"
:precision binary64
(* (* (* (sqrt (* PI 2)) (pow (+ (+ (- z 1) 7) 0.5) (+ (- z 1) 0.5))) (exp (- (+ (+ (- z 1) 7) 0.5)))) (+ (+ (+ (+ (+ (+ (+ (+ 0.99999999999980993 (/ 676.520368121885099 (+ (- z 1) 1))) (/ -1259.13921672240281 (+ (- z 1) 2))) (/ 771.32342877765313 (+ (- z 1) 3))) (/ -176.615029162140587 (+ (- z 1) 4))) (/ 12.5073432786869052 (+ (- z 1) 5))) (/ -0.138571095265720118 (+ (- z 1) 6))) (/ 9.98436957801957158e-6 (+ (- z 1) 7))) (/ 1.50563273514931162e-7 (+ (- z 1) 8)))))
| match: no matching clause for 'binary64 | L | C | |
|---|---|---|---|
| match:error | /opt/racket-7.2/collects/racket/match/runtime.rkt | 24 | 0 |
| make-multiplication-subsubsubnode | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/core/reduce.rkt | 266 | 0 |
| for-loop | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/core/reduce.rkt | 246 | 3 |
| make-multiplication-subnode | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/core/reduce.rkt | 244 | 0 |
| make-multiplication-node | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/core/reduce.rkt | 234 | 0 |
| map | /opt/racket-7.2/collects/racket/private/map.rkt | 35 | 13 |
| f45 | /opt/racket-7.2/collects/racket/match/compiler.rkt | 507 | 40 |
| simplify | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/core/reduce.rkt | 16 | 0 |
| (unnamed) | /opt/racket-7.2/collects/racket/contract/private/arrow-higher-order.rkt | 360 | 33 |
| location-do | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/programs.rkt | 64 | 0 |
| location-do | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/programs.rkt | 64 | 0 |
| location-do | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/programs.rkt | 64 | 0 |
| (unnamed) | /opt/racket-7.2/collects/racket/contract/private/arrow-higher-order.rkt | 360 | 33 |
| for-loop | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt | 151 | 6 |
| taylor-alt | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt | 145 | 0 |
| for-loop | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt | 166 | 7 |
| gen-series! | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt | 159 | 0 |
| run-iter! | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt | 326 | 0 |
| for-loop | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt | 364 | 4 |
| run-improve47 | /data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt | 345 | 0 |
| (unnamed) | /opt/racket-7.2/collects/racket/private/more-scheme.rkt | 261 | 28 |
| run | /opt/racket-7.2/share/pkgs/profile-lib/main.rkt | 39 | 2 |
| profile-thunk16 | /opt/racket-7.2/share/pkgs/profile-lib/main.rkt | 9 | 0 |
| (unnamed) | /opt/racket-7.2/collects/racket/private/more-scheme.rkt | 261 | 28 |