Diagrams.Solve.Polynomial:quadForm from diagrams-solve-0.1, C

Time bar (total: 1.3s)

analyze34.0ms (2.6%)

Algorithm
search
Search
ProbabilityValidUnknownPreconditionInfiniteDomainCan'tIter
0%0%99.9%0.1%0%0%0%0
0%0%99.9%0.1%0%0%0%1
0%0%99.9%0.1%0%0%0%2
50%50%50%0.1%0%0%0%3
50%50%50%0.1%0%0%0%4
75%74.9%25%0.1%0%0%0%5
75%74.9%25%0.1%0%0%0%6
87.5%87.4%12.5%0.1%0%0%0%7
87.5%87.4%12.5%0.1%0%0%0%8
93.8%93.7%6.2%0.1%0%0%0%9
93.8%93.7%6.2%0.1%0%0%0%10
96.9%96.8%3.1%0.1%0%0%0%11
96.9%96.8%3.1%0.1%0%0%0%12
Compiler

Compiled 8 to 6 computations (25% saved)

sample910.0ms (69.3%)

Results
801.0ms7183×body256valid
96.0ms1073×body256infinite
Bogosity

preprocess46.0ms (3.5%)

Algorithm
egg-herbie
Rules
124×fma-def
64×div-sub
50×sub-neg
48×associate-/r*
36×fma-neg
Problems
255×No Errors
(/.f64 x (*.f64 y 2))
Iterations

Useful iterations: 1 (0.0ms)

IterNodesCost
01573
13465
26965
311865
415465
519465
625465
735265
835965
Stop Event
saturated
Calls
Call 1
Inputs
(/.f64 x (*.f64 y 2))
(/.f64 x (*.f64 y 2))
(/.f64 (neg.f64 x) (*.f64 y 2))
(/.f64 x (*.f64 (neg.f64 y) 2))
(/.f64 y (*.f64 x 2))
Outputs
(/.f64 x (*.f64 y 2))
(/.f64 x (*.f64 y 2))
(/.f64 (neg.f64 x) (*.f64 y 2))
(*.f64 -1/2 (/.f64 x y))
(*.f64 (/.f64 x y) -1/2)
(*.f64 x (/.f64 -1/2 y))
(/.f64 x (*.f64 (neg.f64 y) 2))
(/.f64 (neg.f64 x) (*.f64 y 2))
(*.f64 -1/2 (/.f64 x y))
(*.f64 (/.f64 x y) -1/2)
(*.f64 x (/.f64 -1/2 y))
(/.f64 y (*.f64 x 2))
Compiler

Compiled 21 to 14 computations (33.3% saved)

eval0.0ms (0%)

Compiler

Compiled 2 to 2 computations (0% saved)

prune1.0ms (0.1%)

Alt Table
Click to see full alt table
StatusAccuracyProgram
99.6%
(/.f64 x (*.f64 y 2))
99.6%
(/.f64 x (*.f64 y 2))
Compiler

Compiled 14 to 10 computations (28.6% saved)

localize24.0ms (1.9%)

Localize:

Found 1 expressions with local error:

NewAccuracyProgram
99.6%
(/.f64 x (*.f64 y 2))
Compiler

Compiled 14 to 9 computations (35.7% saved)

series3.0ms (0.2%)

Counts
1 → 24
Calls

6 calls:

TimeVariablePointExpression
1.0ms
y
@0
(/.f64 x (*.f64 y 2))
1.0ms
x
@inf
(/.f64 x (*.f64 y 2))
1.0ms
x
@0
(/.f64 x (*.f64 y 2))
0.0ms
y
@inf
(/.f64 x (*.f64 y 2))
0.0ms
x
@-inf
(/.f64 x (*.f64 y 2))

rewrite81.0ms (6.2%)

Algorithm
batch-egg-rewrite
Rules
1558×pow1
1444×add-exp-log
1444×log1p-expm1-u
1444×expm1-log1p-u
1414×add-log-exp
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
0713
114713
2191213
Stop Event
node limit
Counts
1 → 37
Calls
Call 1
Inputs
(/.f64 x (*.f64 y 2))
Outputs
(-.f64 (+.f64 1 (*.f64 x (/.f64 1/2 y))) 1)
(*.f64 x (/.f64 1/2 y))
(*.f64 (*.f64 x (/.f64 1/2 y)) 1)
(*.f64 1 (*.f64 x (/.f64 1/2 y)))
(*.f64 (pow.f64 (cbrt.f64 x) 2) (*.f64 (cbrt.f64 x) (/.f64 1/2 y)))
(*.f64 (cbrt.f64 (*.f64 x (/.f64 1/2 y))) (pow.f64 (cbrt.f64 (*.f64 x (/.f64 1/2 y))) 2))
(*.f64 (pow.f64 (cbrt.f64 (*.f64 x (/.f64 1/2 y))) 2) (cbrt.f64 (*.f64 x (/.f64 1/2 y))))
(*.f64 (sqrt.f64 x) (*.f64 (sqrt.f64 x) (/.f64 1/2 y)))
(*.f64 (sqrt.f64 (*.f64 x (/.f64 1/2 y))) (sqrt.f64 (*.f64 x (/.f64 1/2 y))))
(*.f64 (neg.f64 x) (/.f64 1 (*.f64 y -2)))
(*.f64 (/.f64 1/2 y) x)
(*.f64 (/.f64 x y) 1/2)
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) 1) (/.f64 (cbrt.f64 x) (*.f64 y 2)))
(*.f64 (/.f64 (sqrt.f64 x) 1) (/.f64 (sqrt.f64 x) (*.f64 y 2)))
(*.f64 (/.f64 (sqrt.f64 x) y) (/.f64 (sqrt.f64 x) 2))
(*.f64 (/.f64 1 y) (/.f64 x 2))
(*.f64 (/.f64 1 (pow.f64 (cbrt.f64 (*.f64 y 2)) 2)) (/.f64 x (cbrt.f64 (*.f64 y 2))))
(*.f64 (/.f64 1 (sqrt.f64 (*.f64 y 2))) (/.f64 x (sqrt.f64 (*.f64 y 2))))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) y) (/.f64 (cbrt.f64 x) 2))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) 2) (/.f64 (cbrt.f64 x) y))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) (pow.f64 (cbrt.f64 (*.f64 y 2)) 2)) (cbrt.f64 (*.f64 x (/.f64 1/2 y))))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) (sqrt.f64 (*.f64 y 2))) (/.f64 (cbrt.f64 x) (sqrt.f64 (*.f64 y 2))))
(*.f64 (/.f64 (sqrt.f64 x) 2) (/.f64 (sqrt.f64 x) y))
(*.f64 (/.f64 (sqrt.f64 x) (pow.f64 (cbrt.f64 (*.f64 y 2)) 2)) (/.f64 (sqrt.f64 x) (cbrt.f64 (*.f64 y 2))))
(pow.f64 (*.f64 x (/.f64 1/2 y)) 1)
(pow.f64 (cbrt.f64 (*.f64 x (/.f64 1/2 y))) 3)
(pow.f64 (pow.f64 (*.f64 x (/.f64 1/2 y)) 3) 1/3)
(pow.f64 (sqrt.f64 (*.f64 x (/.f64 1/2 y))) 2)
(pow.f64 (/.f64 y (/.f64 x 2)) -1)
(neg.f64 (/.f64 x (*.f64 y -2)))
(sqrt.f64 (pow.f64 (*.f64 x (/.f64 1/2 y)) 2))
(log.f64 (sqrt.f64 (exp.f64 (/.f64 x y))))
(cbrt.f64 (pow.f64 (*.f64 x (/.f64 1/2 y)) 3))
(expm1.f64 (log1p.f64 (*.f64 x (/.f64 1/2 y))))
(exp.f64 (log.f64 (*.f64 x (/.f64 1/2 y))))
(log1p.f64 (expm1.f64 (*.f64 x (/.f64 1/2 y))))

simplify65.0ms (5%)

Algorithm
egg-herbie
Rules
1064×fma-def
808×log-prod
736×distribute-lft-in
730×log-div
524×cube-div
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
097841
1236841
2789841
33622841
Stop Event
node limit
Counts
61 → 38
Calls
Call 1
Inputs
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(-.f64 (+.f64 1 (*.f64 x (/.f64 1/2 y))) 1)
(*.f64 x (/.f64 1/2 y))
(*.f64 (*.f64 x (/.f64 1/2 y)) 1)
(*.f64 1 (*.f64 x (/.f64 1/2 y)))
(*.f64 (pow.f64 (cbrt.f64 x) 2) (*.f64 (cbrt.f64 x) (/.f64 1/2 y)))
(*.f64 (cbrt.f64 (*.f64 x (/.f64 1/2 y))) (pow.f64 (cbrt.f64 (*.f64 x (/.f64 1/2 y))) 2))
(*.f64 (pow.f64 (cbrt.f64 (*.f64 x (/.f64 1/2 y))) 2) (cbrt.f64 (*.f64 x (/.f64 1/2 y))))
(*.f64 (sqrt.f64 x) (*.f64 (sqrt.f64 x) (/.f64 1/2 y)))
(*.f64 (sqrt.f64 (*.f64 x (/.f64 1/2 y))) (sqrt.f64 (*.f64 x (/.f64 1/2 y))))
(*.f64 (neg.f64 x) (/.f64 1 (*.f64 y -2)))
(*.f64 (/.f64 1/2 y) x)
(*.f64 (/.f64 x y) 1/2)
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) 1) (/.f64 (cbrt.f64 x) (*.f64 y 2)))
(*.f64 (/.f64 (sqrt.f64 x) 1) (/.f64 (sqrt.f64 x) (*.f64 y 2)))
(*.f64 (/.f64 (sqrt.f64 x) y) (/.f64 (sqrt.f64 x) 2))
(*.f64 (/.f64 1 y) (/.f64 x 2))
(*.f64 (/.f64 1 (pow.f64 (cbrt.f64 (*.f64 y 2)) 2)) (/.f64 x (cbrt.f64 (*.f64 y 2))))
(*.f64 (/.f64 1 (sqrt.f64 (*.f64 y 2))) (/.f64 x (sqrt.f64 (*.f64 y 2))))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) y) (/.f64 (cbrt.f64 x) 2))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) 2) (/.f64 (cbrt.f64 x) y))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) (pow.f64 (cbrt.f64 (*.f64 y 2)) 2)) (cbrt.f64 (*.f64 x (/.f64 1/2 y))))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) (sqrt.f64 (*.f64 y 2))) (/.f64 (cbrt.f64 x) (sqrt.f64 (*.f64 y 2))))
(*.f64 (/.f64 (sqrt.f64 x) 2) (/.f64 (sqrt.f64 x) y))
(*.f64 (/.f64 (sqrt.f64 x) (pow.f64 (cbrt.f64 (*.f64 y 2)) 2)) (/.f64 (sqrt.f64 x) (cbrt.f64 (*.f64 y 2))))
(pow.f64 (*.f64 x (/.f64 1/2 y)) 1)
(pow.f64 (cbrt.f64 (*.f64 x (/.f64 1/2 y))) 3)
(pow.f64 (pow.f64 (*.f64 x (/.f64 1/2 y)) 3) 1/3)
(pow.f64 (sqrt.f64 (*.f64 x (/.f64 1/2 y))) 2)
(pow.f64 (/.f64 y (/.f64 x 2)) -1)
(neg.f64 (/.f64 x (*.f64 y -2)))
(sqrt.f64 (pow.f64 (*.f64 x (/.f64 1/2 y)) 2))
(log.f64 (sqrt.f64 (exp.f64 (/.f64 x y))))
(cbrt.f64 (pow.f64 (*.f64 x (/.f64 1/2 y)) 3))
(expm1.f64 (log1p.f64 (*.f64 x (/.f64 1/2 y))))
(exp.f64 (log.f64 (*.f64 x (/.f64 1/2 y))))
(log1p.f64 (expm1.f64 (*.f64 x (/.f64 1/2 y))))
Outputs
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(-.f64 (+.f64 1 (*.f64 x (/.f64 1/2 y))) 1)
(*.f64 1/2 (/.f64 x y))
(*.f64 x (/.f64 1/2 y))
(*.f64 1/2 (/.f64 x y))
(*.f64 (*.f64 x (/.f64 1/2 y)) 1)
(*.f64 1/2 (/.f64 x y))
(*.f64 1 (*.f64 x (/.f64 1/2 y)))
(*.f64 1/2 (/.f64 x y))
(*.f64 (pow.f64 (cbrt.f64 x) 2) (*.f64 (cbrt.f64 x) (/.f64 1/2 y)))
(*.f64 1/2 (/.f64 x y))
(*.f64 (cbrt.f64 (*.f64 x (/.f64 1/2 y))) (pow.f64 (cbrt.f64 (*.f64 x (/.f64 1/2 y))) 2))
(*.f64 1/2 (/.f64 x y))
(*.f64 (pow.f64 (cbrt.f64 (*.f64 x (/.f64 1/2 y))) 2) (cbrt.f64 (*.f64 x (/.f64 1/2 y))))
(*.f64 1/2 (/.f64 x y))
(*.f64 (sqrt.f64 x) (*.f64 (sqrt.f64 x) (/.f64 1/2 y)))
(*.f64 1/2 (/.f64 x y))
(*.f64 (sqrt.f64 (*.f64 x (/.f64 1/2 y))) (sqrt.f64 (*.f64 x (/.f64 1/2 y))))
(*.f64 1/2 (/.f64 x y))
(*.f64 (neg.f64 x) (/.f64 1 (*.f64 y -2)))
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 1/2 y) x)
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 x y) 1/2)
(*.f64 1/2 (/.f64 x y))
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) 1) (/.f64 (cbrt.f64 x) (*.f64 y 2)))
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 (sqrt.f64 x) 1) (/.f64 (sqrt.f64 x) (*.f64 y 2)))
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 (sqrt.f64 x) y) (/.f64 (sqrt.f64 x) 2))
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 1 y) (/.f64 x 2))
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 1 (pow.f64 (cbrt.f64 (*.f64 y 2)) 2)) (/.f64 x (cbrt.f64 (*.f64 y 2))))
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 1 (sqrt.f64 (*.f64 y 2))) (/.f64 x (sqrt.f64 (*.f64 y 2))))
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) y) (/.f64 (cbrt.f64 x) 2))
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) 2) (/.f64 (cbrt.f64 x) y))
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) (pow.f64 (cbrt.f64 (*.f64 y 2)) 2)) (cbrt.f64 (*.f64 x (/.f64 1/2 y))))
(*.f64 (cbrt.f64 (*.f64 1/2 (/.f64 x y))) (/.f64 (pow.f64 (cbrt.f64 x) 2) (pow.f64 (cbrt.f64 (*.f64 y 2)) 2)))
(*.f64 (/.f64 (pow.f64 (cbrt.f64 x) 2) (sqrt.f64 (*.f64 y 2))) (/.f64 (cbrt.f64 x) (sqrt.f64 (*.f64 y 2))))
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 (sqrt.f64 x) 2) (/.f64 (sqrt.f64 x) y))
(*.f64 1/2 (/.f64 x y))
(*.f64 (/.f64 (sqrt.f64 x) (pow.f64 (cbrt.f64 (*.f64 y 2)) 2)) (/.f64 (sqrt.f64 x) (cbrt.f64 (*.f64 y 2))))
(*.f64 1/2 (/.f64 x y))
(pow.f64 (*.f64 x (/.f64 1/2 y)) 1)
(*.f64 1/2 (/.f64 x y))
(pow.f64 (cbrt.f64 (*.f64 x (/.f64 1/2 y))) 3)
(*.f64 1/2 (/.f64 x y))
(pow.f64 (pow.f64 (*.f64 x (/.f64 1/2 y)) 3) 1/3)
(*.f64 1/2 (/.f64 x y))
(pow.f64 (sqrt.f64 (*.f64 x (/.f64 1/2 y))) 2)
(*.f64 1/2 (/.f64 x y))
(pow.f64 (/.f64 y (/.f64 x 2)) -1)
(*.f64 1/2 (/.f64 x y))
(neg.f64 (/.f64 x (*.f64 y -2)))
(*.f64 1/2 (/.f64 x y))
(sqrt.f64 (pow.f64 (*.f64 x (/.f64 1/2 y)) 2))
(*.f64 1/2 (/.f64 x y))
(log.f64 (sqrt.f64 (exp.f64 (/.f64 x y))))
(cbrt.f64 (pow.f64 (*.f64 x (/.f64 1/2 y)) 3))
(*.f64 1/2 (/.f64 x y))
(expm1.f64 (log1p.f64 (*.f64 x (/.f64 1/2 y))))
(*.f64 1/2 (/.f64 x y))
(exp.f64 (log.f64 (*.f64 x (/.f64 1/2 y))))
(*.f64 1/2 (/.f64 x y))
(log1p.f64 (expm1.f64 (*.f64 x (/.f64 1/2 y))))
(*.f64 1/2 (/.f64 x y))

eval9.0ms (0.7%)

Compiler

Compiled 373 to 249 computations (33.2% saved)

prune6.0ms (0.5%)

Pruning

1 alts after pruning (1 fresh and 0 done)

PrunedKeptTotal
New37138
Fresh000
Picked101
Done000
Total38139
Accuracy
100.0%
Counts
39 → 1
Alt Table
Click to see full alt table
StatusAccuracyProgram
100.0%
(*.f64 (/.f64 x y) 1/2)
100.0%
(*.f64 (/.f64 x y) 1/2)
Compiler

Compiled 14 to 10 computations (28.6% saved)

localize27.0ms (2%)

Compiler

Compiled 14 to 7 computations (50% saved)

eval0.0ms (0%)

Compiler

Compiled 2 to 2 computations (0% saved)

prune4.0ms (0.3%)

Pruning

1 alts after pruning (0 fresh and 1 done)

PrunedKeptTotal
New000
Fresh000
Picked011
Done000
Total011
Accuracy
100.0%
Counts
1 → 1
Alt Table
Click to see full alt table
StatusAccuracyProgram
100.0%
(*.f64 (/.f64 x y) 1/2)
Compiler

Compiled 19 to 13 computations (31.6% saved)

regimes12.0ms (0.9%)

Accuracy

Total -0.2b remaining (-∞%)

Threshold costs -0.2b (-∞%)

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

4 calls:

4.0ms
x
3.0ms
(*.f64 y 2)
3.0ms
(/.f64 x (*.f64 y 2))
2.0ms
y
Results
AccuracySegmentsBranch
100.0%1x
100.0%1y
100.0%1(/.f64 x (*.f64 y 2))
100.0%1(*.f64 y 2)
Compiler

Compiled 30 to 21 computations (30% saved)

simplify4.0ms (0.3%)

Algorithm
egg-herbie
Rules
*-commutative
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
0713
1813
Stop Event
done
saturated
Calls
Call 1
Inputs
(*.f64 (/.f64 x y) 1/2)
Outputs
(*.f64 (/.f64 x y) 1/2)
Compiler

Compiled 7 to 5 computations (28.6% saved)

soundness64.0ms (4.8%)

Rules
1558×pow1
1444×add-exp-log
1444×log1p-expm1-u
1444×expm1-log1p-u
1414×add-log-exp
Iterations

Useful iterations: 0 (0.0ms)

IterNodesCost
0713
114713
2191213
Stop Event
node limit
Compiler

Compiled 24 to 14 computations (41.7% saved)

end0.0ms (0%)

preprocess23.0ms (1.8%)

Compiler

Compiled 28 to 20 computations (28.6% saved)

Profiling

Loading profile data...