Average Error: 18.5 → 18.5
Time: 2.6s
Precision: binary64
\[\frac{1}{\sqrt{\left(a - b\right) \cdot \left(a - b\right) + r \cdot r}} - \frac{1}{\sqrt{\left(a + b\right) \cdot \left(a + b\right) + r \cdot r}}\]
↓
\[\frac{1}{\sqrt{\left(a - b\right) \cdot \left(a - b\right) + r \cdot r}} - \frac{1}{\sqrt{\left(a + b\right) \cdot \left(a + b\right) + r \cdot r}}\]
\frac{1}{\sqrt{\left(a - b\right) \cdot \left(a - b\right) + r \cdot r}} - \frac{1}{\sqrt{\left(a + b\right) \cdot \left(a + b\right) + r \cdot r}}
↓
\frac{1}{\sqrt{\left(a - b\right) \cdot \left(a - b\right) + r \cdot r}} - \frac{1}{\sqrt{\left(a + b\right) \cdot \left(a + b\right) + r \cdot r}}
double code(double a, double b, double r) {
	return ((double) (((double) (1.0 / ((double) sqrt(((double) (((double) (((double) (a - b)) * ((double) (a - b)))) + ((double) (r * r)))))))) - ((double) (1.0 / ((double) sqrt(((double) (((double) (((double) (a + b)) * ((double) (a + b)))) + ((double) (r * r))))))))));
}
↓
double code(double a, double b, double r) {
	return ((double) (((double) (1.0 / ((double) sqrt(((double) (((double) (((double) (a - b)) * ((double) (a - b)))) + ((double) (r * r)))))))) - ((double) (1.0 / ((double) sqrt(((double) (((double) (((double) (a + b)) * ((double) (a + b)))) + ((double) (r * r))))))))));
}

Error

Bits error versus a

Bits error versus b

Bits error versus r

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 18.5

    \[\frac{1}{\sqrt{\left(a - b\right) \cdot \left(a - b\right) + r \cdot r}} - \frac{1}{\sqrt{\left(a + b\right) \cdot \left(a + b\right) + r \cdot r}}\]

  2. Final simplification18.5

    \[\leadsto \frac{1}{\sqrt{\left(a - b\right) \cdot \left(a - b\right) + r \cdot r}} - \frac{1}{\sqrt{\left(a + b\right) \cdot \left(a + b\right) + r \cdot r}}\]

Reproduce

herbie shell --seed 2020153 
(FPCore (a b r)
  :name "(- (/ 1 (sqrt (+ (* (- a b) (- a b)) (* r r)))) (/ 1 (sqrt (+ (* (+ a b) (+ a b)) (* r r)))))"
  :precision binary64
  (- (/ 1.0 (sqrt (+ (* (- a b) (- a b)) (* r r)))) (/ 1.0 (sqrt (+ (* (+ a b) (+ a b)) (* r r))))))