Details

Time bar (total: 8.1s)

analyze900.0ms (11.1%)

Algorithm
search
Search
TrueOtherFalseIter
0%99.9%0.1%0
0%99.9%0.1%1
0%99.9%0.1%2
0%99.9%0.1%3
0%99.9%0.1%4
0%99.9%0.1%5
12.5%87.4%0.1%6
20.3%73.3%6.4%7
27.3%66.3%6.4%8
38.2%53%8.7%9
46%40.6%13.4%10
49.9%35.1%15%11
54.9%26.9%18.1%12
57.4%22%20.6%13
59.7%18%22.4%14
Compiler

Compiled 42 to 27 computations (35.7% saved)

sample2.2s (27%)

Results
1.8s8256×body128valid
263.0ms1052×body128invalid
Compiler

Compiled 125 to 80 computations (36% saved)

simplify130.0ms (1.6%)

Algorithm
egg-herbie
Rules
1002×fma-def_binary64
782×fma-neg_binary64
748×distribute-rgt-in_binary64
560×div-sub_binary64
469×associate-/l*_binary64
Counts
1 → 2
Iterations

Useful iterations: 3 (0.0ms)

IterNodesCost
02258
14452
29452
321150
434150
559450
6108250
7170250
8206650
9287050
10349550
11350150
12350350
13419650
14424950
15488150

prune8.0ms (0.1%)

Pruning

3 alts after pruning (3 fresh and 0 done)

PrunedKeptTotal
New022
Fresh011
Picked000
Done000
Total033
Error
20.1b
Counts
3 → 2
Alt Table
StatusErrorProgram
20.2b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a)))
20.3b
(if (>=.f64 b 0) (*.f64 c (/.f64 -2 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))))) (*.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b) (/.f64 1/2 a)))
Compiler

Compiled 231 to 143 computations (38.1% saved)

localize21.0ms (0.3%)

Local error

Found 4 expressions with local error:

NewErrorProgram
0.1b
(/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a))
1.3b
(-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))
1.7b
(+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))
23.5b
(sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))
Compiler

Compiled 1536 to 1280 computations (16.7% saved)

series276.0ms (3.4%)

Counts
4 → 80
Calls

4 calls:

109.0ms
(/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a))
65.0ms
(-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))
57.0ms
(+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))
44.0ms
(sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))

rewrite49.0ms (0.6%)

Algorithm
rewrite-expression-head
Rules
69×*-un-lft-identity_binary64
60×prod-diff_binary64
58×add-sqr-sqrt_binary64
53×add-cube-cbrt_binary64
37×sqrt-prod_binary64
Counts
4 → 169
Calls

4 calls:

9.0ms
(/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a))
8.0ms
(-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))
7.0ms
(+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))
4.0ms
(sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))

simplify219.0ms (2.7%)

Algorithm
egg-herbie
Rules
465×fma-neg_binary64
405×cancel-sign-sub-inv_binary64
186×associate-*r*_binary64
177×sub-neg_binary64
170×unswap-sqr_binary64
Counts
249 → 280
Iterations

Useful iterations: 2 (0.0ms)

IterNodesCost
03108186
18216808
229505670
348255670
449995670
548775670

prune531.0ms (6.6%)

Pruning

12 alts after pruning (12 fresh and 0 done)

PrunedKeptTotal
New31711328
Fresh011
Picked101
Done000
Total31812330
Error
5.0b
Counts
330 → 12
Alt Table
StatusErrorProgram
20.2b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 1 (/.f64 (*.f64 2 a) (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))))))
16.7b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (-.f64 (/.f64 c b) (/.f64 b a)))
34.3b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (/.f64 (*.f64 4 (*.f64 c a)) (-.f64 (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b) b))) (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a)))
24.4b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) b)) (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a)))
20.3b
(if (>=.f64 b 0) (*.f64 c (/.f64 -2 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))))) (*.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b) (/.f64 1/2 a)))
26.1b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (/.f64 (+.f64 (pow.f64 (neg.f64 b) 3) (pow.f64 (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))) 3)) (+.f64 (*.f64 (neg.f64 b) (neg.f64 b)) (-.f64 (*.f64 (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))))) (*.f64 2 a)))
28.4b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (+.f64 (neg.f64 b) (/.f64 (sqrt.f64 (-.f64 (pow.f64 (*.f64 b b) 3) (pow.f64 (*.f64 (*.f64 4 a) c) 3))) (sqrt.f64 (+.f64 (*.f64 (*.f64 b b) (*.f64 b b)) (+.f64 (*.f64 (*.f64 (*.f64 4 a) c) (*.f64 (*.f64 4 a) c)) (*.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))))) (*.f64 2 a)))
20.5b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (*.f64 (cbrt.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (cbrt.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))))) (/.f64 (*.f64 2 a) (cbrt.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))))))
31.4b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (-.f64 (*.f64 (neg.f64 b) (neg.f64 b)) (*.f64 (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (*.f64 (-.f64 (neg.f64 b) (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b)) (*.f64 a 2))))
27.2b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (exp.f64 (log.f64 (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b))))) (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a)))
16.8b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (+.f64 (neg.f64 b) (*.f64 -1 b)) (*.f64 2 a)))
26.4b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (*.f64 (*.f64 (cbrt.f64 (-.f64 (neg.f64 b) (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b))) (cbrt.f64 (-.f64 (neg.f64 b) (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b)))) (cbrt.f64 (-.f64 (neg.f64 b) (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b))))) (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a)))
Compiler

Compiled 19294 to 11038 computations (42.8% saved)

localize23.0ms (0.3%)

Local error

Found 4 expressions with local error:

NewErrorProgram
0.0b
(-.f64 (/.f64 c b) (/.f64 b a))
0.1b
(/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))))
1.3b
(-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))
23.5b
(sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))
Compiler

Compiled 1536 to 1280 computations (16.7% saved)

series293.0ms (3.6%)

Counts
2 → 29
Calls

2 calls:

281.0ms
(/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))))
12.0ms
(-.f64 (/.f64 c b) (/.f64 b a))

rewrite71.0ms (0.9%)

Algorithm
rewrite-expression-head
Rules
265×times-frac_binary64
220×add-sqr-sqrt_binary64
214×*-un-lft-identity_binary64
200×add-cube-cbrt_binary64
169×prod-diff_binary64
Counts
2 → 248
Calls

2 calls:

25.0ms
(-.f64 (/.f64 c b) (/.f64 b a))
9.0ms
(/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))))

simplify209.0ms (2.6%)

Algorithm
egg-herbie
Rules
459×fma-neg_binary64
332×associate-/l/_binary64
175×distribute-rgt-neg-in_binary64
163×sub-neg_binary64
156×fma-udef_binary64
Counts
277 → 308
Iterations

Useful iterations: 2 (0.0ms)

IterNodesCost
04298755
110316629
226215171
348965171
449315171

prune924.0ms (11.4%)

Pruning

13 alts after pruning (13 fresh and 0 done)

PrunedKeptTotal
New4604464
Fresh2911
Picked101
Done000
Total46313476
Error
5.0b
Counts
476 → 13
Alt Table
StatusErrorProgram
22.8b
(if (>=.f64 b 0) (*.f64 (*.f64 (cbrt.f64 (/.f64 (*.f64 c 2) (-.f64 (neg.f64 b) (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))) (cbrt.f64 (/.f64 (*.f64 c 2) (-.f64 (neg.f64 b) (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))))) (cbrt.f64 (/.f64 (*.f64 c 2) (-.f64 (neg.f64 b) (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))))) (-.f64 (/.f64 c b) (/.f64 b a)))
28.4b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (+.f64 (neg.f64 b) (/.f64 (sqrt.f64 (-.f64 (pow.f64 (*.f64 b b) 3) (pow.f64 (*.f64 (*.f64 4 a) c) 3))) (sqrt.f64 (+.f64 (*.f64 (*.f64 b b) (*.f64 b b)) (+.f64 (*.f64 (*.f64 (*.f64 4 a) c) (*.f64 (*.f64 4 a) c)) (*.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))))) (*.f64 2 a)))
16.7b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (+.f64 (-.f64 (pow.f64 (/.f64 (cbrt.f64 c) (cbrt.f64 b)) 3) (/.f64 b a)) 0))
20.2b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 1 (/.f64 (*.f64 2 a) (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))))))
20.9b
(if (>=.f64 b 0) (*.f64 -1 (/.f64 c b)) (-.f64 (/.f64 c b) (/.f64 b a)))
34.3b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (/.f64 (*.f64 4 (*.f64 c a)) (-.f64 (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b) b))) (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a)))
20.3b
(if (>=.f64 b 0) (*.f64 c (/.f64 -2 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))))) (*.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b) (/.f64 1/2 a)))
26.1b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (/.f64 (+.f64 (pow.f64 (neg.f64 b) 3) (pow.f64 (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))) 3)) (+.f64 (*.f64 (neg.f64 b) (neg.f64 b)) (-.f64 (*.f64 (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))))) (*.f64 2 a)))
24.0b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (expm1.f64 (log1p.f64 (-.f64 (/.f64 c b) (/.f64 b a)))))
20.5b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (*.f64 (cbrt.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (cbrt.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))))) (/.f64 (*.f64 2 a) (cbrt.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))))))
31.4b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (-.f64 (*.f64 (neg.f64 b) (neg.f64 b)) (*.f64 (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (*.f64 (-.f64 (neg.f64 b) (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b)) (*.f64 a 2))))
27.2b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (exp.f64 (log.f64 (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b))))) (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a)))
16.8b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (+.f64 (neg.f64 b) (*.f64 -1 b)) (*.f64 2 a)))
Compiler

Compiled 26148 to 10976 computations (58% saved)

localize26.0ms (0.3%)

Local error

Found 4 expressions with local error:

NewErrorProgram
0.6b
(cbrt.f64 b)
0.6b
(cbrt.f64 c)
1.3b
(-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))
23.5b
(sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))
Compiler

Compiled 1536 to 1280 computations (16.7% saved)

series63.0ms (0.8%)

Counts
2 → 24
Calls

2 calls:

34.0ms
(cbrt.f64 c)
29.0ms
(cbrt.f64 b)

rewrite2.0ms (0%)

Algorithm
rewrite-expression-head
Rules
cbrt-prod_binary64
*-un-lft-identity_binary64
add-sqr-sqrt_binary64
add-cube-cbrt_binary64
add-cbrt-cube_binary64
Counts
2 → 26
Calls

2 calls:

1.0ms
(cbrt.f64 c)
1.0ms
(cbrt.f64 b)

simplify33.0ms (0.4%)

Algorithm
egg-herbie
Rules
422×distribute-rgt-in_binary64
414×distribute-lft-in_binary64
409×associate-*l*_binary64
378×unswap-sqr_binary64
274×associate-*r*_binary64
Counts
50 → 36
Iterations

Useful iterations: 4 (0.0ms)

IterNodesCost
042286
164274
2148274
3384274
41209266
54945266

prune448.0ms (5.5%)

Pruning

13 alts after pruning (13 fresh and 0 done)

PrunedKeptTotal
New1902192
Fresh11112
Picked101
Done000
Total19213205
Error
5.0b
Counts
205 → 13
Alt Table
StatusErrorProgram
22.8b
(if (>=.f64 b 0) (*.f64 (*.f64 (cbrt.f64 (/.f64 (*.f64 c 2) (-.f64 (neg.f64 b) (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))) (cbrt.f64 (/.f64 (*.f64 c 2) (-.f64 (neg.f64 b) (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))))) (cbrt.f64 (/.f64 (*.f64 c 2) (-.f64 (neg.f64 b) (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))))) (-.f64 (/.f64 c b) (/.f64 b a)))
28.4b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (+.f64 (neg.f64 b) (/.f64 (sqrt.f64 (-.f64 (pow.f64 (*.f64 b b) 3) (pow.f64 (*.f64 (*.f64 4 a) c) 3))) (sqrt.f64 (+.f64 (*.f64 (*.f64 b b) (*.f64 b b)) (+.f64 (*.f64 (*.f64 (*.f64 4 a) c) (*.f64 (*.f64 4 a) c)) (*.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))))) (*.f64 2 a)))
20.9b
(if (>=.f64 b 0) (*.f64 -1 (/.f64 c b)) (-.f64 (/.f64 c b) (/.f64 b a)))
34.3b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (/.f64 (*.f64 4 (*.f64 c a)) (-.f64 (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b) b))) (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a)))
20.3b
(if (>=.f64 b 0) (*.f64 c (/.f64 -2 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))))) (*.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b) (/.f64 1/2 a)))
26.1b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (/.f64 (+.f64 (pow.f64 (neg.f64 b) 3) (pow.f64 (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))) 3)) (+.f64 (*.f64 (neg.f64 b) (neg.f64 b)) (-.f64 (*.f64 (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))))) (*.f64 2 a)))
27.4b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (+.f64 (-.f64 (pow.f64 (/.f64 (pow.f64 c 1/3) (cbrt.f64 b)) 3) (/.f64 b a)) 0))
20.5b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (*.f64 (cbrt.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (cbrt.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))))) (/.f64 (*.f64 2 a) (cbrt.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))))))
31.4b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (-.f64 (*.f64 (neg.f64 b) (neg.f64 b)) (*.f64 (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (*.f64 (-.f64 (neg.f64 b) (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b)) (*.f64 a 2))))
27.2b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (exp.f64 (log.f64 (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b))))) (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a)))
20.2b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 1 (/.f64 (*.f64 2 a) (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))))))
16.8b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (+.f64 (neg.f64 b) (*.f64 -1 b)) (*.f64 2 a)))
16.7b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (+.f64 (-.f64 (pow.f64 (/.f64 (cbrt.f64 c) (*.f64 (pow.f64 (*.f64 -1 b) 1/3) (cbrt.f64 -1))) 3) (/.f64 b a)) 0))
Compiler

Compiled 11981 to 6604 computations (44.9% saved)

localize32.0ms (0.4%)

Local error

Found 4 expressions with local error:

NewErrorProgram
0.6b
(cbrt.f64 c)
1.3b
(-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))
2.0b
(pow.f64 (*.f64 -1 b) 1/3)
23.5b
(sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))
Compiler

Compiled 1536 to 1280 computations (16.7% saved)

series23.0ms (0.3%)

Counts
1 → 8
Calls

1 calls:

23.0ms
(pow.f64 (*.f64 -1 b) 1/3)

rewrite3.0ms (0%)

Algorithm
rewrite-expression-head
Rules
add-exp-log_binary64
pow1_binary64
pow-unpow_binary64
pow-exp_binary64
pow-pow_binary64
Counts
1 → 20
Calls

1 calls:

2.0ms
(pow.f64 (*.f64 -1 b) 1/3)

simplify55.0ms (0.7%)

Algorithm
egg-herbie
Rules
538×fma-def_binary64
526×log-prod_binary64
432×fma-neg_binary64
346×associate-*l*_binary64
253×associate-*r*_binary64
Counts
28 → 32
Iterations

Useful iterations: 2 (0.0ms)

IterNodesCost
043207
187175
2152171
3437171
41654171
54539171
64847171

prune515.0ms (6.4%)

Pruning

13 alts after pruning (13 fresh and 0 done)

PrunedKeptTotal
New2051206
Fresh01212
Picked101
Done000
Total20613219
Error
5.0b
Counts
219 → 13
Alt Table
StatusErrorProgram
22.8b
(if (>=.f64 b 0) (*.f64 (*.f64 (cbrt.f64 (/.f64 (*.f64 c 2) (-.f64 (neg.f64 b) (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))) (cbrt.f64 (/.f64 (*.f64 c 2) (-.f64 (neg.f64 b) (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))))) (cbrt.f64 (/.f64 (*.f64 c 2) (-.f64 (neg.f64 b) (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))))) (-.f64 (/.f64 c b) (/.f64 b a)))
28.4b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (+.f64 (neg.f64 b) (/.f64 (sqrt.f64 (-.f64 (pow.f64 (*.f64 b b) 3) (pow.f64 (*.f64 (*.f64 4 a) c) 3))) (sqrt.f64 (+.f64 (*.f64 (*.f64 b b) (*.f64 b b)) (+.f64 (*.f64 (*.f64 (*.f64 4 a) c) (*.f64 (*.f64 4 a) c)) (*.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))))) (*.f64 2 a)))
20.9b
(if (>=.f64 b 0) (*.f64 -1 (/.f64 c b)) (-.f64 (/.f64 c b) (/.f64 b a)))
34.3b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (/.f64 (*.f64 4 (*.f64 c a)) (-.f64 (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b) b))) (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a)))
20.3b
(if (>=.f64 b 0) (*.f64 c (/.f64 -2 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))))) (*.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b) (/.f64 1/2 a)))
26.1b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (/.f64 (+.f64 (pow.f64 (neg.f64 b) 3) (pow.f64 (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))) 3)) (+.f64 (*.f64 (neg.f64 b) (neg.f64 b)) (-.f64 (*.f64 (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))))) (*.f64 2 a)))
27.4b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (+.f64 (-.f64 (pow.f64 (/.f64 (pow.f64 c 1/3) (cbrt.f64 b)) 3) (/.f64 b a)) 0))
20.5b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (*.f64 (cbrt.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (cbrt.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))))) (/.f64 (*.f64 2 a) (cbrt.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))))))
31.4b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (-.f64 (*.f64 (neg.f64 b) (neg.f64 b)) (*.f64 (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (*.f64 (-.f64 (neg.f64 b) (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b)) (*.f64 a 2))))
27.2b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (exp.f64 (log.f64 (hypot.f64 (sqrt.f64 (*.f64 (*.f64 c a) -4)) b))))) (/.f64 (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))) (*.f64 2 a)))
20.2b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 1 (/.f64 (*.f64 2 a) (+.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c)))))))
16.8b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (/.f64 (+.f64 (neg.f64 b) (*.f64 -1 b)) (*.f64 2 a)))
16.7b
(if (>=.f64 b 0) (/.f64 (*.f64 2 c) (-.f64 (neg.f64 b) (sqrt.f64 (-.f64 (*.f64 b b) (*.f64 (*.f64 4 a) c))))) (+.f64 (-.f64 (pow.f64 (/.f64 (cbrt.f64 c) (*.f64 (pow.f64 (pow.f64 (neg.f64 b) (*.f64 (cbrt.f64 1/3) (cbrt.f64 1/3))) (cbrt.f64 1/3)) (cbrt.f64 -1))) 3) (/.f64 b a)) 0))
Compiler

Compiled 13069 to 8062 computations (38.3% saved)

regimes605.0ms (7.5%)

Accuracy

Total 1.1b remaining (15.8%)

Threshold costs 0b (0%)

Counts
72 → 3
Compiler

Compiled 12774 to 8094 computations (36.6% saved)

bsearch79.0ms (1%)

Steps
ItersRangePoint
9
9.317915900075717e+54
6.212821628630246e+62
6.108410036495975e+62
8
-1.3016370291905524e+112
-7.914045696574179e+108
-1.8019041569772403e+111
Compiler

Compiled 4 to 3 computations (25% saved)

simplify6.0ms (0.1%)

Algorithm
egg-herbie
Rules
16×*-commutative_binary64
15×+-commutative_binary64
13×sub-neg_binary64
10×neg-sub0_binary64
neg-mul-1_binary64
Iterations

Useful iterations: 1 (0.0ms)

IterNodesCost
043132
165130
295130
3114130
4125130
5128130
6127130

end356.0ms (4.4%)

Compiler

Compiled 851 to 534 computations (37.3% saved)

Profiling

Loading profile data...