Details

sample19.9m (71.6%)

Algorithm
72×intervals
Results
7.5m80056×body10240exit
50.6s32843×body1280valid
26.2s198227×body80valid
22.3s24326×body640valid
20.7s152416×body80nan
14.5s24770×body320valid
11.3s35302×body160valid
9.8s153658×pre80true
3.7s3044×body2560valid
1.8s7939×body160nan
881.0ms4047×pre10240exit
287.0ms256×body5120valid

simplify3.4m (12.1%)

prune3.4m (12.1%)

series34.6s (2.1%)

Calls

399 calls:

Slowest
677.0ms
(/ (- (pow (log (+ (pow 1 3) (pow (exp x) 3))) 3) (pow (+ (* x y) (log (+ 1 (* (- (exp x) 1) (exp x))))) 3)) (+ (* (log (+ (pow 1 3) (pow (exp x) 3))) (log (+ (pow 1 3) (pow (exp x) 3)))) (+ (* (+ (* x y) (log (+ 1 (* (- (exp x) 1) (exp x))))) (+ (* x y) (log (+ 1 (* (- (exp x) 1) (exp x)))))) (* (log (+ (pow 1 3) (pow (exp x) 3))) (+ (* x y) (log (+ 1 (* (- (exp x) 1) (exp x)))))))))
607.0ms
(sqrt (+ (* (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2))) (* (- lambda1 lambda2) (cos (/ (+ phi1 phi2) 2)))) (* (- phi1 phi2) (- phi1 phi2))))
581.0ms
(sqrt (+ (* (* (- lambda1 lambda2) (- lambda1 lambda2)) (* (cos (* 1/2 (+ phi2 phi1))) (cbrt (* (* (cos (* 1/2 (+ phi2 phi1))) (cos (* 1/2 (+ phi2 phi1)))) (cos (* 1/2 (+ phi2 phi1))))))) (* (- phi1 phi2) (- phi1 phi2))))
578.0ms
(/ (* (pow (/ 1 (+ 1 (exp (- s)))) c_p) (pow (- 1 (/ 1 (+ 1 (exp (- s))))) c_n)) (* (pow (/ 1 (+ 1 (exp (- t)))) c_p) (pow (- 1 (/ 1 (+ 1 (exp (- t))))) c_n)))
557.0ms
(sqrt (+ (* (* (- lambda1 lambda2) (- lambda1 lambda2)) (* (cos (* 1/2 (+ phi2 phi1))) (cos (* 1/2 (+ phi2 phi1))))) (* (- phi1 phi2) (- phi1 phi2))))

regimes16.0s (1.0%)

Accuracy

Total 66.8b remaining (38.7%)

Threshold costs 4.8b (2.8%)

12.8b43%Equirectangular approximation to distance on a great circle
7.6b43%Cubic critical, narrow range
7.4b46.4%Quadratic roots, narrow range
6.9b0%Quadratic roots, medium range
5.9b11.6%Cubic critical, medium range

rewrite12.8s (0.8%)

Algorithm
128×rewrite-expression-head
Calls

399 calls:

Slowest
435.0ms
(* (* (+ (/ 1 (- x 1)) (/ x (+ x 1))) (+ (/ 1 (- x 1)) (/ x (+ x 1)))) (+ (/ 1 (- x 1)) (/ x (+ x 1))))
430.0ms
(sqrt (+ (* (* (- lambda1 lambda2) (- lambda1 lambda2)) (* (cos (* 1/2 (+ phi2 phi1))) (cbrt (* (* (cos (* 1/2 (+ phi2 phi1))) (cos (* 1/2 (+ phi2 phi1)))) (cos (* 1/2 (+ phi2 phi1))))))) (* (- phi1 phi2) (- phi1 phi2))))
398.0ms
(* (/ 1 3) (/ (- (sqrt (- (* b b) (* (* 3 a) c))) b) a))
335.0ms
(- (- (/ 1/8 (* (* x x) x)) (/ 1/2 x)) (/ 1/16 (pow x 5)))
216.0ms
(/ (- (pow (log (+ (pow 1 3) (pow (exp x) 3))) 3) (pow (+ (* x y) (log (+ 1 (* (- (exp x) 1) (exp x))))) 3)) (+ (* (log (+ (pow 1 3) (pow (exp x) 3))) (log (+ (pow 1 3) (pow (exp x) 3)))) (+ (* (+ (* x y) (log (+ 1 (* (- (exp x) 1) (exp x))))) (+ (* x y) (log (+ 1 (* (- (exp x) 1) (exp x)))))) (* (log (+ (pow 1 3) (pow (exp x) 3))) (+ (* x y) (log (+ 1 (* (- (exp x) 1) (exp x)))))))))
Rules
5170×times-frac
5087×*-un-lft-identity
5079×add-sqr-sqrt
3234×add-exp-log
2960×add-cube-cbrt
1742×add-cbrt-cube
1285×prod-exp
1215×pow1
1107×div-exp
966×difference-of-squares
939×flip--
938×flip3--
826×frac-times
718×associate-*l/
696×pow-exp
636×add-log-exp
580×associate-*l*
559×associate-/r*
552×sqrt-div
549×associate-*r*
540×frac-add
513×distribute-lft-out--
483×cbrt-unprod
476×associate-*r/
445×associate-/l*
411×cbrt-undiv
399×insert-posit16
387×cbrt-prod
364×distribute-lft-out
351×pow-to-exp
319×div-inv
315×log-prod
291×associate-/r/
281×pow-prod-down
279×difference-of-sqr-1
228×flip-+
227×sqrt-prod flip3-+
219×cos-mult
170×1-exp rec-exp
144×unswap-sqr
128×swap-sqr
119×unpow-prod-down
108×frac-2neg cbrt-div clear-num
102×*-commutative
101×sub-neg
95×associate-/l/
94×pow-prod-up
85×distribute-rgt-in distribute-lft-in
84×diff-log
76×associate-+r+
71×cube-prod
70×log-pow
64×pow-sqr
60×pow-plus tan-quot
52×frac-sub
50×sqr-pow
49×log-div
44×sum-log
39×pow2 associate--l+
35×rem-log-exp
33×pow1/3 distribute-rgt-neg-in
30×exp-prod
29×pow-flip +-commutative distribute-lft-neg-in
28×sqrt-pow1 pow1/2
26×rem-sqrt-square
24×div-sub
23×inv-pow
22×cos-sum
18×neg-mul-1
17×associate-+l+
15×associate--r+
12×distribute-rgt-out--
10×distribute-rgt-neg-out cos-diff associate--r-
associate--l- log-rec acos-asin distribute-rgt-out
exp-sum pow-div
associate-+r- pow-pow
sin-sum pow-unpow sin-mult rem-exp-log
exp-diff distribute-frac-neg
pow3 difference-cubes un-div-inv associate-+l-
neg-sub0 asin-acos acos-neg unsub-neg
cube-mult exp-to-pow distribute-neg-frac sqr-cos sin-diff
rem-cube-cbrt distribute-lft-neg-out distribute-lft1-in sqrt-unprod unpow3 tanh-undef tanh-def cube-unmult sum-cubes rem-cbrt-cube distribute-rgt1-in

localize5.2s (0.3%)

bsearch2.0s (0.1%)

end1.0ms (0.0%)