Average Error: 4.0 → 2.9
Time: 14.2s
Precision: binary64
\[\]
↓
\[\]
↓
double code(double kx, double ky, double th) {
	return ((double) (((double) (((double) sin(ky)) / ((double) sqrt(((double) (((double) pow(((double) sin(kx)), 2.0)) + ((double) pow(((double) sin(ky)), 2.0)))))))) * ((double) sin(th))));
}

↓
double code(double kx, double ky, double th) {
	double VAR;
	if ((((double) (((double) sin(ky)) / ((double) sqrt(((double) (((double) pow(((double) sin(kx)), 2.0)) + ((double) pow(((double) sin(ky)), 2.0)))))))) <= 1.0)) {
		VAR = ((double) (((double) (((double) sin(ky)) / ((double) sqrt(((double) (((double) pow(((double) sin(kx)), 2.0)) + ((double) pow(((double) sin(ky)), 2.0)))))))) * ((double) sin(th))));
	} else {
		VAR = ((double) (((double) sin(th)) * ((double) (((double) sin(ky)) / ((double) (ky + ((double) (((double) (kx * ((double) (kx * ((double) (ky * 0.08333333333333333)))))) - ((double) (0.16666666666666666 * ((double) pow(ky, 3.0))))))))))));
	}
	return VAR;
}

Error
Bits error versus kx
Bits error versus ky
Bits error versus th
Try it out
Your Program's Arguments
Results
Enter valid numbers for all inputs
Derivation
  1. Split input into 2 regimes
  2. if  (/ (sin ky) (sqrt (+ (pow (sin kx) 2.0) (pow (sin ky) 2.0))))  < 1
    1. Initial program 2.2
      \[\]
    if 1 <  (/ (sin ky) (sqrt (+ (pow (sin kx) 2.0) (pow (sin ky) 2.0)))) 
    1. Initial program 62.3
      \[\]
    2. Taylor expanded around 0 24.8
      \[\leadsto \]
    3. Simplified24.8
      \[\leadsto \]
  3. Recombined 2 regimes into one program.
  4. Final simplification2.9
    \[\leadsto \]
Reproduce
herbie shell --seed 2020181 
(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)))