Details

prune3.3m (39.2%)

sample3.3m (39.1%)

Algorithm
78×intervals
Results
1.9m23266×body10240exit
33.7s44024×body1280valid
11.1s30814×body640valid
6.8s202184×body80valid
4.6s27677×body320valid
4.4s152103×body80nan
2.8s35629×body160valid
2.5s160053×pre80true
763.0ms3187×body2560valid
347.0ms7694×body160nan
344.0ms4104×pre10240exit
91.0ms8222×pre80false
82.0ms266×body5120valid

series46.7s (9.2%)

Calls

507 calls:

876.0ms
(/ (- (pow (cos delta) 2) (+ (* (pow (sin phi1) 3) (* (cos phi1) (* (cos delta) (* (sin delta) (cos theta))))) (+ (* (pow (sin phi1) 2) (* (pow (cos phi1) 2) (* (pow (cos theta) 2) (pow (sin delta) 2)))) (+ (* (pow (sin phi1) 4) (pow (cos delta) 2)) (* (pow (sin phi1) 3) (* (cos phi1) (* (sin delta) (* (cos delta) (cos theta))))))))) (+ (cos delta) (+ (* (pow (sin phi1) 2) (cos delta)) (* (sin phi1) (* (cos phi1) (* (cos theta) (sin delta)))))))
863.0ms
(/ (- (* (cos delta) (cos delta)) (* (* (sin phi1) (sin (asin (+ (* (sin phi1) (cos delta)) (* (* (cos phi1) (sin delta)) (cos theta)))))) (* (sin phi1) (sin (asin (+ (* (sin phi1) (cos delta)) (* (* (cos phi1) (sin delta)) (cos theta)))))))) (+ (cos delta) (* (sin phi1) (sin (asin (+ (* (sin phi1) (cos delta)) (* (* (cos phi1) (sin delta)) (cos theta))))))))
814.0ms
(/ (* (pow (/ 1.0 (+ 1.0 (exp (- s)))) c_p) (pow (- 1.0 (/ 1.0 (+ 1.0 (exp (- s))))) c_n)) (* (pow (/ 1.0 (+ 1.0 (exp (- t)))) c_p) (pow (- 1.0 (/ 1.0 (+ 1.0 (exp (- t))))) c_n)))
552.0ms
(/ (- (pow (/ x (+ x 1.0)) 3) (pow (/ (+ x 1.0) (- x 1.0)) 3)) (+ (* (/ (+ x 1.0) (- x 1.0)) (+ (/ (+ x 1.0) (- x 1.0)) (/ x (+ x 1.0)))) (* (/ x (+ x 1.0)) (/ x (+ x 1.0)))))
548.0ms
(/ (- (pow (/ x (+ x 1.0)) 3) (log (exp (pow (/ (+ x 1.0) (- x 1.0)) 3)))) (+ (* (/ (+ x 1.0) (- x 1.0)) (+ (/ (+ x 1.0) (- x 1.0)) (/ x (+ x 1.0)))) (* (/ x (+ x 1.0)) (/ x (+ x 1.0)))))

rewrite28.9s (5.7%)

Algorithm
144×rewrite-expression-head
Calls

507 calls:

2.0s
(/ (- (pow (cos delta) 2) (+ (* (pow (sin phi1) 3) (* (cos phi1) (* (cos delta) (* (sin delta) (cos theta))))) (+ (* (pow (sin phi1) 2) (* (pow (cos phi1) 2) (* (pow (cos theta) 2) (pow (sin delta) 2)))) (+ (* (pow (sin phi1) 4) (pow (cos delta) 2)) (* (pow (sin phi1) 3) (* (cos phi1) (* (sin delta) (* (cos delta) (cos theta))))))))) (+ (cos delta) (+ (* (pow (sin phi1) 2) (cos delta)) (* (sin phi1) (* (cos phi1) (* (cos theta) (sin delta)))))))
1.1s
(/ (- (pow (/ x (+ x 1.0)) 3) (log (exp (pow (/ (+ x 1.0) (- x 1.0)) 3)))) (+ (* (/ (+ x 1.0) (- x 1.0)) (+ (/ (+ x 1.0) (- x 1.0)) (/ x (+ x 1.0)))) (* (/ x (+ x 1.0)) (/ x (+ x 1.0)))))
802.0ms
(sqrt (+ (* (* (exp (log (cbrt (* (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2.0))) (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2.0))))))) (cbrt (* (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2.0))) (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2.0)))))) (cbrt (* (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2.0))) (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2.0)))))) (* (- phi1 phi2) (- phi1 phi2))))
771.0ms
(/ (- (pow (/ x (+ x 1.0)) 3) (pow (/ (+ x 1.0) (- x 1.0)) 3)) (+ (* (/ (+ x 1.0) (- x 1.0)) (+ (/ (+ x 1.0) (- x 1.0)) (/ x (+ x 1.0)))) (* (/ x (+ x 1.0)) (/ x (+ x 1.0)))))
727.0ms
(sqrt (+ (* (* (cbrt (* (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2.0))) (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2.0))))) (cbrt (* (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2.0))) (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2.0)))))) (cbrt (* (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2.0))) (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2.0)))))) (* (- phi1 phi2) (- phi1 phi2))))
Rules
60354×associate-*l/
41172×frac-times
32673×flip--
32672×flip3--
31241×cbrt-div
12246×associate-*r/
11960×sqrt-div
11933×frac-add
7761×times-frac
7146×*-un-lft-identity
6733×add-sqr-sqrt
5235×log-div
5196×exp-diff
4054×add-cube-cbrt
2795×add-exp-log
1837×add-cbrt-cube
1318×sqrt-prod
1176×pow1
1095×associate-*r*
1021×associate-/l* unpow-prod-down
1001×difference-of-squares
968×distribute-lft-out--
870×prod-exp
863×div-exp
816×add-log-exp
650×cube-prod
604×associate-/r*
513×distribute-lft-out
509×associate-*l*
436×cbrt-undiv
426×cbrt-unprod
383×cbrt-prod
368×div-inv
347×associate-/r/
331×pow-exp
285×distribute-rgt-neg-in
283×pow-prod-down
274×sqr-pow
242×pow-to-exp
241×unswap-sqr
203×flip-+
202×flip3-+
197×log-prod
140×associate-/l/
126×clear-num frac-2neg
105×sub-neg
103×*-commutative
98×unpow2
95×diff-log
77×pow-unpow
71×sum-log
70×pow1/3
61×div-sub
58×exp-prod sqrt-pow1
56×unpow3 cube-mult
54×pow1/2 distribute-lft-neg-in swap-sqr
53×log-pow
48×+-commutative
45×pow-pow sin-mult
42×1-exp rec-exp
41×rem-sqrt-square
40×exp-sum
36×pow-prod-up
35×neg-sub0
34×pow-sqr frac-sub
31×distribute-rgt-in distribute-rgt-out distribute-lft-in
30×sin-cos-mult
27×associate-+l-
25×tan-quot
24×associate-+r+
23×neg-log associate--l+
20×rem-log-exp associate--l-
18×pow-plus neg-mul-1
17×sin-diff
16×cos-diff
15×cos-sum
14×rem-cube-cbrt pow2
13×cos-mult associate--r+
11×associate-+l+
10×cube-div acos-asin
distribute-frac-neg
rem-exp-log distribute-lft1-in pow-flip associate--r-
sum-cubes associate-+r-
difference-cubes pow-sub
sin-sum distribute-rgt-neg-out
inv-pow
acos-neg un-div-inv unpow-prod-up
pow-div sqrt-pow2 sin-neg sqrt-unprod sub-div sqrt-undiv
e-exp-1 exp-neg unpow1/2 rem-square-sqrt tanh-undef asin-acos tanh-def rem-cbrt-cube

simplify22.8s (4.5%)

regimes9.1s (1.8%)

Accuracy

Total 20.5b remaining (17.8%)

Threshold costs 6.5b (5.7%)

4.3b83.8%Quadratic roots, full range
3.9b71.6%Equirectangular approximation to distance on a great circle
3.4b70%Quotient of products
1.2b94.6%Cubic critical
1.1b0%2-ancestry mixing, zero discriminant

localize2.4s (0.5%)

bsearch549.0ms (0.1%)

end1.0ms (0.0%)