cos2 (problem 3.4.1)

Average Error: 32.1 → 0.4

Time: 7.9s

Precision: binary64

Math

\[\frac{1 - \cos x}{x \cdot x}\]

↓

\[\begin{array}{l} \mathbf{if}\;x \leq -0.03140225698466425:\\ \;\;\;\;\frac{\sqrt{1 - \cos x}}{x} \cdot \frac{\sqrt{\sqrt{\log \left(e^{1 - \cos x}\right)}}}{\frac{x}{\sqrt{\sqrt{1 - \cos x}}}}\\ \mathbf{elif}\;x \leq 0.03002211543800249:\\ \;\;\;\;{x}^{4} \cdot 0.001388888888888889 + \left(0.5 + x \cdot \left(x \cdot -0.041666666666666664\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{{1}^{3} - {\left(\cos x\right)}^{3}}{x \cdot \left(x \cdot \left(1 \cdot 1 + \cos x \cdot \left(1 + \cos x\right)\right)\right)}\\ \end{array}\]

C

double code(double x) {
	return (((double) (1.0 - ((double) cos(x)))) / ((double) (x * x)));
}

↓

double code(double x) {
	double VAR;
	if ((x <= -0.03140225698466425)) {
		VAR = ((double) ((((double) sqrt(((double) (1.0 - ((double) cos(x)))))) / x) * (((double) sqrt(((double) sqrt(((double) log(((double) exp(((double) (1.0 - ((double) cos(x)))))))))))) / (x / ((double) sqrt(((double) sqrt(((double) (1.0 - ((double) cos(x))))))))))));
	} else {
		double VAR_1;
		if ((x <= 0.03002211543800249)) {
			VAR_1 = ((double) (((double) (((double) pow(x, 4.0)) * 0.001388888888888889)) + ((double) (0.5 + ((double) (x * ((double) (x * -0.041666666666666664))))))));
		} else {
			VAR_1 = (((double) (((double) pow(1.0, 3.0)) - ((double) pow(((double) cos(x)), 3.0)))) / ((double) (x * ((double) (x * ((double) (((double) (1.0 * 1.0)) + ((double) (((double) cos(x)) * ((double) (1.0 + ((double) cos(x)))))))))))));
		}
		VAR = VAR_1;
	}
	return VAR;
}

Error

Bits error versus x

Derivation

1. Split input into 3 regimes

2. if x < -0.0314022569846642494

   1. Initial program 1.0

      \[\frac{1 - \cos x}{x \cdot x}\]
   2. Using strategy rm

   3. Applied add-sqr-sqrt 1.1

      \[\leadsto \frac{\color{blue}{\sqrt{1 - \cos x} \cdot \sqrt{1 - \cos x}}}{x \cdot x}\]

   4. Applied times-frac 0.6

      \[\leadsto \color{blue}{\frac{\sqrt{1 - \cos x}}{x} \cdot \frac{\sqrt{1 - \cos x}}{x}}\]
   5. Using strategy rm

   6. Applied add-log-exp 0.6

      \[\leadsto \frac{\sqrt{1 - \cos x}}{x} \cdot \frac{\sqrt{1 - \color{blue}{\log \left(e^{\cos x}\right)}}}{x}\]

   7. Applied add-log-exp 0.6

      \[\leadsto \frac{\sqrt{1 - \cos x}}{x} \cdot \frac{\sqrt{\color{blue}{\log \left(e^{1}\right)} - \log \left(e^{\cos x}\right)}}{x}\]

   8. Applied diff-log 0.6

      \[\leadsto \frac{\sqrt{1 - \cos x}}{x} \cdot \frac{\sqrt{\color{blue}{\log \left(\frac{e^{1}}{e^{\cos x}}\right)}}}{x}\]

   9. Simplified 0.6

      \[\leadsto \frac{\sqrt{1 - \cos x}}{x} \cdot \frac{\sqrt{\log \color{blue}{\left(e^{1 - \cos x}\right)}}}{x}\]
   10. Using strategy rm

   11. Applied add-sqr-sqrt 0.6

       \[\leadsto \frac{\sqrt{1 - \cos x}}{x} \cdot \frac{\sqrt{\color{blue}{\sqrt{\log \left(e^{1 - \cos x}\right)} \cdot \sqrt{\log \left(e^{1 - \cos x}\right)}}}}{x}\]

   12. Applied sqrt-prod 0.7

       \[\leadsto \frac{\sqrt{1 - \cos x}}{x} \cdot \frac{\color{blue}{\sqrt{\sqrt{\log \left(e^{1 - \cos x}\right)}} \cdot \sqrt{\sqrt{\log \left(e^{1 - \cos x}\right)}}}}{x}\]

   13. Applied associate-/l* 0.7

       \[\leadsto \frac{\sqrt{1 - \cos x}}{x} \cdot \color{blue}{\frac{\sqrt{\sqrt{\log \left(e^{1 - \cos x}\right)}}}{\frac{x}{\sqrt{\sqrt{\log \left(e^{1 - \cos x}\right)}}}}}\]

   14. Simplified 0.7

       \[\leadsto \frac{\sqrt{1 - \cos x}}{x} \cdot \frac{\sqrt{\sqrt{\log \left(e^{1 - \cos x}\right)}}}{\color{blue}{\frac{x}{\sqrt{\sqrt{1 - \cos x}}}}}\]

   if -0.0314022569846642494 < x < 0.03002211543800249

   1. Initial program 62.3

      \[\frac{1 - \cos x}{x \cdot x}\]

   2. Taylor expanded around 0 0.0

      \[\leadsto \color{blue}{\left(0.001388888888888889 \cdot {x}^{4} + 0.5\right) - 0.041666666666666664 \cdot {x}^{2}}\]

   3. Simplified 0.0

      \[\leadsto \color{blue}{{x}^{4} \cdot 0.001388888888888889 + \left(0.5 + x \cdot \left(x \cdot -0.041666666666666664\right)\right)}\]

   if 0.03002211543800249 < x

   1. Initial program 1.1

      \[\frac{1 - \cos x}{x \cdot x}\]
   2. Using strategy rm

   3. Applied flip3-- 1.1

      \[\leadsto \frac{\color{blue}{\frac{{1}^{3} - {\left(\cos x\right)}^{3}}{1 \cdot 1 + \left(\cos x \cdot \cos x + 1 \cdot \cos x\right)}}}{x \cdot x}\]

   4. Applied associate-/l/ 1.1

      \[\leadsto \color{blue}{\frac{{1}^{3} - {\left(\cos x\right)}^{3}}{\left(x \cdot x\right) \cdot \left(1 \cdot 1 + \left(\cos x \cdot \cos x + 1 \cdot \cos x\right)\right)}}\]

   5. Simplified 1.1

      \[\leadsto \frac{{1}^{3} - {\left(\cos x\right)}^{3}}{\color{blue}{x \cdot \left(x \cdot \left(1 \cdot 1 + \cos x \cdot \left(1 + \cos x\right)\right)\right)}}\]
3. Recombined 3 regimes into one program.

4. Final simplification 0.4

   \[\leadsto \begin{array}{l} \mathbf{if}\;x \leq -0.03140225698466425:\\ \;\;\;\;\frac{\sqrt{1 - \cos x}}{x} \cdot \frac{\sqrt{\sqrt{\log \left(e^{1 - \cos x}\right)}}}{\frac{x}{\sqrt{\sqrt{1 - \cos x}}}}\\ \mathbf{elif}\;x \leq 0.03002211543800249:\\ \;\;\;\;{x}^{4} \cdot 0.001388888888888889 + \left(0.5 + x \cdot \left(x \cdot -0.041666666666666664\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{{1}^{3} - {\left(\cos x\right)}^{3}}{x \cdot \left(x \cdot \left(1 \cdot 1 + \cos x \cdot \left(1 + \cos x\right)\right)\right)}\\ \end{array}\]

Reproduce

herbie shell --seed 2020196 
(FPCore (x)
  :name "cos2 (problem 3.4.1)"
  :precision binary64
  (/ (- 1.0 (cos x)) (* x x)))