Average Error: 43.2 → 43.2
Time: 1.9s
Precision: binary64
\[\frac{p \cdot q}{\sqrt{p \cdot p + \left(q - r\right) \cdot \left(\left(q - r\right) + \sqrt{p \cdot p + \left(q - r\right) \cdot \left(q - r\right)}\right)}}\]
↓
\[\frac{p \cdot q}{\sqrt{p \cdot p + \left(q - r\right) \cdot \left(\left(q - r\right) + \sqrt{p \cdot p + \left(q - r\right) \cdot \left(q - r\right)}\right)}}\]
\frac{p \cdot q}{\sqrt{p \cdot p + \left(q - r\right) \cdot \left(\left(q - r\right) + \sqrt{p \cdot p + \left(q - r\right) \cdot \left(q - r\right)}\right)}}
↓
\frac{p \cdot q}{\sqrt{p \cdot p + \left(q - r\right) \cdot \left(\left(q - r\right) + \sqrt{p \cdot p + \left(q - r\right) \cdot \left(q - r\right)}\right)}}
double code(double p, double q, double r) {
	return ((double) (((double) (p * q)) / ((double) sqrt(((double) (((double) (p * p)) + ((double) (((double) (q - r)) * ((double) (((double) (q - r)) + ((double) sqrt(((double) (((double) (p * p)) + ((double) (((double) (q - r)) * ((double) (q - r))))))))))))))))));
}
↓
double code(double p, double q, double r) {
	return ((double) (((double) (p * q)) / ((double) sqrt(((double) (((double) (p * p)) + ((double) (((double) (q - r)) * ((double) (((double) (q - r)) + ((double) sqrt(((double) (((double) (p * p)) + ((double) (((double) (q - r)) * ((double) (q - r))))))))))))))))));
}

Error

Bits error versus p

Bits error versus q

Bits error versus r

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 43.2

    \[\frac{p \cdot q}{\sqrt{p \cdot p + \left(q - r\right) \cdot \left(\left(q - r\right) + \sqrt{p \cdot p + \left(q - r\right) \cdot \left(q - r\right)}\right)}}\]

  2. Final simplification43.2

    \[\leadsto \frac{p \cdot q}{\sqrt{p \cdot p + \left(q - r\right) \cdot \left(\left(q - r\right) + \sqrt{p \cdot p + \left(q - r\right) \cdot \left(q - r\right)}\right)}}\]

Reproduce

herbie shell --seed 2020153 
(FPCore (p q r)
  :name "(/ (* p q) (sqrt (+ (* p p) (* (- q r) (+ (- q r) (sqrt (+ (* p p) (* (- q r) (- q r)))))))))"
  :precision binary64
  (/ (* p q) (sqrt (+ (* p p) (* (- q r) (+ (- q r) (sqrt (+ (* p p) (* (- q r) (- q r))))))))))