Details

Time bar (total: 1.4s)

analyze137.0ms (9.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
12.5%87.4%0.1%5
25%68.7%6.3%6
34.3%53.1%12.6%7
46.8%39%14.1%8
53.1%28.9%18%9
60.1%20.7%19.2%10
63.6%15%21.4%11
67.3%10.6%22%12
69.2%7.7%23.2%13
71.1%5.4%23.5%14
Compiler

Compiled 9 to 6 computations (33.3% saved)

sample858.0ms (60.9%)

Results
834.0ms8256×body128valid
12.0ms164×body128invalid
Compiler

Compiled 18 to 12 computations (33.3% saved)

preprocess16.0ms (1.2%)

Algorithm
egg-herbie
Rules
+-commutative_binary64
fma-def_binary64
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
0914
11314
022
Stop Event
saturated
saturated
Calls
Call 1
Inputs
0
1
Outputs
0
1
Call 2
Inputs
(sqrt.f64 (+.f64 (*.f64 x x) y))
(sqrt.f64 (+.f64 (*.f64 y y) x))
Outputs
(sqrt.f64 (+.f64 (*.f64 x x) y))
(sqrt.f64 (fma.f64 x x y))
(sqrt.f64 (+.f64 (*.f64 y y) x))
(sqrt.f64 (fma.f64 y y x))
Compiler

Compiled 8 to 5 computations (37.5% saved)

simplify6.0ms (0.4%)

Algorithm
egg-herbie
Rules
+-commutative_binary64
fma-def_binary64
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
067
187
Stop Event
saturated
Counts
1 → 2
Calls
Call 1
Inputs
(sqrt.f64 (+.f64 (*.f64 x x) y))
Outputs
(sqrt.f64 (+.f64 (*.f64 x x) y))
(sqrt.f64 (fma.f64 x x y))

prune2.0ms (0.2%)

Pruning

1 alts after pruning (1 fresh and 0 done)

PrunedKeptTotal
New112
Fresh101
Picked000
Done000
Total213
Error
21.0b
Counts
3 → 1
Alt Table
StatusErrorProgram
21.0b
(sqrt.f64 (fma.f64 x x y))
Compiler

Compiled 27 to 15 computations (44.4% saved)

localize3.0ms (0.2%)

Local error

Found 1 expressions with local error:

NewErrorProgram
19.5b
(sqrt.f64 (fma.f64 x x y))

series8.0ms (0.6%)

Counts
1 → 16
Calls

6 calls:

TimeVariablePointExpression
3.0ms
x
@-inf
(sqrt.f64 (fma.f64 x x y))
1.0ms
y
@inf
(sqrt.f64 (fma.f64 x x y))
1.0ms
x
@0
(sqrt.f64 (fma.f64 x x y))
1.0ms
x
@inf
(sqrt.f64 (fma.f64 x x y))
1.0ms
y
@-inf
(sqrt.f64 (fma.f64 x x y))

rewrite32.0ms (2.3%)

Algorithm
batch-egg-rewrite
Rules
439×pow1_binary64
404×add-log-exp_binary64
404×log1p-expm1-u_binary64
404×expm1-log1p-u_binary64
393×add-exp-log_binary64
Iterations

Useful iterations: 1 (0.0ms)

IterNodesCost
057
1986
211216
Stop Event
node limit
Counts
1 → 21
Calls
Call 1
Inputs
(sqrt.f64 (fma.f64 x x y))
Outputs
((#(struct:change #<rule egg-rr> (2) ((x -.f64 (+.f64 1 (hypot.f64 x (sqrt.f64 y))) 1)))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (hypot.f64 x (sqrt.f64 y)) 1)))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 1 (hypot.f64 x (sqrt.f64 y)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (cbrt.f64 (hypot.f64 x (sqrt.f64 y))) (cbrt.f64 (fma.f64 x x y)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (cbrt.f64 (fma.f64 x x y)) (cbrt.f64 (hypot.f64 x (sqrt.f64 y))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (pow.f64 (fma.f64 x x y) 1/4) (pow.f64 (fma.f64 x x y) 1/4))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (sqrt.f64 (pow.f64 (cbrt.f64 (fma.f64 x x y)) 2)) (sqrt.f64 (cbrt.f64 (fma.f64 x x y))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (pow.f64 1 1/2) (hypot.f64 x (sqrt.f64 y)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (pow.f64 (pow.f64 (cbrt.f64 (fma.f64 x x y)) 2) 1/2) (pow.f64 (cbrt.f64 (fma.f64 x x y)) 1/2))))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (fma.f64 x x y) 1/2)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (hypot.f64 x (sqrt.f64 y)) 1)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (pow.f64 (fma.f64 x x y) 3/2) 1/3)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (cbrt.f64 (hypot.f64 x (sqrt.f64 y))) 3)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (pow.f64 (fma.f64 x x y) 1/4) 2)))) (#(struct:change #<rule egg-rr> (2) ((x fabs.f64 (hypot.f64 x (sqrt.f64 y)))))) (#(struct:change #<rule egg-rr> (2) ((x log.f64 (exp.f64 (hypot.f64 x (sqrt.f64 y))))))) (#(struct:change #<rule egg-rr> (2) ((x cbrt.f64 (pow.f64 (fma.f64 x x y) 3/2))))) (#(struct:change #<rule egg-rr> (2) ((x exp.f64 (log.f64 (hypot.f64 x (sqrt.f64 y))))))) (#(struct:change #<rule egg-rr> (2) ((x log1p.f64 (expm1.f64 (hypot.f64 x (sqrt.f64 y))))))) (#(struct:change #<rule egg-rr> (2) ((x expm1.f64 (log1p.f64 (hypot.f64 x (sqrt.f64 y))))))) (#(struct:change #<rule egg-rr> (2) ((x hypot.f64 x (sqrt.f64 y))))))

simplify24.0ms (1.7%)

Algorithm
egg-herbie
Rules
653×times-frac_binary64
423×cancel-sign-sub-inv_binary64
421×associate-/l*_binary64
345×fma-neg_binary64
339×associate-/r*_binary64
Iterations

Useful iterations: 3 (0.0ms)

IterNodesCost
054329
1156325
2453311
31743309
Stop Event
node limit
Counts
37 → 43
Calls
Call 1
Inputs
(sqrt.f64 y)
(+.f64 (sqrt.f64 y) (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))))
(-.f64 (+.f64 (sqrt.f64 y) (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2)))) (*.f64 1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))))
(-.f64 (+.f64 (sqrt.f64 y) (+.f64 (*.f64 1/16 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (pow.f64 x 6))) (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))))) (*.f64 1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(-.f64 (+.f64 (*.f64 1/2 (/.f64 y x)) x) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))
(-.f64 (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) x)) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))
(*.f64 -1 x)
(neg.f64 (+.f64 (*.f64 1/2 (/.f64 y x)) x))
(-.f64 (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) (+.f64 (*.f64 1/2 (/.f64 y x)) x))
(-.f64 (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) x)))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(-.f64 (+.f64 (*.f64 1/2 (/.f64 y x)) x) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))
(-.f64 (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) x)) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))
Outputs
(sqrt.f64 y)
(+.f64 (sqrt.f64 y) (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))))
(+.f64 (sqrt.f64 y) (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x))))
(fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y))
(-.f64 (+.f64 (sqrt.f64 y) (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2)))) (*.f64 1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))))
(+.f64 (sqrt.f64 y) (-.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x))) (*.f64 1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4)))))
(+.f64 (sqrt.f64 y) (fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (*.f64 (*.f64 -1/8 (pow.f64 x 4)) (sqrt.f64 (/.f64 1 (pow.f64 y 3))))))
(fma.f64 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4)) -1/8 (fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y)))
(-.f64 (+.f64 (sqrt.f64 y) (+.f64 (*.f64 1/16 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (pow.f64 x 6))) (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))))) (*.f64 1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))))
(+.f64 (sqrt.f64 y) (-.f64 (fma.f64 1/16 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (pow.f64 x 6)) (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)))) (*.f64 1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4)))))
(+.f64 (+.f64 (fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y)) (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (*.f64 (pow.f64 x 6) 1/16))) (*.f64 (*.f64 -1/8 (pow.f64 x 4)) (sqrt.f64 (/.f64 1 (pow.f64 y 3)))))
(fma.f64 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4)) -1/8 (fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (*.f64 1/16 (pow.f64 x 6)) (sqrt.f64 y))))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(fma.f64 1/2 (/.f64 y x) x)
(-.f64 (+.f64 (*.f64 1/2 (/.f64 y x)) x) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))
(+.f64 (fma.f64 1/2 (/.f64 y x) x) (*.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3))))
(+.f64 (fma.f64 1/2 (/.f64 y x) x) (/.f64 (*.f64 -1/8 (*.f64 y y)) (pow.f64 x 3)))
(fma.f64 (*.f64 (/.f64 y (pow.f64 x 3)) y) -1/8 (fma.f64 1/2 (/.f64 y x) x))
(-.f64 (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) x)) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))
(+.f64 (fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) x)) (*.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3))))
(+.f64 (fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) x)) (/.f64 (*.f64 -1/8 (*.f64 y y)) (pow.f64 x 3)))
(fma.f64 (*.f64 (/.f64 y (pow.f64 x 3)) y) -1/8 (fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) x)))
(*.f64 -1 x)
(neg.f64 x)
(neg.f64 (+.f64 (*.f64 1/2 (/.f64 y x)) x))
(neg.f64 (fma.f64 1/2 (/.f64 y x) x))
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
(-.f64 (/.f64 -1/2 (/.f64 x y)) x)
(-.f64 (*.f64 (/.f64 -1/2 x) y) x)
(-.f64 (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) (+.f64 (*.f64 1/2 (/.f64 y x)) x))
(-.f64 (*.f64 1/8 (/.f64 (*.f64 y y) (pow.f64 x 3))) (fma.f64 1/2 (/.f64 y x) x))
(-.f64 (/.f64 (*.f64 1/8 (*.f64 y y)) (pow.f64 x 3)) (fma.f64 1/2 (/.f64 y x) x))
(-.f64 (*.f64 (/.f64 y x) (-.f64 (/.f64 y (/.f64 x (/.f64 1/8 x))) 1/2)) x)
(-.f64 (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) x)))
(-.f64 (*.f64 1/8 (/.f64 (*.f64 y y) (pow.f64 x 3))) (fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) x)))
(-.f64 (/.f64 (*.f64 1/8 (*.f64 y y)) (pow.f64 x 3)) (fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) x)))
(-.f64 (*.f64 (/.f64 y x) (-.f64 (/.f64 y (/.f64 x (/.f64 1/8 x))) 1/2)) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) x))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(fma.f64 1/2 (/.f64 y x) x)
(-.f64 (+.f64 (*.f64 1/2 (/.f64 y x)) x) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))
(+.f64 (fma.f64 1/2 (/.f64 y x) x) (*.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3))))
(+.f64 (fma.f64 1/2 (/.f64 y x) x) (/.f64 (*.f64 -1/8 (*.f64 y y)) (pow.f64 x 3)))
(fma.f64 (*.f64 (/.f64 y (pow.f64 x 3)) y) -1/8 (fma.f64 1/2 (/.f64 y x) x))
(-.f64 (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) x)) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))
(+.f64 (fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) x)) (*.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3))))
(+.f64 (fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) x)) (/.f64 (*.f64 -1/8 (*.f64 y y)) (pow.f64 x 3)))
(fma.f64 (*.f64 (/.f64 y (pow.f64 x 3)) y) -1/8 (fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) x)))

prune30.0ms (2.2%)

Pruning

4 alts after pruning (3 fresh and 1 done)

PrunedKeptTotal
New40343
Fresh000
Picked011
Done000
Total40444
Error
0b
Counts
44 → 4
Alt Table
StatusErrorProgram
23.0b
(hypot.f64 x (sqrt.f64 y))
38.2b
(neg.f64 x)
21.0b
(sqrt.f64 (fma.f64 x x y))
43.0b
x
Compiler

Compiled 557 to 318 computations (42.9% saved)

localize6.0ms (0.5%)

prune2.0ms (0.1%)

Pruning

4 alts after pruning (2 fresh and 2 done)

PrunedKeptTotal
New000
Fresh022
Picked011
Done011
Total044
Error
0b
Counts
4 → 4
Alt Table
StatusErrorProgram
23.0b
(hypot.f64 x (sqrt.f64 y))
38.2b
(neg.f64 x)
21.0b
(sqrt.f64 (fma.f64 x x y))
43.0b
x
Compiler

Compiled 27 to 18 computations (33.3% saved)

localize1.0ms (0.1%)

prune2.0ms (0.1%)

Pruning

4 alts after pruning (1 fresh and 3 done)

PrunedKeptTotal
New000
Fresh011
Picked011
Done022
Total044
Error
0b
Counts
4 → 4
Alt Table
StatusErrorProgram
23.0b
(hypot.f64 x (sqrt.f64 y))
38.2b
(neg.f64 x)
21.0b
(sqrt.f64 (fma.f64 x x y))
43.0b
x
Compiler

Compiled 23 to 15 computations (34.8% saved)

localize0.0ms (0%)

prune1.0ms (0.1%)

Pruning

4 alts after pruning (0 fresh and 4 done)

PrunedKeptTotal
New000
Fresh000
Picked011
Done033
Total044
Error
0b
Counts
4 → 4
Alt Table
StatusErrorProgram
23.0b
(hypot.f64 x (sqrt.f64 y))
38.2b
(neg.f64 x)
21.0b
(sqrt.f64 (fma.f64 x x y))
43.0b
x
Compiler

Compiled 20 to 13 computations (35% saved)

regimes89.0ms (6.4%)

Accuracy

Total 0.1b remaining (96.7%)

Threshold costs 0.1b (96.7%)

Counts
10 → 3
Compiler

Compiled 75 to 38 computations (49.3% saved)

bsearch12.0ms (0.8%)

Algorithm
binary-search
Calls

2 calls:

8.0ms
x
4.0ms
x
Steps
ItersPointRange
1e+100
∈ [
2.9845752798042047e+93
,
1.6639317211339316e+103
]
-1e+154
∈ [
-1.4748852761592935e+159
,
-2.75761153735336e+153
]
Compiler

Compiled 3 to 2 computations (33.3% saved)

simplify5.0ms (0.4%)

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

Useful iterations: 0 (0.0ms)

IterNodesCost
01323
11723
21923
32023
Stop Event
saturated
Calls
Call 1
Inputs
(if (<=.f64 x -10000000000000000369475456880582265409809179829842688451922778552150543659347219597216513109705408327446511753687232667314337003349573404171046192448274432) (neg.f64 x) (if (<=.f64 x 10000000000000000159028911097599180468360808563945281389781327557747838772170381060813469985856815104) (sqrt.f64 (fma.f64 x x y)) x))
Outputs
(if (<=.f64 x -10000000000000000369475456880582265409809179829842688451922778552150543659347219597216513109705408327446511753687232667314337003349573404171046192448274432) (neg.f64 x) (if (<=.f64 x 10000000000000000159028911097599180468360808563945281389781327557747838772170381060813469985856815104) (sqrt.f64 (fma.f64 x x y)) x))

end171.0ms (12.1%)

Stop Event
done
Compiler

Compiled 56 to 35 computations (37.5% saved)

Profiling

Loading profile data...