Average Error: 30.3 → 15.0
Time: 4.6s
Precision: binary64
\[\]
\[\]
double code(double x) {
	return ((double) sqrt(((double) (((double) pow(x, 2.0)) + ((double) pow(x, 2.0))))));
}
double code(double x) {
	double VAR;
	if ((x <= 3.59275856969827e-310)) {
		VAR = ((double) sqrt(((double) (((double) pow(x, 2.0)) * 2.0))));
	} else {
		VAR = ((double) (((double) (((double) pow(((double) sqrt(2.0)), 0.6666666666666666)) * ((double) pow(x, 1.0)))) * ((double) cbrt(((double) sqrt(2.0))))));
	}
	return VAR;
}

Error

Bits error versus x

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 2 regimes
  2. if x < 3.59275856969827e-310

    1. Initial program 30.0

      \[\]
    2. Simplified30.0

      \[\leadsto \]

    if 3.59275856969827e-310 < x

    1. Initial program 30.6

      \[\]
    2. Simplified30.6

      \[\leadsto \]
    3. Taylor expanded around 0 5.7

      \[\leadsto \]
    4. Simplified0.4

      \[\leadsto \]
    5. Using strategy rm
    6. Applied add-cube-cbrt0.4

      \[\leadsto \]
    7. Applied associate-*l*0.4

      \[\leadsto \]
    8. Simplified0.4

      \[\leadsto \]
    9. Using strategy rm
    10. Applied associate-*r*0.4

      \[\leadsto \]
    11. Simplified0.4

      \[\leadsto \]
  3. Recombined 2 regimes into one program.
  4. Final simplification15.0

    \[\leadsto \]

Reproduce

herbie shell --seed 2020181 
(FPCore (x)
  :name "sqrt E"
  :precision binary64
  (sqrt (+ (pow x 2.0) (pow x 2.0))))