\frac{\sin ky}{\sqrt{{\sin kx}^{2} + {\sin ky}^{2}}} \cdot \sin th\frac{\sin ky}{\sqrt{{\sin kx}^{2} + {\sin ky}^{2}}} \cdot \sin th(FPCore (kx ky th) :precision binary64 (* (/ (sin ky) (sqrt (+ (pow (sin kx) 2.0) (pow (sin ky) 2.0)))) (sin th)))
(FPCore (kx ky th) :precision binary64 (* (/ (sin ky) (sqrt (+ (pow (sin kx) 2.0) (pow (sin ky) 2.0)))) (sin th)))
double code(double kx, double ky, double th) {
return (sin(ky) / sqrt(pow(sin(kx), 2.0) + pow(sin(ky), 2.0))) * sin(th);
}
double code(double kx, double ky, double th) {
return (sin(ky) / sqrt(pow(sin(kx), 2.0) + pow(sin(ky), 2.0))) * sin(th);
}















Bits error versus kx















Bits error versus ky















Bits error versus th
Results
| Alternative 1 | |
|---|---|
| Error | 18.0 |
| Cost | 40067 |
| Alternative 2 | |
|---|---|
| Error | 18.1 |
| Cost | 40067 |
| Alternative 3 | |
|---|---|
| Error | 17.0 |
| Cost | 45768 |
| Alternative 4 | |
|---|---|
| Error | 21.2 |
| Cost | 111622 |
| Alternative 5 | |
|---|---|
| Error | 28.6 |
| Cost | 39874 |
| Alternative 6 | |
|---|---|
| Error | 28.8 |
| Cost | 39560 |
| Alternative 7 | |
|---|---|
| Error | 28.9 |
| Cost | 39048 |
| Alternative 8 | |
|---|---|
| Error | 37.5 |
| Cost | 26241 |
| Alternative 9 | |
|---|---|
| Error | 38.9 |
| Cost | 13448 |
| Alternative 10 | |
|---|---|
| Error | 44.8 |
| Cost | 6792 |
| Alternative 11 | |
|---|---|
| Error | 56.0 |
| Cost | 64 |
| Alternative 12 | |
|---|---|
| Error | 60.3 |
| Cost | 64 |


Initial program 4.0
rmApplied *-un-lft-identity_binary64_784.0
Applied *-un-lft-identity_binary64_784.0
Applied times-frac_binary64_844.0
Applied associate-*l*_binary64_194.0
rmApplied *-un-lft-identity_binary64_784.0
Applied *-un-lft-identity_binary64_784.0
Applied times-frac_binary64_844.0
Applied associate-*l*_binary64_194.0
rmApplied *-un-lft-identity_binary64_784.0
Applied unpow-prod-down_binary64_1574.0
Simplified4.0
Simplified4.0
Final simplification4.0
herbie shell --seed 2021065
(FPCore (kx ky th)
:name "Toniolo and Linder, Equation (3b), real"
:precision binary64
(* (/ (sin ky) (sqrt (+ (pow (sin kx) 2.0) (pow (sin ky) 2.0)))) (sin th)))