Details

Time bar (total: 9.3s)

analyze721.0ms (7.7%)

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
18.7%81.1%0.1%6
26.5%67.1%6.4%7
35.1%58.5%6.4%8
41.3%49.9%8.7%9
46.4%40.2%13.4%10
50.8%34.2%15%11
55.2%26.7%18.1%12
57.4%22%20.6%13
59.8%17.9%22.3%14
Compiler

Compiled 42 to 27 computations (35.7% saved)

sample148.0ms (1.6%)

Algorithm
intervals
Results
31.0ms256×body128valid
2.0ms16×body128nan
1.0ms10×body128invalid
Compiler

Compiled 86 to 56 computations (34.9% saved)

simplify160.0ms (1.7%)

Algorithm
egg-herbie
Rules
1753×fma-neg_binary64
1288×fma-def_binary64
631×div-sub_binary64
512×times-frac_binary64
430×associate-/l*_binary64
Counts
1 → 2
Iterations

Useful iterations: 3 (0.0ms)

IterNodesCost
02258
14752
211852
326950
442950
574350
6131550
7173850
8217650
9285350
10406650
11352350
12357550
13354350
14355750
15451050
16462250
17476150

prune6.0ms (0.1%)

Pruning

3 alts after pruning (3 fresh and 0 done)

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

Compiled 231 to 143 computations (38.1% saved)

localize26.0ms (0.3%)

Local error

Found 4 expressions with local error:

NewErrorProgram
0.1b
(/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)
0.9b
(-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)
1.5b
(+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))
21.6b
(sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))
Compiler

Compiled 1530 to 1275 computations (16.7% saved)

series249.0ms (2.7%)

Counts
4 → 80
Calls

4 calls:

84.0ms
(/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)
58.0ms
(-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)
57.0ms
(+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))
49.0ms
(sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))

rewrite23.0ms (0.2%)

Algorithm
rewrite-expression-head
Rules
42×*-un-lft-identity_binary64
35×add-sqr-sqrt_binary64
30×add-cube-cbrt_binary64
18×prod-diff_binary64
16×sqrt-prod_binary64
Counts
4 → 110
Calls

4 calls:

5.0ms
(/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)
4.0ms
(-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)
3.0ms
(+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))
1.0ms
(sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))

simplify198.0ms (2.1%)

Algorithm
egg-herbie
Rules
878×fma-neg_binary64
423×cancel-sign-sub-inv_binary64
171×distribute-rgt-in_binary64
160×distribute-rgt-neg-in_binary64
153×fma-def_binary64
Counts
190 → 240
Iterations

Useful iterations: 2 (0.0ms)

IterNodesCost
02515471
16384667
221874480
344024480
449474480
549864480
648724480

prune475.0ms (5.1%)

Pruning

13 alts after pruning (13 fresh and 0 done)

PrunedKeptTotal
New27113284
Fresh101
Picked101
Done000
Total27313286
Error
3.6b
Counts
286 → 13
Alt Table
StatusErrorProgram
20.3b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b b) a)) (/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
19.2b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 (fabs.f64 (cbrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (cbrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))) b)))
28.4b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 (/.f64 (sqrt.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))) (sqrt.f64 a)) (/.f64 (sqrt.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))) (sqrt.f64 a)))) (/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
21.0b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (exp.f64 (log.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))) a)) (/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
20.2b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 (/.f64 (*.f64 (cbrt.f64 (+.f64 b (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 (+.f64 b (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))) (*.f64 (cbrt.f64 a) (cbrt.f64 a))) (/.f64 (cbrt.f64 (+.f64 b (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 a)))) (/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
23.4b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (exp.f64 (log.f64 (-.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))) b)))))
22.1b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 (*.f64 (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) b)))
20.1b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 (/.f64 (sqrt.f64 (+.f64 b (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (*.f64 (cbrt.f64 a) (cbrt.f64 a))) (/.f64 (sqrt.f64 (+.f64 b (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 a)))) (/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
19.0b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (sqrt.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))) (/.f64 a (sqrt.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))))))) (/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
19.0b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (/.f64 1 a))) (/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
29.4b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (cbrt.f64 (pow.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))) 3)) b)))
20.1b
(if (>=.f64 b 0) (*.f64 -1/2 (-.f64 (*.f64 2 (/.f64 b a)) (*.f64 2 (/.f64 c b)))) (/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
17.7b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
Compiler

Compiled 14150 to 8300 computations (41.3% saved)

localize24.0ms (0.3%)

Local error

Found 4 expressions with local error:

NewErrorProgram
0.0b
(fma.f64 a (*.f64 c -4) (*.f64 b b))
0.1b
(/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)
1.5b
(+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))
21.6b
(sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))
Compiler

Compiled 1530 to 1275 computations (16.7% saved)

series15.0ms (0.2%)

Counts
1 → 36
Calls

1 calls:

14.0ms
(fma.f64 a (*.f64 c -4) (*.f64 b b))

rewrite1.0ms (0%)

Algorithm
rewrite-expression-head
Rules
*-un-lft-identity_binary64
add-sqr-sqrt_binary64
add-cube-cbrt_binary64
add-cbrt-cube_binary64
add-exp-log_binary64
Counts
1 → 10
Calls

1 calls:

0.0ms
(fma.f64 a (*.f64 c -4) (*.f64 b b))

simplify77.0ms (0.8%)

Algorithm
egg-herbie
Rules
490×*-commutative_binary64
334×neg-sub0_binary64
314×neg-mul-1_binary64
291×associate-*l*_binary64
290×times-frac_binary64
Counts
46 → 24
Iterations

Useful iterations: 1 (0.0ms)

IterNodesCost
025565
142488
290488
3294488
41369488
52959488
64609488
75240488

prune385.0ms (4.1%)

Pruning

12 alts after pruning (12 fresh and 0 done)

PrunedKeptTotal
New1857192
Fresh7512
Picked101
Done000
Total19312205
Error
3.6b
Counts
205 → 12
Alt Table
StatusErrorProgram
19.0b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b b) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
27.1b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 (/.f64 1 (sqrt.f64 a)) (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 a)))) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
23.4b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (exp.f64 (log.f64 (-.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))) b)))))
22.1b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 (*.f64 (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) b)))
21.5b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (sqrt.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))) (/.f64 a (sqrt.f64 (+.f64 b (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))))) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
29.4b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (cbrt.f64 (pow.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))) 3)) b)))
18.8b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 2 (-.f64 (/.f64 b a) (/.f64 c b)))) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
18.8b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (*.f64 (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
19.2b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 (fabs.f64 (cbrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (cbrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))) b)))
19.0b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (/.f64 1 a))) (/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
19.7b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (exp.f64 (log.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
18.8b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 (/.f64 (sqrt.f64 (+.f64 b (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (*.f64 (cbrt.f64 a) (cbrt.f64 a))) (/.f64 (sqrt.f64 (+.f64 b (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 a)))) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
Compiler

Compiled 8290 to 5250 computations (36.7% saved)

localize24.0ms (0.3%)

Local error

Found 4 expressions with local error:

NewErrorProgram
0.3b
(*.f64 (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))
0.4b
(cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))
8.1b
(+.f64 b (*.f64 (*.f64 (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))))
13.6b
(sqrt.f64 (*.f64 a (*.f64 c -4)))
Compiler

Compiled 1530 to 1275 computations (16.7% saved)

series771.0ms (8.2%)

Counts
4 → 42
Calls

4 calls:

318.0ms
(*.f64 (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))
276.0ms
(cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))
155.0ms
(+.f64 b (*.f64 (*.f64 (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))))
22.0ms
(sqrt.f64 (*.f64 a (*.f64 c -4)))

rewrite17.0ms (0.2%)

Algorithm
rewrite-expression-head
Rules
22×add-sqr-sqrt_binary64
19×cbrt-prod_binary64
17×pow1_binary64
14×*-un-lft-identity_binary64
14×add-cube-cbrt_binary64
Counts
4 → 85
Calls

4 calls:

5.0ms
(*.f64 (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))
3.0ms
(sqrt.f64 (*.f64 a (*.f64 c -4)))
2.0ms
(+.f64 b (*.f64 (*.f64 (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))))
1.0ms
(cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))

simplify112.0ms (1.2%)

Algorithm
egg-herbie
Rules
764×fma-def_binary64
535×times-frac_binary64
306×cancel-sign-sub-inv_binary64
248×fma-neg_binary64
215×associate-*r*_binary64
Counts
127 → 193
Iterations

Useful iterations: 3 (0.0ms)

IterNodesCost
02483637
16443413
219893114
347273080
449543080

prune700.0ms (7.5%)

Pruning

11 alts after pruning (10 fresh and 1 done)

PrunedKeptTotal
New3261327
Fresh2911
Picked011
Done000
Total32811339
Error
3.6b
Counts
339 → 11
Alt Table
StatusErrorProgram
19.0b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b b) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
27.1b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 (/.f64 1 (sqrt.f64 a)) (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 a)))) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
23.4b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (exp.f64 (log.f64 (-.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))) b)))))
22.1b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 (*.f64 (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) b)))
21.5b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (sqrt.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))) (/.f64 a (sqrt.f64 (+.f64 b (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))))) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
29.4b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (cbrt.f64 (pow.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))) 3)) b)))
18.8b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 2 (-.f64 (/.f64 b a) (/.f64 c b)))) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
18.8b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (*.f64 (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
19.2b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 (fabs.f64 (cbrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (cbrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))) b)))
19.0b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (/.f64 1 a))) (/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
19.8b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (pow.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))) 2/3) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
Compiler

Compiled 19786 to 13733 computations (30.6% saved)

localize19.0ms (0.2%)

Local error

Found 1 expressions with local error:

NewErrorProgram
0.0b
(-.f64 (/.f64 b a) (/.f64 c b))
Compiler

Compiled 1530 to 1275 computations (16.7% saved)

series13.0ms (0.1%)

Counts
1 → 9
Calls

1 calls:

13.0ms
(-.f64 (/.f64 b a) (/.f64 c b))

rewrite50.0ms (0.5%)

Algorithm
rewrite-expression-head
Rules
258×times-frac_binary64
209×add-sqr-sqrt_binary64
203×*-un-lft-identity_binary64
197×add-cube-cbrt_binary64
169×prod-diff_binary64
Counts
1 → 215
Calls

1 calls:

17.0ms
(-.f64 (/.f64 b a) (/.f64 c b))

simplify315.0ms (3.4%)

Algorithm
egg-herbie
Rules
1366×fma-def_binary64
504×associate-*l/_binary64
361×associate-/l*_binary64
190×exp-prod_binary64
172×fma-neg_binary64
Counts
224 → 230
Iterations

Useful iterations: 3 (0.0ms)

IterNodesCost
03067669
16745128
213242983
337912904
448332904

prune448.0ms (4.8%)

Pruning

12 alts after pruning (11 fresh and 1 done)

PrunedKeptTotal
New2282230
Fresh099
Picked101
Done011
Total22912241
Error
3.6b
Counts
241 → 12
Alt Table
StatusErrorProgram
19.0b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b b) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
27.1b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 (/.f64 1 (sqrt.f64 a)) (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 a)))) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
19.2b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 (fabs.f64 (cbrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (cbrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))) b)))
30.5b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 2 (cbrt.f64 (pow.f64 (-.f64 (/.f64 b a) (/.f64 c b)) 3)))) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
23.4b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (exp.f64 (log.f64 (-.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))) b)))))
22.1b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 (*.f64 (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) b)))
21.5b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (sqrt.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))) (/.f64 a (sqrt.f64 (+.f64 b (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))))) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
18.8b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (*.f64 (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))))) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
29.4b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (cbrt.f64 (pow.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))) 3)) b)))
19.0b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (/.f64 1 a))) (/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
19.8b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (pow.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))) 2/3) (cbrt.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4))))))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
18.8b
(if (>=.f64 b 0) (*.f64 -1/2 (*.f64 2 (+.f64 (/.f64 b a) (neg.f64 (/.f64 c b))))) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
Compiler

Compiled 13013 to 5018 computations (61.4% saved)

regimes2.6s (28.1%)

Accuracy

Total 0.9b remaining (14.4%)

Threshold costs 0.9b (14.4%)

Counts
174 → 3
Compiler

Compiled 28766 to 19010 computations (33.9% saved)

bsearch37.0ms (0.4%)

Steps
ItersRangePoint
1
9.459365096577122e+142
1.4975735133232363e+146
3.8060811703142394e+144
5
-3.235337976882003e+87
-1.52982733063236e+87
-2.795577173814291e+87
Compiler

Compiled 4 to 3 computations (25% saved)

simplify6.0ms (0.1%)

Algorithm
egg-herbie
Rules
14×*-commutative_binary64
13×+-commutative_binary64
12×sub-neg_binary64
neg-sub0_binary64
neg-mul-1_binary64
Iterations

Useful iterations: 1 (0.0ms)

IterNodesCost
038100
15698
27998
39898
411198
511798
611898
711798

end1.0ms (0%)

Compiler

Compiled 68 to 44 computations (35.3% saved)

sample1.7s (18.2%)

Algorithm
intervals
Results
948.0ms8000×body128valid
77.0ms674×body128nan
45.0ms400×body128invalid
Compiler

Compiled 590 to 400 computations (32.2% saved)

Profiling

Loading profile data...