
(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 14 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
(let* ((t_0 (+ (/ 2.0 (* r r)) -1.5))
(t_1 (+ t_0 (* (+ 0.375 (* v -0.25)) (* (* r w) (/ (* r w) v))))))
(if (<= v -2.0)
t_1
(if (<= v 9.5e-6) (+ t_0 (* (* r w) (* w (* r -0.375)))) t_1))))
double code(double v, double w, double r) {
double t_0 = (2.0 / (r * r)) + -1.5;
double t_1 = t_0 + ((0.375 + (v * -0.25)) * ((r * w) * ((r * w) / v)));
double tmp;
if (v <= -2.0) {
tmp = t_1;
} else if (v <= 9.5e-6) {
tmp = t_0 + ((r * w) * (w * (r * -0.375)));
} else {
tmp = t_1;
}
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)) + (-1.5d0)
t_1 = t_0 + ((0.375d0 + (v * (-0.25d0))) * ((r * w) * ((r * w) / v)))
if (v <= (-2.0d0)) then
tmp = t_1
else if (v <= 9.5d-6) then
tmp = t_0 + ((r * w) * (w * (r * (-0.375d0))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double t_0 = (2.0 / (r * r)) + -1.5;
double t_1 = t_0 + ((0.375 + (v * -0.25)) * ((r * w) * ((r * w) / v)));
double tmp;
if (v <= -2.0) {
tmp = t_1;
} else if (v <= 9.5e-6) {
tmp = t_0 + ((r * w) * (w * (r * -0.375)));
} else {
tmp = t_1;
}
return tmp;
}
def code(v, w, r): t_0 = (2.0 / (r * r)) + -1.5 t_1 = t_0 + ((0.375 + (v * -0.25)) * ((r * w) * ((r * w) / v))) tmp = 0 if v <= -2.0: tmp = t_1 elif v <= 9.5e-6: tmp = t_0 + ((r * w) * (w * (r * -0.375))) else: tmp = t_1 return tmp
function code(v, w, r) t_0 = Float64(Float64(2.0 / Float64(r * r)) + -1.5) t_1 = Float64(t_0 + Float64(Float64(0.375 + Float64(v * -0.25)) * Float64(Float64(r * w) * Float64(Float64(r * w) / v)))) tmp = 0.0 if (v <= -2.0) tmp = t_1; elseif (v <= 9.5e-6) tmp = Float64(t_0 + Float64(Float64(r * w) * Float64(w * Float64(r * -0.375)))); else tmp = t_1; end return tmp end
function tmp_2 = code(v, w, r) t_0 = (2.0 / (r * r)) + -1.5; t_1 = t_0 + ((0.375 + (v * -0.25)) * ((r * w) * ((r * w) / v))); tmp = 0.0; if (v <= -2.0) tmp = t_1; elseif (v <= 9.5e-6) tmp = t_0 + ((r * w) * (w * (r * -0.375))); else tmp = t_1; end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 + N[(N[(0.375 + N[(v * -0.25), $MachinePrecision]), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[v, -2.0], t$95$1, If[LessEqual[v, 9.5e-6], N[(t$95$0 + N[(N[(r * w), $MachinePrecision] * N[(w * N[(r * -0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r} + -1.5\\
t_1 := t\_0 + \left(0.375 + v \cdot -0.25\right) \cdot \left(\left(r \cdot w\right) \cdot \frac{r \cdot w}{v}\right)\\
\mathbf{if}\;v \leq -2:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;v \leq 9.5 \cdot 10^{-6}:\\
\;\;\;\;t\_0 + \left(r \cdot w\right) \cdot \left(w \cdot \left(r \cdot -0.375\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if v < -2 or 9.5000000000000005e-6 < v Initial program 85.2%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified89.9%
associate-/l*N/A
*-commutativeN/A
associate-*r*N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6498.5%
Applied egg-rr98.5%
Taylor expanded in v around inf
/-lowering-/.f64N/A
*-lowering-*.f6498.6%
Simplified98.6%
if -2 < v < 9.5000000000000005e-6Initial program 86.8%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified86.8%
Taylor expanded in v around 0
associate--l+N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6483.3%
Simplified83.3%
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.8%
Applied egg-rr99.8%
Final simplification99.2%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (+ (/ 2.0 (* r r)) -1.5))
(t_1 (+ t_0 (* (* v -0.25) (* (* r w) (* r (/ w (+ v -1.0))))))))
(if (<= v -3.2)
t_1
(if (<= v 9.5e-6) (+ t_0 (* (* r w) (* w (* r -0.375)))) t_1))))
double code(double v, double w, double r) {
double t_0 = (2.0 / (r * r)) + -1.5;
double t_1 = t_0 + ((v * -0.25) * ((r * w) * (r * (w / (v + -1.0)))));
double tmp;
if (v <= -3.2) {
tmp = t_1;
} else if (v <= 9.5e-6) {
tmp = t_0 + ((r * w) * (w * (r * -0.375)));
} else {
tmp = t_1;
}
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)) + (-1.5d0)
t_1 = t_0 + ((v * (-0.25d0)) * ((r * w) * (r * (w / (v + (-1.0d0))))))
if (v <= (-3.2d0)) then
tmp = t_1
else if (v <= 9.5d-6) then
tmp = t_0 + ((r * w) * (w * (r * (-0.375d0))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double t_0 = (2.0 / (r * r)) + -1.5;
double t_1 = t_0 + ((v * -0.25) * ((r * w) * (r * (w / (v + -1.0)))));
double tmp;
if (v <= -3.2) {
tmp = t_1;
} else if (v <= 9.5e-6) {
tmp = t_0 + ((r * w) * (w * (r * -0.375)));
} else {
tmp = t_1;
}
return tmp;
}
def code(v, w, r): t_0 = (2.0 / (r * r)) + -1.5 t_1 = t_0 + ((v * -0.25) * ((r * w) * (r * (w / (v + -1.0))))) tmp = 0 if v <= -3.2: tmp = t_1 elif v <= 9.5e-6: tmp = t_0 + ((r * w) * (w * (r * -0.375))) else: tmp = t_1 return tmp
function code(v, w, r) t_0 = Float64(Float64(2.0 / Float64(r * r)) + -1.5) t_1 = Float64(t_0 + Float64(Float64(v * -0.25) * Float64(Float64(r * w) * Float64(r * Float64(w / Float64(v + -1.0)))))) tmp = 0.0 if (v <= -3.2) tmp = t_1; elseif (v <= 9.5e-6) tmp = Float64(t_0 + Float64(Float64(r * w) * Float64(w * Float64(r * -0.375)))); else tmp = t_1; end return tmp end
function tmp_2 = code(v, w, r) t_0 = (2.0 / (r * r)) + -1.5; t_1 = t_0 + ((v * -0.25) * ((r * w) * (r * (w / (v + -1.0))))); tmp = 0.0; if (v <= -3.2) tmp = t_1; elseif (v <= 9.5e-6) tmp = t_0 + ((r * w) * (w * (r * -0.375))); else tmp = t_1; end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 + N[(N[(v * -0.25), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] * N[(r * N[(w / N[(v + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[v, -3.2], t$95$1, If[LessEqual[v, 9.5e-6], N[(t$95$0 + N[(N[(r * w), $MachinePrecision] * N[(w * N[(r * -0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r} + -1.5\\
t_1 := t\_0 + \left(v \cdot -0.25\right) \cdot \left(\left(r \cdot w\right) \cdot \left(r \cdot \frac{w}{v + -1}\right)\right)\\
\mathbf{if}\;v \leq -3.2:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;v \leq 9.5 \cdot 10^{-6}:\\
\;\;\;\;t\_0 + \left(r \cdot w\right) \cdot \left(w \cdot \left(r \cdot -0.375\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if v < -3.2000000000000002 or 9.5000000000000005e-6 < v Initial program 85.2%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified89.9%
associate-/l*N/A
*-commutativeN/A
associate-*r*N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6498.5%
Applied egg-rr98.5%
Taylor expanded in v around inf
*-commutativeN/A
*-lowering-*.f6497.2%
Simplified97.2%
if -3.2000000000000002 < v < 9.5000000000000005e-6Initial program 86.8%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified86.8%
Taylor expanded in v around 0
associate--l+N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6483.3%
Simplified83.3%
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.8%
Applied egg-rr99.8%
Final simplification98.4%
(FPCore (v w r)
:precision binary64
(if (<= r 7e-5)
(+ (/ 2.0 (* r r)) -1.5)
(if (<= r 1e+143)
(* (* r r) (+ (* -0.375 (* w w)) (/ -1.5 (* r r))))
(* (* r w) (* (* r w) (/ (+ 0.375 (* v -0.25)) (+ v -1.0)))))))
double code(double v, double w, double r) {
double tmp;
if (r <= 7e-5) {
tmp = (2.0 / (r * r)) + -1.5;
} else if (r <= 1e+143) {
tmp = (r * r) * ((-0.375 * (w * w)) + (-1.5 / (r * r)));
} else {
tmp = (r * w) * ((r * w) * ((0.375 + (v * -0.25)) / (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) :: tmp
if (r <= 7d-5) then
tmp = (2.0d0 / (r * r)) + (-1.5d0)
else if (r <= 1d+143) then
tmp = (r * r) * (((-0.375d0) * (w * w)) + ((-1.5d0) / (r * r)))
else
tmp = (r * w) * ((r * w) * ((0.375d0 + (v * (-0.25d0))) / (v + (-1.0d0))))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 7e-5) {
tmp = (2.0 / (r * r)) + -1.5;
} else if (r <= 1e+143) {
tmp = (r * r) * ((-0.375 * (w * w)) + (-1.5 / (r * r)));
} else {
tmp = (r * w) * ((r * w) * ((0.375 + (v * -0.25)) / (v + -1.0)));
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 7e-5: tmp = (2.0 / (r * r)) + -1.5 elif r <= 1e+143: tmp = (r * r) * ((-0.375 * (w * w)) + (-1.5 / (r * r))) else: tmp = (r * w) * ((r * w) * ((0.375 + (v * -0.25)) / (v + -1.0))) return tmp
function code(v, w, r) tmp = 0.0 if (r <= 7e-5) tmp = Float64(Float64(2.0 / Float64(r * r)) + -1.5); elseif (r <= 1e+143) tmp = Float64(Float64(r * r) * Float64(Float64(-0.375 * Float64(w * w)) + Float64(-1.5 / Float64(r * r)))); else tmp = Float64(Float64(r * w) * Float64(Float64(r * w) * Float64(Float64(0.375 + Float64(v * -0.25)) / Float64(v + -1.0)))); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 7e-5) tmp = (2.0 / (r * r)) + -1.5; elseif (r <= 1e+143) tmp = (r * r) * ((-0.375 * (w * w)) + (-1.5 / (r * r))); else tmp = (r * w) * ((r * w) * ((0.375 + (v * -0.25)) / (v + -1.0))); end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 7e-5], N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision], If[LessEqual[r, 1e+143], N[(N[(r * r), $MachinePrecision] * N[(N[(-0.375 * N[(w * w), $MachinePrecision]), $MachinePrecision] + N[(-1.5 / N[(r * r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(r * w), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] * N[(N[(0.375 + N[(v * -0.25), $MachinePrecision]), $MachinePrecision] / N[(v + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 7 \cdot 10^{-5}:\\
\;\;\;\;\frac{2}{r \cdot r} + -1.5\\
\mathbf{elif}\;r \leq 10^{+143}:\\
\;\;\;\;\left(r \cdot r\right) \cdot \left(-0.375 \cdot \left(w \cdot w\right) + \frac{-1.5}{r \cdot r}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(r \cdot w\right) \cdot \left(\left(r \cdot w\right) \cdot \frac{0.375 + v \cdot -0.25}{v + -1}\right)\\
\end{array}
\end{array}
if r < 6.9999999999999994e-5Initial program 84.8%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified86.1%
Taylor expanded in w around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6464.8%
Simplified64.8%
if 6.9999999999999994e-5 < r < 1e143Initial program 89.6%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified99.7%
Taylor expanded in v around 0
associate--l+N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6491.0%
Simplified91.0%
Taylor expanded in r around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6491.0%
Simplified91.0%
if 1e143 < r Initial program 90.4%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified93.2%
Taylor expanded in r around inf
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6477.6%
Simplified77.6%
associate-/l*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6485.4%
Applied egg-rr85.4%
Final simplification70.1%
(FPCore (v w r)
:precision binary64
(if (<= r 7.8e-5)
(+ (/ 2.0 (* r r)) -1.5)
(if (<= r 2.7e+143)
(* (* r r) (+ (* -0.375 (* w w)) (/ -1.5 (* r r))))
(* (* r w) (* w (/ (* (+ 0.375 (* v -0.25)) r) (+ v -1.0)))))))
double code(double v, double w, double r) {
double tmp;
if (r <= 7.8e-5) {
tmp = (2.0 / (r * r)) + -1.5;
} else if (r <= 2.7e+143) {
tmp = (r * r) * ((-0.375 * (w * w)) + (-1.5 / (r * r)));
} else {
tmp = (r * w) * (w * (((0.375 + (v * -0.25)) * r) / (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) :: tmp
if (r <= 7.8d-5) then
tmp = (2.0d0 / (r * r)) + (-1.5d0)
else if (r <= 2.7d+143) then
tmp = (r * r) * (((-0.375d0) * (w * w)) + ((-1.5d0) / (r * r)))
else
tmp = (r * w) * (w * (((0.375d0 + (v * (-0.25d0))) * r) / (v + (-1.0d0))))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 7.8e-5) {
tmp = (2.0 / (r * r)) + -1.5;
} else if (r <= 2.7e+143) {
tmp = (r * r) * ((-0.375 * (w * w)) + (-1.5 / (r * r)));
} else {
tmp = (r * w) * (w * (((0.375 + (v * -0.25)) * r) / (v + -1.0)));
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 7.8e-5: tmp = (2.0 / (r * r)) + -1.5 elif r <= 2.7e+143: tmp = (r * r) * ((-0.375 * (w * w)) + (-1.5 / (r * r))) else: tmp = (r * w) * (w * (((0.375 + (v * -0.25)) * r) / (v + -1.0))) return tmp
function code(v, w, r) tmp = 0.0 if (r <= 7.8e-5) tmp = Float64(Float64(2.0 / Float64(r * r)) + -1.5); elseif (r <= 2.7e+143) tmp = Float64(Float64(r * r) * Float64(Float64(-0.375 * Float64(w * w)) + Float64(-1.5 / Float64(r * r)))); else tmp = Float64(Float64(r * w) * Float64(w * Float64(Float64(Float64(0.375 + Float64(v * -0.25)) * r) / Float64(v + -1.0)))); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 7.8e-5) tmp = (2.0 / (r * r)) + -1.5; elseif (r <= 2.7e+143) tmp = (r * r) * ((-0.375 * (w * w)) + (-1.5 / (r * r))); else tmp = (r * w) * (w * (((0.375 + (v * -0.25)) * r) / (v + -1.0))); end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 7.8e-5], N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision], If[LessEqual[r, 2.7e+143], N[(N[(r * r), $MachinePrecision] * N[(N[(-0.375 * N[(w * w), $MachinePrecision]), $MachinePrecision] + N[(-1.5 / N[(r * r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(r * w), $MachinePrecision] * N[(w * N[(N[(N[(0.375 + N[(v * -0.25), $MachinePrecision]), $MachinePrecision] * r), $MachinePrecision] / N[(v + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 7.8 \cdot 10^{-5}:\\
\;\;\;\;\frac{2}{r \cdot r} + -1.5\\
\mathbf{elif}\;r \leq 2.7 \cdot 10^{+143}:\\
\;\;\;\;\left(r \cdot r\right) \cdot \left(-0.375 \cdot \left(w \cdot w\right) + \frac{-1.5}{r \cdot r}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(r \cdot w\right) \cdot \left(w \cdot \frac{\left(0.375 + v \cdot -0.25\right) \cdot r}{v + -1}\right)\\
\end{array}
\end{array}
if r < 7.7999999999999999e-5Initial program 84.8%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified86.1%
Taylor expanded in w around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6464.8%
Simplified64.8%
if 7.7999999999999999e-5 < r < 2.7000000000000002e143Initial program 89.6%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified99.7%
Taylor expanded in v around 0
associate--l+N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6491.0%
Simplified91.0%
Taylor expanded in r around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6491.0%
Simplified91.0%
if 2.7000000000000002e143 < r Initial program 90.4%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified93.2%
Taylor expanded in r around inf
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6477.6%
Simplified77.6%
associate-/l*N/A
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6482.4%
Applied egg-rr82.4%
Final simplification69.8%
(FPCore (v w r)
:precision binary64
(if (<= r 7.4e-5)
(+ (/ 2.0 (* r r)) -1.5)
(if (<= r 7.4e+152)
(* (* r r) (+ (* -0.375 (* w w)) (/ -1.5 (* r r))))
(* (* r w) (* r (* -0.25 w))))))
double code(double v, double w, double r) {
double tmp;
if (r <= 7.4e-5) {
tmp = (2.0 / (r * r)) + -1.5;
} else if (r <= 7.4e+152) {
tmp = (r * r) * ((-0.375 * (w * w)) + (-1.5 / (r * r)));
} else {
tmp = (r * w) * (r * (-0.25 * w));
}
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 <= 7.4d-5) then
tmp = (2.0d0 / (r * r)) + (-1.5d0)
else if (r <= 7.4d+152) then
tmp = (r * r) * (((-0.375d0) * (w * w)) + ((-1.5d0) / (r * r)))
else
tmp = (r * w) * (r * ((-0.25d0) * w))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 7.4e-5) {
tmp = (2.0 / (r * r)) + -1.5;
} else if (r <= 7.4e+152) {
tmp = (r * r) * ((-0.375 * (w * w)) + (-1.5 / (r * r)));
} else {
tmp = (r * w) * (r * (-0.25 * w));
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 7.4e-5: tmp = (2.0 / (r * r)) + -1.5 elif r <= 7.4e+152: tmp = (r * r) * ((-0.375 * (w * w)) + (-1.5 / (r * r))) else: tmp = (r * w) * (r * (-0.25 * w)) return tmp
function code(v, w, r) tmp = 0.0 if (r <= 7.4e-5) tmp = Float64(Float64(2.0 / Float64(r * r)) + -1.5); elseif (r <= 7.4e+152) tmp = Float64(Float64(r * r) * Float64(Float64(-0.375 * Float64(w * w)) + Float64(-1.5 / Float64(r * r)))); else tmp = Float64(Float64(r * w) * Float64(r * Float64(-0.25 * w))); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 7.4e-5) tmp = (2.0 / (r * r)) + -1.5; elseif (r <= 7.4e+152) tmp = (r * r) * ((-0.375 * (w * w)) + (-1.5 / (r * r))); else tmp = (r * w) * (r * (-0.25 * w)); end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 7.4e-5], N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision], If[LessEqual[r, 7.4e+152], N[(N[(r * r), $MachinePrecision] * N[(N[(-0.375 * N[(w * w), $MachinePrecision]), $MachinePrecision] + N[(-1.5 / N[(r * r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(r * w), $MachinePrecision] * N[(r * N[(-0.25 * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 7.4 \cdot 10^{-5}:\\
\;\;\;\;\frac{2}{r \cdot r} + -1.5\\
\mathbf{elif}\;r \leq 7.4 \cdot 10^{+152}:\\
\;\;\;\;\left(r \cdot r\right) \cdot \left(-0.375 \cdot \left(w \cdot w\right) + \frac{-1.5}{r \cdot r}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(r \cdot w\right) \cdot \left(r \cdot \left(-0.25 \cdot w\right)\right)\\
\end{array}
\end{array}
if r < 7.39999999999999962e-5Initial program 84.8%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified86.1%
Taylor expanded in w around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6464.8%
Simplified64.8%
if 7.39999999999999962e-5 < r < 7.39999999999999992e152Initial program 89.6%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified99.7%
Taylor expanded in v around 0
associate--l+N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6491.0%
Simplified91.0%
Taylor expanded in r around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6491.0%
Simplified91.0%
if 7.39999999999999992e152 < r Initial program 90.4%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified93.2%
Taylor expanded in r around inf
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6477.6%
Simplified77.6%
Taylor expanded in v around inf
*-commutativeN/A
associate-*l*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6477.8%
Simplified77.8%
*-commutativeN/A
associate-*l*N/A
associate-*l*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6482.2%
Applied egg-rr82.2%
Final simplification69.7%
(FPCore (v w r) :precision binary64 (+ (* (+ 0.375 (* v -0.25)) (* (* r w) (* r (/ w (+ v -1.0))))) (+ (/ 2.0 (* r r)) -1.5)))
double code(double v, double w, double r) {
return ((0.375 + (v * -0.25)) * ((r * w) * (r * (w / (v + -1.0))))) + ((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 = ((0.375d0 + (v * (-0.25d0))) * ((r * w) * (r * (w / (v + (-1.0d0)))))) + ((2.0d0 / (r * r)) + (-1.5d0))
end function
public static double code(double v, double w, double r) {
return ((0.375 + (v * -0.25)) * ((r * w) * (r * (w / (v + -1.0))))) + ((2.0 / (r * r)) + -1.5);
}
def code(v, w, r): return ((0.375 + (v * -0.25)) * ((r * w) * (r * (w / (v + -1.0))))) + ((2.0 / (r * r)) + -1.5)
function code(v, w, r) return Float64(Float64(Float64(0.375 + Float64(v * -0.25)) * Float64(Float64(r * w) * Float64(r * Float64(w / Float64(v + -1.0))))) + Float64(Float64(2.0 / Float64(r * r)) + -1.5)) end
function tmp = code(v, w, r) tmp = ((0.375 + (v * -0.25)) * ((r * w) * (r * (w / (v + -1.0))))) + ((2.0 / (r * r)) + -1.5); end
code[v_, w_, r_] := N[(N[(N[(0.375 + N[(v * -0.25), $MachinePrecision]), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] * N[(r * N[(w / N[(v + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(0.375 + v \cdot -0.25\right) \cdot \left(\left(r \cdot w\right) \cdot \left(r \cdot \frac{w}{v + -1}\right)\right) + \left(\frac{2}{r \cdot r} + -1.5\right)
\end{array}
Initial program 86.0%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified88.5%
associate-/l*N/A
*-commutativeN/A
associate-*r*N/A
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6499.1%
Applied egg-rr99.1%
Final simplification99.1%
(FPCore (v w r) :precision binary64 (if (<= r 1.25e+222) (+ (+ (/ 2.0 (* r r)) -1.5) (* w (* -0.375 (* r (* r w))))) (* (* r w) (* (* r w) (/ (+ 0.375 (* v -0.25)) (+ v -1.0))))))
double code(double v, double w, double r) {
double tmp;
if (r <= 1.25e+222) {
tmp = ((2.0 / (r * r)) + -1.5) + (w * (-0.375 * (r * (r * w))));
} else {
tmp = (r * w) * ((r * w) * ((0.375 + (v * -0.25)) / (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) :: tmp
if (r <= 1.25d+222) then
tmp = ((2.0d0 / (r * r)) + (-1.5d0)) + (w * ((-0.375d0) * (r * (r * w))))
else
tmp = (r * w) * ((r * w) * ((0.375d0 + (v * (-0.25d0))) / (v + (-1.0d0))))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 1.25e+222) {
tmp = ((2.0 / (r * r)) + -1.5) + (w * (-0.375 * (r * (r * w))));
} else {
tmp = (r * w) * ((r * w) * ((0.375 + (v * -0.25)) / (v + -1.0)));
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 1.25e+222: tmp = ((2.0 / (r * r)) + -1.5) + (w * (-0.375 * (r * (r * w)))) else: tmp = (r * w) * ((r * w) * ((0.375 + (v * -0.25)) / (v + -1.0))) return tmp
function code(v, w, r) tmp = 0.0 if (r <= 1.25e+222) tmp = Float64(Float64(Float64(2.0 / Float64(r * r)) + -1.5) + Float64(w * Float64(-0.375 * Float64(r * Float64(r * w))))); else tmp = Float64(Float64(r * w) * Float64(Float64(r * w) * Float64(Float64(0.375 + Float64(v * -0.25)) / Float64(v + -1.0)))); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 1.25e+222) tmp = ((2.0 / (r * r)) + -1.5) + (w * (-0.375 * (r * (r * w)))); else tmp = (r * w) * ((r * w) * ((0.375 + (v * -0.25)) / (v + -1.0))); end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 1.25e+222], N[(N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision] + N[(w * N[(-0.375 * N[(r * N[(r * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(r * w), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] * N[(N[(0.375 + N[(v * -0.25), $MachinePrecision]), $MachinePrecision] / N[(v + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 1.25 \cdot 10^{+222}:\\
\;\;\;\;\left(\frac{2}{r \cdot r} + -1.5\right) + w \cdot \left(-0.375 \cdot \left(r \cdot \left(r \cdot w\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(r \cdot w\right) \cdot \left(\left(r \cdot w\right) \cdot \frac{0.375 + v \cdot -0.25}{v + -1}\right)\\
\end{array}
\end{array}
if r < 1.25000000000000006e222Initial program 85.8%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified88.1%
Taylor expanded in v around 0
associate--l+N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6480.1%
Simplified80.1%
associate-*l*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6493.3%
Applied egg-rr93.3%
if 1.25000000000000006e222 < r Initial program 87.9%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified92.7%
Taylor expanded in r around inf
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6477.7%
Simplified77.7%
associate-/l*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6489.9%
Applied egg-rr89.9%
Final simplification93.0%
(FPCore (v w r) :precision binary64 (+ (+ (/ 2.0 (* r r)) -1.5) (* (* r w) (* w (* r -0.375)))))
double code(double v, double w, double r) {
return ((2.0 / (r * r)) + -1.5) + ((r * w) * (w * (r * -0.375)));
}
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)) + ((r * w) * (w * (r * (-0.375d0))))
end function
public static double code(double v, double w, double r) {
return ((2.0 / (r * r)) + -1.5) + ((r * w) * (w * (r * -0.375)));
}
def code(v, w, r): return ((2.0 / (r * r)) + -1.5) + ((r * w) * (w * (r * -0.375)))
function code(v, w, r) return Float64(Float64(Float64(2.0 / Float64(r * r)) + -1.5) + Float64(Float64(r * w) * Float64(w * Float64(r * -0.375)))) end
function tmp = code(v, w, r) tmp = ((2.0 / (r * r)) + -1.5) + ((r * w) * (w * (r * -0.375))); end
code[v_, w_, r_] := N[(N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision] + N[(N[(r * w), $MachinePrecision] * N[(w * N[(r * -0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\frac{2}{r \cdot r} + -1.5\right) + \left(r \cdot w\right) \cdot \left(w \cdot \left(r \cdot -0.375\right)\right)
\end{array}
Initial program 86.0%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified88.5%
Taylor expanded in v around 0
associate--l+N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6479.3%
Simplified79.3%
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6493.8%
Applied egg-rr93.8%
Final simplification93.8%
(FPCore (v w r) :precision binary64 (if (<= r 7e+81) (+ (/ 2.0 (* r r)) -1.5) (* (* r w) (* r (* -0.25 w)))))
double code(double v, double w, double r) {
double tmp;
if (r <= 7e+81) {
tmp = (2.0 / (r * r)) + -1.5;
} else {
tmp = (r * w) * (r * (-0.25 * w));
}
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 <= 7d+81) then
tmp = (2.0d0 / (r * r)) + (-1.5d0)
else
tmp = (r * w) * (r * ((-0.25d0) * w))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 7e+81) {
tmp = (2.0 / (r * r)) + -1.5;
} else {
tmp = (r * w) * (r * (-0.25 * w));
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 7e+81: tmp = (2.0 / (r * r)) + -1.5 else: tmp = (r * w) * (r * (-0.25 * w)) return tmp
function code(v, w, r) tmp = 0.0 if (r <= 7e+81) tmp = Float64(Float64(2.0 / Float64(r * r)) + -1.5); else tmp = Float64(Float64(r * w) * Float64(r * Float64(-0.25 * w))); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 7e+81) tmp = (2.0 / (r * r)) + -1.5; else tmp = (r * w) * (r * (-0.25 * w)); end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 7e+81], N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision], N[(N[(r * w), $MachinePrecision] * N[(r * N[(-0.25 * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 7 \cdot 10^{+81}:\\
\;\;\;\;\frac{2}{r \cdot r} + -1.5\\
\mathbf{else}:\\
\;\;\;\;\left(r \cdot w\right) \cdot \left(r \cdot \left(-0.25 \cdot w\right)\right)\\
\end{array}
\end{array}
if r < 7.0000000000000001e81Initial program 85.2%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified87.0%
Taylor expanded in w around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6463.8%
Simplified63.8%
if 7.0000000000000001e81 < r Initial program 89.2%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified95.2%
Taylor expanded in r around inf
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6476.2%
Simplified76.2%
Taylor expanded in v around inf
*-commutativeN/A
associate-*l*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6478.5%
Simplified78.5%
*-commutativeN/A
associate-*l*N/A
associate-*l*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6481.6%
Applied egg-rr81.6%
Final simplification67.0%
(FPCore (v w r) :precision binary64 (if (<= r 5.4e+81) (+ (/ 2.0 (* r r)) -1.5) (* (* r w) (* r (* w -0.375)))))
double code(double v, double w, double r) {
double tmp;
if (r <= 5.4e+81) {
tmp = (2.0 / (r * r)) + -1.5;
} else {
tmp = (r * w) * (r * (w * -0.375));
}
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 <= 5.4d+81) then
tmp = (2.0d0 / (r * r)) + (-1.5d0)
else
tmp = (r * w) * (r * (w * (-0.375d0)))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 5.4e+81) {
tmp = (2.0 / (r * r)) + -1.5;
} else {
tmp = (r * w) * (r * (w * -0.375));
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 5.4e+81: tmp = (2.0 / (r * r)) + -1.5 else: tmp = (r * w) * (r * (w * -0.375)) return tmp
function code(v, w, r) tmp = 0.0 if (r <= 5.4e+81) tmp = Float64(Float64(2.0 / Float64(r * r)) + -1.5); else tmp = Float64(Float64(r * w) * Float64(r * Float64(w * -0.375))); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 5.4e+81) tmp = (2.0 / (r * r)) + -1.5; else tmp = (r * w) * (r * (w * -0.375)); end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 5.4e+81], N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision], N[(N[(r * w), $MachinePrecision] * N[(r * N[(w * -0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 5.4 \cdot 10^{+81}:\\
\;\;\;\;\frac{2}{r \cdot r} + -1.5\\
\mathbf{else}:\\
\;\;\;\;\left(r \cdot w\right) \cdot \left(r \cdot \left(w \cdot -0.375\right)\right)\\
\end{array}
\end{array}
if r < 5.3999999999999999e81Initial program 85.2%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified87.0%
Taylor expanded in w around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6463.8%
Simplified63.8%
if 5.3999999999999999e81 < r Initial program 89.2%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified95.2%
Taylor expanded in r around inf
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6476.2%
Simplified76.2%
associate-/l*N/A
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6483.5%
Applied egg-rr83.5%
Taylor expanded in v around 0
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6473.5%
Simplified73.5%
Final simplification65.5%
(FPCore (v w r) :precision binary64 (if (<= r 8.2e+81) (+ (/ 2.0 (* r r)) -1.5) (* r (* r (* -0.25 (* w w))))))
double code(double v, double w, double r) {
double tmp;
if (r <= 8.2e+81) {
tmp = (2.0 / (r * r)) + -1.5;
} else {
tmp = r * (r * (-0.25 * (w * w)));
}
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 <= 8.2d+81) then
tmp = (2.0d0 / (r * r)) + (-1.5d0)
else
tmp = r * (r * ((-0.25d0) * (w * w)))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 8.2e+81) {
tmp = (2.0 / (r * r)) + -1.5;
} else {
tmp = r * (r * (-0.25 * (w * w)));
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 8.2e+81: tmp = (2.0 / (r * r)) + -1.5 else: tmp = r * (r * (-0.25 * (w * w))) return tmp
function code(v, w, r) tmp = 0.0 if (r <= 8.2e+81) tmp = Float64(Float64(2.0 / Float64(r * r)) + -1.5); else tmp = Float64(r * Float64(r * Float64(-0.25 * Float64(w * w)))); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 8.2e+81) tmp = (2.0 / (r * r)) + -1.5; else tmp = r * (r * (-0.25 * (w * w))); end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 8.2e+81], N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision], N[(r * N[(r * N[(-0.25 * N[(w * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 8.2 \cdot 10^{+81}:\\
\;\;\;\;\frac{2}{r \cdot r} + -1.5\\
\mathbf{else}:\\
\;\;\;\;r \cdot \left(r \cdot \left(-0.25 \cdot \left(w \cdot w\right)\right)\right)\\
\end{array}
\end{array}
if r < 8.20000000000000024e81Initial program 85.2%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified87.0%
Taylor expanded in w around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6463.8%
Simplified63.8%
if 8.20000000000000024e81 < r Initial program 89.2%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified95.2%
Taylor expanded in r around inf
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6476.2%
Simplified76.2%
Taylor expanded in v around inf
*-commutativeN/A
associate-*l*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6478.5%
Simplified78.5%
Final simplification66.4%
(FPCore (v w r) :precision binary64 (if (<= r 5.5e+81) (+ (/ 2.0 (* r r)) -1.5) (* r (* r (* -0.375 (* w w))))))
double code(double v, double w, double r) {
double tmp;
if (r <= 5.5e+81) {
tmp = (2.0 / (r * r)) + -1.5;
} else {
tmp = r * (r * (-0.375 * (w * w)));
}
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 <= 5.5d+81) then
tmp = (2.0d0 / (r * r)) + (-1.5d0)
else
tmp = r * (r * ((-0.375d0) * (w * w)))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 5.5e+81) {
tmp = (2.0 / (r * r)) + -1.5;
} else {
tmp = r * (r * (-0.375 * (w * w)));
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 5.5e+81: tmp = (2.0 / (r * r)) + -1.5 else: tmp = r * (r * (-0.375 * (w * w))) return tmp
function code(v, w, r) tmp = 0.0 if (r <= 5.5e+81) tmp = Float64(Float64(2.0 / Float64(r * r)) + -1.5); else tmp = Float64(r * Float64(r * Float64(-0.375 * Float64(w * w)))); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 5.5e+81) tmp = (2.0 / (r * r)) + -1.5; else tmp = r * (r * (-0.375 * (w * w))); end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 5.5e+81], N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision], N[(r * N[(r * N[(-0.375 * N[(w * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 5.5 \cdot 10^{+81}:\\
\;\;\;\;\frac{2}{r \cdot r} + -1.5\\
\mathbf{else}:\\
\;\;\;\;r \cdot \left(r \cdot \left(-0.375 \cdot \left(w \cdot w\right)\right)\right)\\
\end{array}
\end{array}
if r < 5.5000000000000003e81Initial program 85.2%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified87.0%
Taylor expanded in w around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6463.8%
Simplified63.8%
if 5.5000000000000003e81 < r Initial program 89.2%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified95.2%
Taylor expanded in v around 0
associate--l+N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6473.3%
Simplified73.3%
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6487.7%
Applied egg-rr87.7%
Taylor expanded in r around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6473.0%
Simplified73.0%
(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 86.0%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified88.5%
Taylor expanded in w around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6455.2%
Simplified55.2%
(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 86.0%
sub-negN/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
Simplified88.5%
Taylor expanded in r around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6442.7%
Simplified42.7%
herbie shell --seed 2024164
(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))