Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2, C

Average Error: 26.1 → 0.7

Time: 20.9s

Precision: 64

Math

\[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922227999963610045597306452691555 + 78.69949241540000173245061887428164482117\right) \cdot x + 137.5194164160000127594685181975364685059\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514000013984514225739985704422\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x + 313.3992158940000081202015280723571777344\right) \cdot x + 47.06687660600000100430406746454536914825}\]

↓

\[\begin{array}{l} \mathbf{if}\;x \le -1.486655794493063927988268380417398027792 \cdot 10^{52}:\\ \;\;\;\;\left(4.16438922227999963610045597306452691555 \cdot x - 110.1139242984810948655649553984403610229\right) + \frac{\frac{y}{x}}{x}\\ \mathbf{elif}\;x \le 2.375361050611561796555312227670044766263 \cdot 10^{59}:\\ \;\;\;\;\left(\left(z + x \cdot \left(y + x \cdot \left(x \cdot \left(78.69949241540000173245061887428164482117 + 4.16438922227999963610045597306452691555 \cdot x\right) + 137.5194164160000127594685181975364685059\right)\right)\right) \cdot \frac{1}{\left(\left(\left(x + 43.3400022514000013984514225739985704422\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x + 313.3992158940000081202015280723571777344\right) \cdot x + 47.06687660600000100430406746454536914825}\right) \cdot \left(x - 2\right)\\ \mathbf{else}:\\ \;\;\;\;\left(4.16438922227999963610045597306452691555 + \left(\frac{y}{{x}^{3}} - \frac{101.785145853921093817007204052060842514}{x}\right)\right) \cdot \left(x - 2\right)\\ \end{array}\]

C

double f(double x, double y, double z) {
        double r871234 = x;
        double r871235 = 2.0;
        double r871236 = r871234 - r871235;
        double r871237 = 4.16438922228;
        double r871238 = r871234 * r871237;
        double r871239 = 78.6994924154;
        double r871240 = r871238 + r871239;
        double r871241 = r871240 * r871234;
        double r871242 = 137.519416416;
        double r871243 = r871241 + r871242;
        double r871244 = r871243 * r871234;
        double r871245 = y;
        double r871246 = r871244 + r871245;
        double r871247 = r871246 * r871234;
        double r871248 = z;
        double r871249 = r871247 + r871248;
        double r871250 = r871236 * r871249;
        double r871251 = 43.3400022514;
        double r871252 = r871234 + r871251;
        double r871253 = r871252 * r871234;
        double r871254 = 263.505074721;
        double r871255 = r871253 + r871254;
        double r871256 = r871255 * r871234;
        double r871257 = 313.399215894;
        double r871258 = r871256 + r871257;
        double r871259 = r871258 * r871234;
        double r871260 = 47.066876606;
        double r871261 = r871259 + r871260;
        double r871262 = r871250 / r871261;
        return r871262;
}

↓

double f(double x, double y, double z) {
        double r871263 = x;
        double r871264 = -1.486655794493064e+52;
        bool r871265 = r871263 <= r871264;
        double r871266 = 4.16438922228;
        double r871267 = r871266 * r871263;
        double r871268 = 110.1139242984811;
        double r871269 = r871267 - r871268;
        double r871270 = y;
        double r871271 = r871270 / r871263;
        double r871272 = r871271 / r871263;
        double r871273 = r871269 + r871272;
        double r871274 = 2.375361050611562e+59;
        bool r871275 = r871263 <= r871274;
        double r871276 = z;
        double r871277 = 78.6994924154;
        double r871278 = r871277 + r871267;
        double r871279 = r871263 * r871278;
        double r871280 = 137.519416416;
        double r871281 = r871279 + r871280;
        double r871282 = r871263 * r871281;
        double r871283 = r871270 + r871282;
        double r871284 = r871263 * r871283;
        double r871285 = r871276 + r871284;
        double r871286 = 1.0;
        double r871287 = 43.3400022514;
        double r871288 = r871263 + r871287;
        double r871289 = r871288 * r871263;
        double r871290 = 263.505074721;
        double r871291 = r871289 + r871290;
        double r871292 = r871291 * r871263;
        double r871293 = 313.399215894;
        double r871294 = r871292 + r871293;
        double r871295 = r871294 * r871263;
        double r871296 = 47.066876606;
        double r871297 = r871295 + r871296;
        double r871298 = r871286 / r871297;
        double r871299 = r871285 * r871298;
        double r871300 = 2.0;
        double r871301 = r871263 - r871300;
        double r871302 = r871299 * r871301;
        double r871303 = 3.0;
        double r871304 = pow(r871263, r871303);
        double r871305 = r871270 / r871304;
        double r871306 = 101.7851458539211;
        double r871307 = r871306 / r871263;
        double r871308 = r871305 - r871307;
        double r871309 = r871266 + r871308;
        double r871310 = r871309 * r871301;
        double r871311 = r871275 ? r871302 : r871310;
        double r871312 = r871265 ? r871273 : r871311;
        return r871312;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Target

Original	26.1
Target	0.5
Herbie	0.7

\[\begin{array}{l} \mathbf{if}\;x \lt -3.326128725870004842699683658678411714981 \cdot 10^{62}:\\ \;\;\;\;\left(\frac{y}{x \cdot x} + 4.16438922227999963610045597306452691555 \cdot x\right) - 110.1139242984810948655649553984403610229\\ \mathbf{elif}\;x \lt 9.429991714554672672712552870340896976735 \cdot 10^{55}:\\ \;\;\;\;\frac{x - 2}{1} \cdot \frac{\left(\left(\left(x \cdot 4.16438922227999963610045597306452691555 + 78.69949241540000173245061887428164482117\right) \cdot x + 137.5194164160000127594685181975364685059\right) \cdot x + y\right) \cdot x + z}{\left(\left(263.5050747210000281484099105000495910645 \cdot x + \left(43.3400022514000013984514225739985704422 \cdot \left(x \cdot x\right) + x \cdot \left(x \cdot x\right)\right)\right) + 313.3992158940000081202015280723571777344\right) \cdot x + 47.06687660600000100430406746454536914825}\\ \mathbf{else}:\\ \;\;\;\;\left(\frac{y}{x \cdot x} + 4.16438922227999963610045597306452691555 \cdot x\right) - 110.1139242984810948655649553984403610229\\ \end{array}\]

Derivation

1. Split input into 3 regimes

2. if x < -1.486655794493064e+52

   1. Initial program 62.1

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922227999963610045597306452691555 + 78.69949241540000173245061887428164482117\right) \cdot x + 137.5194164160000127594685181975364685059\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514000013984514225739985704422\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x + 313.3992158940000081202015280723571777344\right) \cdot x + 47.06687660600000100430406746454536914825}\]

   2. Simplified 58.6

      \[\leadsto \color{blue}{\frac{x - 2}{x \cdot \left(313.3992158940000081202015280723571777344 + \left(263.5050747210000281484099105000495910645 + x \cdot \left(43.3400022514000013984514225739985704422 + x\right)\right) \cdot x\right) + 47.06687660600000100430406746454536914825} \cdot \left(\left(y + \left(x \cdot \left(78.69949241540000173245061887428164482117 + 4.16438922227999963610045597306452691555 \cdot x\right) + 137.5194164160000127594685181975364685059\right) \cdot x\right) \cdot x + z\right)}\]

   3. Taylor expanded around inf 0.4

      \[\leadsto \color{blue}{\left(\frac{y}{{x}^{2}} + 4.16438922227999963610045597306452691555 \cdot x\right) - 110.1139242984810948655649553984403610229}\]

   4. Simplified 0.4

      \[\leadsto \color{blue}{\left(4.16438922227999963610045597306452691555 \cdot x - 110.1139242984810948655649553984403610229\right) + \frac{\frac{y}{x}}{x}}\]

   if -1.486655794493064e+52 < x < 2.375361050611562e+59

   1. Initial program 1.4

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922227999963610045597306452691555 + 78.69949241540000173245061887428164482117\right) \cdot x + 137.5194164160000127594685181975364685059\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514000013984514225739985704422\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x + 313.3992158940000081202015280723571777344\right) \cdot x + 47.06687660600000100430406746454536914825}\]

   2. Simplified 0.9

      \[\leadsto \color{blue}{\frac{x - 2}{x \cdot \left(313.3992158940000081202015280723571777344 + \left(263.5050747210000281484099105000495910645 + x \cdot \left(43.3400022514000013984514225739985704422 + x\right)\right) \cdot x\right) + 47.06687660600000100430406746454536914825} \cdot \left(\left(y + \left(x \cdot \left(78.69949241540000173245061887428164482117 + 4.16438922227999963610045597306452691555 \cdot x\right) + 137.5194164160000127594685181975364685059\right) \cdot x\right) \cdot x + z\right)}\]
   3. Using strategy rm

   4. Applied div-inv 0.9

      \[\leadsto \color{blue}{\left(\left(x - 2\right) \cdot \frac{1}{x \cdot \left(313.3992158940000081202015280723571777344 + \left(263.5050747210000281484099105000495910645 + x \cdot \left(43.3400022514000013984514225739985704422 + x\right)\right) \cdot x\right) + 47.06687660600000100430406746454536914825}\right)} \cdot \left(\left(y + \left(x \cdot \left(78.69949241540000173245061887428164482117 + 4.16438922227999963610045597306452691555 \cdot x\right) + 137.5194164160000127594685181975364685059\right) \cdot x\right) \cdot x + z\right)\]

   5. Applied associate-*l* 0.9

      \[\leadsto \color{blue}{\left(x - 2\right) \cdot \left(\frac{1}{x \cdot \left(313.3992158940000081202015280723571777344 + \left(263.5050747210000281484099105000495910645 + x \cdot \left(43.3400022514000013984514225739985704422 + x\right)\right) \cdot x\right) + 47.06687660600000100430406746454536914825} \cdot \left(\left(y + \left(x \cdot \left(78.69949241540000173245061887428164482117 + 4.16438922227999963610045597306452691555 \cdot x\right) + 137.5194164160000127594685181975364685059\right) \cdot x\right) \cdot x + z\right)\right)}\]

   6. Simplified 0.5

      \[\leadsto \left(x - 2\right) \cdot \color{blue}{\frac{\left(\left(\left(4.16438922227999963610045597306452691555 \cdot x + 78.69949241540000173245061887428164482117\right) \cdot x + 137.5194164160000127594685181975364685059\right) \cdot x + y\right) \cdot x + z}{47.06687660600000100430406746454536914825 + \left(313.3992158940000081202015280723571777344 + \left(\left(43.3400022514000013984514225739985704422 + x\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x\right) \cdot x}}\]
   7. Using strategy rm

   8. Applied add-sqr-sqrt 0.7

      \[\leadsto \left(x - 2\right) \cdot \frac{\left(\left(\left(4.16438922227999963610045597306452691555 \cdot x + 78.69949241540000173245061887428164482117\right) \cdot x + 137.5194164160000127594685181975364685059\right) \cdot x + y\right) \cdot x + z}{\color{blue}{\sqrt{47.06687660600000100430406746454536914825 + \left(313.3992158940000081202015280723571777344 + \left(\left(43.3400022514000013984514225739985704422 + x\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x\right) \cdot x} \cdot \sqrt{47.06687660600000100430406746454536914825 + \left(313.3992158940000081202015280723571777344 + \left(\left(43.3400022514000013984514225739985704422 + x\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x\right) \cdot x}}}\]

   9. Applied associate-/r* 0.9

      \[\leadsto \left(x - 2\right) \cdot \color{blue}{\frac{\frac{\left(\left(\left(4.16438922227999963610045597306452691555 \cdot x + 78.69949241540000173245061887428164482117\right) \cdot x + 137.5194164160000127594685181975364685059\right) \cdot x + y\right) \cdot x + z}{\sqrt{47.06687660600000100430406746454536914825 + \left(313.3992158940000081202015280723571777344 + \left(\left(43.3400022514000013984514225739985704422 + x\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x\right) \cdot x}}}{\sqrt{47.06687660600000100430406746454536914825 + \left(313.3992158940000081202015280723571777344 + \left(\left(43.3400022514000013984514225739985704422 + x\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x\right) \cdot x}}}\]

   10. Simplified 0.9

       \[\leadsto \left(x - 2\right) \cdot \frac{\color{blue}{\frac{z + \left(x \cdot \left(x \cdot \left(4.16438922227999963610045597306452691555 \cdot x + 78.69949241540000173245061887428164482117\right) + 137.5194164160000127594685181975364685059\right) + y\right) \cdot x}{\sqrt{x \cdot \left(\left(263.5050747210000281484099105000495910645 + x \cdot \left(x + 43.3400022514000013984514225739985704422\right)\right) \cdot x + 313.3992158940000081202015280723571777344\right) + 47.06687660600000100430406746454536914825}}}}{\sqrt{47.06687660600000100430406746454536914825 + \left(313.3992158940000081202015280723571777344 + \left(\left(43.3400022514000013984514225739985704422 + x\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x\right) \cdot x}}\]
   11. Using strategy rm

   12. Applied *-un-lft-identity 0.9

       \[\leadsto \left(x - 2\right) \cdot \frac{\frac{z + \left(x \cdot \left(x \cdot \left(4.16438922227999963610045597306452691555 \cdot x + 78.69949241540000173245061887428164482117\right) + 137.5194164160000127594685181975364685059\right) + y\right) \cdot x}{\sqrt{x \cdot \left(\left(263.5050747210000281484099105000495910645 + x \cdot \left(x + 43.3400022514000013984514225739985704422\right)\right) \cdot x + 313.3992158940000081202015280723571777344\right) + 47.06687660600000100430406746454536914825}}}{\sqrt{\color{blue}{1 \cdot \left(47.06687660600000100430406746454536914825 + \left(313.3992158940000081202015280723571777344 + \left(\left(43.3400022514000013984514225739985704422 + x\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x\right) \cdot x\right)}}}\]

   13. Applied sqrt-prod 0.9

       \[\leadsto \left(x - 2\right) \cdot \frac{\frac{z + \left(x \cdot \left(x \cdot \left(4.16438922227999963610045597306452691555 \cdot x + 78.69949241540000173245061887428164482117\right) + 137.5194164160000127594685181975364685059\right) + y\right) \cdot x}{\sqrt{x \cdot \left(\left(263.5050747210000281484099105000495910645 + x \cdot \left(x + 43.3400022514000013984514225739985704422\right)\right) \cdot x + 313.3992158940000081202015280723571777344\right) + 47.06687660600000100430406746454536914825}}}{\color{blue}{\sqrt{1} \cdot \sqrt{47.06687660600000100430406746454536914825 + \left(313.3992158940000081202015280723571777344 + \left(\left(43.3400022514000013984514225739985704422 + x\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x\right) \cdot x}}}\]

   14. Applied div-inv 0.9

       \[\leadsto \left(x - 2\right) \cdot \frac{\color{blue}{\left(z + \left(x \cdot \left(x \cdot \left(4.16438922227999963610045597306452691555 \cdot x + 78.69949241540000173245061887428164482117\right) + 137.5194164160000127594685181975364685059\right) + y\right) \cdot x\right) \cdot \frac{1}{\sqrt{x \cdot \left(\left(263.5050747210000281484099105000495910645 + x \cdot \left(x + 43.3400022514000013984514225739985704422\right)\right) \cdot x + 313.3992158940000081202015280723571777344\right) + 47.06687660600000100430406746454536914825}}}}{\sqrt{1} \cdot \sqrt{47.06687660600000100430406746454536914825 + \left(313.3992158940000081202015280723571777344 + \left(\left(43.3400022514000013984514225739985704422 + x\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x\right) \cdot x}}\]

   15. Applied times-frac 1.1

       \[\leadsto \left(x - 2\right) \cdot \color{blue}{\left(\frac{z + \left(x \cdot \left(x \cdot \left(4.16438922227999963610045597306452691555 \cdot x + 78.69949241540000173245061887428164482117\right) + 137.5194164160000127594685181975364685059\right) + y\right) \cdot x}{\sqrt{1}} \cdot \frac{\frac{1}{\sqrt{x \cdot \left(\left(263.5050747210000281484099105000495910645 + x \cdot \left(x + 43.3400022514000013984514225739985704422\right)\right) \cdot x + 313.3992158940000081202015280723571777344\right) + 47.06687660600000100430406746454536914825}}}{\sqrt{47.06687660600000100430406746454536914825 + \left(313.3992158940000081202015280723571777344 + \left(\left(43.3400022514000013984514225739985704422 + x\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x\right) \cdot x}}\right)}\]

   16. Simplified 1.1

       \[\leadsto \left(x - 2\right) \cdot \left(\color{blue}{\left(x \cdot \left(y + \left(137.5194164160000127594685181975364685059 + \left(x \cdot 4.16438922227999963610045597306452691555 + 78.69949241540000173245061887428164482117\right) \cdot x\right) \cdot x\right) + z\right)} \cdot \frac{\frac{1}{\sqrt{x \cdot \left(\left(263.5050747210000281484099105000495910645 + x \cdot \left(x + 43.3400022514000013984514225739985704422\right)\right) \cdot x + 313.3992158940000081202015280723571777344\right) + 47.06687660600000100430406746454536914825}}}{\sqrt{47.06687660600000100430406746454536914825 + \left(313.3992158940000081202015280723571777344 + \left(\left(43.3400022514000013984514225739985704422 + x\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x\right) \cdot x}}\right)\]

   17. Simplified 0.9

       \[\leadsto \left(x - 2\right) \cdot \left(\left(x \cdot \left(y + \left(137.5194164160000127594685181975364685059 + \left(x \cdot 4.16438922227999963610045597306452691555 + 78.69949241540000173245061887428164482117\right) \cdot x\right) \cdot x\right) + z\right) \cdot \color{blue}{\frac{1}{47.06687660600000100430406746454536914825 + x \cdot \left(\left(x \cdot \left(x + 43.3400022514000013984514225739985704422\right) + 263.5050747210000281484099105000495910645\right) \cdot x + 313.3992158940000081202015280723571777344\right)}}\right)\]

   if 2.375361050611562e+59 < x

   1. Initial program 63.9

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922227999963610045597306452691555 + 78.69949241540000173245061887428164482117\right) \cdot x + 137.5194164160000127594685181975364685059\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514000013984514225739985704422\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x + 313.3992158940000081202015280723571777344\right) \cdot x + 47.06687660600000100430406746454536914825}\]

   2. Simplified 60.9

      \[\leadsto \color{blue}{\frac{x - 2}{x \cdot \left(313.3992158940000081202015280723571777344 + \left(263.5050747210000281484099105000495910645 + x \cdot \left(43.3400022514000013984514225739985704422 + x\right)\right) \cdot x\right) + 47.06687660600000100430406746454536914825} \cdot \left(\left(y + \left(x \cdot \left(78.69949241540000173245061887428164482117 + 4.16438922227999963610045597306452691555 \cdot x\right) + 137.5194164160000127594685181975364685059\right) \cdot x\right) \cdot x + z\right)}\]
   3. Using strategy rm

   4. Applied div-inv 60.9

      \[\leadsto \color{blue}{\left(\left(x - 2\right) \cdot \frac{1}{x \cdot \left(313.3992158940000081202015280723571777344 + \left(263.5050747210000281484099105000495910645 + x \cdot \left(43.3400022514000013984514225739985704422 + x\right)\right) \cdot x\right) + 47.06687660600000100430406746454536914825}\right)} \cdot \left(\left(y + \left(x \cdot \left(78.69949241540000173245061887428164482117 + 4.16438922227999963610045597306452691555 \cdot x\right) + 137.5194164160000127594685181975364685059\right) \cdot x\right) \cdot x + z\right)\]

   5. Applied associate-*l* 60.9

      \[\leadsto \color{blue}{\left(x - 2\right) \cdot \left(\frac{1}{x \cdot \left(313.3992158940000081202015280723571777344 + \left(263.5050747210000281484099105000495910645 + x \cdot \left(43.3400022514000013984514225739985704422 + x\right)\right) \cdot x\right) + 47.06687660600000100430406746454536914825} \cdot \left(\left(y + \left(x \cdot \left(78.69949241540000173245061887428164482117 + 4.16438922227999963610045597306452691555 \cdot x\right) + 137.5194164160000127594685181975364685059\right) \cdot x\right) \cdot x + z\right)\right)}\]

   6. Simplified 60.9

      \[\leadsto \left(x - 2\right) \cdot \color{blue}{\frac{\left(\left(\left(4.16438922227999963610045597306452691555 \cdot x + 78.69949241540000173245061887428164482117\right) \cdot x + 137.5194164160000127594685181975364685059\right) \cdot x + y\right) \cdot x + z}{47.06687660600000100430406746454536914825 + \left(313.3992158940000081202015280723571777344 + \left(\left(43.3400022514000013984514225739985704422 + x\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x\right) \cdot x}}\]

   7. Taylor expanded around inf 0.4

      \[\leadsto \left(x - 2\right) \cdot \color{blue}{\left(\left(\frac{y}{{x}^{3}} + 4.16438922227999963610045597306452691555\right) - 101.785145853921093817007204052060842514 \cdot \frac{1}{x}\right)}\]

   8. Simplified 0.4

      \[\leadsto \left(x - 2\right) \cdot \color{blue}{\left(4.16438922227999963610045597306452691555 + \left(\frac{y}{{x}^{3}} - \frac{101.785145853921093817007204052060842514}{x}\right)\right)}\]
3. Recombined 3 regimes into one program.

4. Final simplification 0.7

   \[\leadsto \begin{array}{l} \mathbf{if}\;x \le -1.486655794493063927988268380417398027792 \cdot 10^{52}:\\ \;\;\;\;\left(4.16438922227999963610045597306452691555 \cdot x - 110.1139242984810948655649553984403610229\right) + \frac{\frac{y}{x}}{x}\\ \mathbf{elif}\;x \le 2.375361050611561796555312227670044766263 \cdot 10^{59}:\\ \;\;\;\;\left(\left(z + x \cdot \left(y + x \cdot \left(x \cdot \left(78.69949241540000173245061887428164482117 + 4.16438922227999963610045597306452691555 \cdot x\right) + 137.5194164160000127594685181975364685059\right)\right)\right) \cdot \frac{1}{\left(\left(\left(x + 43.3400022514000013984514225739985704422\right) \cdot x + 263.5050747210000281484099105000495910645\right) \cdot x + 313.3992158940000081202015280723571777344\right) \cdot x + 47.06687660600000100430406746454536914825}\right) \cdot \left(x - 2\right)\\ \mathbf{else}:\\ \;\;\;\;\left(4.16438922227999963610045597306452691555 + \left(\frac{y}{{x}^{3}} - \frac{101.785145853921093817007204052060842514}{x}\right)\right) \cdot \left(x - 2\right)\\ \end{array}\]

Reproduce

herbie shell --seed 2019179 
(FPCore (x y z)
  :name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2, C"

  :herbie-target
  (if (< x -3.326128725870005e+62) (- (+ (/ y (* x x)) (* 4.16438922228 x)) 110.1139242984811) (if (< x 9.429991714554673e+55) (* (/ (- x 2.0) 1.0) (/ (+ (* (+ (* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x) y) x) z) (+ (* (+ (+ (* 263.505074721 x) (+ (* 43.3400022514 (* x x)) (* x (* x x)))) 313.399215894) x) 47.066876606))) (- (+ (/ y (* x x)) (* 4.16438922228 x)) 110.1139242984811)))

  (/ (* (- x 2.0) (+ (* (+ (* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x) y) x) z)) (+ (* (+ (* (+ (* (+ x 43.3400022514) x) 263.505074721) x) 313.399215894) x) 47.066876606)))