Details

Time bar (total: 1.4s)

analyze15.0ms (1.1%)

Algorithm
search
Search
TrueOtherFalseIter
0%37.5%62.5%0
0%37.5%62.5%1
0%18.7%81.3%2
0%9.4%90.6%3
0%9.4%90.6%4
0%7%93%5
1.2%4.7%94.1%6
1.2%3.5%95.3%7
1.8%2.6%95.6%8
2%1.9%96%9
2.4%1.3%96.3%10
2.6%1%96.5%11
2.7%0.7%96.6%12
Compiler

Compiled 26 to 15 computations (42.3% saved)

sample1.0s (72.2%)

Results
915.0ms8256×body64valid
93.0ms784×body64invalid

preprocess43.0ms (3%)

Algorithm
egg-herbie
Rules
746×fma-neg_binary64
359×fabs-div_binary64
352×distribute-neg-frac_binary64
303×associate-/l*_binary64
297×fabs-mul_binary64
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
01340
14734
213730
339430
4128230
5545830
022
Stop Event
saturated
node limit
Calls
Call 1
Inputs
0
1
Outputs
0
1
Call 2
Inputs
(sqrt.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 a a) (*.f64 b b)) (*.f64 a a))))
(sqrt.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 b b) (*.f64 a a)) (*.f64 b b))))
Outputs
(sqrt.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 a a) (*.f64 b b)) (*.f64 a a))))
(sqrt.f64 (fabs.f64 (-.f64 1 (*.f64 (/.f64 b a) (/.f64 b a)))))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b (*.f64 a a)) b -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 b b) (*.f64 a a)) (*.f64 b b))))
(sqrt.f64 (fabs.f64 (-.f64 1 (/.f64 (*.f64 a a) (*.f64 b b)))))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 a b) (/.f64 a b) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 a (/.f64 a (*.f64 b b)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 a (/.f64 (/.f64 a b) b) -1)))
Compiler

Compiled 15 to 8 computations (46.7% saved)

simplify40.0ms (2.8%)

Algorithm
egg-herbie
Rules
704×fma-def_binary64
600×associate-/r*_binary64
558×associate-*l*_binary64
489×associate-*r*_binary64
362×distribute-rgt-in_binary64
Iterations

Useful iterations: 2 (0.0ms)

IterNodesCost
0920
12617
27215
320115
454915
5197915
6766915
Stop Event
node limit
Counts
1 → 2
Calls
Call 1
Inputs
(sqrt.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 a a) (*.f64 b b)) (*.f64 a a))))
Outputs
(sqrt.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 a a) (*.f64 b b)) (*.f64 a a))))
(sqrt.f64 (fabs.f64 (-.f64 1 (*.f64 (/.f64 b a) (/.f64 b a)))))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))

prune3.0ms (0.2%)

Pruning

1 alts after pruning (1 fresh and 0 done)

PrunedKeptTotal
New112
Fresh101
Picked000
Done000
Total213
Error
14.9b
Counts
3 → 1
Alt Table
StatusErrorProgram
14.9b
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
Compiler

Compiled 49 to 29 computations (40.8% saved)

localize8.0ms (0.6%)

Local error

Found 2 expressions with local error:

NewErrorProgram
0.0b
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
14.7b
(/.f64 b (*.f64 a a))

series7.0ms (0.5%)

Counts
2 → 48
Calls

12 calls:

TimeVariablePointExpression
1.0ms
b
@0
(/.f64 b (*.f64 a a))
1.0ms
a
@0
(/.f64 b (*.f64 a a))
1.0ms
b
@inf
(/.f64 b (*.f64 a a))
1.0ms
b
@-inf
(/.f64 b (*.f64 a a))
0.0ms
a
@inf
(/.f64 b (*.f64 a a))

rewrite48.0ms (3.4%)

Algorithm
batch-egg-rewrite
Rules
832×pow1_binary64
772×add-log-exp_binary64
772×log1p-expm1-u_binary64
772×expm1-log1p-u_binary64
752×add-exp-log_binary64
Iterations

Useful iterations: 1 (0.0ms)

IterNodesCost
0922
118320
2234420
Stop Event
node limit
Counts
2 → 52
Calls
Call 1
Inputs
(/.f64 b (*.f64 a a))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
Outputs
((#(struct:change #<rule egg-rr> (2) ((x -.f64 (+.f64 1 (*.f64 b (pow.f64 a -2))) 1)))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 b (pow.f64 a -2))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (*.f64 b (pow.f64 a -2)) 1)))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 1 (*.f64 b (pow.f64 a -2)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (cbrt.f64 (*.f64 b (pow.f64 a -2))) (cbrt.f64 (/.f64 (*.f64 b b) (pow.f64 a 4))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (cbrt.f64 (/.f64 (*.f64 b b) (pow.f64 a 4))) (cbrt.f64 (*.f64 b (pow.f64 a -2))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (pow.f64 (cbrt.f64 b) 2) (*.f64 (cbrt.f64 b) (pow.f64 a -2)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (/.f64 (sqrt.f64 b) a) (/.f64 (sqrt.f64 b) a))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (sqrt.f64 b) (*.f64 (sqrt.f64 b) (pow.f64 a -2)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (neg.f64 b) (/.f64 1 (*.f64 a (neg.f64 a))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (pow.f64 a -2) b)))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (/.f64 b a) (/.f64 1 a))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (/.f64 (pow.f64 (cbrt.f64 b) 2) 1) (/.f64 (cbrt.f64 b) (*.f64 a a)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (/.f64 (sqrt.f64 b) 1) (/.f64 (sqrt.f64 b) (*.f64 a a)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (/.f64 1 a) (/.f64 b a))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (/.f64 (pow.f64 (cbrt.f64 b) 2) a) (/.f64 (cbrt.f64 b) a))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (/.f64 1 (cbrt.f64 (pow.f64 a 4))) (/.f64 b (cbrt.f64 (*.f64 a a))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (/.f64 (pow.f64 (cbrt.f64 b) 2) (cbrt.f64 (pow.f64 a 4))) (cbrt.f64 (*.f64 b (pow.f64 a -2))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (/.f64 (sqrt.f64 b) (cbrt.f64 (pow.f64 a 4))) (/.f64 (sqrt.f64 b) (cbrt.f64 (*.f64 a a))))))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (*.f64 b (pow.f64 a -2)) 1)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (pow.f64 (*.f64 b (pow.f64 a -2)) 3) 1/3)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (cbrt.f64 (*.f64 b (pow.f64 a -2))) 3)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (/.f64 (sqrt.f64 b) a) 2)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (/.f64 a (/.f64 b a)) -1)))) (#(struct:change #<rule egg-rr> (2) ((x neg.f64 (/.f64 b (*.f64 a (neg.f64 a))))))) (#(struct:change #<rule egg-rr> (2) ((x sqrt.f64 (/.f64 (*.f64 b b) (pow.f64 a 4)))))) (#(struct:change #<rule egg-rr> (2) ((x log.f64 (pow.f64 (exp.f64 b) (pow.f64 a -2)))))) (#(struct:change #<rule egg-rr> (2) ((x cbrt.f64 (pow.f64 (*.f64 b (pow.f64 a -2)) 3))))) (#(struct:change #<rule egg-rr> (2) ((x cbrt.f64 (/.f64 (pow.f64 b 3) (pow.f64 (*.f64 a a) 3)))))) (#(struct:change #<rule egg-rr> (2) ((x exp.f64 (log.f64 (*.f64 b (pow.f64 a -2))))))) (#(struct:change #<rule egg-rr> (2) ((x log1p.f64 (expm1.f64 (*.f64 b (pow.f64 a -2))))))) (#(struct:change #<rule egg-rr> (2) ((x expm1.f64 (log1p.f64 (*.f64 b (pow.f64 a -2))))))))
((#(struct:change #<rule egg-rr> (2) ((x -.f64 (+.f64 1 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1))) 1)))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1)) 1)))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 1 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (cbrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1)) (cbrt.f64 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (cbrt.f64 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1))) (cbrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (pow.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1) 1/4) (pow.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1) 1/4))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (sqrt.f64 (pow.f64 (cbrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1)) 2)) (sqrt.f64 (cbrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (pow.f64 1 1/2) (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1)))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (pow.f64 (pow.f64 (cbrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1)) 2) 1/2) (pow.f64 (cbrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1)) 1/2))))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1) 1/2)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1)) 1)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (pow.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1) 3/2) 1/3)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (cbrt.f64 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1))) 3)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (pow.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1) 1/4) 2)))) (#(struct:change #<rule egg-rr> (2) ((x fabs.f64 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1)))))) (#(struct:change #<rule egg-rr> (2) ((x log.f64 (exp.f64 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1))))))) (#(struct:change #<rule egg-rr> (2) ((x cbrt.f64 (pow.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1) 3/2))))) (#(struct:change #<rule egg-rr> (2) ((x exp.f64 (log.f64 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1))))))) (#(struct:change #<rule egg-rr> (2) ((x log1p.f64 (expm1.f64 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1))))))) (#(struct:change #<rule egg-rr> (2) ((x expm1.f64 (log1p.f64 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1))))))))

simplify64.0ms (4.6%)

Algorithm
egg-herbie
Rules
966×fma-def_binary64
693×associate-/r/_binary64
434×unswap-sqr_binary64
239×cancel-sign-sub-inv_binary64
232×associate-/r*_binary64
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
010528
117528
230528
355528
4136528
5444528
61879528
73242528
83964528
94573528
105561528
Stop Event
node limit
Counts
100 → 56
Calls
Call 1
Inputs
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(/.f64 b (pow.f64 a 2))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
Outputs
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(/.f64 b (pow.f64 a 2))
(/.f64 b (*.f64 a a))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (pow.f64 a 2)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))

prune61.0ms (4.3%)

Pruning

1 alts after pruning (1 fresh and 0 done)

PrunedKeptTotal
New55156
Fresh000
Picked101
Done000
Total56157
Error
0.0b
Counts
57 → 1
Alt Table
StatusErrorProgram
0.0b
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
Compiler

Compiled 797 to 542 computations (32% saved)

localize5.0ms (0.4%)

Local error

Found 1 expressions with local error:

NewErrorProgram
0.0b
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))

series3.0ms (0.2%)

Counts
1 → 0
Calls

6 calls:

TimeVariablePointExpression
1.0ms
b
@inf
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
1.0ms
b
@0
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
0.0ms
b
@-inf
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
0.0ms
a
@-inf
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
0.0ms
a
@inf
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))

rewrite41.0ms (2.9%)

Algorithm
batch-egg-rewrite
Rules
710×pow1_binary64
656×add-log-exp_binary64
656×log1p-expm1-u_binary64
656×expm1-log1p-u_binary64
641×add-exp-log_binary64
Iterations

Useful iterations: 1 (0.0ms)

IterNodesCost
0815
116012
2198612
Stop Event
node limit
Counts
1 → 20
Calls
Call 1
Inputs
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
Outputs
((#(struct:change #<rule egg-rr> (2) ((x -.f64 (+.f64 1 (sqrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)))) 1)))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (sqrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2))) 1)))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 1 (sqrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (cbrt.f64 (sqrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)))) (cbrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (cbrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2))) (cbrt.f64 (sqrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (pow.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)) 1/4) (pow.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)) 1/4))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (sqrt.f64 (pow.f64 (cbrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2))) 2)) (sqrt.f64 (cbrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (pow.f64 1 1/2) (sqrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2))))))) (#(struct:change #<rule egg-rr> (2) ((x *.f64 (pow.f64 (pow.f64 (cbrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2))) 2) 1/2) (pow.f64 (cbrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2))) 1/2))))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)) 1/2)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (sqrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2))) 1)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (pow.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)) 3/2) 1/3)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (cbrt.f64 (sqrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)))) 3)))) (#(struct:change #<rule egg-rr> (2) ((x pow.f64 (pow.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)) 1/4) 2)))) (#(struct:change #<rule egg-rr> (2) ((x fabs.f64 (sqrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2))))))) (#(struct:change #<rule egg-rr> (2) ((x log.f64 (exp.f64 (sqrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)))))))) (#(struct:change #<rule egg-rr> (2) ((x cbrt.f64 (pow.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)) 3/2))))) (#(struct:change #<rule egg-rr> (2) ((x exp.f64 (log.f64 (sqrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)))))))) (#(struct:change #<rule egg-rr> (2) ((x log1p.f64 (expm1.f64 (sqrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)))))))) (#(struct:change #<rule egg-rr> (2) ((x expm1.f64 (log1p.f64 (sqrt.f64 (+.f64 -1 (pow.f64 (/.f64 b a) 2)))))))))

simplify6.0ms (0.4%)

Algorithm
egg-herbie
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
000
Stop Event
saturated
Counts
20 → 20
Calls
Call 1
Inputs
Outputs

prune10.0ms (0.7%)

Pruning

1 alts after pruning (0 fresh and 1 done)

PrunedKeptTotal
New20020
Fresh000
Picked011
Done000
Total20121
Error
0.0b
Counts
21 → 1
Alt Table
StatusErrorProgram
0.0b
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
Compiler

Compiled 268 to 189 computations (29.5% saved)

regimes22.0ms (1.6%)

Accuracy

Total -14.9b remaining (-127547.4%)

Threshold costs -14.9b (-127547.4%)

Counts
3 → 1
Calls
Call 1
Inputs
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
(sqrt.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 a a) (*.f64 b b)) (*.f64 a a))))
Outputs
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
Results
5.0ms0.0ba
3.0ms0.0bb
2.0ms0.0b(sqrt.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 a a) (*.f64 b b)) (*.f64 a a))))
2.0ms0.0b(fabs.f64 (/.f64 (-.f64 (*.f64 a a) (*.f64 b b)) (*.f64 a a)))
2.0ms0.0b(/.f64 (-.f64 (*.f64 a a) (*.f64 b b)) (*.f64 a a))
3.0ms0.0b(*.f64 a a)
3.0ms0.0b(*.f64 a a)
2.0ms0.0b(*.f64 b b)
Compiler

Compiled 98 to 52 computations (46.9% saved)

simplify2.0ms (0.1%)

Algorithm
egg-herbie
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
0815
Stop Event
saturated
Calls
Call 1
Inputs
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))
Outputs
(sqrt.f64 (fabs.f64 (fma.f64 (/.f64 b a) (/.f64 b a) -1)))

end15.0ms (1.1%)

Stop Event
done
Compiler

Compiled 24 to 14 computations (41.7% saved)

Profiling

Loading profile data...