Average Error: 7.0 → 6.9
Time: 12.2s
Precision: 64
Internal Precision: 128
\[\left(x.re \cdot x.re - x.im \cdot x.im\right) \cdot x.re - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im\]
\[(\left(x.re \cdot -3\right) \cdot \left(x.im \cdot x.im\right) + \left({x.re}^{3}\right))_*\]

Error

Bits error versus x.re

Bits error versus x.im

Target

Original7.0
Target0.3
Herbie6.9
\[\left(x.re \cdot x.re\right) \cdot \left(x.re - x.im\right) + \left(x.re \cdot x.im\right) \cdot \left(x.re - 3 \cdot x.im\right)\]

Derivation

  1. Initial program 7.0

    \[\left(x.re \cdot x.re - x.im \cdot x.im\right) \cdot x.re - \left(x.re \cdot x.im + x.im \cdot x.re\right) \cdot x.im\]
  2. Initial simplification7.0

    \[\leadsto (x.re \cdot \left((\left(x.im + x.im\right) \cdot \left(-x.im\right) + \left(x.re \cdot x.re\right))_*\right) + \left(\left(x.re \cdot x.im\right) \cdot \left(-x.im\right)\right))_*\]
  3. Taylor expanded around 0 6.9

    \[\leadsto \color{blue}{{x.re}^{3} - 3 \cdot \left({x.im}^{2} \cdot x.re\right)}\]
  4. Simplified6.9

    \[\leadsto \color{blue}{(\left(x.re \cdot -3\right) \cdot \left(x.im \cdot x.im\right) + \left({x.re}^{3}\right))_*}\]
  5. Final simplification6.9

    \[\leadsto (\left(x.re \cdot -3\right) \cdot \left(x.im \cdot x.im\right) + \left({x.re}^{3}\right))_*\]

Reproduce

herbie shell --seed 2018360 +o rules:numerics
(FPCore (x.re x.im)
  :name "math.cube on complex, real part"

  :herbie-target
  (+ (* (* x.re x.re) (- x.re x.im)) (* (* x.re x.im) (- x.re (* 3 x.im))))

  (- (* (- (* x.re x.re) (* x.im x.im)) x.re) (* (+ (* x.re x.im) (* x.im x.re)) x.im)))

Details

Time bar (total: 8.2s)Debug log

start116.0ms

Algorithm
intervals

setup335.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 8.8b

localize27.0ms

Local error

Found 2 expressions with local error:

8.8b
(fma x.re (fma (+ x.im x.im) (- x.im) (* x.re x.re)) (* (* x.re x.im) (- x.im)))
0.2b
(* (* x.re x.im) (- x.im))

rewrite4.0ms

Algorithm
rewrite-expression-head
Counts
2 → 29
Calls

2 calls. Slowest were:

4.0ms
(* (* x.re x.im) (- x.im))
0.0ms
(fma x.re (fma (+ x.im x.im) (- x.im) (* x.re x.re)) (* (* x.re x.im) (- x.im)))

series65.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

44.0ms
(fma x.re (fma (+ x.im x.im) (- x.im) (* x.re x.re)) (* (* x.re x.im) (- x.im)))
20.0ms
(* (* x.re x.im) (- x.im))

simplify933.0ms

Counts
12 → 35
Calls

12 calls. Slowest were:

493.0ms
(* (* (* (* x.re x.im) (* x.re x.im)) (* x.re x.im)) (* (* (- x.im) (- x.im)) (- x.im)))
120.0ms
(- (pow x.re 3) (* 3 (* (pow x.im 2) x.re)))
114.0ms
(- (pow x.re 3) (* 3 (* (pow x.im 2) x.re)))

prune440.0ms

Pruning

6 alts after pruning (6 fresh and 0 done)

Merged error: 8.7b

localize18.0ms

Local error

Found 1 expressions with local error:

8.8b
(fma (* x.re -3) (* x.im x.im) (pow x.re 3))

rewrite1.0ms

Algorithm
rewrite-expression-head
Counts
1 → 10
Calls

1 calls. Slowest were:

0.0ms
(fma (* x.re -3) (* x.im x.im) (pow x.re 3))

series40.0ms

Counts
1 → 3
Calls

1 calls. Slowest were:

40.0ms
(fma (* x.re -3) (* x.im x.im) (pow x.re 3))

simplify291.0ms

Counts
3 → 13
Calls

3 calls. Slowest were:

142.0ms
(- (pow x.re 3) (* 3 (* (pow x.im 2) x.re)))
78.0ms
(- (pow x.re 3) (* 3 (* (pow x.im 2) x.re)))
70.0ms
(- (pow x.re 3) (* 3 (* (pow x.im 2) x.re)))

prune157.0ms

Pruning

6 alts after pruning (5 fresh and 1 done)

Merged error: 8.7b

localize5.0ms

Local error

Found 3 expressions with local error:

8.8b
(* x.re (fma (+ x.im x.im) (- x.im) (* x.re x.re)))
0.2b
(* (* x.re x.im) (- x.im))
0.2b
(+ (* x.re (fma (+ x.im x.im) (- x.im) (* x.re x.re))) (* (* x.re x.im) (- x.im)))

rewrite9.0ms

Algorithm
rewrite-expression-head
Counts
3 → 52
Calls

3 calls. Slowest were:

3.0ms
(+ (* x.re (fma (+ x.im x.im) (- x.im) (* x.re x.re))) (* (* x.re x.im) (- x.im)))
3.0ms
(* (* x.re x.im) (- x.im))
1.0ms
(* x.re (fma (+ x.im x.im) (- x.im) (* x.re x.re)))

series121.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

48.0ms
(* x.re (fma (+ x.im x.im) (- x.im) (* x.re x.re)))
45.0ms
(+ (* x.re (fma (+ x.im x.im) (- x.im) (* x.re x.re))) (* (* x.re x.im) (- x.im)))
27.0ms
(* (* x.re x.im) (- x.im))

simplify1.7s

Counts
25 → 61
Calls

25 calls. Slowest were:

503.0ms
(* (* (* (* x.re x.im) (* x.re x.im)) (* x.re x.im)) (* (* (- x.im) (- x.im)) (- x.im)))
218.0ms
(* (exp (* x.re (fma (+ x.im x.im) (- x.im) (* x.re x.re)))) (exp (* (* x.re x.im) (- x.im))))
127.0ms
(- (pow x.re 3) (* 2 (* (pow x.im 2) x.re)))

prune545.0ms

Pruning

6 alts after pruning (4 fresh and 2 done)

Merged error: 8.7b

localize12.0ms

Local error

Found 4 expressions with local error:

22.5b
(cbrt (* (* x.re x.im) (- x.im)))
22.5b
(cbrt (* (* x.re x.im) (- x.im)))
22.5b
(cbrt (* (* x.re x.im) (- x.im)))
8.8b
(fma x.re (fma (+ x.im x.im) (- x.im) (* x.re x.re)) (* (* (cbrt (* (* x.re x.im) (- x.im))) (cbrt (* (* x.re x.im) (- x.im)))) (cbrt (* (* x.re x.im) (- x.im)))))

rewrite13.0ms

Algorithm
rewrite-expression-head
Counts
4 → 43
Calls

4 calls. Slowest were:

5.0ms
(cbrt (* (* x.re x.im) (- x.im)))
4.0ms
(cbrt (* (* x.re x.im) (- x.im)))
4.0ms
(cbrt (* (* x.re x.im) (- x.im)))

series479.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

144.0ms
(cbrt (* (* x.re x.im) (- x.im)))
125.0ms
(cbrt (* (* x.re x.im) (- x.im)))
111.0ms
(cbrt (* (* x.re x.im) (- x.im)))
99.0ms
(fma x.re (fma (+ x.im x.im) (- x.im) (* x.re x.re)) (* (* (cbrt (* (* x.re x.im) (- x.im))) (cbrt (* (* x.re x.im) (- x.im)))) (cbrt (* (* x.re x.im) (- x.im)))))

simplify2.1s

Counts
15 → 55
Calls

15 calls. Slowest were:

346.0ms
(* (cbrt -1) (exp (* -1/3 (+ (* 2 (log (/ 1 x.im))) (log (/ 1 x.re))))))
310.0ms
(* (cbrt -1) (exp (* -1/3 (+ (* 2 (log (/ 1 x.im))) (log (/ 1 x.re))))))
284.0ms
(* (cbrt -1) (exp (* -1/3 (+ (* 2 (log (/ 1 x.im))) (log (/ 1 x.re))))))

prune687.0ms

Pruning

6 alts after pruning (4 fresh and 2 done)

Merged error: 8.6b

regimes158.0ms

Accuracy

0% (0.1b remaining)

Error of 6.9b against oracle of 6.8b and baseline of 6.9b

bsearch1.0ms