Average Error: 32.7 → 29.3
Time: 1.2m
Precision: 64
Internal Precision: 128
\[\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell \cdot \ell}{Om}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}\]
\[\begin{array}{l} \mathbf{if}\;t \le 8.152146896984202 \cdot 10^{-243}:\\ \;\;\;\;\sqrt{\sqrt{\left(n \cdot 0\right) \cdot U + \left(U \cdot \left(n \cdot 2\right)\right) \cdot (\left(\sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}} \cdot \sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}}\right) \cdot \left(\sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}}\right) + \left(\left({\left(\frac{\ell}{Om}\right)}^{2} \cdot \left(-n\right)\right) \cdot \left(U - U*\right)\right))_*}} \cdot \sqrt{\sqrt{\left(U \cdot \left(n \cdot 2\right)\right) \cdot \left(\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - \left(U - U*\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right)\right)}}\\ \mathbf{else}:\\ \;\;\;\;\sqrt{\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - \left(U - U*\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right)} \cdot \sqrt{U \cdot \left(n \cdot 2\right)}\\ \end{array}\]

Error

Bits error versus n

Bits error versus U

Bits error versus t

Bits error versus l

Bits error versus Om

Bits error versus U*

Derivation

  1. Split input into 2 regimes
  2. if t < 8.152146896984202e-243

    1. Initial program 33.7

      \[\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell \cdot \ell}{Om}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}\]
    2. Using strategy rm
    3. Applied associate-/l*30.7

      \[\leadsto \sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \color{blue}{\frac{\ell}{\frac{Om}{\ell}}}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}\]
    4. Using strategy rm
    5. Applied add-sqr-sqrt30.9

      \[\leadsto \color{blue}{\sqrt{\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}} \cdot \sqrt{\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}}}\]
    6. Using strategy rm
    7. Applied add-cube-cbrt31.0

      \[\leadsto \sqrt{\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\color{blue}{\left(\sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}} \cdot \sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}}\right) \cdot \sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}}} - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}} \cdot \sqrt{\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}}\]
    8. Applied prod-diff31.0

      \[\leadsto \sqrt{\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \color{blue}{\left((\left(\sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}} \cdot \sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}}\right) \cdot \left(\sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}}\right) + \left(-\left(U - U*\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right)\right))_* + (\left(-\left(U - U*\right)\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) + \left(\left(U - U*\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right)\right))_*\right)}}} \cdot \sqrt{\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}}\]
    9. Applied distribute-rgt-in31.0

      \[\leadsto \sqrt{\sqrt{\color{blue}{(\left(\sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}} \cdot \sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}}\right) \cdot \left(\sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}}\right) + \left(-\left(U - U*\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right)\right))_* \cdot \left(\left(2 \cdot n\right) \cdot U\right) + (\left(-\left(U - U*\right)\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) + \left(\left(U - U*\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right)\right))_* \cdot \left(\left(2 \cdot n\right) \cdot U\right)}}} \cdot \sqrt{\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}}\]
    10. Simplified31.0

      \[\leadsto \sqrt{\sqrt{(\left(\sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}} \cdot \sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}}\right) \cdot \left(\sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}}\right) + \left(-\left(U - U*\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right)\right))_* \cdot \left(\left(2 \cdot n\right) \cdot U\right) + \color{blue}{\left(0 \cdot n\right) \cdot U}}} \cdot \sqrt{\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}}\]

    if 8.152146896984202e-243 < t

    1. Initial program 31.6

      \[\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell \cdot \ell}{Om}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}\]
    2. Using strategy rm
    3. Applied associate-/l*28.9

      \[\leadsto \sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \color{blue}{\frac{\ell}{\frac{Om}{\ell}}}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}\]
    4. Using strategy rm
    5. Applied sqrt-prod27.2

      \[\leadsto \color{blue}{\sqrt{\left(2 \cdot n\right) \cdot U} \cdot \sqrt{\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)}}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification29.3

    \[\leadsto \begin{array}{l} \mathbf{if}\;t \le 8.152146896984202 \cdot 10^{-243}:\\ \;\;\;\;\sqrt{\sqrt{\left(n \cdot 0\right) \cdot U + \left(U \cdot \left(n \cdot 2\right)\right) \cdot (\left(\sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}} \cdot \sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}}\right) \cdot \left(\sqrt[3]{t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}}\right) + \left(\left({\left(\frac{\ell}{Om}\right)}^{2} \cdot \left(-n\right)\right) \cdot \left(U - U*\right)\right))_*}} \cdot \sqrt{\sqrt{\left(U \cdot \left(n \cdot 2\right)\right) \cdot \left(\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - \left(U - U*\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right)\right)}}\\ \mathbf{else}:\\ \;\;\;\;\sqrt{\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - \left(U - U*\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right)} \cdot \sqrt{U \cdot \left(n \cdot 2\right)}\\ \end{array}\]

Reproduce

herbie shell --seed 2019010 +o rules:numerics
(FPCore (n U t l Om U*)
  :name "Toniolo and Linder, Equation (13)"
  (sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))))

Details

Time bar (total: 1.2m)Debug log

sample387.0ms

Algorithm
intervals

simplify132.0ms

Counts
1 → 1
Calls
1 calls:
Slowest
131.0ms
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))

prune28.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 31.7b

localize76.0ms

Local error

Found 4 expressions with local error:

18.1b
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
10.3b
(* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))
7.7b
(/ (* l l) Om)
7.3b
(* n (pow (/ l Om) 2))

rewrite84.0ms

Algorithm
rewrite-expression-head
Rules
20×associate-*r/
11×associate-*r*
10×flip--
10×flip3--
add-cube-cbrt
*-un-lft-identity
add-sqr-sqrt
frac-sub
distribute-lft-in
distribute-rgt-in
add-exp-log
prod-diff
sqrt-div
add-cbrt-cube
pow1
add-log-exp
log1p-expm1-u
associate--l+
unpow-prod-down
sub-neg
expm1-log1p-u
div-inv
*-commutative
unpow2
associate-/l*
cbrt-unprod
sqrt-prod
prod-exp
pow-prod-down
associate-*l*
pow1/2
frac-2neg
clear-num
Counts
4 → 85
Calls
4 calls:
Slowest
37.0ms
(* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))
36.0ms
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
4.0ms
(/ (* l l) Om)
3.0ms
(* n (pow (/ l Om) 2))

series1.5s

Counts
4 → 12
Calls
4 calls:
Slowest
1.2s
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
290.0ms
(* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))
28.0ms
(* n (pow (/ l Om) 2))
12.0ms
(/ (* l l) Om)

simplify7.8s

Counts
54 → 97
Calls
54 calls:
Slowest
614.0ms
(* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U))
588.0ms
(* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U))
524.0ms
(* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U))
358.0ms
(sqrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))
353.0ms
(cbrt (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))

prune2.0s

Pruning

11 alts after pruning (11 fresh and 0 done)

Merged error: 23.3b

localize18.0ms

Local error

Found 4 expressions with local error:

18.1b
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
10.3b
(* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))
7.3b
(* n (pow (/ l Om) 2))
6.2b
(* (* n (pow (/ l Om) 2)) (- U U*))

rewrite70.0ms

Algorithm
rewrite-expression-head
Rules
22×associate-*r/
14×associate-*r*
11×flip--
11×flip3--
10×add-cube-cbrt
10×*-un-lft-identity
10×add-sqr-sqrt
distribute-lft-in
distribute-rgt-in
frac-sub
add-exp-log
add-cbrt-cube
pow1
prod-diff
sqrt-div
sub-neg
add-log-exp
log1p-expm1-u
associate--l+
unpow-prod-down
expm1-log1p-u
*-commutative
cbrt-unprod
prod-exp
pow-prod-down
associate-*l*
unpow2
div-inv
sqrt-prod
pow1/2
Counts
4 → 95
Calls
4 calls:
Slowest
40.0ms
(* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))
12.0ms
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
11.0ms
(* (* n (pow (/ l Om) 2)) (- U U*))
4.0ms
(* n (pow (/ l Om) 2))

series1.6s

Counts
4 → 12
Calls
4 calls:
Slowest
1.2s
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
248.0ms
(* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))
95.0ms
(* (* n (pow (/ l Om) 2)) (- U U*))
27.0ms
(* n (pow (/ l Om) 2))

simplify10.9s

Counts
66 → 107
Calls
66 calls:
Slowest
570.0ms
(sqrt (* (* (* 2 n) U) (- (* (- (* t t) (* (* 2 (/ l (/ Om l))) (* 2 (/ l (/ Om l))))) (+ (* U U) (+ (* U* U*) (* U U*)))) (* (+ t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- (pow U 3) (pow U* 3)))))))
552.0ms
(* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U))
529.0ms
(* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U))
512.0ms
(* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U))
459.0ms
(+ (log (* (* 2 n) U)) (log (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))

prune2.0s

Pruning

9 alts after pruning (9 fresh and 0 done)

Merged error: 21.9b

localize17.0ms

Local error

Found 4 expressions with local error:

18.1b
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
18.1b
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
10.3b
(* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))
10.3b
(* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))

rewrite87.0ms

Algorithm
rewrite-expression-head
Rules
40×associate-*r/
20×flip--
20×flip3--
16×frac-sub
14×distribute-lft-in
14×distribute-rgt-in
12×prod-diff
12×sqrt-div
10×add-cube-cbrt
10×*-un-lft-identity
10×add-sqr-sqrt
add-exp-log
associate--l+
add-cbrt-cube
sub-neg
pow1
associate-*r*
add-log-exp
log1p-expm1-u
expm1-log1p-u
cbrt-unprod
sqrt-prod
*-commutative
prod-exp
pow-prod-down
associate-*l*
pow1/2
Counts
4 → 108
Calls
4 calls:
Slowest
22.0ms
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
21.0ms
(* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))
21.0ms
(* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))
20.0ms
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))

series2.9s

Counts
4 → 12
Calls
4 calls:
Slowest
1.2s
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
1.1s
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
326.0ms
(* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))
276.0ms
(* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*))))

simplify16.5s

Counts
80 → 120
Calls
80 calls:
Slowest
976.0ms
(* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U))
663.0ms
(* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U))
591.0ms
(* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U))
575.0ms
(* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U))
528.0ms
(* (fma (- (- U U*)) (* n (pow (/ l Om) 2)) (* (- U U*) (* n (pow (/ l Om) 2)))) (* (* 2 n) U))

prune3.1s

Pruning

11 alts after pruning (10 fresh and 1 done)

Merged error: 21.9b

localize56.0ms

Local error

Found 4 expressions with local error:

18.1b
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
18.1b
(sqrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U)))
10.5b
(cbrt (- t (* 2 (/ l (/ Om l)))))
10.5b
(cbrt (- t (* 2 (/ l (/ Om l)))))

rewrite25.0ms

Algorithm
rewrite-expression-head
Rules
10×associate-*r/
sqrt-div
add-sqr-sqrt
flip--
add-cube-cbrt
flip3--
*-un-lft-identity
cbrt-prod
add-log-exp
log1p-expm1-u
sqrt-prod
frac-sub
add-exp-log
add-cbrt-cube
cbrt-div
pow1
expm1-log1p-u
pow1/3
pow1/2
flip-+
flip3-+
rem-sqrt-square
Counts
4 → 63
Calls
4 calls:
Slowest
16.0ms
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
7.0ms
(sqrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U)))
1.0ms
(cbrt (- t (* 2 (/ l (/ Om l)))))
1.0ms
(cbrt (- t (* 2 (/ l (/ Om l)))))

series4.3s

Counts
4 → 12
Calls
4 calls:
Slowest
2.6s
(sqrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U)))
1.2s
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ l (/ Om l)))) (* (* n (pow (/ l Om) 2)) (- U U*)))))
239.0ms
(cbrt (- t (* 2 (/ l (/ Om l)))))
204.0ms
(cbrt (- t (* 2 (/ l (/ Om l)))))

simplify7.6s

Counts
35 → 75
Calls
35 calls:
Slowest
829.0ms
(sqrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U)))
637.0ms
(cbrt (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3)))
608.0ms
(cbrt (- (pow t 3) (pow (* 2 (/ l (/ Om l))) 3)))
535.0ms
(sqrt (sqrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U))))
500.0ms
(sqrt (* (cbrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U))) (cbrt (+ (* (fma (* (cbrt (- t (* 2 (/ l (/ Om l))))) (cbrt (- t (* 2 (/ l (/ Om l)))))) (cbrt (- t (* 2 (/ l (/ Om l))))) (- (* (- U U*) (* n (pow (/ l Om) 2))))) (* (* 2 n) U)) (* (* 0 n) U)))))

prune2.4s

Pruning

11 alts after pruning (9 fresh and 2 done)

Merged error: 21.9b

regimes920.0ms

Accuracy

9.3% (8.0b remaining)

Error of 29.3b against oracle of 21.3b and baseline of 30.1b

bsearch562.0ms

end0.0ms

sample6.9s

Algorithm
intervals