Average Error: 9.1 → 0.5
Time: 8.3s
Precision: 64
\[\left(x \cdot \log y + z \cdot \log \left(1 - y\right)\right) - t\]
↓
\[\left(\log \left(\sqrt[3]{y} \cdot \sqrt[3]{y}\right) \cdot x + \left(\log \left(\sqrt[3]{y}\right) \cdot x + z \cdot \left(\log 1 - \left(1 \cdot y + \frac{1}{2} \cdot \frac{{y}^{2}}{{1}^{2}}\right)\right)\right)\right) - t\]
\left(x \cdot \log y + z \cdot \log \left(1 - y\right)\right) - t
↓
\left(\log \left(\sqrt[3]{y} \cdot \sqrt[3]{y}\right) \cdot x + \left(\log \left(\sqrt[3]{y}\right) \cdot x + z \cdot \left(\log 1 - \left(1 \cdot y + \frac{1}{2} \cdot \frac{{y}^{2}}{{1}^{2}}\right)\right)\right)\right) - t
double f(double x, double y, double z, double t) {
        double r448798 = x;
        double r448799 = y;
        double r448800 = log(r448799);
        double r448801 = r448798 * r448800;
        double r448802 = z;
        double r448803 = 1.0;
        double r448804 = r448803 - r448799;
        double r448805 = log(r448804);
        double r448806 = r448802 * r448805;
        double r448807 = r448801 + r448806;
        double r448808 = t;
        double r448809 = r448807 - r448808;
        return r448809;
}
↓
double f(double x, double y, double z, double t) {
        double r448810 = y;
        double r448811 = cbrt(r448810);
        double r448812 = r448811 * r448811;
        double r448813 = log(r448812);
        double r448814 = x;
        double r448815 = r448813 * r448814;
        double r448816 = log(r448811);
        double r448817 = r448816 * r448814;
        double r448818 = z;
        double r448819 = 1.0;
        double r448820 = log(r448819);
        double r448821 = r448819 * r448810;
        double r448822 = 0.5;
        double r448823 = 2.0;
        double r448824 = pow(r448810, r448823);
        double r448825 = pow(r448819, r448823);
        double r448826 = r448824 / r448825;
        double r448827 = r448822 * r448826;
        double r448828 = r448821 + r448827;
        double r448829 = r448820 - r448828;
        double r448830 = r448818 * r448829;
        double r448831 = r448817 + r448830;
        double r448832 = r448815 + r448831;
        double r448833 = t;
        double r448834 = r448832 - r448833;
        return r448834;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Bits error versus t

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Target

Original9.1
Target0.3
Herbie0.5
\[\left(-z\right) \cdot \left(\left(0.5 \cdot \left(y \cdot y\right) + y\right) + \frac{0.333333333333333315}{1 \cdot \left(1 \cdot 1\right)} \cdot \left(y \cdot \left(y \cdot y\right)\right)\right) - \left(t - x \cdot \log y\right)\]

Derivation

  1. Initial program 9.1

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

  2. Taylor expanded around 0 0.4

    \[\leadsto \left(x \cdot \log y + z \cdot \color{blue}{\left(\log 1 - \left(1 \cdot y + \frac{1}{2} \cdot \frac{{y}^{2}}{{1}^{2}}\right)\right)}\right) - t\]
  3. Using strategy rm

  4. Applied add-cube-cbrt0.4

    \[\leadsto \left(x \cdot \log \color{blue}{\left(\left(\sqrt[3]{y} \cdot \sqrt[3]{y}\right) \cdot \sqrt[3]{y}\right)} + z \cdot \left(\log 1 - \left(1 \cdot y + \frac{1}{2} \cdot \frac{{y}^{2}}{{1}^{2}}\right)\right)\right) - t\]

  5. Applied log-prod0.5

    \[\leadsto \left(x \cdot \color{blue}{\left(\log \left(\sqrt[3]{y} \cdot \sqrt[3]{y}\right) + \log \left(\sqrt[3]{y}\right)\right)} + z \cdot \left(\log 1 - \left(1 \cdot y + \frac{1}{2} \cdot \frac{{y}^{2}}{{1}^{2}}\right)\right)\right) - t\]

  6. Applied distribute-rgt-in0.5

    \[\leadsto \left(\color{blue}{\left(\log \left(\sqrt[3]{y} \cdot \sqrt[3]{y}\right) \cdot x + \log \left(\sqrt[3]{y}\right) \cdot x\right)} + z \cdot \left(\log 1 - \left(1 \cdot y + \frac{1}{2} \cdot \frac{{y}^{2}}{{1}^{2}}\right)\right)\right) - t\]

  7. Applied associate-+l+0.5

    \[\leadsto \color{blue}{\left(\log \left(\sqrt[3]{y} \cdot \sqrt[3]{y}\right) \cdot x + \left(\log \left(\sqrt[3]{y}\right) \cdot x + z \cdot \left(\log 1 - \left(1 \cdot y + \frac{1}{2} \cdot \frac{{y}^{2}}{{1}^{2}}\right)\right)\right)\right)} - t\]

  8. Final simplification0.5

    \[\leadsto \left(\log \left(\sqrt[3]{y} \cdot \sqrt[3]{y}\right) \cdot x + \left(\log \left(\sqrt[3]{y}\right) \cdot x + z \cdot \left(\log 1 - \left(1 \cdot y + \frac{1}{2} \cdot \frac{{y}^{2}}{{1}^{2}}\right)\right)\right)\right) - t\]

Reproduce

herbie shell --seed 2020046 
(FPCore (x y z t)
  :name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, B"
  :precision binary64

  :herbie-target
  (- (* (- z) (+ (+ (* 0.5 (* y y)) y) (* (/ 0.3333333333333333 (* 1 (* 1 1))) (* y (* y y))))) (- t (* x (log y))))

  (- (+ (* x (log y)) (* z (log (- 1 y)))) t))