Average Error: 2.7 → 1.6
Time: 3.5s
Precision: binary64
\[\]
\[\]
double code(double x, double y, double z) {
	return (x * (sin(y) / y)) / z;
}
double code(double x, double y, double z) {
	double tmp;
	if (x <= -4.4107549280265686e+38) {
		tmp = ((x * sin(y)) / y) / z;
	} else {
		tmp = x / (z / (sin(y) / y));
	}
	return tmp;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Target

Original2.7
Target0.3
Herbie1.6
\[\]

Derivation

  1. Split input into 2 regimes
  2. if x < -4.4107549280265686e38

    1. Initial program 0.2

      \[\]
    2. Using strategy rm
    3. Applied associate-*r/0.3

      \[\leadsto \]

    if -4.4107549280265686e38 < x

    1. Initial program 3.3

      \[\]
    2. Using strategy rm
    3. Applied associate-/l*1.9

      \[\leadsto \]
  3. Recombined 2 regimes into one program.
  4. Final simplification1.6

    \[\leadsto \]

Reproduce

herbie shell --seed 2020338 
(FPCore (x y z)
  :name "Linear.Quaternion:$ctanh from linear-1.19.1.3"
  :precision binary64

  :herbie-target
  (if (< z -4.2173720203427147e-29) (/ (* x (/ 1.0 (/ y (sin y)))) z) (if (< z 4.446702369113811e+64) (/ x (* z (/ y (sin y)))) (/ (* x (/ 1.0 (/ y (sin y)))) z)))

  (/ (* x (/ (sin y) y)) z))