Average Error: 0.0 → 0.0
Time: 6.3s
Precision: binary64
Cost: 14592
\[\left(\frac{\sqrt{2}}{4} \cdot \sqrt{1 - 3 \cdot \left(v \cdot v\right)}\right) \cdot \left(1 - v \cdot v\right) \]
\[\sqrt{\frac{1 + {v}^{6} \cdot -27}{1 + \left(v \cdot v\right) \cdot \left(3 + \left(v \cdot v\right) \cdot 9\right)} \cdot 0.125} \cdot \left(1 - v \cdot v\right) \]
(FPCore (v)
 :precision binary64
 (* (* (/ (sqrt 2.0) 4.0) (sqrt (- 1.0 (* 3.0 (* v v))))) (- 1.0 (* v v))))
(FPCore (v)
 :precision binary64
 (*
  (sqrt
   (*
    (/
     (+ 1.0 (* (pow v 6.0) -27.0))
     (+ 1.0 (* (* v v) (+ 3.0 (* (* v v) 9.0)))))
    0.125))
  (- 1.0 (* v v))))
double code(double v) {
	return ((sqrt(2.0) / 4.0) * sqrt((1.0 - (3.0 * (v * v))))) * (1.0 - (v * v));
}
double code(double v) {
	return sqrt((((1.0 + (pow(v, 6.0) * -27.0)) / (1.0 + ((v * v) * (3.0 + ((v * v) * 9.0))))) * 0.125)) * (1.0 - (v * v));
}
real(8) function code(v)
    real(8), intent (in) :: v
    code = ((sqrt(2.0d0) / 4.0d0) * sqrt((1.0d0 - (3.0d0 * (v * v))))) * (1.0d0 - (v * v))
end function
real(8) function code(v)
    real(8), intent (in) :: v
    code = sqrt((((1.0d0 + ((v ** 6.0d0) * (-27.0d0))) / (1.0d0 + ((v * v) * (3.0d0 + ((v * v) * 9.0d0))))) * 0.125d0)) * (1.0d0 - (v * v))
end function
public static double code(double v) {
	return ((Math.sqrt(2.0) / 4.0) * Math.sqrt((1.0 - (3.0 * (v * v))))) * (1.0 - (v * v));
}
public static double code(double v) {
	return Math.sqrt((((1.0 + (Math.pow(v, 6.0) * -27.0)) / (1.0 + ((v * v) * (3.0 + ((v * v) * 9.0))))) * 0.125)) * (1.0 - (v * v));
}
def code(v):
	return ((math.sqrt(2.0) / 4.0) * math.sqrt((1.0 - (3.0 * (v * v))))) * (1.0 - (v * v))
def code(v):
	return math.sqrt((((1.0 + (math.pow(v, 6.0) * -27.0)) / (1.0 + ((v * v) * (3.0 + ((v * v) * 9.0))))) * 0.125)) * (1.0 - (v * v))
function code(v)
	return Float64(Float64(Float64(sqrt(2.0) / 4.0) * sqrt(Float64(1.0 - Float64(3.0 * Float64(v * v))))) * Float64(1.0 - Float64(v * v)))
end
function code(v)
	return Float64(sqrt(Float64(Float64(Float64(1.0 + Float64((v ^ 6.0) * -27.0)) / Float64(1.0 + Float64(Float64(v * v) * Float64(3.0 + Float64(Float64(v * v) * 9.0))))) * 0.125)) * Float64(1.0 - Float64(v * v)))
end
function tmp = code(v)
	tmp = ((sqrt(2.0) / 4.0) * sqrt((1.0 - (3.0 * (v * v))))) * (1.0 - (v * v));
end
function tmp = code(v)
	tmp = sqrt((((1.0 + ((v ^ 6.0) * -27.0)) / (1.0 + ((v * v) * (3.0 + ((v * v) * 9.0))))) * 0.125)) * (1.0 - (v * v));
end
code[v_] := N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] / 4.0), $MachinePrecision] * N[Sqrt[N[(1.0 - N[(3.0 * N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(1.0 - N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[v_] := N[(N[Sqrt[N[(N[(N[(1.0 + N[(N[Power[v, 6.0], $MachinePrecision] * -27.0), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(N[(v * v), $MachinePrecision] * N[(3.0 + N[(N[(v * v), $MachinePrecision] * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.125), $MachinePrecision]], $MachinePrecision] * N[(1.0 - N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\left(\frac{\sqrt{2}}{4} \cdot \sqrt{1 - 3 \cdot \left(v \cdot v\right)}\right) \cdot \left(1 - v \cdot v\right)
\sqrt{\frac{1 + {v}^{6} \cdot -27}{1 + \left(v \cdot v\right) \cdot \left(3 + \left(v \cdot v\right) \cdot 9\right)} \cdot 0.125} \cdot \left(1 - v \cdot v\right)

Error

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 0.0

    \[\left(\frac{\sqrt{2}}{4} \cdot \sqrt{1 - 3 \cdot \left(v \cdot v\right)}\right) \cdot \left(1 - v \cdot v\right) \]
  2. Applied egg-rr0.0

    \[\leadsto \color{blue}{\sqrt{\mathsf{fma}\left(v \cdot v, -3, 1\right) \cdot 0.125}} \cdot \left(1 - v \cdot v\right) \]
  3. Applied egg-rr0.0

    \[\leadsto \sqrt{\color{blue}{\frac{{\left(v \cdot \left(v \cdot -3\right)\right)}^{3} + 1}{{\left(v \cdot \left(v \cdot -3\right)\right)}^{2} + \left(1 - \left(v \cdot \left(v \cdot -3\right)\right) \cdot 1\right)}} \cdot 0.125} \cdot \left(1 - v \cdot v\right) \]
  4. Simplified0.0

    \[\leadsto \sqrt{\color{blue}{\frac{1 + {v}^{6} \cdot -27}{1 + \left(v \cdot v\right) \cdot \left(3 + \left(v \cdot v\right) \cdot 9\right)}} \cdot 0.125} \cdot \left(1 - v \cdot v\right) \]
    Proof
    (/.f64 (+.f64 1 (*.f64 (pow.f64 v 6) -27)) (+.f64 1 (*.f64 (*.f64 v v) (+.f64 3 (*.f64 (*.f64 v v) 9))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 1 (*.f64 (pow.f64 v (Rewrite<= metadata-eval (*.f64 2 3))) -27)) (+.f64 1 (*.f64 (*.f64 v v) (+.f64 3 (*.f64 (*.f64 v v) 9))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 1 (*.f64 (Rewrite<= pow-sqr_binary64 (*.f64 (pow.f64 v 3) (pow.f64 v 3))) -27)) (+.f64 1 (*.f64 (*.f64 v v) (+.f64 3 (*.f64 (*.f64 v v) 9))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 1 (*.f64 (Rewrite<= cube-prod_binary64 (pow.f64 (*.f64 v v) 3)) -27)) (+.f64 1 (*.f64 (*.f64 v v) (+.f64 3 (*.f64 (*.f64 v v) 9))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 1 (*.f64 (pow.f64 (Rewrite<= unpow2_binary64 (pow.f64 v 2)) 3) -27)) (+.f64 1 (*.f64 (*.f64 v v) (+.f64 3 (*.f64 (*.f64 v v) 9))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 1 (*.f64 (pow.f64 (pow.f64 v 2) 3) (Rewrite<= metadata-eval (pow.f64 -3 3)))) (+.f64 1 (*.f64 (*.f64 v v) (+.f64 3 (*.f64 (*.f64 v v) 9))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 1 (Rewrite<= cube-prod_binary64 (pow.f64 (*.f64 (pow.f64 v 2) -3) 3))) (+.f64 1 (*.f64 (*.f64 v v) (+.f64 3 (*.f64 (*.f64 v v) 9))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 1 (pow.f64 (*.f64 (Rewrite=> unpow2_binary64 (*.f64 v v)) -3) 3)) (+.f64 1 (*.f64 (*.f64 v v) (+.f64 3 (*.f64 (*.f64 v v) 9))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 1 (pow.f64 (Rewrite<= associate-*r*_binary64 (*.f64 v (*.f64 v -3))) 3)) (+.f64 1 (*.f64 (*.f64 v v) (+.f64 3 (*.f64 (*.f64 v v) 9))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (Rewrite<= +-commutative_binary64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1)) (+.f64 1 (*.f64 (*.f64 v v) (+.f64 3 (*.f64 (*.f64 v v) 9))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (+.f64 1 (*.f64 (Rewrite<= unpow2_binary64 (pow.f64 v 2)) (+.f64 3 (*.f64 (*.f64 v v) 9))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (+.f64 1 (*.f64 (pow.f64 v 2) (+.f64 3 (*.f64 (*.f64 v v) (Rewrite<= metadata-eval (*.f64 -3 -3))))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (+.f64 1 (*.f64 (pow.f64 v 2) (+.f64 3 (Rewrite<= swap-sqr_binary64 (*.f64 (*.f64 v -3) (*.f64 v -3))))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (+.f64 1 (Rewrite<= distribute-lft-out_binary64 (+.f64 (*.f64 (pow.f64 v 2) 3) (*.f64 (pow.f64 v 2) (*.f64 (*.f64 v -3) (*.f64 v -3))))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (+.f64 1 (+.f64 (*.f64 (pow.f64 v 2) (Rewrite<= metadata-eval (neg.f64 -3))) (*.f64 (pow.f64 v 2) (*.f64 (*.f64 v -3) (*.f64 v -3)))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (+.f64 1 (+.f64 (Rewrite<= distribute-rgt-neg-in_binary64 (neg.f64 (*.f64 (pow.f64 v 2) -3))) (*.f64 (pow.f64 v 2) (*.f64 (*.f64 v -3) (*.f64 v -3)))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (+.f64 1 (+.f64 (neg.f64 (*.f64 (Rewrite=> unpow2_binary64 (*.f64 v v)) -3)) (*.f64 (pow.f64 v 2) (*.f64 (*.f64 v -3) (*.f64 v -3)))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (+.f64 1 (+.f64 (neg.f64 (Rewrite<= associate-*r*_binary64 (*.f64 v (*.f64 v -3)))) (*.f64 (pow.f64 v 2) (*.f64 (*.f64 v -3) (*.f64 v -3)))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (+.f64 1 (+.f64 (neg.f64 (*.f64 v (*.f64 v -3))) (*.f64 (Rewrite=> unpow2_binary64 (*.f64 v v)) (*.f64 (*.f64 v -3) (*.f64 v -3)))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (+.f64 1 (+.f64 (neg.f64 (*.f64 v (*.f64 v -3))) (Rewrite<= swap-sqr_binary64 (*.f64 (*.f64 v (*.f64 v -3)) (*.f64 v (*.f64 v -3))))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (+.f64 1 (+.f64 (neg.f64 (*.f64 v (*.f64 v -3))) (Rewrite<= unpow2_binary64 (pow.f64 (*.f64 v (*.f64 v -3)) 2))))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (Rewrite<= associate-+l+_binary64 (+.f64 (+.f64 1 (neg.f64 (*.f64 v (*.f64 v -3)))) (pow.f64 (*.f64 v (*.f64 v -3)) 2)))): 2 points increase in error, 2 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (+.f64 (+.f64 1 (neg.f64 (Rewrite<= *-rgt-identity_binary64 (*.f64 (*.f64 v (*.f64 v -3)) 1)))) (pow.f64 (*.f64 v (*.f64 v -3)) 2))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (+.f64 (Rewrite<= sub-neg_binary64 (-.f64 1 (*.f64 (*.f64 v (*.f64 v -3)) 1))) (pow.f64 (*.f64 v (*.f64 v -3)) 2))): 0 points increase in error, 0 points decrease in error
    (/.f64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 3) 1) (Rewrite<= +-commutative_binary64 (+.f64 (pow.f64 (*.f64 v (*.f64 v -3)) 2) (-.f64 1 (*.f64 (*.f64 v (*.f64 v -3)) 1))))): 0 points increase in error, 0 points decrease in error
  5. Final simplification0.0

    \[\leadsto \sqrt{\frac{1 + {v}^{6} \cdot -27}{1 + \left(v \cdot v\right) \cdot \left(3 + \left(v \cdot v\right) \cdot 9\right)} \cdot 0.125} \cdot \left(1 - v \cdot v\right) \]

Alternatives

Alternative 1
Error0.0
Cost7360
\[\left(1 - v \cdot v\right) \cdot \sqrt{0.125 \cdot \left(1 + v \cdot \left(v \cdot -3\right)\right)} \]
Alternative 2
Error0.3
Cost6976
\[\sqrt{2} \cdot \left(\left(v \cdot v\right) \cdot -0.625 + 0.25\right) \]
Alternative 3
Error0.6
Cost6464
\[\sqrt{0.125} \]

Error

Reproduce

herbie shell --seed 2022330 
(FPCore (v)
  :name "Falkner and Boettcher, Appendix B, 2"
  :precision binary64
  (* (* (/ (sqrt 2.0) 4.0) (sqrt (- 1.0 (* 3.0 (* v v))))) (- 1.0 (* v v))))