
(FPCore (v w r) :precision binary64 (- (- (+ 3.0 (/ 2.0 (* r r))) (/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v))) 4.5))
double code(double v, double w, double r) {
return ((3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) - 4.5;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
code = ((3.0d0 + (2.0d0 / (r * r))) - (((0.125d0 * (3.0d0 - (2.0d0 * v))) * (((w * w) * r) * r)) / (1.0d0 - v))) - 4.5d0
end function
public static double code(double v, double w, double r) {
return ((3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) - 4.5;
}
def code(v, w, r): return ((3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) - 4.5
function code(v, w, r) return Float64(Float64(Float64(3.0 + Float64(2.0 / Float64(r * r))) - Float64(Float64(Float64(0.125 * Float64(3.0 - Float64(2.0 * v))) * Float64(Float64(Float64(w * w) * r) * r)) / Float64(1.0 - v))) - 4.5) end
function tmp = code(v, w, r) tmp = ((3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) - 4.5; end
code[v_, w_, r_] := N[(N[(N[(3.0 + N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(0.125 * N[(3.0 - N[(2.0 * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(w * w), $MachinePrecision] * r), $MachinePrecision] * r), $MachinePrecision]), $MachinePrecision] / N[(1.0 - v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(3 + \frac{2}{r \cdot r}\right) - \frac{\left(0.125 \cdot \left(3 - 2 \cdot v\right)\right) \cdot \left(\left(\left(w \cdot w\right) \cdot r\right) \cdot r\right)}{1 - v}\right) - 4.5
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 10 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (v w r) :precision binary64 (- (- (+ 3.0 (/ 2.0 (* r r))) (/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v))) 4.5))
double code(double v, double w, double r) {
return ((3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) - 4.5;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
code = ((3.0d0 + (2.0d0 / (r * r))) - (((0.125d0 * (3.0d0 - (2.0d0 * v))) * (((w * w) * r) * r)) / (1.0d0 - v))) - 4.5d0
end function
public static double code(double v, double w, double r) {
return ((3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) - 4.5;
}
def code(v, w, r): return ((3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) - 4.5
function code(v, w, r) return Float64(Float64(Float64(3.0 + Float64(2.0 / Float64(r * r))) - Float64(Float64(Float64(0.125 * Float64(3.0 - Float64(2.0 * v))) * Float64(Float64(Float64(w * w) * r) * r)) / Float64(1.0 - v))) - 4.5) end
function tmp = code(v, w, r) tmp = ((3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) - 4.5; end
code[v_, w_, r_] := N[(N[(N[(3.0 + N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(0.125 * N[(3.0 - N[(2.0 * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(w * w), $MachinePrecision] * r), $MachinePrecision] * r), $MachinePrecision]), $MachinePrecision] / N[(1.0 - v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(3 + \frac{2}{r \cdot r}\right) - \frac{\left(0.125 \cdot \left(3 - 2 \cdot v\right)\right) \cdot \left(\left(\left(w \cdot w\right) \cdot r\right) \cdot r\right)}{1 - v}\right) - 4.5
\end{array}
(FPCore (v w r) :precision binary64 (- (+ 3.0 (* 2.0 (pow r -2.0))) (+ (* (* 0.125 (+ 3.0 (* -2.0 v))) (* (* r w) (/ (* r w) (- 1.0 v)))) 4.5)))
double code(double v, double w, double r) {
return (3.0 + (2.0 * pow(r, -2.0))) - (((0.125 * (3.0 + (-2.0 * v))) * ((r * w) * ((r * w) / (1.0 - v)))) + 4.5);
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
code = (3.0d0 + (2.0d0 * (r ** (-2.0d0)))) - (((0.125d0 * (3.0d0 + ((-2.0d0) * v))) * ((r * w) * ((r * w) / (1.0d0 - v)))) + 4.5d0)
end function
public static double code(double v, double w, double r) {
return (3.0 + (2.0 * Math.pow(r, -2.0))) - (((0.125 * (3.0 + (-2.0 * v))) * ((r * w) * ((r * w) / (1.0 - v)))) + 4.5);
}
def code(v, w, r): return (3.0 + (2.0 * math.pow(r, -2.0))) - (((0.125 * (3.0 + (-2.0 * v))) * ((r * w) * ((r * w) / (1.0 - v)))) + 4.5)
function code(v, w, r) return Float64(Float64(3.0 + Float64(2.0 * (r ^ -2.0))) - Float64(Float64(Float64(0.125 * Float64(3.0 + Float64(-2.0 * v))) * Float64(Float64(r * w) * Float64(Float64(r * w) / Float64(1.0 - v)))) + 4.5)) end
function tmp = code(v, w, r) tmp = (3.0 + (2.0 * (r ^ -2.0))) - (((0.125 * (3.0 + (-2.0 * v))) * ((r * w) * ((r * w) / (1.0 - v)))) + 4.5); end
code[v_, w_, r_] := N[(N[(3.0 + N[(2.0 * N[Power[r, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(0.125 * N[(3.0 + N[(-2.0 * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] / N[(1.0 - v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 4.5), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(3 + 2 \cdot {r}^{-2}\right) - \left(\left(0.125 \cdot \left(3 + -2 \cdot v\right)\right) \cdot \left(\left(r \cdot w\right) \cdot \frac{r \cdot w}{1 - v}\right) + 4.5\right)
\end{array}
Initial program 85.3%
associate--l-85.3%
associate-*l*78.7%
sqr-neg78.7%
associate-*l*85.3%
+-commutative85.3%
+-commutative85.3%
associate-/l*87.4%
fma-define87.5%
Simplified87.4%
*-un-lft-identity87.4%
add-sqr-sqrt87.4%
times-frac87.4%
associate-*r*80.9%
sqrt-prod80.9%
sqrt-prod41.9%
add-sqr-sqrt72.5%
sqrt-prod38.4%
add-sqr-sqrt75.6%
associate-*r*65.8%
sqrt-prod65.8%
sqrt-prod37.5%
add-sqr-sqrt75.6%
sqrt-prod50.2%
add-sqr-sqrt99.8%
Applied egg-rr99.8%
*-un-lft-identity99.8%
div-inv99.8%
pow299.8%
pow-flip99.9%
metadata-eval99.9%
Applied egg-rr99.9%
*-lft-identity99.9%
Simplified99.9%
Final simplification99.9%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (or (<= v -0.5) (not (<= v 2.2e-24)))
(- (+ 3.0 t_0) (+ 4.5 (/ (* w (* r 0.25)) (/ 1.0 (* r w)))))
(+ t_0 (+ -1.5 (* 0.375 (* (* r w) (/ (* r w) (+ v -1.0)))))))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if ((v <= -0.5) || !(v <= 2.2e-24)) {
tmp = (3.0 + t_0) - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w))));
} else {
tmp = t_0 + (-1.5 + (0.375 * ((r * w) * ((r * w) / (v + -1.0)))));
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: t_0
real(8) :: tmp
t_0 = 2.0d0 / (r * r)
if ((v <= (-0.5d0)) .or. (.not. (v <= 2.2d-24))) then
tmp = (3.0d0 + t_0) - (4.5d0 + ((w * (r * 0.25d0)) / (1.0d0 / (r * w))))
else
tmp = t_0 + ((-1.5d0) + (0.375d0 * ((r * w) * ((r * w) / (v + (-1.0d0))))))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if ((v <= -0.5) || !(v <= 2.2e-24)) {
tmp = (3.0 + t_0) - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w))));
} else {
tmp = t_0 + (-1.5 + (0.375 * ((r * w) * ((r * w) / (v + -1.0)))));
}
return tmp;
}
def code(v, w, r): t_0 = 2.0 / (r * r) tmp = 0 if (v <= -0.5) or not (v <= 2.2e-24): tmp = (3.0 + t_0) - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w)))) else: tmp = t_0 + (-1.5 + (0.375 * ((r * w) * ((r * w) / (v + -1.0))))) return tmp
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if ((v <= -0.5) || !(v <= 2.2e-24)) tmp = Float64(Float64(3.0 + t_0) - Float64(4.5 + Float64(Float64(w * Float64(r * 0.25)) / Float64(1.0 / Float64(r * w))))); else tmp = Float64(t_0 + Float64(-1.5 + Float64(0.375 * Float64(Float64(r * w) * Float64(Float64(r * w) / Float64(v + -1.0)))))); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 2.0 / (r * r); tmp = 0.0; if ((v <= -0.5) || ~((v <= 2.2e-24))) tmp = (3.0 + t_0) - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w)))); else tmp = t_0 + (-1.5 + (0.375 * ((r * w) * ((r * w) / (v + -1.0))))); end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[v, -0.5], N[Not[LessEqual[v, 2.2e-24]], $MachinePrecision]], N[(N[(3.0 + t$95$0), $MachinePrecision] - N[(4.5 + N[(N[(w * N[(r * 0.25), $MachinePrecision]), $MachinePrecision] / N[(1.0 / N[(r * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$0 + N[(-1.5 + N[(0.375 * N[(N[(r * w), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] / N[(v + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;v \leq -0.5 \lor \neg \left(v \leq 2.2 \cdot 10^{-24}\right):\\
\;\;\;\;\left(3 + t\_0\right) - \left(4.5 + \frac{w \cdot \left(r \cdot 0.25\right)}{\frac{1}{r \cdot w}}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0 + \left(-1.5 + 0.375 \cdot \left(\left(r \cdot w\right) \cdot \frac{r \cdot w}{v + -1}\right)\right)\\
\end{array}
\end{array}
if v < -0.5 or 2.20000000000000002e-24 < v Initial program 88.1%
associate--l-88.1%
associate-*l*84.3%
sqr-neg84.3%
associate-*l*88.1%
+-commutative88.1%
+-commutative88.1%
associate-/l*92.4%
fma-define92.4%
Simplified92.4%
*-un-lft-identity92.4%
add-sqr-sqrt92.3%
times-frac92.3%
associate-*r*88.6%
sqrt-prod88.6%
sqrt-prod47.6%
add-sqr-sqrt73.8%
sqrt-prod33.0%
add-sqr-sqrt71.6%
associate-*r*67.0%
sqrt-prod66.9%
sqrt-prod38.4%
add-sqr-sqrt76.2%
sqrt-prod47.5%
add-sqr-sqrt99.7%
Applied egg-rr99.7%
/-rgt-identity99.7%
*-commutative99.7%
clear-num99.7%
/-rgt-identity99.7%
clear-num99.7%
frac-times99.7%
metadata-eval99.7%
Applied egg-rr99.7%
*-commutative99.7%
associate-/r*99.7%
associate-*l/99.7%
clear-num99.7%
distribute-rgt-in99.7%
metadata-eval99.7%
*-commutative99.7%
associate-*l*99.7%
metadata-eval99.7%
Applied egg-rr99.7%
Taylor expanded in v around inf 98.8%
associate-*r*98.8%
Simplified98.8%
if -0.5 < v < 2.20000000000000002e-24Initial program 82.3%
Simplified82.3%
Taylor expanded in v around 0 81.8%
add-sqr-sqrt82.3%
*-un-lft-identity82.3%
times-frac82.3%
associate-*r*72.9%
sqrt-prod72.9%
sqrt-prod36.0%
add-sqr-sqrt71.1%
sqrt-prod43.9%
add-sqr-sqrt79.8%
associate-*r*64.7%
sqrt-prod64.7%
sqrt-prod36.6%
add-sqr-sqrt74.9%
sqrt-prod53.1%
add-sqr-sqrt99.8%
Applied egg-rr99.3%
Final simplification99.1%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (+ 3.0 (/ 2.0 (* r r)))))
(if (or (<= v -1.2) (not (<= v 2.2e-24)))
(- t_0 (+ 4.5 (/ (* w (* r 0.25)) (/ 1.0 (* r w)))))
(+ t_0 (- (/ (* (* r w) 0.375) (/ -1.0 (* r w))) 4.5)))))
double code(double v, double w, double r) {
double t_0 = 3.0 + (2.0 / (r * r));
double tmp;
if ((v <= -1.2) || !(v <= 2.2e-24)) {
tmp = t_0 - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w))));
} else {
tmp = t_0 + ((((r * w) * 0.375) / (-1.0 / (r * w))) - 4.5);
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: t_0
real(8) :: tmp
t_0 = 3.0d0 + (2.0d0 / (r * r))
if ((v <= (-1.2d0)) .or. (.not. (v <= 2.2d-24))) then
tmp = t_0 - (4.5d0 + ((w * (r * 0.25d0)) / (1.0d0 / (r * w))))
else
tmp = t_0 + ((((r * w) * 0.375d0) / ((-1.0d0) / (r * w))) - 4.5d0)
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double t_0 = 3.0 + (2.0 / (r * r));
double tmp;
if ((v <= -1.2) || !(v <= 2.2e-24)) {
tmp = t_0 - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w))));
} else {
tmp = t_0 + ((((r * w) * 0.375) / (-1.0 / (r * w))) - 4.5);
}
return tmp;
}
def code(v, w, r): t_0 = 3.0 + (2.0 / (r * r)) tmp = 0 if (v <= -1.2) or not (v <= 2.2e-24): tmp = t_0 - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w)))) else: tmp = t_0 + ((((r * w) * 0.375) / (-1.0 / (r * w))) - 4.5) return tmp
function code(v, w, r) t_0 = Float64(3.0 + Float64(2.0 / Float64(r * r))) tmp = 0.0 if ((v <= -1.2) || !(v <= 2.2e-24)) tmp = Float64(t_0 - Float64(4.5 + Float64(Float64(w * Float64(r * 0.25)) / Float64(1.0 / Float64(r * w))))); else tmp = Float64(t_0 + Float64(Float64(Float64(Float64(r * w) * 0.375) / Float64(-1.0 / Float64(r * w))) - 4.5)); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 3.0 + (2.0 / (r * r)); tmp = 0.0; if ((v <= -1.2) || ~((v <= 2.2e-24))) tmp = t_0 - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w)))); else tmp = t_0 + ((((r * w) * 0.375) / (-1.0 / (r * w))) - 4.5); end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(3.0 + N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[v, -1.2], N[Not[LessEqual[v, 2.2e-24]], $MachinePrecision]], N[(t$95$0 - N[(4.5 + N[(N[(w * N[(r * 0.25), $MachinePrecision]), $MachinePrecision] / N[(1.0 / N[(r * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$0 + N[(N[(N[(N[(r * w), $MachinePrecision] * 0.375), $MachinePrecision] / N[(-1.0 / N[(r * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 3 + \frac{2}{r \cdot r}\\
\mathbf{if}\;v \leq -1.2 \lor \neg \left(v \leq 2.2 \cdot 10^{-24}\right):\\
\;\;\;\;t\_0 - \left(4.5 + \frac{w \cdot \left(r \cdot 0.25\right)}{\frac{1}{r \cdot w}}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0 + \left(\frac{\left(r \cdot w\right) \cdot 0.375}{\frac{-1}{r \cdot w}} - 4.5\right)\\
\end{array}
\end{array}
if v < -1.19999999999999996 or 2.20000000000000002e-24 < v Initial program 88.1%
associate--l-88.1%
associate-*l*84.3%
sqr-neg84.3%
associate-*l*88.1%
+-commutative88.1%
+-commutative88.1%
associate-/l*92.4%
fma-define92.4%
Simplified92.4%
*-un-lft-identity92.4%
add-sqr-sqrt92.3%
times-frac92.3%
associate-*r*88.6%
sqrt-prod88.6%
sqrt-prod47.6%
add-sqr-sqrt73.8%
sqrt-prod33.0%
add-sqr-sqrt71.6%
associate-*r*67.0%
sqrt-prod66.9%
sqrt-prod38.4%
add-sqr-sqrt76.2%
sqrt-prod47.5%
add-sqr-sqrt99.7%
Applied egg-rr99.7%
/-rgt-identity99.7%
*-commutative99.7%
clear-num99.7%
/-rgt-identity99.7%
clear-num99.7%
frac-times99.7%
metadata-eval99.7%
Applied egg-rr99.7%
*-commutative99.7%
associate-/r*99.7%
associate-*l/99.7%
clear-num99.7%
distribute-rgt-in99.7%
metadata-eval99.7%
*-commutative99.7%
associate-*l*99.7%
metadata-eval99.7%
Applied egg-rr99.7%
Taylor expanded in v around inf 98.8%
associate-*r*98.8%
Simplified98.8%
if -1.19999999999999996 < v < 2.20000000000000002e-24Initial program 82.3%
associate--l-82.3%
associate-*l*72.9%
sqr-neg72.9%
associate-*l*82.3%
+-commutative82.3%
+-commutative82.3%
associate-/l*82.3%
fma-define82.3%
Simplified82.3%
*-un-lft-identity82.3%
add-sqr-sqrt82.3%
times-frac82.3%
associate-*r*72.9%
sqrt-prod72.9%
sqrt-prod36.0%
add-sqr-sqrt71.1%
sqrt-prod43.9%
add-sqr-sqrt79.8%
associate-*r*64.7%
sqrt-prod64.7%
sqrt-prod36.6%
add-sqr-sqrt74.9%
sqrt-prod53.1%
add-sqr-sqrt99.8%
Applied egg-rr99.8%
/-rgt-identity99.8%
*-commutative99.8%
clear-num99.8%
/-rgt-identity99.8%
clear-num99.8%
frac-times99.8%
metadata-eval99.8%
Applied egg-rr99.8%
*-commutative99.8%
associate-/r*99.8%
associate-*l/99.8%
clear-num99.8%
distribute-rgt-in99.8%
metadata-eval99.8%
*-commutative99.8%
associate-*l*99.8%
metadata-eval99.8%
Applied egg-rr99.8%
Taylor expanded in v around 0 99.3%
Final simplification99.1%
(FPCore (v w r)
:precision binary64
(if (<= r 0.22)
(- (+ 3.0 (/ 2.0 (* r r))) (+ 4.5 (/ (* w (* r 0.25)) (/ 1.0 (* r w)))))
(+
3.0
(-
(* (* 0.125 (+ 3.0 (* -2.0 v))) (* (* r w) (/ (* r w) (+ v -1.0))))
4.5))))
double code(double v, double w, double r) {
double tmp;
if (r <= 0.22) {
tmp = (3.0 + (2.0 / (r * r))) - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w))));
} else {
tmp = 3.0 + (((0.125 * (3.0 + (-2.0 * v))) * ((r * w) * ((r * w) / (v + -1.0)))) - 4.5);
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: tmp
if (r <= 0.22d0) then
tmp = (3.0d0 + (2.0d0 / (r * r))) - (4.5d0 + ((w * (r * 0.25d0)) / (1.0d0 / (r * w))))
else
tmp = 3.0d0 + (((0.125d0 * (3.0d0 + ((-2.0d0) * v))) * ((r * w) * ((r * w) / (v + (-1.0d0))))) - 4.5d0)
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 0.22) {
tmp = (3.0 + (2.0 / (r * r))) - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w))));
} else {
tmp = 3.0 + (((0.125 * (3.0 + (-2.0 * v))) * ((r * w) * ((r * w) / (v + -1.0)))) - 4.5);
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 0.22: tmp = (3.0 + (2.0 / (r * r))) - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w)))) else: tmp = 3.0 + (((0.125 * (3.0 + (-2.0 * v))) * ((r * w) * ((r * w) / (v + -1.0)))) - 4.5) return tmp
function code(v, w, r) tmp = 0.0 if (r <= 0.22) tmp = Float64(Float64(3.0 + Float64(2.0 / Float64(r * r))) - Float64(4.5 + Float64(Float64(w * Float64(r * 0.25)) / Float64(1.0 / Float64(r * w))))); else tmp = Float64(3.0 + Float64(Float64(Float64(0.125 * Float64(3.0 + Float64(-2.0 * v))) * Float64(Float64(r * w) * Float64(Float64(r * w) / Float64(v + -1.0)))) - 4.5)); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 0.22) tmp = (3.0 + (2.0 / (r * r))) - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w)))); else tmp = 3.0 + (((0.125 * (3.0 + (-2.0 * v))) * ((r * w) * ((r * w) / (v + -1.0)))) - 4.5); end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 0.22], N[(N[(3.0 + N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(4.5 + N[(N[(w * N[(r * 0.25), $MachinePrecision]), $MachinePrecision] / N[(1.0 / N[(r * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(3.0 + N[(N[(N[(0.125 * N[(3.0 + N[(-2.0 * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] / N[(v + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 0.22:\\
\;\;\;\;\left(3 + \frac{2}{r \cdot r}\right) - \left(4.5 + \frac{w \cdot \left(r \cdot 0.25\right)}{\frac{1}{r \cdot w}}\right)\\
\mathbf{else}:\\
\;\;\;\;3 + \left(\left(0.125 \cdot \left(3 + -2 \cdot v\right)\right) \cdot \left(\left(r \cdot w\right) \cdot \frac{r \cdot w}{v + -1}\right) - 4.5\right)\\
\end{array}
\end{array}
if r < 0.220000000000000001Initial program 85.0%
associate--l-85.0%
associate-*l*80.5%
sqr-neg80.5%
associate-*l*85.0%
+-commutative85.0%
+-commutative85.0%
associate-/l*84.9%
fma-define85.0%
Simplified84.9%
*-un-lft-identity84.9%
add-sqr-sqrt84.9%
times-frac84.9%
associate-*r*80.5%
sqrt-prod80.5%
sqrt-prod26.6%
add-sqr-sqrt65.8%
sqrt-prod34.6%
add-sqr-sqrt78.8%
associate-*r*69.7%
sqrt-prod69.7%
sqrt-prod30.0%
add-sqr-sqrt78.7%
sqrt-prod48.9%
add-sqr-sqrt99.8%
Applied egg-rr99.8%
/-rgt-identity99.8%
*-commutative99.8%
clear-num99.8%
/-rgt-identity99.8%
clear-num99.8%
frac-times99.8%
metadata-eval99.8%
Applied egg-rr99.8%
*-commutative99.8%
associate-/r*99.8%
associate-*l/99.8%
clear-num99.8%
distribute-rgt-in99.8%
metadata-eval99.8%
*-commutative99.8%
associate-*l*99.8%
metadata-eval99.8%
Applied egg-rr99.8%
Taylor expanded in v around inf 93.3%
associate-*r*93.3%
Simplified93.3%
if 0.220000000000000001 < r Initial program 86.3%
associate--l-86.3%
associate-*l*72.4%
sqr-neg72.4%
associate-*l*86.3%
+-commutative86.3%
+-commutative86.3%
associate-/l*96.4%
fma-define96.4%
Simplified96.4%
*-un-lft-identity96.4%
add-sqr-sqrt96.2%
times-frac96.1%
associate-*r*82.3%
sqrt-prod82.4%
sqrt-prod96.2%
add-sqr-sqrt96.3%
sqrt-prod51.7%
add-sqr-sqrt64.4%
associate-*r*51.9%
sqrt-prod51.9%
sqrt-prod64.2%
add-sqr-sqrt64.3%
sqrt-prod55.0%
add-sqr-sqrt99.7%
Applied egg-rr99.7%
Taylor expanded in r around inf 99.7%
Final simplification94.7%
(FPCore (v w r)
:precision binary64
(if (<= r 0.22)
(- (+ 3.0 (/ 2.0 (* r r))) (+ 4.5 (/ (* w (* r 0.25)) (/ 1.0 (* r w)))))
(+
3.0
(-
(/ (* (/ (* r w) (- 1.0 v)) (+ 0.375 (* v -0.25))) (/ -1.0 (* r w)))
4.5))))
double code(double v, double w, double r) {
double tmp;
if (r <= 0.22) {
tmp = (3.0 + (2.0 / (r * r))) - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w))));
} else {
tmp = 3.0 + (((((r * w) / (1.0 - v)) * (0.375 + (v * -0.25))) / (-1.0 / (r * w))) - 4.5);
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: tmp
if (r <= 0.22d0) then
tmp = (3.0d0 + (2.0d0 / (r * r))) - (4.5d0 + ((w * (r * 0.25d0)) / (1.0d0 / (r * w))))
else
tmp = 3.0d0 + (((((r * w) / (1.0d0 - v)) * (0.375d0 + (v * (-0.25d0)))) / ((-1.0d0) / (r * w))) - 4.5d0)
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 0.22) {
tmp = (3.0 + (2.0 / (r * r))) - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w))));
} else {
tmp = 3.0 + (((((r * w) / (1.0 - v)) * (0.375 + (v * -0.25))) / (-1.0 / (r * w))) - 4.5);
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 0.22: tmp = (3.0 + (2.0 / (r * r))) - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w)))) else: tmp = 3.0 + (((((r * w) / (1.0 - v)) * (0.375 + (v * -0.25))) / (-1.0 / (r * w))) - 4.5) return tmp
function code(v, w, r) tmp = 0.0 if (r <= 0.22) tmp = Float64(Float64(3.0 + Float64(2.0 / Float64(r * r))) - Float64(4.5 + Float64(Float64(w * Float64(r * 0.25)) / Float64(1.0 / Float64(r * w))))); else tmp = Float64(3.0 + Float64(Float64(Float64(Float64(Float64(r * w) / Float64(1.0 - v)) * Float64(0.375 + Float64(v * -0.25))) / Float64(-1.0 / Float64(r * w))) - 4.5)); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 0.22) tmp = (3.0 + (2.0 / (r * r))) - (4.5 + ((w * (r * 0.25)) / (1.0 / (r * w)))); else tmp = 3.0 + (((((r * w) / (1.0 - v)) * (0.375 + (v * -0.25))) / (-1.0 / (r * w))) - 4.5); end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 0.22], N[(N[(3.0 + N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(4.5 + N[(N[(w * N[(r * 0.25), $MachinePrecision]), $MachinePrecision] / N[(1.0 / N[(r * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(3.0 + N[(N[(N[(N[(N[(r * w), $MachinePrecision] / N[(1.0 - v), $MachinePrecision]), $MachinePrecision] * N[(0.375 + N[(v * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 / N[(r * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 0.22:\\
\;\;\;\;\left(3 + \frac{2}{r \cdot r}\right) - \left(4.5 + \frac{w \cdot \left(r \cdot 0.25\right)}{\frac{1}{r \cdot w}}\right)\\
\mathbf{else}:\\
\;\;\;\;3 + \left(\frac{\frac{r \cdot w}{1 - v} \cdot \left(0.375 + v \cdot -0.25\right)}{\frac{-1}{r \cdot w}} - 4.5\right)\\
\end{array}
\end{array}
if r < 0.220000000000000001Initial program 85.0%
associate--l-85.0%
associate-*l*80.5%
sqr-neg80.5%
associate-*l*85.0%
+-commutative85.0%
+-commutative85.0%
associate-/l*84.9%
fma-define85.0%
Simplified84.9%
*-un-lft-identity84.9%
add-sqr-sqrt84.9%
times-frac84.9%
associate-*r*80.5%
sqrt-prod80.5%
sqrt-prod26.6%
add-sqr-sqrt65.8%
sqrt-prod34.6%
add-sqr-sqrt78.8%
associate-*r*69.7%
sqrt-prod69.7%
sqrt-prod30.0%
add-sqr-sqrt78.7%
sqrt-prod48.9%
add-sqr-sqrt99.8%
Applied egg-rr99.8%
/-rgt-identity99.8%
*-commutative99.8%
clear-num99.8%
/-rgt-identity99.8%
clear-num99.8%
frac-times99.8%
metadata-eval99.8%
Applied egg-rr99.8%
*-commutative99.8%
associate-/r*99.8%
associate-*l/99.8%
clear-num99.8%
distribute-rgt-in99.8%
metadata-eval99.8%
*-commutative99.8%
associate-*l*99.8%
metadata-eval99.8%
Applied egg-rr99.8%
Taylor expanded in v around inf 93.3%
associate-*r*93.3%
Simplified93.3%
if 0.220000000000000001 < r Initial program 86.3%
associate--l-86.3%
associate-*l*72.4%
sqr-neg72.4%
associate-*l*86.3%
+-commutative86.3%
+-commutative86.3%
associate-/l*96.4%
fma-define96.4%
Simplified96.4%
*-un-lft-identity96.4%
add-sqr-sqrt96.2%
times-frac96.1%
associate-*r*82.3%
sqrt-prod82.4%
sqrt-prod96.2%
add-sqr-sqrt96.3%
sqrt-prod51.7%
add-sqr-sqrt64.4%
associate-*r*51.9%
sqrt-prod51.9%
sqrt-prod64.2%
add-sqr-sqrt64.3%
sqrt-prod55.0%
add-sqr-sqrt99.7%
Applied egg-rr99.7%
/-rgt-identity99.7%
*-commutative99.7%
clear-num99.6%
/-rgt-identity99.6%
clear-num99.6%
frac-times99.6%
metadata-eval99.6%
Applied egg-rr99.6%
*-commutative99.6%
associate-/r*99.6%
associate-*l/99.5%
clear-num99.6%
distribute-rgt-in99.6%
metadata-eval99.6%
*-commutative99.6%
associate-*l*99.6%
metadata-eval99.6%
Applied egg-rr99.6%
Taylor expanded in r around inf 99.6%
Final simplification94.7%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<= r 5e-157)
(- (+ 3.0 t_0) 4.5)
(+ t_0 (+ -1.5 (* 0.375 (/ (* r (* r (* w w))) (+ v -1.0))))))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (r <= 5e-157) {
tmp = (3.0 + t_0) - 4.5;
} else {
tmp = t_0 + (-1.5 + (0.375 * ((r * (r * (w * w))) / (v + -1.0))));
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: t_0
real(8) :: tmp
t_0 = 2.0d0 / (r * r)
if (r <= 5d-157) then
tmp = (3.0d0 + t_0) - 4.5d0
else
tmp = t_0 + ((-1.5d0) + (0.375d0 * ((r * (r * (w * w))) / (v + (-1.0d0)))))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (r <= 5e-157) {
tmp = (3.0 + t_0) - 4.5;
} else {
tmp = t_0 + (-1.5 + (0.375 * ((r * (r * (w * w))) / (v + -1.0))));
}
return tmp;
}
def code(v, w, r): t_0 = 2.0 / (r * r) tmp = 0 if r <= 5e-157: tmp = (3.0 + t_0) - 4.5 else: tmp = t_0 + (-1.5 + (0.375 * ((r * (r * (w * w))) / (v + -1.0)))) return tmp
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (r <= 5e-157) tmp = Float64(Float64(3.0 + t_0) - 4.5); else tmp = Float64(t_0 + Float64(-1.5 + Float64(0.375 * Float64(Float64(r * Float64(r * Float64(w * w))) / Float64(v + -1.0))))); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 2.0 / (r * r); tmp = 0.0; if (r <= 5e-157) tmp = (3.0 + t_0) - 4.5; else tmp = t_0 + (-1.5 + (0.375 * ((r * (r * (w * w))) / (v + -1.0)))); end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[r, 5e-157], N[(N[(3.0 + t$95$0), $MachinePrecision] - 4.5), $MachinePrecision], N[(t$95$0 + N[(-1.5 + N[(0.375 * N[(N[(r * N[(r * N[(w * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(v + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;r \leq 5 \cdot 10^{-157}:\\
\;\;\;\;\left(3 + t\_0\right) - 4.5\\
\mathbf{else}:\\
\;\;\;\;t\_0 + \left(-1.5 + 0.375 \cdot \frac{r \cdot \left(r \cdot \left(w \cdot w\right)\right)}{v + -1}\right)\\
\end{array}
\end{array}
if r < 5.0000000000000002e-157Initial program 84.7%
Simplified79.3%
Taylor expanded in r around 0 68.4%
if 5.0000000000000002e-157 < r Initial program 86.2%
Simplified92.5%
Taylor expanded in v around 0 69.7%
Final simplification68.9%
(FPCore (v w r) :precision binary64 (+ (/ 2.0 (* r r)) (+ (* (+ 0.375 (* v -0.25)) (* (* r w) (/ (* r w) (+ v -1.0)))) -1.5)))
double code(double v, double w, double r) {
return (2.0 / (r * r)) + (((0.375 + (v * -0.25)) * ((r * w) * ((r * w) / (v + -1.0)))) + -1.5);
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
code = (2.0d0 / (r * r)) + (((0.375d0 + (v * (-0.25d0))) * ((r * w) * ((r * w) / (v + (-1.0d0))))) + (-1.5d0))
end function
public static double code(double v, double w, double r) {
return (2.0 / (r * r)) + (((0.375 + (v * -0.25)) * ((r * w) * ((r * w) / (v + -1.0)))) + -1.5);
}
def code(v, w, r): return (2.0 / (r * r)) + (((0.375 + (v * -0.25)) * ((r * w) * ((r * w) / (v + -1.0)))) + -1.5)
function code(v, w, r) return Float64(Float64(2.0 / Float64(r * r)) + Float64(Float64(Float64(0.375 + Float64(v * -0.25)) * Float64(Float64(r * w) * Float64(Float64(r * w) / Float64(v + -1.0)))) + -1.5)) end
function tmp = code(v, w, r) tmp = (2.0 / (r * r)) + (((0.375 + (v * -0.25)) * ((r * w) * ((r * w) / (v + -1.0)))) + -1.5); end
code[v_, w_, r_] := N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(0.375 + N[(v * -0.25), $MachinePrecision]), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] / N[(v + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{r \cdot r} + \left(\left(0.375 + v \cdot -0.25\right) \cdot \left(\left(r \cdot w\right) \cdot \frac{r \cdot w}{v + -1}\right) + -1.5\right)
\end{array}
Initial program 85.3%
Simplified87.5%
Taylor expanded in v around 0 87.5%
*-commutative87.5%
Simplified87.5%
add-sqr-sqrt87.4%
*-un-lft-identity87.4%
times-frac87.4%
associate-*r*80.9%
sqrt-prod80.9%
sqrt-prod41.9%
add-sqr-sqrt72.5%
sqrt-prod38.4%
add-sqr-sqrt75.6%
associate-*r*65.8%
sqrt-prod65.8%
sqrt-prod37.5%
add-sqr-sqrt75.6%
sqrt-prod50.2%
add-sqr-sqrt99.8%
Applied egg-rr99.8%
Final simplification99.8%
(FPCore (v w r) :precision binary64 (+ (+ 3.0 (/ 2.0 (* r r))) (- (/ (* (* r w) 0.375) (/ -1.0 (* r w))) 4.5)))
double code(double v, double w, double r) {
return (3.0 + (2.0 / (r * r))) + ((((r * w) * 0.375) / (-1.0 / (r * w))) - 4.5);
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
code = (3.0d0 + (2.0d0 / (r * r))) + ((((r * w) * 0.375d0) / ((-1.0d0) / (r * w))) - 4.5d0)
end function
public static double code(double v, double w, double r) {
return (3.0 + (2.0 / (r * r))) + ((((r * w) * 0.375) / (-1.0 / (r * w))) - 4.5);
}
def code(v, w, r): return (3.0 + (2.0 / (r * r))) + ((((r * w) * 0.375) / (-1.0 / (r * w))) - 4.5)
function code(v, w, r) return Float64(Float64(3.0 + Float64(2.0 / Float64(r * r))) + Float64(Float64(Float64(Float64(r * w) * 0.375) / Float64(-1.0 / Float64(r * w))) - 4.5)) end
function tmp = code(v, w, r) tmp = (3.0 + (2.0 / (r * r))) + ((((r * w) * 0.375) / (-1.0 / (r * w))) - 4.5); end
code[v_, w_, r_] := N[(N[(3.0 + N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(N[(r * w), $MachinePrecision] * 0.375), $MachinePrecision] / N[(-1.0 / N[(r * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(3 + \frac{2}{r \cdot r}\right) + \left(\frac{\left(r \cdot w\right) \cdot 0.375}{\frac{-1}{r \cdot w}} - 4.5\right)
\end{array}
Initial program 85.3%
associate--l-85.3%
associate-*l*78.7%
sqr-neg78.7%
associate-*l*85.3%
+-commutative85.3%
+-commutative85.3%
associate-/l*87.4%
fma-define87.5%
Simplified87.4%
*-un-lft-identity87.4%
add-sqr-sqrt87.4%
times-frac87.4%
associate-*r*80.9%
sqrt-prod80.9%
sqrt-prod41.9%
add-sqr-sqrt72.5%
sqrt-prod38.4%
add-sqr-sqrt75.6%
associate-*r*65.8%
sqrt-prod65.8%
sqrt-prod37.5%
add-sqr-sqrt75.6%
sqrt-prod50.2%
add-sqr-sqrt99.8%
Applied egg-rr99.8%
/-rgt-identity99.8%
*-commutative99.8%
clear-num99.7%
/-rgt-identity99.7%
clear-num99.7%
frac-times99.7%
metadata-eval99.7%
Applied egg-rr99.7%
*-commutative99.7%
associate-/r*99.7%
associate-*l/99.7%
clear-num99.8%
distribute-rgt-in99.8%
metadata-eval99.8%
*-commutative99.8%
associate-*l*99.8%
metadata-eval99.8%
Applied egg-rr99.8%
Taylor expanded in v around 0 91.7%
Final simplification91.7%
(FPCore (v w r) :precision binary64 (- (+ 3.0 (/ 2.0 (* r r))) 4.5))
double code(double v, double w, double r) {
return (3.0 + (2.0 / (r * r))) - 4.5;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
code = (3.0d0 + (2.0d0 / (r * r))) - 4.5d0
end function
public static double code(double v, double w, double r) {
return (3.0 + (2.0 / (r * r))) - 4.5;
}
def code(v, w, r): return (3.0 + (2.0 / (r * r))) - 4.5
function code(v, w, r) return Float64(Float64(3.0 + Float64(2.0 / Float64(r * r))) - 4.5) end
function tmp = code(v, w, r) tmp = (3.0 + (2.0 / (r * r))) - 4.5; end
code[v_, w_, r_] := N[(N[(3.0 + N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]
\begin{array}{l}
\\
\left(3 + \frac{2}{r \cdot r}\right) - 4.5
\end{array}
Initial program 85.3%
Simplified80.3%
Taylor expanded in r around 0 61.7%
(FPCore (v w r) :precision binary64 -1.5)
double code(double v, double w, double r) {
return -1.5;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
code = -1.5d0
end function
public static double code(double v, double w, double r) {
return -1.5;
}
def code(v, w, r): return -1.5
function code(v, w, r) return -1.5 end
function tmp = code(v, w, r) tmp = -1.5; end
code[v_, w_, r_] := -1.5
\begin{array}{l}
\\
-1.5
\end{array}
Initial program 85.3%
Simplified80.3%
Taylor expanded in r around 0 61.7%
Taylor expanded in r around inf 13.8%
herbie shell --seed 2024139
(FPCore (v w r)
:name "Rosa's TurbineBenchmark"
:precision binary64
(- (- (+ 3.0 (/ 2.0 (* r r))) (/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v))) 4.5))