?

Average Error: 0.3 → 0.3
Time: 23.1s
Precision: binary64
Cost: 26304

?

\[\left(\left(\log \left(x + y\right) + \log z\right) - t\right) + \left(a - 0.5\right) \cdot \log t \]
↓
\[\log \left(x + y\right) + \left(\mathsf{fma}\left(a + -0.5, \log t, \log z\right) - t\right) \]
(FPCore (x y z t a)
 :precision binary64
 (+ (- (+ (log (+ x y)) (log z)) t) (* (- a 0.5) (log t))))
↓
(FPCore (x y z t a)
 :precision binary64
 (+ (log (+ x y)) (- (fma (+ a -0.5) (log t) (log z)) t)))
double code(double x, double y, double z, double t, double a) {
	return ((log((x + y)) + log(z)) - t) + ((a - 0.5) * log(t));
}
↓
double code(double x, double y, double z, double t, double a) {
	return log((x + y)) + (fma((a + -0.5), log(t), log(z)) - t);
}

function code(x, y, z, t, a)
	return Float64(Float64(Float64(log(Float64(x + y)) + log(z)) - t) + Float64(Float64(a - 0.5) * log(t)))
end
↓
function code(x, y, z, t, a)
	return Float64(log(Float64(x + y)) + Float64(fma(Float64(a + -0.5), log(t), log(z)) - t))
end

code[x_, y_, z_, t_, a_] := N[(N[(N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[Log[z], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision] + N[(N[(a - 0.5), $MachinePrecision] * N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
↓
code[x_, y_, z_, t_, a_] := N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[(N[(a + -0.5), $MachinePrecision] * N[Log[t], $MachinePrecision] + N[Log[z], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]

\left(\left(\log \left(x + y\right) + \log z\right) - t\right) + \left(a - 0.5\right) \cdot \log t
↓
\log \left(x + y\right) + \left(\mathsf{fma}\left(a + -0.5, \log t, \log z\right) - t\right)

Error?

Target

Original0.3
Target0.3
Herbie0.3
\[\log \left(x + y\right) + \left(\left(\log z - t\right) + \left(a - 0.5\right) \cdot \log t\right) \]

Derivation?

  1. Initial program 0.3

    \[\left(\left(\log \left(x + y\right) + \log z\right) - t\right) + \left(a - 0.5\right) \cdot \log t \]

  2. Simplified0.3

    \[\leadsto \color{blue}{\log \left(x + y\right) + \left(\mathsf{fma}\left(a + -0.5, \log t, \log z\right) - t\right)} \]
    Proof

    [Start]0.3

    \[ \left(\left(\log \left(x + y\right) + \log z\right) - t\right) + \left(a - 0.5\right) \cdot \log t \]

    associate-+l- [=>]0.3

    \[ \color{blue}{\left(\log \left(x + y\right) + \log z\right) - \left(t - \left(a - 0.5\right) \cdot \log t\right)} \]

    associate--l+ [=>]0.3

    \[ \color{blue}{\log \left(x + y\right) + \left(\log z - \left(t - \left(a - 0.5\right) \cdot \log t\right)\right)} \]

    associate-+l- [<=]0.3

    \[ \log \left(x + y\right) + \color{blue}{\left(\left(\log z - t\right) + \left(a - 0.5\right) \cdot \log t\right)} \]

    +-commutative [=>]0.3

    \[ \log \left(x + y\right) + \color{blue}{\left(\left(a - 0.5\right) \cdot \log t + \left(\log z - t\right)\right)} \]

    associate-+r- [=>]0.3

    \[ \log \left(x + y\right) + \color{blue}{\left(\left(\left(a - 0.5\right) \cdot \log t + \log z\right) - t\right)} \]

    fma-def [=>]0.3

    \[ \log \left(x + y\right) + \left(\color{blue}{\mathsf{fma}\left(a - 0.5, \log t, \log z\right)} - t\right) \]

    sub-neg [=>]0.3

    \[ \log \left(x + y\right) + \left(\mathsf{fma}\left(\color{blue}{a + \left(-0.5\right)}, \log t, \log z\right) - t\right) \]

    metadata-eval [=>]0.3

    \[ \log \left(x + y\right) + \left(\mathsf{fma}\left(a + \color{blue}{-0.5}, \log t, \log z\right) - t\right) \]

  3. Final simplification0.3

    \[\leadsto \log \left(x + y\right) + \left(\mathsf{fma}\left(a + -0.5, \log t, \log z\right) - t\right) \]

Alternatives

Alternative 1
Error0.3
Cost20032
\[\left(\log \left(x + y\right) + \left(\log z - t\right)\right) + \left(a + -0.5\right) \cdot \log t \]
Alternative 2
Error12.6
Cost19908
\[\begin{array}{l} \mathbf{if}\;t \leq 0.00095:\\ \;\;\;\;\log z + \left(\left(a + -0.5\right) \cdot \log t + \log y\right)\\ \mathbf{else}:\\ \;\;\;\;\log \left(x + y\right) + \left(a \cdot \log t - t\right)\\ \end{array} \]
Alternative 3
Error20.0
Cost19904
\[\left(\log z - t\right) + \left(\left(a + -0.5\right) \cdot \log t + \log y\right) \]
Alternative 4
Error14.8
Cost13905
\[\begin{array}{l} \mathbf{if}\;t \leq 3.8 \cdot 10^{-222}:\\ \;\;\;\;\left(a + -0.5\right) \cdot \log t - t\\ \mathbf{elif}\;t \leq 4.6 \cdot 10^{-216}:\\ \;\;\;\;\log \left(\left(x + y\right) \cdot \frac{z}{\sqrt{t}}\right)\\ \mathbf{elif}\;t \leq 2.4 \cdot 10^{-150} \lor \neg \left(t \leq 5.6 \cdot 10^{-143}\right):\\ \;\;\;\;\log \left(x + y\right) + \left(a \cdot \log t - t\right)\\ \mathbf{else}:\\ \;\;\;\;\log \left(\left(x + y\right) \cdot z\right) + -0.5 \cdot \log t\\ \end{array} \]
Alternative 5
Error8.8
Cost13764
\[\begin{array}{l} \mathbf{if}\;t \leq 0.00092:\\ \;\;\;\;\left(a + -0.5\right) \cdot \log t + \left(\log \left(\left(x + y\right) \cdot z\right) - t\right)\\ \mathbf{else}:\\ \;\;\;\;\log \left(x + y\right) + \left(a \cdot \log t - t\right)\\ \end{array} \]
Alternative 6
Error14.8
Cost13640
\[\begin{array}{l} \mathbf{if}\;t \leq 3.6 \cdot 10^{-222}:\\ \;\;\;\;\left(a + -0.5\right) \cdot \log t - t\\ \mathbf{elif}\;t \leq 6.2 \cdot 10^{-216}:\\ \;\;\;\;\log \left(\left(x + y\right) \cdot \frac{z}{\sqrt{t}}\right)\\ \mathbf{else}:\\ \;\;\;\;\log \left(x + y\right) + \left(a \cdot \log t - t\right)\\ \end{array} \]
Alternative 7
Error8.9
Cost13636
\[\begin{array}{l} \mathbf{if}\;t \leq 0.0008:\\ \;\;\;\;\left(a + -0.5\right) \cdot \log t + \log \left(\left(x + y\right) \cdot z\right)\\ \mathbf{else}:\\ \;\;\;\;\log \left(x + y\right) + \left(a \cdot \log t - t\right)\\ \end{array} \]
Alternative 8
Error15.0
Cost13513
\[\begin{array}{l} \mathbf{if}\;t \leq 3.8 \cdot 10^{-222} \lor \neg \left(t \leq 1.55 \cdot 10^{-215}\right):\\ \;\;\;\;\left(a + -0.5\right) \cdot \log t - t\\ \mathbf{else}:\\ \;\;\;\;\log \left(\left(x + y\right) \cdot \frac{z}{\sqrt{t}}\right)\\ \end{array} \]
Alternative 9
Error17.5
Cost13508
\[\begin{array}{l} \mathbf{if}\;t \leq 0.00068:\\ \;\;\;\;\left(a + -0.5\right) \cdot \log t + \log \left(y \cdot z\right)\\ \mathbf{else}:\\ \;\;\;\;\log \left(x + y\right) + \left(a \cdot \log t - t\right)\\ \end{array} \]
Alternative 10
Error14.8
Cost6848
\[\left(a + -0.5\right) \cdot \log t - t \]
Alternative 11
Error24.1
Cost6724
\[\begin{array}{l} \mathbf{if}\;t \leq 10^{+25}:\\ \;\;\;\;a \cdot \log t\\ \mathbf{else}:\\ \;\;\;\;-t\\ \end{array} \]
Alternative 12
Error40.0
Cost128
\[-t \]

Error

Reproduce?

herbie shell --seed 2023023 
(FPCore (x y z t a)
  :name "Numeric.SpecFunctions:logGammaL from math-functions-0.1.5.2"
  :precision binary64

  :herbie-target
  (+ (log (+ x y)) (+ (- (log z) t) (* (- a 0.5) (log t))))

  (+ (- (+ (log (+ x y)) (log z)) t) (* (- a 0.5) (log t))))