
(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 21 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 (+ (/ 2.0 (* r r)) (- 3.0 (fma (/ (pow (* w r) 2.0) (- 1.0 v)) (* (fma -2.0 v 3.0) 0.125) 4.5))))
double code(double v, double w, double r) {
return (2.0 / (r * r)) + (3.0 - fma((pow((w * r), 2.0) / (1.0 - v)), (fma(-2.0, v, 3.0) * 0.125), 4.5));
}
function code(v, w, r) return Float64(Float64(2.0 / Float64(r * r)) + Float64(3.0 - fma(Float64((Float64(w * r) ^ 2.0) / Float64(1.0 - v)), Float64(fma(-2.0, v, 3.0) * 0.125), 4.5))) end
code[v_, w_, r_] := N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + N[(3.0 - N[(N[(N[Power[N[(w * r), $MachinePrecision], 2.0], $MachinePrecision] / N[(1.0 - v), $MachinePrecision]), $MachinePrecision] * N[(N[(-2.0 * v + 3.0), $MachinePrecision] * 0.125), $MachinePrecision] + 4.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{r \cdot r} + \left(3 - \mathsf{fma}\left(\frac{{\left(w \cdot r\right)}^{2}}{1 - v}, \mathsf{fma}\left(-2, v, 3\right) \cdot 0.125, 4.5\right)\right)
\end{array}
Initial program 85.0%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.8%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r)))
(t_1
(-
(+ 3.0 t_0)
(/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v)))))
(if (<= t_1 (- INFINITY))
(fma (* (* -0.25 (* r r)) w) w (- t_0 1.5))
(if (<= t_1 -2000000.0)
(* (* (* (/ w (- 1.0 v)) (* (fma -2.0 v 3.0) w)) (* -0.125 r)) r)
(+ (/ (/ 2.0 r) r) -1.5)))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double t_1 = (3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = fma(((-0.25 * (r * r)) * w), w, (t_0 - 1.5));
} else if (t_1 <= -2000000.0) {
tmp = (((w / (1.0 - v)) * (fma(-2.0, v, 3.0) * w)) * (-0.125 * r)) * r;
} else {
tmp = ((2.0 / r) / r) + -1.5;
}
return tmp;
}
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) t_1 = Float64(Float64(3.0 + t_0) - Float64(Float64(Float64(0.125 * Float64(3.0 - Float64(2.0 * v))) * Float64(Float64(Float64(w * w) * r) * r)) / Float64(1.0 - v))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = fma(Float64(Float64(-0.25 * Float64(r * r)) * w), w, Float64(t_0 - 1.5)); elseif (t_1 <= -2000000.0) tmp = Float64(Float64(Float64(Float64(w / Float64(1.0 - v)) * Float64(fma(-2.0, v, 3.0) * w)) * Float64(-0.125 * r)) * r); else tmp = Float64(Float64(Float64(2.0 / r) / r) + -1.5); end return tmp end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(3.0 + t$95$0), $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]}, If[LessEqual[t$95$1, (-Infinity)], N[(N[(N[(-0.25 * N[(r * r), $MachinePrecision]), $MachinePrecision] * w), $MachinePrecision] * w + N[(t$95$0 - 1.5), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -2000000.0], N[(N[(N[(N[(w / N[(1.0 - v), $MachinePrecision]), $MachinePrecision] * N[(N[(-2.0 * v + 3.0), $MachinePrecision] * w), $MachinePrecision]), $MachinePrecision] * N[(-0.125 * r), $MachinePrecision]), $MachinePrecision] * r), $MachinePrecision], N[(N[(N[(2.0 / r), $MachinePrecision] / r), $MachinePrecision] + -1.5), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
t_1 := \left(3 + t\_0\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}\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;\mathsf{fma}\left(\left(-0.25 \cdot \left(r \cdot r\right)\right) \cdot w, w, t\_0 - 1.5\right)\\
\mathbf{elif}\;t\_1 \leq -2000000:\\
\;\;\;\;\left(\left(\frac{w}{1 - v} \cdot \left(\mathsf{fma}\left(-2, v, 3\right) \cdot w\right)\right) \cdot \left(-0.125 \cdot r\right)\right) \cdot r\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{r}}{r} + -1.5\\
\end{array}
\end{array}
if (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -inf.0Initial program 85.9%
Taylor expanded in v around inf
sub-negN/A
+-commutativeN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
associate-+l+N/A
distribute-lft-neg-inN/A
metadata-evalN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower-fma.f64N/A
Applied rewrites94.9%
if -inf.0 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -2e6Initial program 98.3%
Taylor expanded in w around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
lower--.f6479.0
Applied rewrites79.0%
Applied rewrites96.5%
if -2e6 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) Initial program 82.6%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.8%
Taylor expanded in w around 0
Applied rewrites94.8%
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6494.8
Applied rewrites94.8%
(FPCore (v w r)
:precision binary64
(let* ((t_0
(-
(+ 3.0 (/ 2.0 (* r r)))
(/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v)))))
(if (<= t_0 -5e+247)
(* (* (* -0.25 (* w w)) r) r)
(if (<= t_0 -2000000.0)
(* (* (fma (* -0.375 w) w (/ -1.5 (* r r))) r) r)
(+ (/ (/ 2.0 r) r) -1.5)))))
double code(double v, double w, double r) {
double t_0 = (3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v));
double tmp;
if (t_0 <= -5e+247) {
tmp = ((-0.25 * (w * w)) * r) * r;
} else if (t_0 <= -2000000.0) {
tmp = (fma((-0.375 * w), w, (-1.5 / (r * r))) * r) * r;
} else {
tmp = ((2.0 / r) / r) + -1.5;
}
return tmp;
}
function code(v, w, r) t_0 = 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))) tmp = 0.0 if (t_0 <= -5e+247) tmp = Float64(Float64(Float64(-0.25 * Float64(w * w)) * r) * r); elseif (t_0 <= -2000000.0) tmp = Float64(Float64(fma(Float64(-0.375 * w), w, Float64(-1.5 / Float64(r * r))) * r) * r); else tmp = Float64(Float64(Float64(2.0 / r) / r) + -1.5); end return tmp end
code[v_, w_, r_] := Block[{t$95$0 = 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]}, If[LessEqual[t$95$0, -5e+247], N[(N[(N[(-0.25 * N[(w * w), $MachinePrecision]), $MachinePrecision] * r), $MachinePrecision] * r), $MachinePrecision], If[LessEqual[t$95$0, -2000000.0], N[(N[(N[(N[(-0.375 * w), $MachinePrecision] * w + N[(-1.5 / N[(r * r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * r), $MachinePrecision] * r), $MachinePrecision], N[(N[(N[(2.0 / r), $MachinePrecision] / r), $MachinePrecision] + -1.5), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \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}\\
\mathbf{if}\;t\_0 \leq -5 \cdot 10^{+247}:\\
\;\;\;\;\left(\left(-0.25 \cdot \left(w \cdot w\right)\right) \cdot r\right) \cdot r\\
\mathbf{elif}\;t\_0 \leq -2000000:\\
\;\;\;\;\left(\mathsf{fma}\left(-0.375 \cdot w, w, \frac{-1.5}{r \cdot r}\right) \cdot r\right) \cdot r\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{r}}{r} + -1.5\\
\end{array}
\end{array}
if (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -5.00000000000000023e247Initial program 86.0%
Taylor expanded in w around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
lower--.f6488.1
Applied rewrites88.1%
Applied rewrites89.0%
Taylor expanded in v around inf
Applied rewrites91.2%
if -5.00000000000000023e247 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -2e6Initial program 98.3%
Taylor expanded in v around 0
sub-negN/A
+-commutativeN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-inN/A
metadata-evalN/A
associate-+l+N/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower-fma.f64N/A
Applied rewrites73.1%
Taylor expanded in r around inf
Applied rewrites73.1%
Applied rewrites83.8%
if -2e6 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) Initial program 82.6%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.8%
Taylor expanded in w around 0
Applied rewrites94.8%
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6494.8
Applied rewrites94.8%
(FPCore (v w r)
:precision binary64
(let* ((t_0
(-
(+ 3.0 (/ 2.0 (* r r)))
(/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v)))))
(if (<= t_0 -5e+247)
(* (* (* -0.25 (* w w)) r) r)
(if (<= t_0 -2000000.0)
(* (* (* (* w w) 3.0) (* -0.125 r)) r)
(+ (/ (/ 2.0 r) r) -1.5)))))
double code(double v, double w, double r) {
double t_0 = (3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v));
double tmp;
if (t_0 <= -5e+247) {
tmp = ((-0.25 * (w * w)) * r) * r;
} else if (t_0 <= -2000000.0) {
tmp = (((w * w) * 3.0) * (-0.125 * r)) * r;
} else {
tmp = ((2.0 / r) / r) + -1.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))) - (((0.125d0 * (3.0d0 - (2.0d0 * v))) * (((w * w) * r) * r)) / (1.0d0 - v))
if (t_0 <= (-5d+247)) then
tmp = (((-0.25d0) * (w * w)) * r) * r
else if (t_0 <= (-2000000.0d0)) then
tmp = (((w * w) * 3.0d0) * ((-0.125d0) * r)) * r
else
tmp = ((2.0d0 / r) / r) + (-1.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))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v));
double tmp;
if (t_0 <= -5e+247) {
tmp = ((-0.25 * (w * w)) * r) * r;
} else if (t_0 <= -2000000.0) {
tmp = (((w * w) * 3.0) * (-0.125 * r)) * r;
} else {
tmp = ((2.0 / r) / r) + -1.5;
}
return tmp;
}
def code(v, w, r): t_0 = (3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v)) tmp = 0 if t_0 <= -5e+247: tmp = ((-0.25 * (w * w)) * r) * r elif t_0 <= -2000000.0: tmp = (((w * w) * 3.0) * (-0.125 * r)) * r else: tmp = ((2.0 / r) / r) + -1.5 return tmp
function code(v, w, r) t_0 = 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))) tmp = 0.0 if (t_0 <= -5e+247) tmp = Float64(Float64(Float64(-0.25 * Float64(w * w)) * r) * r); elseif (t_0 <= -2000000.0) tmp = Float64(Float64(Float64(Float64(w * w) * 3.0) * Float64(-0.125 * r)) * r); else tmp = Float64(Float64(Float64(2.0 / r) / r) + -1.5); end return tmp end
function tmp_2 = code(v, w, r) t_0 = (3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v)); tmp = 0.0; if (t_0 <= -5e+247) tmp = ((-0.25 * (w * w)) * r) * r; elseif (t_0 <= -2000000.0) tmp = (((w * w) * 3.0) * (-0.125 * r)) * r; else tmp = ((2.0 / r) / r) + -1.5; end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = 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]}, If[LessEqual[t$95$0, -5e+247], N[(N[(N[(-0.25 * N[(w * w), $MachinePrecision]), $MachinePrecision] * r), $MachinePrecision] * r), $MachinePrecision], If[LessEqual[t$95$0, -2000000.0], N[(N[(N[(N[(w * w), $MachinePrecision] * 3.0), $MachinePrecision] * N[(-0.125 * r), $MachinePrecision]), $MachinePrecision] * r), $MachinePrecision], N[(N[(N[(2.0 / r), $MachinePrecision] / r), $MachinePrecision] + -1.5), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \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}\\
\mathbf{if}\;t\_0 \leq -5 \cdot 10^{+247}:\\
\;\;\;\;\left(\left(-0.25 \cdot \left(w \cdot w\right)\right) \cdot r\right) \cdot r\\
\mathbf{elif}\;t\_0 \leq -2000000:\\
\;\;\;\;\left(\left(\left(w \cdot w\right) \cdot 3\right) \cdot \left(-0.125 \cdot r\right)\right) \cdot r\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{r}}{r} + -1.5\\
\end{array}
\end{array}
if (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -5.00000000000000023e247Initial program 86.0%
Taylor expanded in w around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
lower--.f6488.1
Applied rewrites88.1%
Applied rewrites89.0%
Taylor expanded in v around inf
Applied rewrites91.2%
if -5.00000000000000023e247 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -2e6Initial program 98.3%
Taylor expanded in w around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
lower--.f6482.6
Applied rewrites82.6%
Applied rewrites96.3%
Taylor expanded in v around 0
Applied rewrites81.2%
if -2e6 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) Initial program 82.6%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.8%
Taylor expanded in w around 0
Applied rewrites94.8%
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6494.8
Applied rewrites94.8%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r)))
(t_1
(-
(+ 3.0 t_0)
(/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v)))))
(if (<= t_1 -5e+247)
(* (* (* -0.25 (* w w)) r) r)
(if (<= t_1 -2000000.0)
(* (* (* (* w w) 3.0) (* -0.125 r)) r)
(- t_0 1.5)))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double t_1 = (3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v));
double tmp;
if (t_1 <= -5e+247) {
tmp = ((-0.25 * (w * w)) * r) * r;
} else if (t_1 <= -2000000.0) {
tmp = (((w * w) * 3.0) * (-0.125 * r)) * r;
} else {
tmp = t_0 - 1.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) :: t_1
real(8) :: tmp
t_0 = 2.0d0 / (r * r)
t_1 = (3.0d0 + t_0) - (((0.125d0 * (3.0d0 - (2.0d0 * v))) * (((w * w) * r) * r)) / (1.0d0 - v))
if (t_1 <= (-5d+247)) then
tmp = (((-0.25d0) * (w * w)) * r) * r
else if (t_1 <= (-2000000.0d0)) then
tmp = (((w * w) * 3.0d0) * ((-0.125d0) * r)) * r
else
tmp = t_0 - 1.5d0
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double t_1 = (3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v));
double tmp;
if (t_1 <= -5e+247) {
tmp = ((-0.25 * (w * w)) * r) * r;
} else if (t_1 <= -2000000.0) {
tmp = (((w * w) * 3.0) * (-0.125 * r)) * r;
} else {
tmp = t_0 - 1.5;
}
return tmp;
}
def code(v, w, r): t_0 = 2.0 / (r * r) t_1 = (3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v)) tmp = 0 if t_1 <= -5e+247: tmp = ((-0.25 * (w * w)) * r) * r elif t_1 <= -2000000.0: tmp = (((w * w) * 3.0) * (-0.125 * r)) * r else: tmp = t_0 - 1.5 return tmp
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) t_1 = Float64(Float64(3.0 + t_0) - Float64(Float64(Float64(0.125 * Float64(3.0 - Float64(2.0 * v))) * Float64(Float64(Float64(w * w) * r) * r)) / Float64(1.0 - v))) tmp = 0.0 if (t_1 <= -5e+247) tmp = Float64(Float64(Float64(-0.25 * Float64(w * w)) * r) * r); elseif (t_1 <= -2000000.0) tmp = Float64(Float64(Float64(Float64(w * w) * 3.0) * Float64(-0.125 * r)) * r); else tmp = Float64(t_0 - 1.5); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 2.0 / (r * r); t_1 = (3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v)); tmp = 0.0; if (t_1 <= -5e+247) tmp = ((-0.25 * (w * w)) * r) * r; elseif (t_1 <= -2000000.0) tmp = (((w * w) * 3.0) * (-0.125 * r)) * r; else tmp = t_0 - 1.5; end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(3.0 + t$95$0), $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]}, If[LessEqual[t$95$1, -5e+247], N[(N[(N[(-0.25 * N[(w * w), $MachinePrecision]), $MachinePrecision] * r), $MachinePrecision] * r), $MachinePrecision], If[LessEqual[t$95$1, -2000000.0], N[(N[(N[(N[(w * w), $MachinePrecision] * 3.0), $MachinePrecision] * N[(-0.125 * r), $MachinePrecision]), $MachinePrecision] * r), $MachinePrecision], N[(t$95$0 - 1.5), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
t_1 := \left(3 + t\_0\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}\\
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{+247}:\\
\;\;\;\;\left(\left(-0.25 \cdot \left(w \cdot w\right)\right) \cdot r\right) \cdot r\\
\mathbf{elif}\;t\_1 \leq -2000000:\\
\;\;\;\;\left(\left(\left(w \cdot w\right) \cdot 3\right) \cdot \left(-0.125 \cdot r\right)\right) \cdot r\\
\mathbf{else}:\\
\;\;\;\;t\_0 - 1.5\\
\end{array}
\end{array}
if (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -5.00000000000000023e247Initial program 86.0%
Taylor expanded in w around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
lower--.f6488.1
Applied rewrites88.1%
Applied rewrites89.0%
Taylor expanded in v around inf
Applied rewrites91.2%
if -5.00000000000000023e247 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -2e6Initial program 98.3%
Taylor expanded in w around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
lower--.f6482.6
Applied rewrites82.6%
Applied rewrites96.3%
Taylor expanded in v around 0
Applied rewrites81.2%
if -2e6 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) Initial program 82.6%
Taylor expanded in w around 0
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f6494.8
Applied rewrites94.8%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (* (* w w) r))
(t_1 (/ 2.0 (* r r)))
(t_2
(-
(+ 3.0 t_1)
(/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* t_0 r)) (- 1.0 v)))))
(if (<= t_2 -5e+247)
(* (* (* -0.25 (* w w)) r) r)
(if (<= t_2 -2000000.0) (* (* t_0 -0.375) r) (- t_1 1.5)))))
double code(double v, double w, double r) {
double t_0 = (w * w) * r;
double t_1 = 2.0 / (r * r);
double t_2 = (3.0 + t_1) - (((0.125 * (3.0 - (2.0 * v))) * (t_0 * r)) / (1.0 - v));
double tmp;
if (t_2 <= -5e+247) {
tmp = ((-0.25 * (w * w)) * r) * r;
} else if (t_2 <= -2000000.0) {
tmp = (t_0 * -0.375) * r;
} else {
tmp = t_1 - 1.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) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = (w * w) * r
t_1 = 2.0d0 / (r * r)
t_2 = (3.0d0 + t_1) - (((0.125d0 * (3.0d0 - (2.0d0 * v))) * (t_0 * r)) / (1.0d0 - v))
if (t_2 <= (-5d+247)) then
tmp = (((-0.25d0) * (w * w)) * r) * r
else if (t_2 <= (-2000000.0d0)) then
tmp = (t_0 * (-0.375d0)) * r
else
tmp = t_1 - 1.5d0
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double t_0 = (w * w) * r;
double t_1 = 2.0 / (r * r);
double t_2 = (3.0 + t_1) - (((0.125 * (3.0 - (2.0 * v))) * (t_0 * r)) / (1.0 - v));
double tmp;
if (t_2 <= -5e+247) {
tmp = ((-0.25 * (w * w)) * r) * r;
} else if (t_2 <= -2000000.0) {
tmp = (t_0 * -0.375) * r;
} else {
tmp = t_1 - 1.5;
}
return tmp;
}
def code(v, w, r): t_0 = (w * w) * r t_1 = 2.0 / (r * r) t_2 = (3.0 + t_1) - (((0.125 * (3.0 - (2.0 * v))) * (t_0 * r)) / (1.0 - v)) tmp = 0 if t_2 <= -5e+247: tmp = ((-0.25 * (w * w)) * r) * r elif t_2 <= -2000000.0: tmp = (t_0 * -0.375) * r else: tmp = t_1 - 1.5 return tmp
function code(v, w, r) t_0 = Float64(Float64(w * w) * r) t_1 = Float64(2.0 / Float64(r * r)) t_2 = Float64(Float64(3.0 + t_1) - Float64(Float64(Float64(0.125 * Float64(3.0 - Float64(2.0 * v))) * Float64(t_0 * r)) / Float64(1.0 - v))) tmp = 0.0 if (t_2 <= -5e+247) tmp = Float64(Float64(Float64(-0.25 * Float64(w * w)) * r) * r); elseif (t_2 <= -2000000.0) tmp = Float64(Float64(t_0 * -0.375) * r); else tmp = Float64(t_1 - 1.5); end return tmp end
function tmp_2 = code(v, w, r) t_0 = (w * w) * r; t_1 = 2.0 / (r * r); t_2 = (3.0 + t_1) - (((0.125 * (3.0 - (2.0 * v))) * (t_0 * r)) / (1.0 - v)); tmp = 0.0; if (t_2 <= -5e+247) tmp = ((-0.25 * (w * w)) * r) * r; elseif (t_2 <= -2000000.0) tmp = (t_0 * -0.375) * r; else tmp = t_1 - 1.5; end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(N[(w * w), $MachinePrecision] * r), $MachinePrecision]}, Block[{t$95$1 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(3.0 + t$95$1), $MachinePrecision] - N[(N[(N[(0.125 * N[(3.0 - N[(2.0 * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t$95$0 * r), $MachinePrecision]), $MachinePrecision] / N[(1.0 - v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -5e+247], N[(N[(N[(-0.25 * N[(w * w), $MachinePrecision]), $MachinePrecision] * r), $MachinePrecision] * r), $MachinePrecision], If[LessEqual[t$95$2, -2000000.0], N[(N[(t$95$0 * -0.375), $MachinePrecision] * r), $MachinePrecision], N[(t$95$1 - 1.5), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(w \cdot w\right) \cdot r\\
t_1 := \frac{2}{r \cdot r}\\
t_2 := \left(3 + t\_1\right) - \frac{\left(0.125 \cdot \left(3 - 2 \cdot v\right)\right) \cdot \left(t\_0 \cdot r\right)}{1 - v}\\
\mathbf{if}\;t\_2 \leq -5 \cdot 10^{+247}:\\
\;\;\;\;\left(\left(-0.25 \cdot \left(w \cdot w\right)\right) \cdot r\right) \cdot r\\
\mathbf{elif}\;t\_2 \leq -2000000:\\
\;\;\;\;\left(t\_0 \cdot -0.375\right) \cdot r\\
\mathbf{else}:\\
\;\;\;\;t\_1 - 1.5\\
\end{array}
\end{array}
if (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -5.00000000000000023e247Initial program 86.0%
Taylor expanded in w around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
lower--.f6488.1
Applied rewrites88.1%
Applied rewrites89.0%
Taylor expanded in v around inf
Applied rewrites91.2%
if -5.00000000000000023e247 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -2e6Initial program 98.3%
Taylor expanded in w around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
lower--.f6482.6
Applied rewrites82.6%
Applied rewrites96.3%
Taylor expanded in v around 0
Applied rewrites81.1%
if -2e6 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) Initial program 82.6%
Taylor expanded in w around 0
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f6494.8
Applied rewrites94.8%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<=
(-
(+ 3.0 t_0)
(/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v)))
-2000000.0)
(+ t_0 (- 3.0 (fma (* (* (* r r) 0.375) w) w 4.5)))
(+ (/ (/ 2.0 r) r) -1.5))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) {
tmp = t_0 + (3.0 - fma((((r * r) * 0.375) * w), w, 4.5));
} else {
tmp = ((2.0 / r) / r) + -1.5;
}
return tmp;
}
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (Float64(Float64(3.0 + t_0) - Float64(Float64(Float64(0.125 * Float64(3.0 - Float64(2.0 * v))) * Float64(Float64(Float64(w * w) * r) * r)) / Float64(1.0 - v))) <= -2000000.0) tmp = Float64(t_0 + Float64(3.0 - fma(Float64(Float64(Float64(r * r) * 0.375) * w), w, 4.5))); else tmp = Float64(Float64(Float64(2.0 / r) / r) + -1.5); end return tmp end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(3.0 + t$95$0), $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], -2000000.0], N[(t$95$0 + N[(3.0 - N[(N[(N[(N[(r * r), $MachinePrecision] * 0.375), $MachinePrecision] * w), $MachinePrecision] * w + 4.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 / r), $MachinePrecision] / r), $MachinePrecision] + -1.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;\left(3 + t\_0\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} \leq -2000000:\\
\;\;\;\;t\_0 + \left(3 - \mathsf{fma}\left(\left(\left(r \cdot r\right) \cdot 0.375\right) \cdot w, w, 4.5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{r}}{r} + -1.5\\
\end{array}
\end{array}
if (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -2e6Initial program 88.3%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.7%
Taylor expanded in v around 0
+-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6488.9
Applied rewrites88.9%
if -2e6 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) Initial program 82.6%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.8%
Taylor expanded in w around 0
Applied rewrites94.8%
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6494.8
Applied rewrites94.8%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<=
(-
(+ 3.0 t_0)
(/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v)))
-2000000.0)
(+ t_0 (fma (* (* (* -0.375 r) r) w) w -1.5))
(+ (/ (/ 2.0 r) r) -1.5))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) {
tmp = t_0 + fma((((-0.375 * r) * r) * w), w, -1.5);
} else {
tmp = ((2.0 / r) / r) + -1.5;
}
return tmp;
}
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (Float64(Float64(3.0 + t_0) - Float64(Float64(Float64(0.125 * Float64(3.0 - Float64(2.0 * v))) * Float64(Float64(Float64(w * w) * r) * r)) / Float64(1.0 - v))) <= -2000000.0) tmp = Float64(t_0 + fma(Float64(Float64(Float64(-0.375 * r) * r) * w), w, -1.5)); else tmp = Float64(Float64(Float64(2.0 / r) / r) + -1.5); end return tmp end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(3.0 + t$95$0), $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], -2000000.0], N[(t$95$0 + N[(N[(N[(N[(-0.375 * r), $MachinePrecision] * r), $MachinePrecision] * w), $MachinePrecision] * w + -1.5), $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 / r), $MachinePrecision] / r), $MachinePrecision] + -1.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;\left(3 + t\_0\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} \leq -2000000:\\
\;\;\;\;t\_0 + \mathsf{fma}\left(\left(\left(-0.375 \cdot r\right) \cdot r\right) \cdot w, w, -1.5\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{r}}{r} + -1.5\\
\end{array}
\end{array}
if (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -2e6Initial program 88.3%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.7%
lift-fma.f64N/A
Applied rewrites93.8%
Taylor expanded in v around 0
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6488.9
Applied rewrites88.9%
if -2e6 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) Initial program 82.6%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.8%
Taylor expanded in w around 0
Applied rewrites94.8%
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6494.8
Applied rewrites94.8%
Final simplification92.3%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<=
(-
(+ 3.0 t_0)
(/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v)))
-2000000.0)
(fma (* (* -0.375 (* r r)) w) w (- t_0 1.5))
(+ (/ (/ 2.0 r) r) -1.5))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) {
tmp = fma(((-0.375 * (r * r)) * w), w, (t_0 - 1.5));
} else {
tmp = ((2.0 / r) / r) + -1.5;
}
return tmp;
}
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (Float64(Float64(3.0 + t_0) - Float64(Float64(Float64(0.125 * Float64(3.0 - Float64(2.0 * v))) * Float64(Float64(Float64(w * w) * r) * r)) / Float64(1.0 - v))) <= -2000000.0) tmp = fma(Float64(Float64(-0.375 * Float64(r * r)) * w), w, Float64(t_0 - 1.5)); else tmp = Float64(Float64(Float64(2.0 / r) / r) + -1.5); end return tmp end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(3.0 + t$95$0), $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], -2000000.0], N[(N[(N[(-0.375 * N[(r * r), $MachinePrecision]), $MachinePrecision] * w), $MachinePrecision] * w + N[(t$95$0 - 1.5), $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 / r), $MachinePrecision] / r), $MachinePrecision] + -1.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;\left(3 + t\_0\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} \leq -2000000:\\
\;\;\;\;\mathsf{fma}\left(\left(-0.375 \cdot \left(r \cdot r\right)\right) \cdot w, w, t\_0 - 1.5\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{r}}{r} + -1.5\\
\end{array}
\end{array}
if (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -2e6Initial program 88.3%
Taylor expanded in v around 0
associate--l+N/A
sub-negN/A
+-commutativeN/A
associate-+r+N/A
sub-negN/A
lower-+.f64N/A
Applied rewrites64.0%
Taylor expanded in v around 0
Applied rewrites88.9%
if -2e6 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) Initial program 82.6%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.8%
Taylor expanded in w around 0
Applied rewrites94.8%
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6494.8
Applied rewrites94.8%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<=
(-
(+ 3.0 t_0)
(/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v)))
-2000000.0)
(* (* (* -0.375 (* r r)) w) w)
(- t_0 1.5))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) {
tmp = ((-0.375 * (r * r)) * w) * w;
} else {
tmp = t_0 - 1.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 = 2.0d0 / (r * r)
if (((3.0d0 + t_0) - (((0.125d0 * (3.0d0 - (2.0d0 * v))) * (((w * w) * r) * r)) / (1.0d0 - v))) <= (-2000000.0d0)) then
tmp = (((-0.375d0) * (r * r)) * w) * w
else
tmp = t_0 - 1.5d0
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 (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) {
tmp = ((-0.375 * (r * r)) * w) * w;
} else {
tmp = t_0 - 1.5;
}
return tmp;
}
def code(v, w, r): t_0 = 2.0 / (r * r) tmp = 0 if ((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0: tmp = ((-0.375 * (r * r)) * w) * w else: tmp = t_0 - 1.5 return tmp
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (Float64(Float64(3.0 + t_0) - Float64(Float64(Float64(0.125 * Float64(3.0 - Float64(2.0 * v))) * Float64(Float64(Float64(w * w) * r) * r)) / Float64(1.0 - v))) <= -2000000.0) tmp = Float64(Float64(Float64(-0.375 * Float64(r * r)) * w) * w); else tmp = Float64(t_0 - 1.5); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 2.0 / (r * r); tmp = 0.0; if (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) tmp = ((-0.375 * (r * r)) * w) * w; else tmp = t_0 - 1.5; end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(3.0 + t$95$0), $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], -2000000.0], N[(N[(N[(-0.375 * N[(r * r), $MachinePrecision]), $MachinePrecision] * w), $MachinePrecision] * w), $MachinePrecision], N[(t$95$0 - 1.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;\left(3 + t\_0\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} \leq -2000000:\\
\;\;\;\;\left(\left(-0.375 \cdot \left(r \cdot r\right)\right) \cdot w\right) \cdot w\\
\mathbf{else}:\\
\;\;\;\;t\_0 - 1.5\\
\end{array}
\end{array}
if (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -2e6Initial program 88.3%
Taylor expanded in v around 0
sub-negN/A
+-commutativeN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-inN/A
metadata-evalN/A
associate-+l+N/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower-fma.f64N/A
Applied rewrites83.3%
Taylor expanded in w around inf
Applied rewrites84.9%
if -2e6 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) Initial program 82.6%
Taylor expanded in w around 0
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f6494.8
Applied rewrites94.8%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<=
(-
(+ 3.0 t_0)
(/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v)))
-2000000.0)
(* (* (* (* -0.375 r) r) w) w)
(- t_0 1.5))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) {
tmp = (((-0.375 * r) * r) * w) * w;
} else {
tmp = t_0 - 1.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 = 2.0d0 / (r * r)
if (((3.0d0 + t_0) - (((0.125d0 * (3.0d0 - (2.0d0 * v))) * (((w * w) * r) * r)) / (1.0d0 - v))) <= (-2000000.0d0)) then
tmp = ((((-0.375d0) * r) * r) * w) * w
else
tmp = t_0 - 1.5d0
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 (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) {
tmp = (((-0.375 * r) * r) * w) * w;
} else {
tmp = t_0 - 1.5;
}
return tmp;
}
def code(v, w, r): t_0 = 2.0 / (r * r) tmp = 0 if ((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0: tmp = (((-0.375 * r) * r) * w) * w else: tmp = t_0 - 1.5 return tmp
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (Float64(Float64(3.0 + t_0) - Float64(Float64(Float64(0.125 * Float64(3.0 - Float64(2.0 * v))) * Float64(Float64(Float64(w * w) * r) * r)) / Float64(1.0 - v))) <= -2000000.0) tmp = Float64(Float64(Float64(Float64(-0.375 * r) * r) * w) * w); else tmp = Float64(t_0 - 1.5); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 2.0 / (r * r); tmp = 0.0; if (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) tmp = (((-0.375 * r) * r) * w) * w; else tmp = t_0 - 1.5; end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(3.0 + t$95$0), $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], -2000000.0], N[(N[(N[(N[(-0.375 * r), $MachinePrecision] * r), $MachinePrecision] * w), $MachinePrecision] * w), $MachinePrecision], N[(t$95$0 - 1.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;\left(3 + t\_0\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} \leq -2000000:\\
\;\;\;\;\left(\left(\left(-0.375 \cdot r\right) \cdot r\right) \cdot w\right) \cdot w\\
\mathbf{else}:\\
\;\;\;\;t\_0 - 1.5\\
\end{array}
\end{array}
if (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -2e6Initial program 88.3%
Taylor expanded in w around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
lower--.f6487.1
Applied rewrites87.1%
Taylor expanded in v around inf
Applied rewrites79.6%
Taylor expanded in v around 0
Applied rewrites84.9%
if -2e6 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) Initial program 82.6%
Taylor expanded in w around 0
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f6494.8
Applied rewrites94.8%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<=
(-
(+ 3.0 t_0)
(/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v)))
-2000000.0)
(* (* (* w (* -0.25 r)) r) w)
(- t_0 1.5))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) {
tmp = ((w * (-0.25 * r)) * r) * w;
} else {
tmp = t_0 - 1.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 = 2.0d0 / (r * r)
if (((3.0d0 + t_0) - (((0.125d0 * (3.0d0 - (2.0d0 * v))) * (((w * w) * r) * r)) / (1.0d0 - v))) <= (-2000000.0d0)) then
tmp = ((w * ((-0.25d0) * r)) * r) * w
else
tmp = t_0 - 1.5d0
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 (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) {
tmp = ((w * (-0.25 * r)) * r) * w;
} else {
tmp = t_0 - 1.5;
}
return tmp;
}
def code(v, w, r): t_0 = 2.0 / (r * r) tmp = 0 if ((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0: tmp = ((w * (-0.25 * r)) * r) * w else: tmp = t_0 - 1.5 return tmp
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (Float64(Float64(3.0 + t_0) - Float64(Float64(Float64(0.125 * Float64(3.0 - Float64(2.0 * v))) * Float64(Float64(Float64(w * w) * r) * r)) / Float64(1.0 - v))) <= -2000000.0) tmp = Float64(Float64(Float64(w * Float64(-0.25 * r)) * r) * w); else tmp = Float64(t_0 - 1.5); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 2.0 / (r * r); tmp = 0.0; if (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) tmp = ((w * (-0.25 * r)) * r) * w; else tmp = t_0 - 1.5; end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(3.0 + t$95$0), $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], -2000000.0], N[(N[(N[(w * N[(-0.25 * r), $MachinePrecision]), $MachinePrecision] * r), $MachinePrecision] * w), $MachinePrecision], N[(t$95$0 - 1.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;\left(3 + t\_0\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} \leq -2000000:\\
\;\;\;\;\left(\left(w \cdot \left(-0.25 \cdot r\right)\right) \cdot r\right) \cdot w\\
\mathbf{else}:\\
\;\;\;\;t\_0 - 1.5\\
\end{array}
\end{array}
if (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -2e6Initial program 88.3%
Taylor expanded in w around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
lower--.f6487.1
Applied rewrites87.1%
Taylor expanded in v around inf
Applied rewrites79.6%
Applied rewrites80.8%
if -2e6 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) Initial program 82.6%
Taylor expanded in w around 0
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f6494.8
Applied rewrites94.8%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<=
(-
(+ 3.0 t_0)
(/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v)))
-2000000.0)
(* (* (* -0.25 (* r r)) w) w)
(- t_0 1.5))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) {
tmp = ((-0.25 * (r * r)) * w) * w;
} else {
tmp = t_0 - 1.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 = 2.0d0 / (r * r)
if (((3.0d0 + t_0) - (((0.125d0 * (3.0d0 - (2.0d0 * v))) * (((w * w) * r) * r)) / (1.0d0 - v))) <= (-2000000.0d0)) then
tmp = (((-0.25d0) * (r * r)) * w) * w
else
tmp = t_0 - 1.5d0
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 (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) {
tmp = ((-0.25 * (r * r)) * w) * w;
} else {
tmp = t_0 - 1.5;
}
return tmp;
}
def code(v, w, r): t_0 = 2.0 / (r * r) tmp = 0 if ((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0: tmp = ((-0.25 * (r * r)) * w) * w else: tmp = t_0 - 1.5 return tmp
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (Float64(Float64(3.0 + t_0) - Float64(Float64(Float64(0.125 * Float64(3.0 - Float64(2.0 * v))) * Float64(Float64(Float64(w * w) * r) * r)) / Float64(1.0 - v))) <= -2000000.0) tmp = Float64(Float64(Float64(-0.25 * Float64(r * r)) * w) * w); else tmp = Float64(t_0 - 1.5); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 2.0 / (r * r); tmp = 0.0; if (((3.0 + t_0) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) <= -2000000.0) tmp = ((-0.25 * (r * r)) * w) * w; else tmp = t_0 - 1.5; end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(3.0 + t$95$0), $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], -2000000.0], N[(N[(N[(-0.25 * N[(r * r), $MachinePrecision]), $MachinePrecision] * w), $MachinePrecision] * w), $MachinePrecision], N[(t$95$0 - 1.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;\left(3 + t\_0\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} \leq -2000000:\\
\;\;\;\;\left(\left(-0.25 \cdot \left(r \cdot r\right)\right) \cdot w\right) \cdot w\\
\mathbf{else}:\\
\;\;\;\;t\_0 - 1.5\\
\end{array}
\end{array}
if (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) < -2e6Initial program 88.3%
Taylor expanded in w around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
lower--.f6487.1
Applied rewrites87.1%
Taylor expanded in v around inf
Applied rewrites79.6%
if -2e6 < (-.f64 (+.f64 #s(literal 3 binary64) (/.f64 #s(literal 2 binary64) (*.f64 r r))) (/.f64 (*.f64 (*.f64 #s(literal 1/8 binary64) (-.f64 #s(literal 3 binary64) (*.f64 #s(literal 2 binary64) v))) (*.f64 (*.f64 (*.f64 w w) r) r)) (-.f64 #s(literal 1 binary64) v))) Initial program 82.6%
Taylor expanded in w around 0
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f6494.8
Applied rewrites94.8%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))) (t_1 (* (* w r) (* w r))))
(if (<= v -1720.0)
(+ t_0 (fma t_1 (- (/ 0.125 v) 0.25) -1.5))
(- (- (+ 3.0 t_0) (/ (* t_1 (fma -0.25 v 0.375)) (- 1.0 v))) 4.5))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double t_1 = (w * r) * (w * r);
double tmp;
if (v <= -1720.0) {
tmp = t_0 + fma(t_1, ((0.125 / v) - 0.25), -1.5);
} else {
tmp = ((3.0 + t_0) - ((t_1 * fma(-0.25, v, 0.375)) / (1.0 - v))) - 4.5;
}
return tmp;
}
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) t_1 = Float64(Float64(w * r) * Float64(w * r)) tmp = 0.0 if (v <= -1720.0) tmp = Float64(t_0 + fma(t_1, Float64(Float64(0.125 / v) - 0.25), -1.5)); else tmp = Float64(Float64(Float64(3.0 + t_0) - Float64(Float64(t_1 * fma(-0.25, v, 0.375)) / Float64(1.0 - v))) - 4.5); end return tmp end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(w * r), $MachinePrecision] * N[(w * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[v, -1720.0], N[(t$95$0 + N[(t$95$1 * N[(N[(0.125 / v), $MachinePrecision] - 0.25), $MachinePrecision] + -1.5), $MachinePrecision]), $MachinePrecision], N[(N[(N[(3.0 + t$95$0), $MachinePrecision] - N[(N[(t$95$1 * N[(-0.25 * v + 0.375), $MachinePrecision]), $MachinePrecision] / N[(1.0 - v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
t_1 := \left(w \cdot r\right) \cdot \left(w \cdot r\right)\\
\mathbf{if}\;v \leq -1720:\\
\;\;\;\;t\_0 + \mathsf{fma}\left(t\_1, \frac{0.125}{v} - 0.25, -1.5\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(3 + t\_0\right) - \frac{t\_1 \cdot \mathsf{fma}\left(-0.25, v, 0.375\right)}{1 - v}\right) - 4.5\\
\end{array}
\end{array}
if v < -1720Initial program 80.5%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.7%
Taylor expanded in v around inf
Applied rewrites99.2%
if -1720 < v Initial program 87.5%
Taylor expanded in v around inf
lower-*.f6471.6
Applied rewrites71.6%
Taylor expanded in v around 0
+-commutativeN/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
distribute-lft-outN/A
lower-*.f64N/A
*-commutativeN/A
unpow2N/A
unpow2N/A
unswap-sqrN/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f6498.6
Applied rewrites98.6%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<= v -5.8e+184)
(+ t_0 (fma (* (* w r) (* w r)) (- (/ 0.125 v) 0.25) -1.5))
(+
t_0
(-
3.0
(fma (* r (* w (fma -0.25 v 0.375))) (* w (/ r (- 1.0 v))) 4.5))))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (v <= -5.8e+184) {
tmp = t_0 + fma(((w * r) * (w * r)), ((0.125 / v) - 0.25), -1.5);
} else {
tmp = t_0 + (3.0 - fma((r * (w * fma(-0.25, v, 0.375))), (w * (r / (1.0 - v))), 4.5));
}
return tmp;
}
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (v <= -5.8e+184) tmp = Float64(t_0 + fma(Float64(Float64(w * r) * Float64(w * r)), Float64(Float64(0.125 / v) - 0.25), -1.5)); else tmp = Float64(t_0 + Float64(3.0 - fma(Float64(r * Float64(w * fma(-0.25, v, 0.375))), Float64(w * Float64(r / Float64(1.0 - v))), 4.5))); end return tmp end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[v, -5.8e+184], N[(t$95$0 + N[(N[(N[(w * r), $MachinePrecision] * N[(w * r), $MachinePrecision]), $MachinePrecision] * N[(N[(0.125 / v), $MachinePrecision] - 0.25), $MachinePrecision] + -1.5), $MachinePrecision]), $MachinePrecision], N[(t$95$0 + N[(3.0 - N[(N[(r * N[(w * N[(-0.25 * v + 0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(w * N[(r / N[(1.0 - v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 4.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;v \leq -5.8 \cdot 10^{+184}:\\
\;\;\;\;t\_0 + \mathsf{fma}\left(\left(w \cdot r\right) \cdot \left(w \cdot r\right), \frac{0.125}{v} - 0.25, -1.5\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0 + \left(3 - \mathsf{fma}\left(r \cdot \left(w \cdot \mathsf{fma}\left(-0.25, v, 0.375\right)\right), w \cdot \frac{r}{1 - v}, 4.5\right)\right)\\
\end{array}
\end{array}
if v < -5.7999999999999998e184Initial program 73.7%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.7%
Taylor expanded in v around inf
Applied rewrites99.7%
if -5.7999999999999998e184 < v Initial program 86.9%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.8%
lift-fma.f64N/A
Applied rewrites96.7%
lift-fma.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
associate-*l*N/A
lower-fma.f64N/A
Applied rewrites98.0%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (or (<= v -0.72) (not (<= v 0.98)))
(+ t_0 (fma (* (* w r) (* w r)) (- (/ 0.125 v) 0.25) -1.5))
(+ (fma (* (* (* r w) r) w) (- (* -0.125 v) 0.375) -1.5) t_0))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if ((v <= -0.72) || !(v <= 0.98)) {
tmp = t_0 + fma(((w * r) * (w * r)), ((0.125 / v) - 0.25), -1.5);
} else {
tmp = fma((((r * w) * r) * w), ((-0.125 * v) - 0.375), -1.5) + t_0;
}
return tmp;
}
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if ((v <= -0.72) || !(v <= 0.98)) tmp = Float64(t_0 + fma(Float64(Float64(w * r) * Float64(w * r)), Float64(Float64(0.125 / v) - 0.25), -1.5)); else tmp = Float64(fma(Float64(Float64(Float64(r * w) * r) * w), Float64(Float64(-0.125 * v) - 0.375), -1.5) + t_0); end return tmp end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[v, -0.72], N[Not[LessEqual[v, 0.98]], $MachinePrecision]], N[(t$95$0 + N[(N[(N[(w * r), $MachinePrecision] * N[(w * r), $MachinePrecision]), $MachinePrecision] * N[(N[(0.125 / v), $MachinePrecision] - 0.25), $MachinePrecision] + -1.5), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[(r * w), $MachinePrecision] * r), $MachinePrecision] * w), $MachinePrecision] * N[(N[(-0.125 * v), $MachinePrecision] - 0.375), $MachinePrecision] + -1.5), $MachinePrecision] + t$95$0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;v \leq -0.72 \lor \neg \left(v \leq 0.98\right):\\
\;\;\;\;t\_0 + \mathsf{fma}\left(\left(w \cdot r\right) \cdot \left(w \cdot r\right), \frac{0.125}{v} - 0.25, -1.5\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(\left(r \cdot w\right) \cdot r\right) \cdot w, -0.125 \cdot v - 0.375, -1.5\right) + t\_0\\
\end{array}
\end{array}
if v < -0.71999999999999997 or 0.97999999999999998 < v Initial program 82.4%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.8%
Taylor expanded in v around inf
Applied rewrites99.5%
if -0.71999999999999997 < v < 0.97999999999999998Initial program 88.1%
Taylor expanded in v around 0
associate--l+N/A
sub-negN/A
+-commutativeN/A
associate-+r+N/A
sub-negN/A
lower-+.f64N/A
Applied rewrites93.2%
Applied rewrites99.8%
Final simplification99.6%
(FPCore (v w r)
:precision binary64
(if (<= r 2.8e+151)
(+ (/ 2.0 (* r r)) (- 3.0 (fma (* (* (* r r) 0.375) w) w 4.5)))
(if (<= r 2.55e+167)
(- (- 3.0 (/ (* (* -0.25 v) (* (* (* w w) r) r)) (- 1.0 v))) 4.5)
(* (/ (* w (* (* (* w (fma v -2.0 3.0)) -0.125) r)) (- 1.0 v)) r))))
double code(double v, double w, double r) {
double tmp;
if (r <= 2.8e+151) {
tmp = (2.0 / (r * r)) + (3.0 - fma((((r * r) * 0.375) * w), w, 4.5));
} else if (r <= 2.55e+167) {
tmp = (3.0 - (((-0.25 * v) * (((w * w) * r) * r)) / (1.0 - v))) - 4.5;
} else {
tmp = ((w * (((w * fma(v, -2.0, 3.0)) * -0.125) * r)) / (1.0 - v)) * r;
}
return tmp;
}
function code(v, w, r) tmp = 0.0 if (r <= 2.8e+151) tmp = Float64(Float64(2.0 / Float64(r * r)) + Float64(3.0 - fma(Float64(Float64(Float64(r * r) * 0.375) * w), w, 4.5))); elseif (r <= 2.55e+167) tmp = Float64(Float64(3.0 - Float64(Float64(Float64(-0.25 * v) * Float64(Float64(Float64(w * w) * r) * r)) / Float64(1.0 - v))) - 4.5); else tmp = Float64(Float64(Float64(w * Float64(Float64(Float64(w * fma(v, -2.0, 3.0)) * -0.125) * r)) / Float64(1.0 - v)) * r); end return tmp end
code[v_, w_, r_] := If[LessEqual[r, 2.8e+151], N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + N[(3.0 - N[(N[(N[(N[(r * r), $MachinePrecision] * 0.375), $MachinePrecision] * w), $MachinePrecision] * w + 4.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[r, 2.55e+167], N[(N[(3.0 - N[(N[(N[(-0.25 * v), $MachinePrecision] * N[(N[(N[(w * w), $MachinePrecision] * r), $MachinePrecision] * r), $MachinePrecision]), $MachinePrecision] / N[(1.0 - v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision], N[(N[(N[(w * N[(N[(N[(w * N[(v * -2.0 + 3.0), $MachinePrecision]), $MachinePrecision] * -0.125), $MachinePrecision] * r), $MachinePrecision]), $MachinePrecision] / N[(1.0 - v), $MachinePrecision]), $MachinePrecision] * r), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 2.8 \cdot 10^{+151}:\\
\;\;\;\;\frac{2}{r \cdot r} + \left(3 - \mathsf{fma}\left(\left(\left(r \cdot r\right) \cdot 0.375\right) \cdot w, w, 4.5\right)\right)\\
\mathbf{elif}\;r \leq 2.55 \cdot 10^{+167}:\\
\;\;\;\;\left(3 - \frac{\left(-0.25 \cdot v\right) \cdot \left(\left(\left(w \cdot w\right) \cdot r\right) \cdot r\right)}{1 - v}\right) - 4.5\\
\mathbf{else}:\\
\;\;\;\;\frac{w \cdot \left(\left(\left(w \cdot \mathsf{fma}\left(v, -2, 3\right)\right) \cdot -0.125\right) \cdot r\right)}{1 - v} \cdot r\\
\end{array}
\end{array}
if r < 2.79999999999999987e151Initial program 85.9%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.8%
Taylor expanded in v around 0
+-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6492.7
Applied rewrites92.7%
if 2.79999999999999987e151 < r < 2.55000000000000002e167Initial program 100.0%
Taylor expanded in v around inf
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in r around inf
Applied rewrites100.0%
if 2.55000000000000002e167 < r Initial program 74.9%
Taylor expanded in w around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
lower--.f6463.8
Applied rewrites63.8%
Applied rewrites74.0%
Applied rewrites97.7%
(FPCore (v w r) :precision binary64 (if (<= r 1.3e+31) (+ (/ 2.0 (* r r)) (fma (* (* (* -0.375 r) r) w) w -1.5)) (- (- 3.0 (/ (* (* (* w r) (* w r)) (fma -0.25 v 0.375)) (- 1.0 v))) 4.5)))
double code(double v, double w, double r) {
double tmp;
if (r <= 1.3e+31) {
tmp = (2.0 / (r * r)) + fma((((-0.375 * r) * r) * w), w, -1.5);
} else {
tmp = (3.0 - ((((w * r) * (w * r)) * fma(-0.25, v, 0.375)) / (1.0 - v))) - 4.5;
}
return tmp;
}
function code(v, w, r) tmp = 0.0 if (r <= 1.3e+31) tmp = Float64(Float64(2.0 / Float64(r * r)) + fma(Float64(Float64(Float64(-0.375 * r) * r) * w), w, -1.5)); else tmp = Float64(Float64(3.0 - Float64(Float64(Float64(Float64(w * r) * Float64(w * r)) * fma(-0.25, v, 0.375)) / Float64(1.0 - v))) - 4.5); end return tmp end
code[v_, w_, r_] := If[LessEqual[r, 1.3e+31], N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(N[(-0.375 * r), $MachinePrecision] * r), $MachinePrecision] * w), $MachinePrecision] * w + -1.5), $MachinePrecision]), $MachinePrecision], N[(N[(3.0 - N[(N[(N[(N[(w * r), $MachinePrecision] * N[(w * r), $MachinePrecision]), $MachinePrecision] * N[(-0.25 * v + 0.375), $MachinePrecision]), $MachinePrecision] / N[(1.0 - v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 1.3 \cdot 10^{+31}:\\
\;\;\;\;\frac{2}{r \cdot r} + \mathsf{fma}\left(\left(\left(-0.375 \cdot r\right) \cdot r\right) \cdot w, w, -1.5\right)\\
\mathbf{else}:\\
\;\;\;\;\left(3 - \frac{\left(\left(w \cdot r\right) \cdot \left(w \cdot r\right)\right) \cdot \mathsf{fma}\left(-0.25, v, 0.375\right)}{1 - v}\right) - 4.5\\
\end{array}
\end{array}
if r < 1.3e31Initial program 84.0%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.8%
lift-fma.f64N/A
Applied rewrites90.9%
Taylor expanded in v around 0
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6491.7
Applied rewrites91.7%
if 1.3e31 < r Initial program 88.4%
Taylor expanded in v around inf
lower-*.f6466.6
Applied rewrites66.6%
Taylor expanded in v around 0
+-commutativeN/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
distribute-lft-outN/A
lower-*.f64N/A
*-commutativeN/A
unpow2N/A
unpow2N/A
unswap-sqrN/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f6498.1
Applied rewrites98.1%
Taylor expanded in r around inf
Applied rewrites98.1%
Final simplification93.2%
(FPCore (v w r) :precision binary64 (if (<= r 9e+146) (+ (/ 2.0 (* r r)) (fma (* (* (* -0.375 r) r) w) w -1.5)) (* (/ (* w (* (* (* w (fma v -2.0 3.0)) -0.125) r)) (- 1.0 v)) r)))
double code(double v, double w, double r) {
double tmp;
if (r <= 9e+146) {
tmp = (2.0 / (r * r)) + fma((((-0.375 * r) * r) * w), w, -1.5);
} else {
tmp = ((w * (((w * fma(v, -2.0, 3.0)) * -0.125) * r)) / (1.0 - v)) * r;
}
return tmp;
}
function code(v, w, r) tmp = 0.0 if (r <= 9e+146) tmp = Float64(Float64(2.0 / Float64(r * r)) + fma(Float64(Float64(Float64(-0.375 * r) * r) * w), w, -1.5)); else tmp = Float64(Float64(Float64(w * Float64(Float64(Float64(w * fma(v, -2.0, 3.0)) * -0.125) * r)) / Float64(1.0 - v)) * r); end return tmp end
code[v_, w_, r_] := If[LessEqual[r, 9e+146], N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(N[(-0.375 * r), $MachinePrecision] * r), $MachinePrecision] * w), $MachinePrecision] * w + -1.5), $MachinePrecision]), $MachinePrecision], N[(N[(N[(w * N[(N[(N[(w * N[(v * -2.0 + 3.0), $MachinePrecision]), $MachinePrecision] * -0.125), $MachinePrecision] * r), $MachinePrecision]), $MachinePrecision] / N[(1.0 - v), $MachinePrecision]), $MachinePrecision] * r), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 9 \cdot 10^{+146}:\\
\;\;\;\;\frac{2}{r \cdot r} + \mathsf{fma}\left(\left(\left(-0.375 \cdot r\right) \cdot r\right) \cdot w, w, -1.5\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{w \cdot \left(\left(\left(w \cdot \mathsf{fma}\left(v, -2, 3\right)\right) \cdot -0.125\right) \cdot r\right)}{1 - v} \cdot r\\
\end{array}
\end{array}
if r < 9.00000000000000051e146Initial program 85.7%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
+-commutativeN/A
associate--l+N/A
lower-+.f64N/A
lower--.f64N/A
Applied rewrites99.8%
lift-fma.f64N/A
Applied rewrites91.9%
Taylor expanded in v around 0
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6492.6
Applied rewrites92.6%
if 9.00000000000000051e146 < r Initial program 80.6%
Taylor expanded in w around inf
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
lower--.f6466.3
Applied rewrites66.3%
Applied rewrites74.3%
Applied rewrites92.8%
Final simplification92.6%
(FPCore (v w r) :precision binary64 (- (/ 2.0 (* r r)) 1.5))
double code(double v, double w, double r) {
return (2.0 / (r * r)) - 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)) - 1.5d0
end function
public static double code(double v, double w, double r) {
return (2.0 / (r * r)) - 1.5;
}
def code(v, w, r): return (2.0 / (r * r)) - 1.5
function code(v, w, r) return Float64(Float64(2.0 / Float64(r * r)) - 1.5) end
function tmp = code(v, w, r) tmp = (2.0 / (r * r)) - 1.5; end
code[v_, w_, r_] := N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] - 1.5), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{r \cdot r} - 1.5
\end{array}
Initial program 85.0%
Taylor expanded in w around 0
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f6457.1
Applied rewrites57.1%
(FPCore (v w r) :precision binary64 (/ 2.0 (* r r)))
double code(double v, double w, double r) {
return 2.0 / (r * r);
}
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)
end function
public static double code(double v, double w, double r) {
return 2.0 / (r * r);
}
def code(v, w, r): return 2.0 / (r * r)
function code(v, w, r) return Float64(2.0 / Float64(r * r)) end
function tmp = code(v, w, r) tmp = 2.0 / (r * r); end
code[v_, w_, r_] := N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{r \cdot r}
\end{array}
Initial program 85.0%
Taylor expanded in r around 0
lower-/.f64N/A
unpow2N/A
lower-*.f6443.1
Applied rewrites43.1%
herbie shell --seed 2024318
(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))