


Bits error versus R



Bits error versus lambda1



Bits error versus lambda2



Bits error versus phi1



Bits error versus phi2
Initial program 16.9
Simplified16.9
rmApplied cos-diff3.9
rmApplied add-log-exp3.9
rmApplied add-exp-log3.9
Final simplification3.9
herbie shell --seed 2019005 +o rules:numerics
(FPCore (R lambda1 lambda2 phi1 phi2)
:name "Spherical law of cosines"
(* (acos (+ (* (sin phi1) (sin phi2)) (* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2))))) R))
Time bar (total: 53.3s)Debug log
1 calls. Slowest were:
| 50.0ms | (* (acos (+ (* (sin phi1) (sin phi2)) (* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2))))) R) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 15.8b
Found 4 expressions with local error:
| 2.4b | (cos (- lambda1 lambda2)) |
| 0.5b | (acos (fma (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)) (* (sin phi2) (sin phi1)))) |
| 0.3b | (fma (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)) (* (sin phi2) (sin phi1))) |
| 0.3b | (* R (acos (fma (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)) (* (sin phi2) (sin phi1))))) |
4 calls. Slowest were:
| 6.0ms | (cos (- lambda1 lambda2)) |
| 2.0ms | (* R (acos (fma (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)) (* (sin phi2) (sin phi1))))) |
| 1.0ms | (acos (fma (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)) (* (sin phi2) (sin phi1)))) |
4 calls. Slowest were:
| 210.0ms | (fma (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)) (* (sin phi2) (sin phi1))) |
| 35.0ms | (cos (- lambda1 lambda2)) |
| 25.0ms | (* R (acos (fma (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)) (* (sin phi2) (sin phi1))))) |
| 18.0ms | (acos (fma (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)) (* (sin phi2) (sin phi1)))) |
18 calls. Slowest were:
| 83.0ms | (- (+ 1 (* lambda2 lambda1)) (* 1/2 (pow lambda1 2))) |
| 82.0ms | (- (+ (* phi1 phi2) 1) (* 1/2 (pow phi1 2))) |
| 44.0ms | (* R (acos (fma (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)) (* (sin phi1) (sin phi2))))) |
15 alts after pruning (15 fresh and 0 done)
Merged error: 3.3b
Found 4 expressions with local error:
| 0.4b | (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1)))) |
| 0.2b | (* R (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1))))) |
| 0.2b | (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1))) |
| 0.2b | (* (sin lambda1) (sin lambda2)) |
4 calls. Slowest were:
| 7.0ms | (* (sin lambda1) (sin lambda2)) |
| 3.0ms | (* R (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1))))) |
| 1.0ms | (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1)))) |
4 calls. Slowest were:
| 302.0ms | (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1))) |
| 62.0ms | (* (sin lambda1) (sin lambda2)) |
| 49.0ms | (* R (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1))))) |
| 20.0ms | (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1)))) |
24 calls. Slowest were:
| 283.0ms | (* R (acos (fma (* (cos phi1) (cos phi2)) (+ (* (sin lambda2) (sin lambda1)) (* (cos lambda1) (cos lambda2))) (* (sin phi1) (sin phi2))))) |
| 281.0ms | (acos (fma (* (cos phi1) (cos phi2)) (+ (* (sin lambda2) (sin lambda1)) (* (cos lambda1) (cos lambda2))) (* (sin phi1) (sin phi2)))) |
| 277.0ms | (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1)))) |
15 alts after pruning (15 fresh and 0 done)
Merged error: 3.2b
Found 4 expressions with local error:
| 0.4b | (acos (+ (* (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2)))) (* (sin phi2) (sin phi1)))) |
| 0.2b | (* R (acos (+ (* (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2)))) (* (sin phi2) (sin phi1))))) |
| 0.2b | (* (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2)))) |
| 0.2b | (* (sin lambda1) (sin lambda2)) |
4 calls. Slowest were:
| 27.0ms | (* (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2)))) |
| 4.0ms | (* (sin lambda1) (sin lambda2)) |
| 2.0ms | (acos (+ (* (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2)))) (* (sin phi2) (sin phi1)))) |
4 calls. Slowest were:
| 224.0ms | (* (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2)))) |
| 73.0ms | (* (sin lambda1) (sin lambda2)) |
| 42.0ms | (* R (acos (+ (* (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2)))) (* (sin phi2) (sin phi1))))) |
| 26.0ms | (acos (+ (* (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2)))) (* (sin phi2) (sin phi1)))) |
46 calls. Slowest were:
| 480.0ms | (* (+ (cos (+ phi1 phi2)) (cos (- phi1 phi2))) (+ (* (+ (cos (+ lambda1 lambda2)) (cos (- lambda1 lambda2))) 2) (* 2 (- (cos (- lambda1 lambda2)) (cos (+ lambda1 lambda2)))))) |
| 420.0ms | (* 2 (+ (* (* (cos lambda1) (cos lambda2)) (* (cos lambda1) (cos lambda2))) (- (* (* (sin lambda1) (sin lambda2)) (* (sin lambda1) (sin lambda2))) (* (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2)))))) |
| 323.0ms | (cbrt (acos (+ (* (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2)))) (* (sin phi2) (sin phi1))))) |
14 alts after pruning (14 fresh and 0 done)
Merged error: 3.2b
Found 4 expressions with local error:
| 0.4b | (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1)))) |
| 0.4b | (log (exp (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1)))))) |
| 0.3b | (exp (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1))))) |
| 0.2b | (* R (log (exp (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1))))))) |
4 calls. Slowest were:
| 7.0ms | (* R (log (exp (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1))))))) |
| 3.0ms | (log (exp (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1)))))) |
| 3.0ms | (exp (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1))))) |
4 calls. Slowest were:
| 66.0ms | (exp (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1))))) |
| 28.0ms | (* R (log (exp (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1))))))) |
| 22.0ms | (log (exp (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1)))))) |
| 21.0ms | (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1)))) |
38 calls. Slowest were:
| 571.0ms | (log (exp (sqrt (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1))))))) |
| 333.0ms | (acos (fma (* (cos phi1) (cos phi2)) (+ (* (cos lambda1) (cos lambda2)) (* (sin lambda1) (sin lambda2))) (* (sin phi2) (sin phi1)))) |
| 285.0ms | (acos (fma (* (cos phi1) (cos phi2)) (+ (* (sin lambda2) (sin lambda1)) (* (cos lambda1) (cos lambda2))) (* (sin phi1) (sin phi2)))) |
17 alts after pruning (17 fresh and 0 done)
Merged error: 3.2b
0% (0.4b remaining)
Error of 3.9b against oracle of 3.5b and baseline of 3.9b