Details

Time bar (total: 1.3s)

analyze33.0ms (2.5%)

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
2.8%0.5%96.7%13
2.9%0.3%96.7%14
Compiler

Compiled 26 to 15 computations (42.3% saved)

sample993.0ms (75.1%)

Results
936.0ms8256×body128valid
44.0ms357×body128invalid

preprocess22.0ms (1.7%)

Algorithm
egg-herbie
Rules
303×associate-/l*_binary64
282×associate-/r/_binary64
258×times-frac_binary64
232×fma-def_binary64
190×distribute-rgt-in_binary64
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
01340
14734
213730
339430
4128230
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)))
Compiler

Compiled 15 to 8 computations (46.7% saved)

simplify21.0ms (1.6%)

Algorithm
egg-herbie
Rules
704×fma-def_binary64
600×associate-/r*_binary64
362×distribute-rgt-in_binary64
331×distribute-lft-in_binary64
206×div-sub_binary64
Iterations

Useful iterations: 2 (0.0ms)

IterNodesCost
0920
12617
27215
320115
454915
5197915
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.3b
Counts
3 → 1
Alt Table
StatusErrorProgram
14.3b
(sqrt.f64 (fabs.f64 (fma.f64 b (/.f64 b (*.f64 a a)) -1)))
Compiler

Compiled 49 to 29 computations (40.8% saved)

localize7.0ms (0.5%)

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
a
@0
(/.f64 b (*.f64 a a))
1.0ms
b
@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))

rewrite40.0ms (3%)

Algorithm
batch-egg-rewrite
Rules
772×log1p-expm1-u_binary64
772×expm1-log1p-u_binary64
92×add-sqr-sqrt_binary64
89×*-un-lft-identity_binary64
89×pow1_binary64
Iterations

Useful iterations: 1 (0.0ms)

IterNodesCost
0922
118320
2234420
Stop Event
node limit
Counts
2 → 57
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 (exp.f64 (log1p.f64 (*.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 log.f64 (+.f64 1 (expm1.f64 (*.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 expm1.f64 (log1p.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 exp.f64 (log.f64 (*.f64 b (pow.f64 a -2))))))) (#(struct:change #<rule egg-rr> (2) ((x exp.f64 (*.f64 (log.f64 (*.f64 b (pow.f64 a -2))) 1))))))
((#(struct:change #<rule egg-rr> (2) ((x -.f64 (exp.f64 (log1p.f64 (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 log.f64 (+.f64 1 (expm1.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 expm1.f64 (log1p.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 exp.f64 (log.f64 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1))))))) (#(struct:change #<rule egg-rr> (2) ((x exp.f64 (*.f64 (log.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1)) 1/2))))) (#(struct:change #<rule egg-rr> (2) ((x exp.f64 (*.f64 (log.f64 (sqrt.f64 (fma.f64 b (*.f64 b (pow.f64 a -2)) -1))) 1))))))

simplify52.0ms (4%)

Algorithm
egg-herbie
Rules
693×associate-/r/_binary64
497×fma-def_binary64
434×unswap-sqr_binary64
232×associate-/r*_binary64
211×times-frac_binary64
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
010528
117528
230528
355528
4136528
5444528
61879528
73242528
83964528
94573528
Stop Event
node limit
Counts
105 → 61
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)))

prune29.0ms (2.2%)

Pruning

1 alts after pruning (1 fresh and 0 done)

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

Compiled 863 to 593 computations (31.3% 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.3%)

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
a
@-inf
(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)))

rewrite35.0ms (2.6%)

Algorithm
batch-egg-rewrite
Rules
710×pow1_binary64
656×add-log-exp_binary64
656×log1p-expm1-u_binary64
656×expm1-log1p-u_binary64
80×add-sqr-sqrt_binary64
Iterations

Useful iterations: 1 (0.0ms)

IterNodesCost
0815
116012
2198612
Stop Event
node limit
Counts
1 → 21
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 (exp.f64 (log1p.f64 (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 log.f64 (+.f64 1 (expm1.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 expm1.f64 (log1p.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 exp.f64 (log.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
21 → 21
Calls
Call 1
Inputs
Outputs

prune7.0ms (0.5%)

Pruning

1 alts after pruning (0 fresh and 1 done)

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

Compiled 280 to 198 computations (29.3% saved)

regimes22.0ms (1.7%)

Accuracy

Total 0.0b remaining (0%)

Threshold costs 0.0b (0%)

Counts
3 → 1
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)))

end35.0ms (2.6%)

Stop Event
done
Compiler

Compiled 36 to 21 computations (41.7% saved)

Profiling

Loading profile data...