


Bits error versus J



Bits error versus K



Bits error versus U
Results
if J < -7.384923856485112e-264 or 7.035784431961942e-228 < J Initial program 14.1
Initial simplification5.3
rmApplied associate-/l/5.3
if -7.384923856485112e-264 < J < 7.035784431961942e-228Initial program 40.9
Initial simplification27.6
rmApplied add-sqr-sqrt27.8
Applied associate-*l*27.7
Taylor expanded around inf 33.5
Simplified33.5
Final simplification8.1
herbie shell --seed 2018362 +o rules:numerics
(FPCore (J K U)
:name "Maksimov and Kolovsky, Equation (3)"
(* (* (* -2 J) (cos (/ K 2))) (sqrt (+ 1 (pow (/ U (* (* 2 J) (cos (/ K 2)))) 2)))))
Time bar (total: 21.5s)Debug log
2 alts after pruning (2 fresh and 0 done)
Merged error: 9.2b
Found 4 expressions with local error:
| 9.2b | (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (* (cos (/ K 2)) (* -2 J))) |
| 0.1b | (* (cos (/ K 2)) (* -2 J)) |
| 0.1b | (/ (/ (/ U 2) J) (cos (/ K 2))) |
| 0.0b | (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) |
4 calls. Slowest were:
| 18.0ms | (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (* (cos (/ K 2)) (* -2 J))) |
| 6.0ms | (/ (/ (/ U 2) J) (cos (/ K 2))) |
| 5.0ms | (* (cos (/ K 2)) (* -2 J)) |
4 calls. Slowest were:
| 270.0ms | (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (* (cos (/ K 2)) (* -2 J))) |
| 192.0ms | (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) |
| 99.0ms | (/ (/ (/ U 2) J) (cos (/ K 2))) |
| 59.0ms | (* (cos (/ K 2)) (* -2 J)) |
59 calls. Slowest were:
| 611.0ms | (/ (* (* (/ (/ U 2) J) (/ (/ U 2) J)) (/ (/ U 2) J)) (* (* (cos (/ K 2)) (cos (/ K 2))) (cos (/ K 2)))) |
| 219.0ms | (- (log (/ (/ U 2) J)) (log (cos (/ K 2)))) |
| 209.0ms | (* (* (* (cos (/ K 2)) (cos (/ K 2))) (cos (/ K 2))) (* (* (* -2 J) (* -2 J)) (* -2 J))) |
11 alts after pruning (11 fresh and 0 done)
Merged error: 5.1b
Found 4 expressions with local error:
| 9.2b | (* (hypot 1 (/ (/ U 2) (* (cos (/ K 2)) J))) (* (cos (/ K 2)) (* -2 J))) |
| 0.1b | (* (cos (/ K 2)) J) |
| 0.1b | (* (cos (/ K 2)) (* -2 J)) |
| 0.1b | (/ (/ U 2) (* (cos (/ K 2)) J)) |
4 calls. Slowest were:
| 10.0ms | (* (hypot 1 (/ (/ U 2) (* (cos (/ K 2)) J))) (* (cos (/ K 2)) (* -2 J))) |
| 9.0ms | (/ (/ U 2) (* (cos (/ K 2)) J)) |
| 7.0ms | (* (cos (/ K 2)) (* -2 J)) |
4 calls. Slowest were:
| 222.0ms | (* (hypot 1 (/ (/ U 2) (* (cos (/ K 2)) J))) (* (cos (/ K 2)) (* -2 J))) |
| 69.0ms | (/ (/ U 2) (* (cos (/ K 2)) J)) |
| 55.0ms | (* (cos (/ K 2)) (* -2 J)) |
| 51.0ms | (* (cos (/ K 2)) J) |
40 calls. Slowest were:
| 418.0ms | (/ (* (* (/ U 2) (/ U 2)) (/ U 2)) (* (* (* (cos (/ K 2)) J) (* (cos (/ K 2)) J)) (* (cos (/ K 2)) J))) |
| 268.0ms | (+ (* 1/16 (/ (* (pow K 2) U) J)) (* 1/2 (/ U J))) |
| 226.0ms | (* (* (* (cos (/ K 2)) (cos (/ K 2))) (cos (/ K 2))) (* (* (* -2 J) (* -2 J)) (* -2 J))) |
12 alts after pruning (11 fresh and 1 done)
Merged error: 5.1b
Found 4 expressions with local error:
| 9.2b | (* (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (cos (/ K 2))) (* -2 J)) |
| 0.1b | (/ (/ (/ U 2) J) (cos (/ K 2))) |
| 0.1b | (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (cos (/ K 2))) |
| 0.0b | (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) |
4 calls. Slowest were:
| 12.0ms | (* (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (cos (/ K 2))) (* -2 J)) |
| 11.0ms | (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (cos (/ K 2))) |
| 11.0ms | (/ (/ (/ U 2) J) (cos (/ K 2))) |
4 calls. Slowest were:
| 259.0ms | (* (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (cos (/ K 2))) (* -2 J)) |
| 202.0ms | (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (cos (/ K 2))) |
| 182.0ms | (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) |
| 102.0ms | (/ (/ (/ U 2) J) (cos (/ K 2))) |
60 calls. Slowest were:
| 726.0ms | (* (* (* (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (cos (/ K 2))) (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (cos (/ K 2)))) (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (cos (/ K 2)))) (* (* (* -2 J) (* -2 J)) (* -2 J))) |
| 614.0ms | (/ (* (* (/ (/ U 2) J) (/ (/ U 2) J)) (/ (/ U 2) J)) (* (* (cos (/ K 2)) (cos (/ K 2))) (cos (/ K 2)))) |
| 199.0ms | (- (log (/ (/ U 2) J)) (log (cos (/ K 2)))) |
11 alts after pruning (10 fresh and 1 done)
Merged error: 5.1b
Found 4 expressions with local error:
| 9.2b | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
| 9.2b | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
| 0.2b | (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (cos (/ K 2)) (* -2 J))) |
| 0.2b | (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (cos (/ K 2)) (* -2 J)))) |
4 calls. Slowest were:
| 37.0ms | (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (cos (/ K 2)) (* -2 J)))) |
| 11.0ms | (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (cos (/ K 2)) (* -2 J))) |
| 1.0ms | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
4 calls. Slowest were:
| 243.0ms | (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (cos (/ K 2)) (* -2 J)))) |
| 196.0ms | (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (cos (/ K 2)) (* -2 J))) |
| 176.0ms | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
| 168.0ms | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
47 calls. Slowest were:
| 169.0ms | (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (cos (/ K 2)) (* -2 J)))) |
| 167.0ms | (+ (log (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))))) (log (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (cos (/ K 2)) (* -2 J))))) |
| 163.0ms | (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (cos (/ K 2)) (* -2 J)))) |
11 alts after pruning (9 fresh and 2 done)
Merged error: 5.1b
-16.1% (4.3b remaining)
Error of 8.1b against oracle of 3.8b and baseline of 7.5b