
(FPCore (v) :precision binary64 (* (* (/ (sqrt 2.0) 4.0) (sqrt (- 1.0 (* 3.0 (* v v))))) (- 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));
}
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
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));
}
def code(v): return ((math.sqrt(2.0) / 4.0) * math.sqrt((1.0 - (3.0 * (v * v))))) * (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 tmp = code(v) tmp = ((sqrt(2.0) / 4.0) * sqrt((1.0 - (3.0 * (v * v))))) * (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]
\begin{array}{l}
\\
\left(\frac{\sqrt{2}}{4} \cdot \sqrt{1 - 3 \cdot \left(v \cdot v\right)}\right) \cdot \left(1 - v \cdot v\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 6 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (v) :precision binary64 (* (* (/ (sqrt 2.0) 4.0) (sqrt (- 1.0 (* 3.0 (* v v))))) (- 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));
}
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
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));
}
def code(v): return ((math.sqrt(2.0) / 4.0) * math.sqrt((1.0 - (3.0 * (v * v))))) * (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 tmp = code(v) tmp = ((sqrt(2.0) / 4.0) * sqrt((1.0 - (3.0 * (v * v))))) * (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]
\begin{array}{l}
\\
\left(\frac{\sqrt{2}}{4} \cdot \sqrt{1 - 3 \cdot \left(v \cdot v\right)}\right) \cdot \left(1 - v \cdot v\right)
\end{array}
(FPCore (v) :precision binary64 (/ (* (sqrt (* 2.0 (fma (* v v) -3.0 1.0))) (- 1.0 (* v v))) 4.0))
double code(double v) {
return (sqrt((2.0 * fma((v * v), -3.0, 1.0))) * (1.0 - (v * v))) / 4.0;
}
function code(v) return Float64(Float64(sqrt(Float64(2.0 * fma(Float64(v * v), -3.0, 1.0))) * Float64(1.0 - Float64(v * v))) / 4.0) end
code[v_] := N[(N[(N[Sqrt[N[(2.0 * N[(N[(v * v), $MachinePrecision] * -3.0 + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(1.0 - N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 4.0), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sqrt{2 \cdot \mathsf{fma}\left(v \cdot v, -3, 1\right)} \cdot \left(1 - v \cdot v\right)}{4}
\end{array}
Initial program 100.0%
associate-*l/100.0%
associate-/r/100.0%
associate-*r/100.0%
sub-neg100.0%
+-commutative100.0%
*-commutative100.0%
distribute-rgt-neg-in100.0%
fma-def100.0%
metadata-eval100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
associate-/r*100.0%
metadata-eval100.0%
fma-neg100.0%
metadata-eval100.0%
Simplified100.0%
clear-num100.0%
un-div-inv100.0%
frac-2neg100.0%
metadata-eval100.0%
associate-/l/100.0%
fma-udef100.0%
metadata-eval100.0%
distribute-rgt-neg-in100.0%
*-commutative100.0%
+-commutative100.0%
sub-neg100.0%
fma-udef100.0%
distribute-neg-in100.0%
metadata-eval100.0%
+-commutative100.0%
sub-neg100.0%
Applied egg-rr100.0%
Final simplification100.0%
(FPCore (v) :precision binary64 (/ (sqrt (* 2.0 (fma (* v v) -3.0 1.0))) (/ 4.0 (- 1.0 (* v v)))))
double code(double v) {
return sqrt((2.0 * fma((v * v), -3.0, 1.0))) / (4.0 / (1.0 - (v * v)));
}
function code(v) return Float64(sqrt(Float64(2.0 * fma(Float64(v * v), -3.0, 1.0))) / Float64(4.0 / Float64(1.0 - Float64(v * v)))) end
code[v_] := N[(N[Sqrt[N[(2.0 * N[(N[(v * v), $MachinePrecision] * -3.0 + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(4.0 / N[(1.0 - N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sqrt{2 \cdot \mathsf{fma}\left(v \cdot v, -3, 1\right)}}{\frac{4}{1 - v \cdot v}}
\end{array}
Initial program 100.0%
associate-*l/100.0%
associate-/r/100.0%
associate-*r/100.0%
sub-neg100.0%
+-commutative100.0%
*-commutative100.0%
distribute-rgt-neg-in100.0%
fma-def100.0%
metadata-eval100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
associate-/r*100.0%
metadata-eval100.0%
fma-neg100.0%
metadata-eval100.0%
Simplified100.0%
clear-num100.0%
un-div-inv100.0%
frac-2neg100.0%
metadata-eval100.0%
associate-/l/100.0%
fma-udef100.0%
metadata-eval100.0%
distribute-rgt-neg-in100.0%
*-commutative100.0%
+-commutative100.0%
sub-neg100.0%
fma-udef100.0%
distribute-neg-in100.0%
metadata-eval100.0%
+-commutative100.0%
sub-neg100.0%
Applied egg-rr100.0%
associate-/l*100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (v) :precision binary64 (+ (sqrt 0.125) (* (* v v) (* -2.5 (sqrt 0.125)))))
double code(double v) {
return sqrt(0.125) + ((v * v) * (-2.5 * sqrt(0.125)));
}
real(8) function code(v)
real(8), intent (in) :: v
code = sqrt(0.125d0) + ((v * v) * ((-2.5d0) * sqrt(0.125d0)))
end function
public static double code(double v) {
return Math.sqrt(0.125) + ((v * v) * (-2.5 * Math.sqrt(0.125)));
}
def code(v): return math.sqrt(0.125) + ((v * v) * (-2.5 * math.sqrt(0.125)))
function code(v) return Float64(sqrt(0.125) + Float64(Float64(v * v) * Float64(-2.5 * sqrt(0.125)))) end
function tmp = code(v) tmp = sqrt(0.125) + ((v * v) * (-2.5 * sqrt(0.125))); end
code[v_] := N[(N[Sqrt[0.125], $MachinePrecision] + N[(N[(v * v), $MachinePrecision] * N[(-2.5 * N[Sqrt[0.125], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\sqrt{0.125} + \left(v \cdot v\right) \cdot \left(-2.5 \cdot \sqrt{0.125}\right)
\end{array}
Initial program 100.0%
associate-*l*100.0%
Simplified100.0%
Taylor expanded in v around 0 99.4%
*-commutative99.4%
unpow299.4%
Simplified99.4%
+-commutative99.4%
distribute-rgt-in99.4%
associate-*l*99.4%
add-sqr-sqrt99.4%
sqrt-unprod99.4%
frac-times99.4%
add-sqr-sqrt99.4%
metadata-eval99.4%
metadata-eval99.4%
*-un-lft-identity99.4%
add-sqr-sqrt97.9%
sqrt-unprod99.4%
frac-times99.4%
add-sqr-sqrt99.4%
metadata-eval99.4%
metadata-eval99.4%
Applied egg-rr99.4%
Final simplification99.4%
(FPCore (v) :precision binary64 (/ (sqrt 2.0) (/ 4.0 (fma v (* v -2.5) 1.0))))
double code(double v) {
return sqrt(2.0) / (4.0 / fma(v, (v * -2.5), 1.0));
}
function code(v) return Float64(sqrt(2.0) / Float64(4.0 / fma(v, Float64(v * -2.5), 1.0))) end
code[v_] := N[(N[Sqrt[2.0], $MachinePrecision] / N[(4.0 / N[(v * N[(v * -2.5), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sqrt{2}}{\frac{4}{\mathsf{fma}\left(v, v \cdot -2.5, 1\right)}}
\end{array}
Initial program 100.0%
associate-*l*100.0%
Simplified100.0%
Taylor expanded in v around 0 99.4%
*-commutative99.4%
unpow299.4%
Simplified99.4%
associate-*l/99.4%
associate-/l*99.4%
+-commutative99.4%
associate-*l*99.4%
fma-def99.4%
Applied egg-rr99.4%
Final simplification99.4%
(FPCore (v) :precision binary64 (* (/ (sqrt 2.0) 4.0) (+ 1.0 (* (* v v) -2.5))))
double code(double v) {
return (sqrt(2.0) / 4.0) * (1.0 + ((v * v) * -2.5));
}
real(8) function code(v)
real(8), intent (in) :: v
code = (sqrt(2.0d0) / 4.0d0) * (1.0d0 + ((v * v) * (-2.5d0)))
end function
public static double code(double v) {
return (Math.sqrt(2.0) / 4.0) * (1.0 + ((v * v) * -2.5));
}
def code(v): return (math.sqrt(2.0) / 4.0) * (1.0 + ((v * v) * -2.5))
function code(v) return Float64(Float64(sqrt(2.0) / 4.0) * Float64(1.0 + Float64(Float64(v * v) * -2.5))) end
function tmp = code(v) tmp = (sqrt(2.0) / 4.0) * (1.0 + ((v * v) * -2.5)); end
code[v_] := N[(N[(N[Sqrt[2.0], $MachinePrecision] / 4.0), $MachinePrecision] * N[(1.0 + N[(N[(v * v), $MachinePrecision] * -2.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sqrt{2}}{4} \cdot \left(1 + \left(v \cdot v\right) \cdot -2.5\right)
\end{array}
Initial program 100.0%
associate-*l*100.0%
Simplified100.0%
Taylor expanded in v around 0 99.4%
*-commutative99.4%
unpow299.4%
Simplified99.4%
Final simplification99.4%
(FPCore (v) :precision binary64 (sqrt 0.125))
double code(double v) {
return sqrt(0.125);
}
real(8) function code(v)
real(8), intent (in) :: v
code = sqrt(0.125d0)
end function
public static double code(double v) {
return Math.sqrt(0.125);
}
def code(v): return math.sqrt(0.125)
function code(v) return sqrt(0.125) end
function tmp = code(v) tmp = sqrt(0.125); end
code[v_] := N[Sqrt[0.125], $MachinePrecision]
\begin{array}{l}
\\
\sqrt{0.125}
\end{array}
Initial program 100.0%
associate-*l*100.0%
Simplified100.0%
Taylor expanded in v around 0 99.4%
*-commutative99.4%
unpow299.4%
Simplified99.4%
+-commutative99.4%
distribute-rgt-in99.4%
associate-*l*99.4%
add-sqr-sqrt99.4%
sqrt-unprod99.4%
frac-times99.4%
add-sqr-sqrt99.4%
metadata-eval99.4%
metadata-eval99.4%
*-un-lft-identity99.4%
add-sqr-sqrt97.9%
sqrt-unprod99.4%
frac-times99.4%
add-sqr-sqrt99.4%
metadata-eval99.4%
metadata-eval99.4%
Applied egg-rr99.4%
Taylor expanded in v around 0 98.6%
Final simplification98.6%
herbie shell --seed 2023187
(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))))