


Bits error versus J



Bits error versus K



Bits error versus U
Results
if J < 2.205439540951895e-236 or 7.706939204047402e-195 < J Initial program 16.4
Simplified7.1
rmApplied associate-*r*7.0
if 2.205439540951895e-236 < J < 7.706939204047402e-195Initial program 34.6
Simplified21.2
rmApplied add-sqr-sqrt21.4
Applied associate-*l*21.4
rmApplied associate-*r*21.4
Taylor expanded around inf 37.0
Simplified37.0
Final simplification8.0
herbie shell --seed 2019005 +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: 23.0s)Debug log
1 calls. Slowest were:
| 272.0ms | (* (* (* -2 J) (cos (/ K 2))) (sqrt (+ 1 (pow (/ U (* (* 2 J) (cos (/ K 2)))) 2)))) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 9.7b
Found 3 expressions with local error:
| 9.6b | (* (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))) |
3 calls. Slowest were:
| 11.0ms | (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (* (cos (/ K 2)) (* -2 J))) |
| 5.0ms | (/ (/ (/ U 2) J) (cos (/ K 2))) |
| 4.0ms | (* (cos (/ K 2)) (* -2 J)) |
3 calls. Slowest were:
| 208.0ms | (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (* (cos (/ K 2)) (* -2 J))) |
| 80.0ms | (/ (/ (/ U 2) J) (cos (/ K 2))) |
| 39.0ms | (* (cos (/ K 2)) (* -2 J)) |
56 calls. Slowest were:
| 488.0ms | (/ (* (* (/ (/ U 2) J) (/ (/ U 2) J)) (/ (/ U 2) J)) (* (* (cos (/ K 2)) (cos (/ K 2))) (cos (/ K 2)))) |
| 191.0ms | (* (* (* (cos (/ K 2)) (cos (/ K 2))) (cos (/ K 2))) (* (* (* -2 J) (* -2 J)) (* -2 J))) |
| 179.0ms | (- (log (/ (/ U 2) J)) (log (cos (/ K 2)))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 5.4b
Found 3 expressions with local error:
| 9.7b | (* (* (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))) |
3 calls. Slowest were:
| 11.0ms | (* (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (cos (/ K 2))) (* -2 J)) |
| 6.0ms | (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (cos (/ K 2))) |
| 5.0ms | (/ (/ (/ U 2) J) (cos (/ K 2))) |
3 calls. Slowest were:
| 254.0ms | (* (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (cos (/ K 2))) (* -2 J)) |
| 155.0ms | (* (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))) (cos (/ K 2))) |
| 80.0ms | (/ (/ (/ U 2) J) (cos (/ K 2))) |
57 calls. Slowest were:
| 616.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))) |
| 465.0ms | (/ (* (* (/ (/ U 2) J) (/ (/ U 2) J)) (/ (/ U 2) J)) (* (* (cos (/ K 2)) (cos (/ K 2))) (cos (/ K 2)))) |
| 222.0ms | (+ (* 1/16 (/ (* (pow K 2) U) J)) (* 1/2 (/ U J))) |
10 alts after pruning (9 fresh and 1 done)
Merged error: 5.4b
Found 4 expressions with local error:
| 9.7b | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
| 9.7b | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
| 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)))) |
| 0.1b | (* (cos (/ K 2)) (* -2 J)) |
4 calls. Slowest were:
| 49.0ms | (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (cos (/ K 2)) (* -2 J)))) |
| 6.0ms | (* (cos (/ K 2)) (* -2 J)) |
| 2.0ms | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
4 calls. Slowest were:
| 239.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))))) |
| 160.0ms | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
| 55.0ms | (* (cos (/ K 2)) (* -2 J)) |
40 calls. Slowest were:
| 252.0ms | (* (* (* (cos (/ K 2)) (cos (/ K 2))) (cos (/ K 2))) (* (* (* -2 J) (* -2 J)) (* -2 J))) |
| 162.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))))) |
| 149.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 (10 fresh and 1 done)
Merged error: 5.4b
Found 4 expressions with local error:
| 9.7b | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
| 9.7b | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
| 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))) |
| 0.1b | (* (* (cos (/ K 2)) -2) J) |
4 calls. Slowest were:
| 35.0ms | (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (* (cos (/ K 2)) -2) J))) |
| 7.0ms | (* (* (cos (/ K 2)) -2) J) |
| 1.0ms | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
4 calls. Slowest were:
| 285.0ms | (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (* (cos (/ K 2)) -2) J))) |
| 189.0ms | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
| 169.0ms | (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) |
| 47.0ms | (* (* (cos (/ K 2)) -2) J) |
33 calls. Slowest were:
| 168.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)))) |
| 117.0ms | (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (* (* (cos (/ K 2)) -2) J))) |
| 113.0ms | (* (* (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))))) (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))))) (* (* (* (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2))))) (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))))) (sqrt (hypot 1 (/ (/ (/ U 2) J) (cos (/ K 2)))))) (* (* (* (* (cos (/ K 2)) -2) J) (* (* (cos (/ K 2)) -2) J)) (* (* (cos (/ K 2)) -2) J)))) |
11 alts after pruning (9 fresh and 2 done)
Merged error: 5.4b
-13.3% (4.2b remaining)
Error of 8.0b against oracle of 3.8b and baseline of 7.5b