Details

Time bar (total: 9.7s)

analyze750.0ms (7.8%)

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)

sample1.9s (19.3%)

Results
1.6s8256×body128valid
191.0ms1074×body128invalid
Compiler

Compiled 125 to 80 computations (36% saved)

simplify157.0ms (1.6%)

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

prune4.0ms (0%)

Pruning

1 alts after pruning (1 fresh and 0 done)

PrunedKeptTotal
New112
Fresh101
Picked000
Done000
Total213
Error
18.1b
Counts
3 → 1
Alt Table
StatusErrorProgram
18.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 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
Compiler

Compiled 149 to 91 computations (38.9% saved)

localize27.0ms (0.3%)

Local error

Found 4 expressions with local error:

NewErrorProgram
0.1b
(/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b))
1.0b
(+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))
1.7b
(-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)
20.8b
(sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))
Compiler

Compiled 1536 to 1280 computations (16.7% saved)

series426.0ms (4.4%)

Counts
4 → 80
Calls

4 calls:

268.0ms
(/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b))
58.0ms
(-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)
51.0ms
(+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))
48.0ms
(sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))

rewrite22.0ms (0.2%)

Algorithm
rewrite-expression-head
Rules
36×add-sqr-sqrt_binary64
31×*-un-lft-identity_binary64
23×add-cube-cbrt_binary64
18×prod-diff_binary64
18×sqrt-prod_binary64
Counts
4 → 106
Calls

4 calls:

6.0ms
(/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b))
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)))

simplify146.0ms (1.5%)

Algorithm
egg-herbie
Rules
885×fma-neg_binary64
381×cancel-sign-sub-inv_binary64
181×distribute-rgt-neg-in_binary64
159×distribute-rgt-in_binary64
157×distribute-lft-neg-in_binary64
Counts
186 → 234
Iterations

Useful iterations: 3 (0.0ms)

IterNodesCost
02615384
16864516
225464321
345004319
448654319

prune431.0ms (4.5%)

Pruning

11 alts after pruning (11 fresh and 0 done)

PrunedKeptTotal
New26711278
Fresh000
Picked101
Done000
Total26811279
Error
4.5b
Counts
279 → 11
Alt Table
StatusErrorProgram
18.4b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (*.f64 (cbrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (cbrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))) (cbrt.f64 (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)))
21.0b
(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)))))
31.5b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (/.f64 (neg.f64 (+.f64 (*.f64 a (*.f64 c -4)) 0)) (-.f64 b (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)))
22.9b
(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.0b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (*.f64 (*.f64 (cbrt.f64 (/.f64 (*.f64 2 c) (-.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))) b))) (cbrt.f64 (/.f64 (*.f64 2 c) (-.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))) b)))) (cbrt.f64 (/.f64 (*.f64 2 c) (-.f64 (hypot.f64 b (sqrt.f64 (*.f64 a (*.f64 c -4)))) b)))))
18.3b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.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)))))) a)) (/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
29.8b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (*.f64 (/.f64 (*.f64 2 c) (+.f64 (*.f64 a (*.f64 c -4)) 0)) (+.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 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
21.8b
(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)))
18.2b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (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)))
18.1b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 1 (/.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b) (*.f64 c 2))))
Compiler

Compiled 13587 to 7973 computations (41.3% saved)

localize20.0ms (0.2%)

Local error

Found 4 expressions with local error:

NewErrorProgram
0.5b
(/.f64 1 (/.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b) (*.f64 c 2)))
1.0b
(+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))
1.7b
(-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)
20.8b
(sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))
Compiler

Compiled 1536 to 1280 computations (16.7% saved)

series271.0ms (2.8%)

Counts
1 → 20
Calls

1 calls:

270.0ms
(/.f64 1 (/.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b) (*.f64 c 2)))

rewrite16.0ms (0.2%)

Algorithm
rewrite-expression-head
Rules
54×times-frac_binary64
36×add-sqr-sqrt_binary64
28×*-un-lft-identity_binary64
22×add-exp-log_binary64
20×add-cube-cbrt_binary64
Counts
1 → 70
Calls

1 calls:

8.0ms
(/.f64 1 (/.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b) (*.f64 c 2)))

simplify97.0ms (1%)

Algorithm
egg-herbie
Rules
731×fma-def_binary64
663×associate-/l*_binary64
592×associate-/r*_binary64
349×associate-*l/_binary64
325×associate-/r/_binary64
Counts
90 → 110
Iterations

Useful iterations: 3 (0.0ms)

IterNodesCost
01553096
13852703
213512577
348092576
453382576

prune603.0ms (6.2%)

Pruning

12 alts after pruning (11 fresh and 1 done)

PrunedKeptTotal
New3152317
Fresh1910
Picked011
Done000
Total31612328
Error
4.5b
Counts
328 → 12
Alt Table
StatusErrorProgram
18.4b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (*.f64 (cbrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (cbrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))) (cbrt.f64 (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)))
21.0b
(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)))))
21.8b
(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)))
20.2b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 1 (/.f64 (-.f64 (neg.f64 b) b) (*.f64 c 2))))
18.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 1 (/.f64 (sqrt.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)) c)) (/.f64 1 (/.f64 (sqrt.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)) 2))))
18.1b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 1 (/.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b) (*.f64 c 2))))
31.5b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (/.f64 (neg.f64 (+.f64 (*.f64 a (*.f64 c -4)) 0)) (-.f64 b (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)))
22.9b
(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)))
18.3b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.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)))))) a)) (/.f64 (*.f64 c 2) (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
29.8b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (*.f64 (/.f64 (*.f64 2 c) (+.f64 (*.f64 a (*.f64 c -4)) 0)) (+.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 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
18.2b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (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 15989 to 9673 computations (39.5% saved)

localize21.0ms (0.2%)

Local error

Found 4 expressions with local error:

NewErrorProgram
0.5b
(*.f64 (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))))
1.7b
(-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)
12.0b
(+.f64 b (*.f64 (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))))
20.8b
(sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))
Compiler

Compiled 1536 to 1280 computations (16.7% saved)

series94.0ms (1%)

Counts
2 → 40
Calls

2 calls:

52.0ms
(+.f64 b (*.f64 (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))))
42.0ms
(*.f64 (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))))

rewrite22.0ms (0.2%)

Algorithm
rewrite-expression-head
Rules
74×sqrt-prod_binary64
61×add-sqr-sqrt_binary64
32×sqrt-pow1_binary64
25×pow1_binary64
16×rem-sqrt-square_binary64
Counts
2 → 101
Calls

2 calls:

8.0ms
(*.f64 (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))))
2.0ms
(+.f64 b (*.f64 (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))))

simplify122.0ms (1.3%)

Algorithm
egg-herbie
Rules
616×fma-neg_binary64
455×times-frac_binary64
367×fma-def_binary64
278×associate-/l*_binary64
255×associate-/r*_binary64
Counts
141 → 150
Iterations

Useful iterations: 3 (0.0ms)

IterNodesCost
01553099
13892882
211552664
341562624
449582624
550082624

prune824.0ms (8.5%)

Pruning

13 alts after pruning (11 fresh and 2 done)

PrunedKeptTotal
New3531354
Fresh01010
Picked011
Done011
Total35313366
Error
4.5b
Counts
366 → 13
Alt Table
StatusErrorProgram
18.4b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (*.f64 (cbrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (cbrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))) (cbrt.f64 (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)))
21.0b
(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)))))
24.8b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (pow.f64 (cbrt.f64 (sqrt.f64 (hypot.f64 (sqrt.f64 (*.f64 c (*.f64 a -4))) b))) 5) (cbrt.f64 (sqrt.f64 (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)))
29.8b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (*.f64 (/.f64 (*.f64 2 c) (+.f64 (*.f64 a (*.f64 c -4)) 0)) (+.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
21.8b
(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)))
20.2b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 1 (/.f64 (-.f64 (neg.f64 b) b) (*.f64 c 2))))
18.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 1 (/.f64 (sqrt.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)) c)) (/.f64 1 (/.f64 (sqrt.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)) 2))))
18.1b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 1 (/.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b) (*.f64 c 2))))
31.5b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (/.f64 (neg.f64 (+.f64 (*.f64 a (*.f64 c -4)) 0)) (-.f64 b (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)))
22.9b
(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)))
18.3b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.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)))))) 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 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
18.2b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (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 20329 to 12012 computations (40.9% saved)

localize22.0ms (0.2%)

Local error

Found 4 expressions with local error:

NewErrorProgram
0.4b
(/.f64 1 (/.f64 (sqrt.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)) c))
1.0b
(+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))))
1.7b
(-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)
20.8b
(sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))
Compiler

Compiled 1536 to 1280 computations (16.7% saved)

series829.0ms (8.6%)

Counts
1 → 4
Calls

1 calls:

828.0ms
(/.f64 1 (/.f64 (sqrt.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)) c))

rewrite40.0ms (0.4%)

Algorithm
rewrite-expression-head
Rules
201×times-frac_binary64
145×add-sqr-sqrt_binary64
121×*-un-lft-identity_binary64
97×add-cube-cbrt_binary64
84×sqrt-prod_binary64
Counts
1 → 150
Calls

1 calls:

12.0ms
(/.f64 1 (/.f64 (sqrt.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)) c))

simplify233.0ms (2.4%)

Algorithm
egg-herbie
Rules
934×fma-def_binary64
409×associate-/l*_binary64
356×cancel-sign-sub-inv_binary64
317×fma-neg_binary64
299×associate-*r*_binary64
Counts
154 → 167
Iterations

Useful iterations: 3 (0.0ms)

IterNodesCost
01686297
13305785
28545484
323145471
444835471
549755471

prune1.1s (11.4%)

Pruning

12 alts after pruning (10 fresh and 2 done)

PrunedKeptTotal
New4891490
Fresh1910
Picked101
Done022
Total49112503
Error
4.5b
Counts
503 → 12
Alt Table
StatusErrorProgram
18.2b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (fma.f64 (sqrt.f64 b) (sqrt.f64 b) (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (*.f64 (/.f64 1 (/.f64 (sqrt.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)) c)) (/.f64 1 (/.f64 (sqrt.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)) 2))))
21.0b
(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)))))
24.8b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (pow.f64 (cbrt.f64 (sqrt.f64 (hypot.f64 (sqrt.f64 (*.f64 c (*.f64 a -4))) b))) 5) (cbrt.f64 (sqrt.f64 (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)))
29.8b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (*.f64 (/.f64 (*.f64 2 c) (+.f64 (*.f64 a (*.f64 c -4)) 0)) (+.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b)))
21.8b
(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)))
20.2b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 1 (/.f64 (-.f64 (neg.f64 b) b) (*.f64 c 2))))
31.5b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (/.f64 (neg.f64 (+.f64 (*.f64 a (*.f64 c -4)) 0)) (-.f64 b (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)))
22.9b
(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)))
18.3b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.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)))))) 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 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 (*.f64 c 2) (-.f64 (*.f64 -1 b) b)))
18.2b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (*.f64 (sqrt.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) (sqrt.f64 (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)))
18.1b
(if (>=.f64 b 0) (*.f64 -1/2 (/.f64 (+.f64 b (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b)))) a)) (/.f64 1 (/.f64 (-.f64 (sqrt.f64 (fma.f64 a (*.f64 c -4) (*.f64 b b))) b) (*.f64 c 2))))
Compiler

Compiled 32801 to 19820 computations (39.6% saved)

regimes1.2s (12.4%)

Accuracy

Total 1.0b remaining (14.3%)

Threshold costs 1.0b (14.3%)

Counts
102 → 3
Compiler

Compiled 19302 to 12550 computations (35% saved)

bsearch36.0ms (0.4%)

Steps
ItersRangePoint
5
2.1665003818594198e+121
7.477859093117982e+121
5.545049901184398e+121
2
-6.23817191079231e+155
-5.018533530183444e+152
-2.9527445303397572e+153
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
neg-sub0_binary64
neg-mul-1_binary64
Iterations

Useful iterations: 1 (0.0ms)

IterNodesCost
048148
167146
290146
3109146
4122146
5128146
6129146
7128146

end288.0ms (3%)

Compiler

Compiled 707 to 454 computations (35.8% saved)

Profiling

Loading profile data...