Details

Time bar (total: 3.6s)

analyze69.0ms (1.9%)

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
Compiler

Compiled 9 to 6 computations (33.3% saved)

sample1.1s (31.1%)

Results
1.1s8256×body256valid
26.0ms320×body256invalid

preprocess20.0ms (0.6%)

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 10 to 7 computations (30% saved)

simplify8.0ms (0.2%)

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))

eval1.0ms (0%)

Compiler

Compiled 13 to 7 computations (46.2% saved)

prune2.0ms (0.1%)

Pruning

3 alts after pruning (3 fresh and 0 done)

PrunedKeptTotal
New022
Fresh011
Picked000
Done000
Total033
Error
20.0b
Counts
3 → 2
Alt Table
Click to see full alt table
StatusErrorProgram
20.0b
(sqrt.f64 (fma.f64 x x y))
20.0b
(sqrt.f64 (+.f64 (*.f64 x x) y))
Compiler

Compiled 15 to 9 computations (40% saved)

localize7.0ms (0.2%)

Local error

Found 2 expressions with local error:

NewErrorProgram
0.0b
(+.f64 (*.f64 x x) y)
18.6b
(sqrt.f64 (+.f64 (*.f64 x x) y))
Compiler

Compiled 18 to 5 computations (72.2% saved)

series13.0ms (0.4%)

Counts
2 → 40
Calls

12 calls:

TimeVariablePointExpression
7.0ms
y
@inf
(sqrt.f64 (+.f64 (*.f64 x x) y))
1.0ms
x
@0
(sqrt.f64 (+.f64 (*.f64 x x) y))
1.0ms
x
@inf
(sqrt.f64 (+.f64 (*.f64 x x) y))
1.0ms
y
@-inf
(sqrt.f64 (+.f64 (*.f64 x x) y))
1.0ms
y
@0
(sqrt.f64 (+.f64 (*.f64 x x) y))

rewrite64.0ms (1.8%)

Algorithm
batch-egg-rewrite
Rules
662×pow1_binary64
619×add-log-exp_binary64
619×log1p-expm1-u_binary64
619×expm1-log1p-u_binary64
602×add-exp-log_binary64
Iterations

Useful iterations: 1 (0.0ms)

IterNodesCost
0612
113011
2175111
Stop Event
node limit
Counts
2 → 62
Calls
Call 1
Inputs
(sqrt.f64 (+.f64 (*.f64 x x) y))
(+.f64 (*.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 /.f64 (sqrt.f64 (-.f64 (pow.f64 x 4) (*.f64 y y))) (sqrt.f64 (-.f64 (*.f64 x x) y)))))) (#(struct:change #<rule egg-rr> (2) ((x /.f64 (sqrt.f64 (+.f64 (pow.f64 x 6) (pow.f64 y 3))) (sqrt.f64 (+.f64 (pow.f64 x 4) (*.f64 y (-.f64 y (*.f64 x x))))))))) (#(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 hypot.f64 x (sqrt.f64 y))))) (#(struct:change #<rule egg-rr> (2) ((x hypot.f64 (sqrt.f64 y) x)))) (#(struct:change #<rule egg-rr> (2) ((x expm1.f64 (log1p.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 exp.f64 (log.f64 (hypot.f64 x (sqrt.f64 y))))))))
((#(struct:change #<rule egg-rr> (2) ((x -.f64 (+.f64 1 (fma.f64 x x y)) 1)))) (#(struct:change #<rule egg-rr> (2) ((x -.f64 (/.f64 (pow.f64 x 4) (-.f64 (*.f64 x x) y)) (/.f64 (*.f64 y y) (-.f64 (*.f64 x x) y)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (fma.f64 x x y) 1)))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (hypot.f64 x (sqrt.f64 y)) (hypot.f64 x (sqrt.f64 y)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 1 (fma.f64 x x y))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (cbrt.f64 (fma.f64 x x y)) (pow.f64 (cbrt.f64 (fma.f64 x x y)) 2))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (cbrt.f64 (fma.f64 x x y)) (*.f64 (cbrt.f64 (hypot.f64 x (sqrt.f64 y))) (hypot.f64 x (sqrt.f64 y))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (pow.f64 (cbrt.f64 (fma.f64 x x y)) 2) (cbrt.f64 (fma.f64 x x y)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (pow.f64 (fma.f64 x x y) 1/4) (*.f64 (pow.f64 (fma.f64 x x y) 1/4) (hypot.f64 x (sqrt.f64 y))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (+.f64 (pow.f64 x 6) (pow.f64 y 3)) (/.f64 1 (+.f64 (pow.f64 x 4) (*.f64 y (-.f64 y (*.f64 x x))))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (-.f64 (pow.f64 x 4) (*.f64 y y)) (/.f64 1 (-.f64 (*.f64 x x) y)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (*.f64 (hypot.f64 x (sqrt.f64 y)) (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 (*.f64 (hypot.f64 x (sqrt.f64 y)) (cbrt.f64 (fma.f64 x x y))) (cbrt.f64 (hypot.f64 x (sqrt.f64 y))))))) (#(struct:change #<rule egg-rr> (2) ((x /.f64 1 (/.f64 (+.f64 (pow.f64 x 4) (*.f64 y (-.f64 y (*.f64 x x)))) (+.f64 (pow.f64 x 6) (pow.f64 y 3))))))) (#(struct:change #<rule egg-rr> (2) ((x /.f64 1 (/.f64 (-.f64 (*.f64 x x) y) (-.f64 (pow.f64 x 4) (*.f64 y y))))))) (#(struct:change #<rule egg-rr> (2) ((x /.f64 (+.f64 (pow.f64 x 6) (pow.f64 y 3)) (+.f64 (pow.f64 x 4) (*.f64 y (-.f64 y (*.f64 x x)))))))) (#(struct:change #<rule egg-rr> (2) ((x /.f64 (+.f64 (pow.f64 x 6) (pow.f64 y 3)) (+.f64 (*.f64 y y) (-.f64 (pow.f64 x 4) (*.f64 (*.f64 x x) y))))))) (#(struct:change #<rule egg-rr> (2) ((x /.f64 (-.f64 (pow.f64 x 4) (*.f64 y y)) (-.f64 (*.f64 x x) y))))) (#(struct:change #<rule egg-rr> (2) ((x /.f64 (neg.f64 (+.f64 (pow.f64 x 6) (pow.f64 y 3))) (neg.f64 (+.f64 (pow.f64 x 4) (*.f64 y (-.f64 y (*.f64 x x))))))))) (#(struct:change #<rule egg-rr> (2) ((x /.f64 (neg.f64 (-.f64 (pow.f64 x 4) (*.f64 y y))) (neg.f64 (-.f64 (*.f64 x x) y)))))) (#(struct:change #<rule egg-rr> (2) ((x /.f64 (-.f64 (*.f64 y y) (pow.f64 x 4)) (-.f64 y (*.f64 x x)))))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (fma.f64 x x y) 1)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (hypot.f64 x (sqrt.f64 y)) 2)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (pow.f64 (fma.f64 x x y) 2) 1/2)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (pow.f64 (fma.f64 x x y) 3) 1/3)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (cbrt.f64 (fma.f64 x x y)) 3)))) (#(struct:change #<rule egg-rr> (2) ((x sqrt.f64 (pow.f64 (fma.f64 x x y) 2))))) (#(struct:change #<rule egg-rr> (2) ((x log.f64 (exp.f64 (fma.f64 x x y)))))) (#(struct:change #<rule egg-rr> (2) ((x cbrt.f64 (pow.f64 (fma.f64 x x y) 3))))) (#(struct:change #<rule egg-rr> (2) ((x expm1.f64 (log1p.f64 (fma.f64 x x y)))))) (#(struct:change #<rule egg-rr> (2) ((x log1p.f64 (expm1.f64 (fma.f64 x x y)))))) (#(struct:change #<rule egg-rr> (2) ((x exp.f64 (log.f64 (fma.f64 x x y)))))) (#(struct:change #<rule egg-rr> (2) ((x fma.f64 x x y)))) (#(struct:change #<rule egg-rr> (2) ((x fma.f64 1 (*.f64 x x) y)))) (#(struct:change #<rule egg-rr> (2) ((x fma.f64 1 y (*.f64 x x))))) (#(struct:change #<rule egg-rr> (2) ((x fma.f64 (cbrt.f64 (pow.f64 x 4)) (pow.f64 (cbrt.f64 x) 2) y)))) (#(struct:change #<rule egg-rr> (2) ((x fma.f64 (pow.f64 (cbrt.f64 y) 2) (cbrt.f64 y) (*.f64 x x))))) (#(struct:change #<rule egg-rr> (2) ((x fma.f64 (sqrt.f64 y) (sqrt.f64 y) (*.f64 x x))))))

simplify56.0ms (1.6%)

Algorithm
egg-herbie
Rules
759×times-frac_binary64
599×associate-/l*_binary64
502×unswap-sqr_binary64
498×fma-neg_binary64
427×cancel-sign-sub-inv_binary64
Iterations

Useful iterations: 3 (0.0ms)

IterNodesCost
060454
1161416
2444416
31586414
45763414
Stop Event
node limit
Counts
102 → 88
Calls
Call 1
Inputs
(sqrt.f64 y)
(+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y))
(+.f64 (*.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))) (+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y)))
(+.f64 (*.f64 1/16 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (pow.f64 x 6))) (+.f64 (*.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))) (+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y))))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x))
(+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x)))
(*.f64 -1 x)
(+.f64 (*.f64 -1/2 (/.f64 y x)) (*.f64 -1 x))
(+.f64 (*.f64 -1/2 (/.f64 y x)) (+.f64 (*.f64 -1 x) (*.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))) (+.f64 (*.f64 -1 x) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x))
(+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x)))
y
(+.f64 y (pow.f64 x 2))
(+.f64 y (pow.f64 x 2))
(+.f64 y (pow.f64 x 2))
(pow.f64 x 2)
(+.f64 y (pow.f64 x 2))
(+.f64 y (pow.f64 x 2))
(+.f64 y (pow.f64 x 2))
(pow.f64 x 2)
(+.f64 y (pow.f64 x 2))
(+.f64 y (pow.f64 x 2))
(+.f64 y (pow.f64 x 2))
(pow.f64 x 2)
(+.f64 y (pow.f64 x 2))
(+.f64 y (pow.f64 x 2))
(+.f64 y (pow.f64 x 2))
y
(+.f64 y (pow.f64 x 2))
(+.f64 y (pow.f64 x 2))
(+.f64 y (pow.f64 x 2))
y
(+.f64 y (pow.f64 x 2))
(+.f64 y (pow.f64 x 2))
(+.f64 y (pow.f64 x 2))
Outputs
(sqrt.f64 y)
(+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y))
(fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y))
(+.f64 (*.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))) (+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y)))
(fma.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4)) (fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y)))
(fma.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 1/2 (*.f64 x x)) (fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (*.f64 -1/8 (pow.f64 x 4)) (sqrt.f64 y)))
(+.f64 (*.f64 1/16 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (pow.f64 x 6))) (+.f64 (*.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))) (+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y))))
(fma.f64 1/16 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (pow.f64 x 6)) (fma.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4)) (fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y))))
(fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (*.f64 1/16 (pow.f64 x 6)) (fma.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 1/2 (*.f64 x x)) (fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (*.f64 -1/8 (pow.f64 x 4)) (sqrt.f64 y))))
(fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (*.f64 -1/8 (pow.f64 x 4)) (fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (*.f64 1/16 (pow.f64 x 6)) (fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y))))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(fma.f64 1/2 (/.f64 y x) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3)) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 y (/.f64 (pow.f64 x 3) y)) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) x))
(fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) (fma.f64 1/2 (/.f64 y x) x))
(+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x)))
(fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3)) x)))
(fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 y (/.f64 (pow.f64 x 3) y)) x)))
(fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1/8 (/.f64 y (/.f64 (pow.f64 x 3) y)) x)))
(fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) x)))
(fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) (fma.f64 1/2 (/.f64 y x) x)))
(*.f64 -1 x)
(neg.f64 x)
(+.f64 (*.f64 -1/2 (/.f64 y x)) (*.f64 -1 x))
(fma.f64 -1/2 (/.f64 y x) (neg.f64 x))
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
(+.f64 (*.f64 -1/2 (/.f64 y x)) (+.f64 (*.f64 -1 x) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3)))))
(fma.f64 -1/2 (/.f64 y x) (fma.f64 -1 x (*.f64 (/.f64 (*.f64 y y) (pow.f64 x 3)) 1/8)))
(fma.f64 (/.f64 y x) -1/2 (-.f64 (*.f64 (/.f64 y (/.f64 (pow.f64 x 3) y)) 1/8) x))
(-.f64 (fma.f64 (/.f64 y x) -1/2 (*.f64 (/.f64 y (/.f64 (pow.f64 x 3) y)) 1/8)) x)
(-.f64 (*.f64 (/.f64 y x) (+.f64 -1/2 (*.f64 (/.f64 y x) (/.f64 1/8 x)))) x)
(fma.f64 (/.f64 y x) (+.f64 -1/2 (*.f64 (/.f64 y x) (/.f64 1/8 x))) (neg.f64 x))
(+.f64 (*.f64 -1/2 (/.f64 y x)) (+.f64 (*.f64 -1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) (+.f64 (*.f64 -1 x) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))))
(fma.f64 -1/2 (/.f64 y x) (fma.f64 -1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1 x (*.f64 (/.f64 (*.f64 y y) (pow.f64 x 3)) 1/8))))
(fma.f64 (/.f64 y x) -1/2 (fma.f64 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) -1/16 (-.f64 (*.f64 (/.f64 y (/.f64 (pow.f64 x 3) y)) 1/8) x)))
(fma.f64 (/.f64 y x) -1/2 (fma.f64 (/.f64 y (/.f64 (pow.f64 x 3) y)) 1/8 (-.f64 (/.f64 (*.f64 (pow.f64 y 3) -1/16) (pow.f64 x 5)) x)))
(fma.f64 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) -1/16 (-.f64 (*.f64 (/.f64 y x) (+.f64 -1/2 (*.f64 (/.f64 y x) (/.f64 1/8 x)))) x))
(fma.f64 (pow.f64 y 3) (/.f64 -1/16 (pow.f64 x 5)) (fma.f64 (/.f64 y x) (+.f64 -1/2 (*.f64 (/.f64 y x) (/.f64 1/8 x))) (neg.f64 x)))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(fma.f64 1/2 (/.f64 y x) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3)) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 y (/.f64 (pow.f64 x 3) y)) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) x))
(fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) (fma.f64 1/2 (/.f64 y x) x))
(+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x)))
(fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3)) x)))
(fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 y (/.f64 (pow.f64 x 3) y)) x)))
(fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1/8 (/.f64 y (/.f64 (pow.f64 x 3) y)) x)))
(fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) x)))
(fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) (fma.f64 1/2 (/.f64 y x) x)))
y
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(pow.f64 x 2)
(*.f64 x x)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(pow.f64 x 2)
(*.f64 x x)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(pow.f64 x 2)
(*.f64 x x)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
y
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
y
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)
(+.f64 y (pow.f64 x 2))
(+.f64 y (*.f64 x x))
(fma.f64 x x y)

localize41.0ms (1.2%)

Local error

Found 1 expressions with local error:

NewErrorProgram
18.6b
(sqrt.f64 (fma.f64 x x y))
Compiler

Compiled 13 to 4 computations (69.2% saved)

series3.0ms (0.1%)

Counts
1 → 16
Calls

6 calls:

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

rewrite59.0ms (1.6%)

Algorithm
batch-egg-rewrite
Rules
439×pow1_binary64
421×add-sqr-sqrt_binary64
414×*-un-lft-identity_binary64
404×add-log-exp_binary64
404×log1p-expm1-u_binary64
Iterations

Useful iterations: 1 (0.0ms)

IterNodesCost
057
1986
211216
Stop Event
node limit
Counts
1 → 25
Calls
Call 1
Inputs
(sqrt.f64 (fma.f64 x x y))
Outputs
((#(struct:change #<rule egg-rr> (2) ((x +.f64 0 (hypot.f64 x (sqrt.f64 y)))))) (#(struct:change #<rule egg-rr> (2) ((x +.f64 (log.f64 (sqrt.f64 (exp.f64 (hypot.f64 x (sqrt.f64 y))))) (log.f64 (sqrt.f64 (exp.f64 (hypot.f64 x (sqrt.f64 y))))))))) (#(struct:change #<rule egg-rr> (2) ((x +.f64 (log.f64 (*.f64 (cbrt.f64 (exp.f64 (hypot.f64 x (sqrt.f64 y)))) (cbrt.f64 (exp.f64 (hypot.f64 x (sqrt.f64 y)))))) (log.f64 (cbrt.f64 (exp.f64 (hypot.f64 x (sqrt.f64 y))))))))) (#(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 hypot.f64 x (sqrt.f64 y))))) (#(struct:change #<rule egg-rr> (2) ((x hypot.f64 (sqrt.f64 y) x)))) (#(struct:change #<rule egg-rr> (2) ((x expm1.f64 (log1p.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 exp.f64 (log.f64 (hypot.f64 x (sqrt.f64 y))))))))

simplify105.0ms (2.9%)

Algorithm
egg-herbie
Rules
759×times-frac_binary64
598×associate-/l*_binary64
502×unswap-sqr_binary64
498×fma-neg_binary64
427×cancel-sign-sub-inv_binary64
Iterations

Useful iterations: 3 (0.0ms)

IterNodesCost
059313
1159311
2440311
31580309
45737309
Stop Event
node limit
Counts
41 → 47
Calls
Call 1
Inputs
(sqrt.f64 y)
(+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y))
(+.f64 (*.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))) (+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y)))
(+.f64 (*.f64 1/16 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (pow.f64 x 6))) (+.f64 (*.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))) (+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y))))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x))
(+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x)))
(*.f64 -1 x)
(+.f64 (*.f64 -1/2 (/.f64 y x)) (*.f64 -1 x))
(+.f64 (*.f64 -1/2 (/.f64 y x)) (+.f64 (*.f64 -1 x) (*.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))) (+.f64 (*.f64 -1 x) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x))
(+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x)))
Outputs
(sqrt.f64 y)
(+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y))
(fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y))
(+.f64 (*.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))) (+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y)))
(fma.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4)) (fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y)))
(fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (*.f64 -1/8 (pow.f64 x 4)) (sqrt.f64 y)))
(fma.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 1/2 (*.f64 x x)) (fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (*.f64 -1/8 (pow.f64 x 4)) (sqrt.f64 y)))
(+.f64 (*.f64 1/16 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (pow.f64 x 6))) (+.f64 (*.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))) (+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y))))
(fma.f64 1/16 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (pow.f64 x 6)) (fma.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4)) (fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y))))
(fma.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4)) (fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (*.f64 1/16 (pow.f64 x 6)) (fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y))))
(fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (*.f64 -1/8 (pow.f64 x 4)) (fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (*.f64 1/16 (pow.f64 x 6)) (fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y))))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(fma.f64 1/2 (/.f64 y x) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3)) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (*.f64 (/.f64 y (pow.f64 x 3)) y) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) x))
(+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x)))
(fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3)) x)))
(fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1/8 (*.f64 (/.f64 y (pow.f64 x 3)) y) x)))
(fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) x)))
(*.f64 -1 x)
(neg.f64 x)
(+.f64 (*.f64 -1/2 (/.f64 y x)) (*.f64 -1 x))
(fma.f64 -1/2 (/.f64 y x) (neg.f64 x))
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
(fma.f64 (/.f64 y x) -1/2 (neg.f64 x))
(+.f64 (*.f64 -1/2 (/.f64 y x)) (+.f64 (*.f64 -1 x) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3)))))
(fma.f64 -1/2 (/.f64 y x) (fma.f64 -1 x (*.f64 (/.f64 (*.f64 y y) (pow.f64 x 3)) 1/8)))
(fma.f64 (/.f64 y x) -1/2 (-.f64 (*.f64 (/.f64 (*.f64 y y) (pow.f64 x 3)) 1/8) x))
(-.f64 (fma.f64 (/.f64 y x) -1/2 (*.f64 (*.f64 (/.f64 y (pow.f64 x 3)) y) 1/8)) x)
(-.f64 (*.f64 (/.f64 y x) (+.f64 -1/2 (/.f64 1/8 (/.f64 x (/.f64 y x))))) x)
(fma.f64 (/.f64 y x) (+.f64 -1/2 (/.f64 1/8 (/.f64 x (/.f64 y x)))) (neg.f64 x))
(+.f64 (*.f64 -1/2 (/.f64 y x)) (+.f64 (*.f64 -1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) (+.f64 (*.f64 -1 x) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))))
(fma.f64 -1/2 (/.f64 y x) (fma.f64 -1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1 x (*.f64 (/.f64 (*.f64 y y) (pow.f64 x 3)) 1/8))))
(fma.f64 (/.f64 y x) -1/2 (fma.f64 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) -1/16 (-.f64 (*.f64 (/.f64 (*.f64 y y) (pow.f64 x 3)) 1/8) x)))
(fma.f64 (/.f64 y x) -1/2 (-.f64 (fma.f64 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) -1/16 (*.f64 (*.f64 (/.f64 y (pow.f64 x 3)) y) 1/8)) x))
(fma.f64 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) -1/16 (-.f64 (*.f64 (/.f64 y x) (+.f64 -1/2 (/.f64 1/8 (/.f64 x (/.f64 y x))))) x))
(fma.f64 (pow.f64 y 3) (/.f64 -1/16 (pow.f64 x 5)) (fma.f64 (/.f64 y x) (+.f64 -1/2 (/.f64 1/8 (/.f64 x (/.f64 y x)))) (neg.f64 x)))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(fma.f64 1/2 (/.f64 y x) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3)) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (*.f64 (/.f64 y (pow.f64 x 3)) y) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) x))
(+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x)))
(fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3)) x)))
(fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1/8 (*.f64 (/.f64 y (pow.f64 x 3)) y) x)))
(fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) x)))

eval47.0ms (1.3%)

Compiler

Compiled 1630 to 805 computations (50.6% saved)

prune54.0ms (1.5%)

Pruning

7 alts after pruning (6 fresh and 1 done)

PrunedKeptTotal
New1296135
Fresh000
Picked101
Done011
Total1307137
Error
0b
Counts
137 → 7
Alt Table
Click to see full alt table
StatusErrorProgram
40.0b
(neg.f64 x)
43.7b
x
20.0b
(sqrt.f64 (+.f64 (*.f64 x x) y))
40.2b
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
20.5b
(hypot.f64 (sqrt.f64 y) x)
43.8b
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
41.7b
(sqrt.f64 y)
Compiler

Compiled 75 to 51 computations (32% saved)

localize103.0ms (2.9%)

Local error

Found 1 expressions with local error:

NewErrorProgram
0.0b
(hypot.f64 (sqrt.f64 y) x)
Compiler

Compiled 10 to 4 computations (60% saved)

series7.0ms (0.2%)

Counts
1 → 13
Calls

6 calls:

TimeVariablePointExpression
2.0ms
y
@-inf
(hypot.f64 (sqrt.f64 y) x)
2.0ms
y
@inf
(hypot.f64 (sqrt.f64 y) x)
1.0ms
y
@0
(hypot.f64 (sqrt.f64 y) x)
1.0ms
x
@0
(hypot.f64 (sqrt.f64 y) x)
0.0ms
x
@inf
(hypot.f64 (sqrt.f64 y) x)

rewrite78.0ms (2.2%)

Algorithm
batch-egg-rewrite
Rules
432×pow1_binary64
413×add-sqr-sqrt_binary64
406×*-un-lft-identity_binary64
397×add-log-exp_binary64
397×log1p-expm1-u_binary64
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
056
1996
211296
Stop Event
node limit
Counts
1 → 20
Calls
Call 1
Inputs
(hypot.f64 (sqrt.f64 y) x)
Outputs
((#(struct:change #<rule egg-rr> (2) ((x +.f64 0 (hypot.f64 (sqrt.f64 y) x))))) (#(struct:change #<rule egg-rr> (2) ((x +.f64 (log.f64 (*.f64 (cbrt.f64 (exp.f64 (hypot.f64 (sqrt.f64 y) x))) (cbrt.f64 (exp.f64 (hypot.f64 (sqrt.f64 y) x))))) (log.f64 (cbrt.f64 (exp.f64 (hypot.f64 (sqrt.f64 y) x)))))))) (#(struct:change #<rule egg-rr> (2) ((x +.f64 (log.f64 (sqrt.f64 (exp.f64 (hypot.f64 (sqrt.f64 y) x)))) (log.f64 (sqrt.f64 (exp.f64 (hypot.f64 (sqrt.f64 y) x)))))))) (#(struct:change #<rule egg-rr> (2) ((x -.f64 (+.f64 1 (hypot.f64 (sqrt.f64 y) x)) 1)))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (hypot.f64 (sqrt.f64 y) x) 1)))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 1 (hypot.f64 (sqrt.f64 y) x))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (cbrt.f64 (hypot.f64 (sqrt.f64 y) x)) (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 (sqrt.f64 y) x)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (sqrt.f64 (hypot.f64 (sqrt.f64 y) x)) (sqrt.f64 (hypot.f64 (sqrt.f64 y) x)))))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (hypot.f64 (sqrt.f64 y) x) 1)))) (#(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 (pow.f64 (hypot.f64 (sqrt.f64 y) x) 3) 1/3)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (cbrt.f64 (hypot.f64 (sqrt.f64 y) x)) 3)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (sqrt.f64 (hypot.f64 (sqrt.f64 y) x)) 2)))) (#(struct:change #<rule egg-rr> (2) ((x sqrt.f64 (fma.f64 x x y))))) (#(struct:change #<rule egg-rr> (2) ((x log.f64 (exp.f64 (hypot.f64 (sqrt.f64 y) x)))))) (#(struct:change #<rule egg-rr> (2) ((x cbrt.f64 (pow.f64 (hypot.f64 (sqrt.f64 y) x) 3))))) (#(struct:change #<rule egg-rr> (2) ((x expm1.f64 (log1p.f64 (hypot.f64 (sqrt.f64 y) x)))))) (#(struct:change #<rule egg-rr> (2) ((x log1p.f64 (expm1.f64 (hypot.f64 (sqrt.f64 y) x)))))) (#(struct:change #<rule egg-rr> (2) ((x exp.f64 (log.f64 (hypot.f64 (sqrt.f64 y) x)))))))

simplify129.0ms (3.6%)

Algorithm
egg-herbie
Rules
759×times-frac_binary64
598×associate-/l*_binary64
502×unswap-sqr_binary64
498×fma-neg_binary64
427×cancel-sign-sub-inv_binary64
Iterations

Useful iterations: 3 (0.0ms)

IterNodesCost
059250
1159248
2440248
31580246
45736246
Stop Event
node limit
Counts
33 → 42
Calls
Call 1
Inputs
x
(sqrt.f64 y)
(+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y))
(+.f64 (*.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))) (+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y)))
(+.f64 (*.f64 1/16 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (pow.f64 x 6))) (+.f64 (*.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))) (+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y))))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x))
(+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x)))
(*.f64 -1 x)
(+.f64 (*.f64 -1/2 (/.f64 y x)) (*.f64 -1 x))
(+.f64 (*.f64 -1/2 (/.f64 y x)) (+.f64 (*.f64 -1 x) (*.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))) (+.f64 (*.f64 -1 x) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))))
Outputs
x
(sqrt.f64 y)
(+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y))
(fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y))
(fma.f64 1/2 (*.f64 x (*.f64 x (sqrt.f64 (/.f64 1 y)))) (sqrt.f64 y))
(fma.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x (*.f64 x 1/2)) (sqrt.f64 y))
(+.f64 (*.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))) (+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y)))
(fma.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4)) (fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y)))
(fma.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4)) (fma.f64 1/2 (*.f64 x (*.f64 x (sqrt.f64 (/.f64 1 y)))) (sqrt.f64 y)))
(fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (*.f64 -1/8 (pow.f64 x 4)) (fma.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x (*.f64 x 1/2)) (sqrt.f64 y)))
(fma.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x (*.f64 x 1/2)) (fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (*.f64 -1/8 (pow.f64 x 4)) (sqrt.f64 y)))
(+.f64 (*.f64 1/16 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (pow.f64 x 6))) (+.f64 (*.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4))) (+.f64 (*.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (pow.f64 x 2))) (sqrt.f64 y))))
(fma.f64 1/16 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (pow.f64 x 6)) (fma.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4)) (fma.f64 1/2 (*.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x x)) (sqrt.f64 y))))
(fma.f64 1/16 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (pow.f64 x 6)) (fma.f64 -1/8 (*.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (pow.f64 x 4)) (fma.f64 1/2 (*.f64 x (*.f64 x (sqrt.f64 (/.f64 1 y)))) (sqrt.f64 y))))
(fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (*.f64 1/16 (pow.f64 x 6)) (fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (*.f64 -1/8 (pow.f64 x 4)) (fma.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x (*.f64 x 1/2)) (sqrt.f64 y))))
(fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 3))) (*.f64 -1/8 (pow.f64 x 4)) (fma.f64 (sqrt.f64 (/.f64 1 (pow.f64 y 5))) (*.f64 1/16 (pow.f64 x 6)) (fma.f64 (sqrt.f64 (/.f64 1 y)) (*.f64 x (*.f64 x 1/2)) (sqrt.f64 y))))
x
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(fma.f64 1/2 (/.f64 y x) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3)) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (*.f64 (/.f64 y (pow.f64 x 3)) y) x))
(fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) x))
(+.f64 (*.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) (+.f64 (*.f64 1/2 (/.f64 y x)) (+.f64 (*.f64 -1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))) x)))
(fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 1/2 (/.f64 y x) (fma.f64 -1/8 (/.f64 (*.f64 y y) (pow.f64 x 3)) x)))
(fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1/8 (*.f64 (/.f64 y (pow.f64 x 3)) y) x)))
(fma.f64 1/2 (/.f64 y x) (fma.f64 1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1/8 (*.f64 y (/.f64 y (pow.f64 x 3))) x)))
(*.f64 -1 x)
(neg.f64 x)
(+.f64 (*.f64 -1/2 (/.f64 y x)) (*.f64 -1 x))
(fma.f64 -1/2 (/.f64 y x) (neg.f64 x))
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
(-.f64 (*.f64 y (/.f64 -1/2 x)) x)
(fma.f64 (/.f64 y x) -1/2 (neg.f64 x))
(+.f64 (*.f64 -1/2 (/.f64 y x)) (+.f64 (*.f64 -1 x) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3)))))
(fma.f64 -1/2 (/.f64 y x) (fma.f64 -1 x (*.f64 (/.f64 (*.f64 y y) (pow.f64 x 3)) 1/8)))
(fma.f64 (/.f64 y x) -1/2 (-.f64 (*.f64 (/.f64 (*.f64 y y) (pow.f64 x 3)) 1/8) x))
(-.f64 (fma.f64 (/.f64 y x) -1/2 (*.f64 (*.f64 (/.f64 y (pow.f64 x 3)) y) 1/8)) x)
(-.f64 (*.f64 (/.f64 y x) (+.f64 -1/2 (*.f64 (/.f64 y x) (/.f64 1/8 x)))) x)
(fma.f64 (/.f64 y x) (+.f64 -1/2 (/.f64 1/8 (/.f64 x (/.f64 y x)))) (neg.f64 x))
(+.f64 (*.f64 -1/2 (/.f64 y x)) (+.f64 (*.f64 -1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5))) (+.f64 (*.f64 -1 x) (*.f64 1/8 (/.f64 (pow.f64 y 2) (pow.f64 x 3))))))
(fma.f64 -1/2 (/.f64 y x) (fma.f64 -1/16 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) (fma.f64 -1 x (*.f64 (/.f64 (*.f64 y y) (pow.f64 x 3)) 1/8))))
(fma.f64 (/.f64 y x) -1/2 (fma.f64 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) -1/16 (-.f64 (*.f64 (/.f64 (*.f64 y y) (pow.f64 x 3)) 1/8) x)))
(fma.f64 (/.f64 y x) -1/2 (-.f64 (fma.f64 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) -1/16 (*.f64 (*.f64 (/.f64 y (pow.f64 x 3)) y) 1/8)) x))
(fma.f64 (/.f64 (pow.f64 y 3) (pow.f64 x 5)) -1/16 (-.f64 (*.f64 (/.f64 y x) (+.f64 -1/2 (*.f64 (/.f64 y x) (/.f64 1/8 x)))) x))
(fma.f64 (pow.f64 y 3) (/.f64 -1/16 (pow.f64 x 5)) (fma.f64 (/.f64 y x) (+.f64 -1/2 (/.f64 1/8 (/.f64 x (/.f64 y x)))) (neg.f64 x)))

localize3.0ms (0.1%)

Compiler

Compiled 5 to 3 computations (40% saved)

localize9.0ms (0.2%)

Compiler

Compiled 20 to 9 computations (55% saved)

localize8.0ms (0.2%)

Compiler

Compiled 20 to 9 computations (55% saved)

localize4.0ms (0.1%)

Compiler

Compiled 5 to 3 computations (40% saved)

eval10.0ms (0.3%)

Compiler

Compiled 525 to 285 computations (45.7% saved)

prune33.0ms (0.9%)

Pruning

7 alts after pruning (1 fresh and 6 done)

PrunedKeptTotal
New41142
Fresh101
Picked011
Done055
Total42749
Error
0b
Counts
49 → 7
Alt Table
Click to see full alt table
StatusErrorProgram
40.0b
(neg.f64 x)
20.0b
(sqrt.f64 (+.f64 (*.f64 x x) y))
40.2b
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
43.7b
x
20.5b
(hypot.f64 (sqrt.f64 y) x)
43.8b
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
41.7b
(sqrt.f64 y)
Compiler

Compiled 43 to 29 computations (32.6% saved)

eval0.0ms (0%)

Compiler

Compiled 2 to 2 computations (0% saved)

prune3.0ms (0.1%)

Pruning

7 alts after pruning (1 fresh and 6 done)

PrunedKeptTotal
New000
Fresh011
Picked000
Done066
Total077
Error
0b
Counts
7 → 7
Alt Table
Click to see full alt table
StatusErrorProgram
40.0b
(neg.f64 x)
20.0b
(sqrt.f64 (+.f64 (*.f64 x x) y))
40.2b
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
43.7b
x
20.5b
(hypot.f64 (sqrt.f64 y) x)
43.8b
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
41.7b
(sqrt.f64 y)
Compiler

Compiled 43 to 29 computations (32.6% saved)

eval0.0ms (0%)

Compiler

Compiled 2 to 2 computations (0% saved)

prune3.0ms (0.1%)

Pruning

7 alts after pruning (1 fresh and 6 done)

PrunedKeptTotal
New000
Fresh011
Picked000
Done066
Total077
Error
0b
Counts
7 → 7
Alt Table
Click to see full alt table
StatusErrorProgram
40.0b
(neg.f64 x)
20.0b
(sqrt.f64 (+.f64 (*.f64 x x) y))
40.2b
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
43.7b
x
20.5b
(hypot.f64 (sqrt.f64 y) x)
43.8b
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
41.7b
(sqrt.f64 y)
Compiler

Compiled 43 to 29 computations (32.6% saved)

regimes236.0ms (6.6%)

Counts
9 → 3
Calls
Call 1
Inputs
x
x
(neg.f64 x)
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(sqrt.f64 y)
(sqrt.f64 (+.f64 (*.f64 x x) y))
(hypot.f64 (sqrt.f64 y) x)
(sqrt.f64 (fma.f64 x x y))
Outputs
(neg.f64 x)
(sqrt.f64 (fma.f64 x x y))
x
Calls

6 calls:

112.0ms
(*.f64 x x)
61.0ms
y
22.0ms
x
15.0ms
(*.f64 x x)
13.0ms
(sqrt.f64 (+.f64 (*.f64 x x) y))
Results
ErrorSegmentsBranch
0b3x
5.2b3y
8.8b2(sqrt.f64 (+.f64 (*.f64 x x) y))
8.8b2(+.f64 (*.f64 x x) y)
8.8b2(*.f64 x x)
8.8b2(*.f64 x x)
Compiler

Compiled 68 to 36 computations (47.1% saved)

bsearch54.0ms (1.5%)

Algorithm
binary-search
Steps
TimeLeftRight
32.0ms
4.2405064940455596e+86
3.780321451447804e+87
22.0ms
-4.618988821803745e+155
-1.205036194840004e+152
Compiler

Compiled 13 to 9 computations (30.8% saved)

regimes195.0ms (5.4%)

Counts
8 → 3
Calls
Call 1
Inputs
x
x
(neg.f64 x)
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(sqrt.f64 y)
(sqrt.f64 (+.f64 (*.f64 x x) y))
(hypot.f64 (sqrt.f64 y) x)
Outputs
(neg.f64 x)
(sqrt.f64 (+.f64 (*.f64 x x) y))
x
Calls

6 calls:

103.0ms
x
38.0ms
y
15.0ms
(sqrt.f64 (+.f64 (*.f64 x x) y))
13.0ms
(+.f64 (*.f64 x x) y)
12.0ms
(*.f64 x x)
Results
ErrorSegmentsBranch
0.0b3x
5.2b3y
8.8b2(sqrt.f64 (+.f64 (*.f64 x x) y))
8.8b2(+.f64 (*.f64 x x) y)
8.8b2(*.f64 x x)
8.8b2(*.f64 x x)
Compiler

Compiled 63 to 34 computations (46% saved)

bsearch97.0ms (2.7%)

Algorithm
binary-search
Steps
TimeLeftRight
32.0ms
4.2405064940455596e+86
3.780321451447804e+87
64.0ms
-4.618988821803745e+155
-1.205036194840004e+152
Compiler

Compiled 13 to 9 computations (30.8% saved)

regimes134.0ms (3.7%)

Counts
6 → 5
Calls
Call 1
Inputs
x
x
(neg.f64 x)
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(sqrt.f64 y)
Outputs
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
(sqrt.f64 y)
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
(sqrt.f64 y)
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
Calls

2 calls:

71.0ms
y
58.0ms
x
Results
ErrorSegmentsBranch
5.6b5x
25.0b6y
Compiler

Compiled 28 to 15 computations (46.4% saved)

bsearch392.0ms (10.9%)

Algorithm
binary-search
Steps
TimeLeftRight
96.0ms
3.9753815869082376e+36
5.029145804959393e+37
134.0ms
1.2288514348048841e+20
7.770298779051505e+21
88.0ms
6.331995661450691e-57
1.056269867795918e-51
73.0ms
-2.303120609406091e-94
-2.8048402946088303e-99
Compiler

Compiled 478 to 350 computations (26.8% saved)

regimes30.0ms (0.8%)

Counts
5 → 2
Calls
Call 1
Inputs
x
x
(neg.f64 x)
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
Outputs
(neg.f64 x)
(+.f64 (*.f64 1/2 (/.f64 y x)) x)
Calls

2 calls:

17.0ms
y
12.0ms
x
Results
ErrorSegmentsBranch
20.4b2x
36.2b2y
Compiler

Compiled 26 to 14 computations (46.2% saved)

bsearch50.0ms (1.4%)

Algorithm
binary-search
Steps
TimeLeftRight
50.0ms
4.6026073125008665e-269
7.231123420675729e-255
Compiler

Compiled 163 to 119 computations (27% saved)

regimes41.0ms (1.1%)

Counts
4 → 2
Calls
Call 1
Inputs
x
x
(neg.f64 x)
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
Outputs
(-.f64 (*.f64 (/.f64 y x) -1/2) x)
x
Calls

2 calls:

22.0ms
x
17.0ms
y
Results
ErrorSegmentsBranch
20.5b2x
36.2b2y
Compiler

Compiled 19 to 11 computations (42.1% saved)

bsearch57.0ms (1.6%)

Algorithm
binary-search
Steps
TimeLeftRight
57.0ms
-7.743155243996523e-153
-1.4479470988250633e-157
Compiler

Compiled 125 to 89 computations (28.8% saved)

regimes25.0ms (0.7%)

Counts
3 → 2
Calls
Call 1
Inputs
x
x
(neg.f64 x)
Outputs
(neg.f64 x)
x
Calls

2 calls:

15.0ms
y
10.0ms
x
Results
ErrorSegmentsBranch
20.5b2x
36.3b2y
Compiler

Compiled 12 to 7 computations (41.7% saved)

bsearch36.0ms (1%)

Algorithm
binary-search
Steps
TimeLeftRight
36.0ms
-2.5994692718385035e-304
6.99145857544507e-304
Compiler

Compiled 94 to 72 computations (23.4% saved)

regimes10.0ms (0.3%)

Accuracy

Total 0.0b remaining (0%)

Threshold costs 0.0b (0%)

Counts
2 → 1
Calls
Call 1
Inputs
x
x
Outputs
x
Calls

2 calls:

5.0ms
x
4.0ms
y
Results
ErrorSegmentsBranch
43.7b1x
43.7b1y
Compiler

Compiled 10 to 6 computations (40% saved)

simplify11.0ms (0.3%)

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
047155
155155
258155
359155
Stop Event
fuel
saturated
Calls
Call 1
Inputs
(if (<=.f64 x -461898882180374491305081930954147208907610703403346851104600231989952547873995861176039816674140037742753948221034974536095772696598478181570329021331275776) (neg.f64 x) (if (<=.f64 x 424050649404555962079885449986528389786092231668027316455562713318732517490200009506816) (sqrt.f64 (fma.f64 x x y)) x))
(if (<=.f64 x -461898882180374491305081930954147208907610703403346851104600231989952547873995861176039816674140037742753948221034974536095772696598478181570329021331275776) (neg.f64 x) (if (<=.f64 x 424050649404555962079885449986528389786092231668027316455562713318732517490200009506816) (sqrt.f64 (+.f64 (*.f64 x x) y)) x))
(if (<=.f64 x -112295717659724272/39402006196394479212279040100143613805079739270465446667948293404245721771497210611414266254884915640806627990306816) (-.f64 (*.f64 (/.f64 y x) -1/2) x) (if (<=.f64 x 5742252960529749/883423532389192164791648750371459257913741948437809479060803100646309888) (sqrt.f64 y) (if (<=.f64 x 7600000000000000000000) (+.f64 (*.f64 1/2 (/.f64 y x)) x) (if (<=.f64 x 4000000000000000169682549496071847936) (sqrt.f64 y) (+.f64 (*.f64 1/2 (/.f64 y x)) x)))))
(if (<=.f64 x 3803720624176790/528294531135665246352339784916516606518847326036121522127960709026673902556724859474417255887657187894674394993257128678882347559502685537250538978462939576908386683999005084168731517676426441053024232908211188404148028292751561738838396898767036476489538580897737998336) (neg.f64 x) (+.f64 (*.f64 1/2 (/.f64 y x)) x))
(if (<=.f64 x -3538089773270195/471745303102692669030194322623533611364615152552055053547025787530703020299018606510765840647207345942683546878845480612471323078740599079070179281429087005456174415872) (-.f64 (*.f64 (/.f64 y x) -1/2) x) x)
(if (<=.f64 x 0) (neg.f64 x) x)
x
Outputs
(if (<=.f64 x -461898882180374491305081930954147208907610703403346851104600231989952547873995861176039816674140037742753948221034974536095772696598478181570329021331275776) (neg.f64 x) (if (<=.f64 x 424050649404555962079885449986528389786092231668027316455562713318732517490200009506816) (sqrt.f64 (fma.f64 x x y)) x))
(if (<=.f64 x -461898882180374491305081930954147208907610703403346851104600231989952547873995861176039816674140037742753948221034974536095772696598478181570329021331275776) (neg.f64 x) (if (<=.f64 x 424050649404555962079885449986528389786092231668027316455562713318732517490200009506816) (sqrt.f64 (+.f64 (*.f64 x x) y)) x))
(if (<=.f64 x -461898882180374491305081930954147208907610703403346851104600231989952547873995861176039816674140037742753948221034974536095772696598478181570329021331275776) (neg.f64 x) (if (<=.f64 x 424050649404555962079885449986528389786092231668027316455562713318732517490200009506816) (sqrt.f64 (+.f64 y (*.f64 x x))) x))
(if (<=.f64 x -112295717659724272/39402006196394479212279040100143613805079739270465446667948293404245721771497210611414266254884915640806627990306816) (-.f64 (*.f64 (/.f64 y x) -1/2) x) (if (<=.f64 x 5742252960529749/883423532389192164791648750371459257913741948437809479060803100646309888) (sqrt.f64 y) (if (<=.f64 x 7600000000000000000000) (+.f64 (*.f64 1/2 (/.f64 y x)) x) (if (<=.f64 x 4000000000000000169682549496071847936) (sqrt.f64 y) (+.f64 (*.f64 1/2 (/.f64 y x)) x)))))
(if (<=.f64 x -7018482353732767/2462625387274654950767440006258975862817483704404090416746768337765357610718575663213391640930307227550414249394176) (-.f64 (*.f64 (/.f64 y x) -1/2) x) (if (<=.f64 x 5742252960529749/883423532389192164791648750371459257913741948437809479060803100646309888) (sqrt.f64 y) (if (<=.f64 x 7600000000000000000000) (+.f64 x (*.f64 (/.f64 y x) 1/2)) (if (<=.f64 x 4000000000000000169682549496071847936) (sqrt.f64 y) (+.f64 x (*.f64 (/.f64 y x) 1/2))))))
(if (<=.f64 x 3803720624176790/528294531135665246352339784916516606518847326036121522127960709026673902556724859474417255887657187894674394993257128678882347559502685537250538978462939576908386683999005084168731517676426441053024232908211188404148028292751561738838396898767036476489538580897737998336) (neg.f64 x) (+.f64 (*.f64 1/2 (/.f64 y x)) x))
(if (<=.f64 x 1901860312088395/264147265567832623176169892458258303259423663018060761063980354513336951278362429737208627943828593947337197496628564339441173779751342768625269489231469788454193341999502542084365758838213220526512116454105594202074014146375780869419198449383518238244769290448868999168) (neg.f64 x) (+.f64 x (*.f64 (/.f64 y x) 1/2)))
(if (<=.f64 x -3538089773270195/471745303102692669030194322623533611364615152552055053547025787530703020299018606510765840647207345942683546878845480612471323078740599079070179281429087005456174415872) (-.f64 (*.f64 (/.f64 y x) -1/2) x) x)
(if (<=.f64 x 0) (neg.f64 x) x)
x
Compiler

Compiled 121 to 75 computations (38% saved)

soundness105.0ms (2.9%)

Algorithm
egg-herbie
Rules
759×times-frac_binary64
598×associate-/l*_binary64
502×unswap-sqr_binary64
498×fma-neg_binary64
427×cancel-sign-sub-inv_binary64
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
067
187
059313
1159311
2440311
31580309
45737309
067
187
067
187
Stop Event
saturated
saturated
node limit
saturated
Compiler

Compiled 54 to 34 computations (37% saved)

end62.0ms (1.7%)

Compiler

Compiled 109 to 58 computations (46.8% saved)

Profiling

Loading profile data...