
(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 18 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (v w r) :precision binary64 (- (- (+ 3.0 (/ 2.0 (* r r))) (/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v))) 4.5))
double code(double v, double w, double r) {
return ((3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) - 4.5;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
code = ((3.0d0 + (2.0d0 / (r * r))) - (((0.125d0 * (3.0d0 - (2.0d0 * v))) * (((w * w) * r) * r)) / (1.0d0 - v))) - 4.5d0
end function
public static double code(double v, double w, double r) {
return ((3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) - 4.5;
}
def code(v, w, r): return ((3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) - 4.5
function code(v, w, r) return Float64(Float64(Float64(3.0 + Float64(2.0 / Float64(r * r))) - Float64(Float64(Float64(0.125 * Float64(3.0 - Float64(2.0 * v))) * Float64(Float64(Float64(w * w) * r) * r)) / Float64(1.0 - v))) - 4.5) end
function tmp = code(v, w, r) tmp = ((3.0 + (2.0 / (r * r))) - (((0.125 * (3.0 - (2.0 * v))) * (((w * w) * r) * r)) / (1.0 - v))) - 4.5; end
code[v_, w_, r_] := N[(N[(N[(3.0 + N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(0.125 * N[(3.0 - N[(2.0 * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(w * w), $MachinePrecision] * r), $MachinePrecision] * r), $MachinePrecision]), $MachinePrecision] / N[(1.0 - v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(3 + \frac{2}{r \cdot r}\right) - \frac{\left(0.125 \cdot \left(3 - 2 \cdot v\right)\right) \cdot \left(\left(\left(w \cdot w\right) \cdot r\right) \cdot r\right)}{1 - v}\right) - 4.5
\end{array}
(FPCore (v w r) :precision binary64 (- (- (+ 3.0 (/ 2.0 (* r r))) (* (* r w) (* w (/ (+ 0.375 (* v -0.25)) (/ (- 1.0 v) r))))) 4.5))
double code(double v, double w, double r) {
return ((3.0 + (2.0 / (r * r))) - ((r * w) * (w * ((0.375 + (v * -0.25)) / ((1.0 - v) / r))))) - 4.5;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
code = ((3.0d0 + (2.0d0 / (r * r))) - ((r * w) * (w * ((0.375d0 + (v * (-0.25d0))) / ((1.0d0 - v) / r))))) - 4.5d0
end function
public static double code(double v, double w, double r) {
return ((3.0 + (2.0 / (r * r))) - ((r * w) * (w * ((0.375 + (v * -0.25)) / ((1.0 - v) / r))))) - 4.5;
}
def code(v, w, r): return ((3.0 + (2.0 / (r * r))) - ((r * w) * (w * ((0.375 + (v * -0.25)) / ((1.0 - v) / r))))) - 4.5
function code(v, w, r) return Float64(Float64(Float64(3.0 + Float64(2.0 / Float64(r * r))) - Float64(Float64(r * w) * Float64(w * Float64(Float64(0.375 + Float64(v * -0.25)) / Float64(Float64(1.0 - v) / r))))) - 4.5) end
function tmp = code(v, w, r) tmp = ((3.0 + (2.0 / (r * r))) - ((r * w) * (w * ((0.375 + (v * -0.25)) / ((1.0 - v) / r))))) - 4.5; end
code[v_, w_, r_] := N[(N[(N[(3.0 + N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(r * w), $MachinePrecision] * N[(w * N[(N[(0.375 + N[(v * -0.25), $MachinePrecision]), $MachinePrecision] / N[(N[(1.0 - v), $MachinePrecision] / r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(3 + \frac{2}{r \cdot r}\right) - \left(r \cdot w\right) \cdot \left(w \cdot \frac{0.375 + v \cdot -0.25}{\frac{1 - v}{r}}\right)\right) - 4.5
\end{array}
Initial program 91.4%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr92.1%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<= r 2.6e-143)
t_0
(if (<= r 1e+86)
(+ t_0 (+ -1.5 (* (* r r) (* (* w w) -0.375))))
(- (- 3.0 (* (* r w) (* w (* r 0.375)))) 4.5)))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (r <= 2.6e-143) {
tmp = t_0;
} else if (r <= 1e+86) {
tmp = t_0 + (-1.5 + ((r * r) * ((w * w) * -0.375)));
} else {
tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5;
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: t_0
real(8) :: tmp
t_0 = 2.0d0 / (r * r)
if (r <= 2.6d-143) then
tmp = t_0
else if (r <= 1d+86) then
tmp = t_0 + ((-1.5d0) + ((r * r) * ((w * w) * (-0.375d0))))
else
tmp = (3.0d0 - ((r * w) * (w * (r * 0.375d0)))) - 4.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 (r <= 2.6e-143) {
tmp = t_0;
} else if (r <= 1e+86) {
tmp = t_0 + (-1.5 + ((r * r) * ((w * w) * -0.375)));
} else {
tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5;
}
return tmp;
}
def code(v, w, r): t_0 = 2.0 / (r * r) tmp = 0 if r <= 2.6e-143: tmp = t_0 elif r <= 1e+86: tmp = t_0 + (-1.5 + ((r * r) * ((w * w) * -0.375))) else: tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5 return tmp
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (r <= 2.6e-143) tmp = t_0; elseif (r <= 1e+86) tmp = Float64(t_0 + Float64(-1.5 + Float64(Float64(r * r) * Float64(Float64(w * w) * -0.375)))); else tmp = Float64(Float64(3.0 - Float64(Float64(r * w) * Float64(w * Float64(r * 0.375)))) - 4.5); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 2.0 / (r * r); tmp = 0.0; if (r <= 2.6e-143) tmp = t_0; elseif (r <= 1e+86) tmp = t_0 + (-1.5 + ((r * r) * ((w * w) * -0.375))); else tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5; end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[r, 2.6e-143], t$95$0, If[LessEqual[r, 1e+86], N[(t$95$0 + N[(-1.5 + N[(N[(r * r), $MachinePrecision] * N[(N[(w * w), $MachinePrecision] * -0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(3.0 - N[(N[(r * w), $MachinePrecision] * N[(w * N[(r * 0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;r \leq 2.6 \cdot 10^{-143}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;r \leq 10^{+86}:\\
\;\;\;\;t\_0 + \left(-1.5 + \left(r \cdot r\right) \cdot \left(\left(w \cdot w\right) \cdot -0.375\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(3 - \left(r \cdot w\right) \cdot \left(w \cdot \left(r \cdot 0.375\right)\right)\right) - 4.5\\
\end{array}
\end{array}
if r < 2.59999999999999987e-143Initial program 90.0%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified98.8%
Taylor expanded in r around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6454.5%
Simplified54.5%
if 2.59999999999999987e-143 < r < 1e86Initial program 97.6%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified99.7%
Taylor expanded in v around 0
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6495.6%
Simplified95.6%
if 1e86 < r Initial program 90.7%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr95.2%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr100.0%
Taylor expanded in v around 0
*-commutativeN/A
*-lowering-*.f6491.1%
Simplified91.1%
Taylor expanded in r around inf
Simplified91.1%
Final simplification67.1%
(FPCore (v w r)
:precision binary64
(if (<= r 1.6e-5)
(- (- (/ 2.0 (* r r)) (* (* r w) (* w (* r 0.375)))) 4.5)
(-
(- 3.0 (* (* r w) (* w (/ (+ 0.375 (* v -0.25)) (/ (- 1.0 v) r)))))
4.5)))
double code(double v, double w, double r) {
double tmp;
if (r <= 1.6e-5) {
tmp = ((2.0 / (r * r)) - ((r * w) * (w * (r * 0.375)))) - 4.5;
} else {
tmp = (3.0 - ((r * w) * (w * ((0.375 + (v * -0.25)) / ((1.0 - v) / r))))) - 4.5;
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: tmp
if (r <= 1.6d-5) then
tmp = ((2.0d0 / (r * r)) - ((r * w) * (w * (r * 0.375d0)))) - 4.5d0
else
tmp = (3.0d0 - ((r * w) * (w * ((0.375d0 + (v * (-0.25d0))) / ((1.0d0 - v) / r))))) - 4.5d0
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 1.6e-5) {
tmp = ((2.0 / (r * r)) - ((r * w) * (w * (r * 0.375)))) - 4.5;
} else {
tmp = (3.0 - ((r * w) * (w * ((0.375 + (v * -0.25)) / ((1.0 - v) / r))))) - 4.5;
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 1.6e-5: tmp = ((2.0 / (r * r)) - ((r * w) * (w * (r * 0.375)))) - 4.5 else: tmp = (3.0 - ((r * w) * (w * ((0.375 + (v * -0.25)) / ((1.0 - v) / r))))) - 4.5 return tmp
function code(v, w, r) tmp = 0.0 if (r <= 1.6e-5) tmp = Float64(Float64(Float64(2.0 / Float64(r * r)) - Float64(Float64(r * w) * Float64(w * Float64(r * 0.375)))) - 4.5); else tmp = Float64(Float64(3.0 - Float64(Float64(r * w) * Float64(w * Float64(Float64(0.375 + Float64(v * -0.25)) / Float64(Float64(1.0 - v) / r))))) - 4.5); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 1.6e-5) tmp = ((2.0 / (r * r)) - ((r * w) * (w * (r * 0.375)))) - 4.5; else tmp = (3.0 - ((r * w) * (w * ((0.375 + (v * -0.25)) / ((1.0 - v) / r))))) - 4.5; end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 1.6e-5], N[(N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] - N[(N[(r * w), $MachinePrecision] * N[(w * N[(r * 0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision], N[(N[(3.0 - N[(N[(r * w), $MachinePrecision] * N[(w * N[(N[(0.375 + N[(v * -0.25), $MachinePrecision]), $MachinePrecision] / N[(N[(1.0 - v), $MachinePrecision] / r), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 1.6 \cdot 10^{-5}:\\
\;\;\;\;\left(\frac{2}{r \cdot r} - \left(r \cdot w\right) \cdot \left(w \cdot \left(r \cdot 0.375\right)\right)\right) - 4.5\\
\mathbf{else}:\\
\;\;\;\;\left(3 - \left(r \cdot w\right) \cdot \left(w \cdot \frac{0.375 + v \cdot -0.25}{\frac{1 - v}{r}}\right)\right) - 4.5\\
\end{array}
\end{array}
if r < 1.59999999999999993e-5Initial program 91.3%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr91.2%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in v around 0
*-commutativeN/A
*-lowering-*.f6497.4%
Simplified97.4%
Taylor expanded in r around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6489.6%
Simplified89.6%
if 1.59999999999999993e-5 < r Initial program 91.9%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr95.0%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in r around inf
Simplified99.9%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<= r 2.5e-143)
t_0
(if (<= r 2.4e-6)
(+ t_0 (* (* r r) (* -0.25 (* w w))))
(- (- 3.0 (* (* r w) (* w (* r 0.375)))) 4.5)))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (r <= 2.5e-143) {
tmp = t_0;
} else if (r <= 2.4e-6) {
tmp = t_0 + ((r * r) * (-0.25 * (w * w)));
} else {
tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5;
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: t_0
real(8) :: tmp
t_0 = 2.0d0 / (r * r)
if (r <= 2.5d-143) then
tmp = t_0
else if (r <= 2.4d-6) then
tmp = t_0 + ((r * r) * ((-0.25d0) * (w * w)))
else
tmp = (3.0d0 - ((r * w) * (w * (r * 0.375d0)))) - 4.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 (r <= 2.5e-143) {
tmp = t_0;
} else if (r <= 2.4e-6) {
tmp = t_0 + ((r * r) * (-0.25 * (w * w)));
} else {
tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5;
}
return tmp;
}
def code(v, w, r): t_0 = 2.0 / (r * r) tmp = 0 if r <= 2.5e-143: tmp = t_0 elif r <= 2.4e-6: tmp = t_0 + ((r * r) * (-0.25 * (w * w))) else: tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5 return tmp
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (r <= 2.5e-143) tmp = t_0; elseif (r <= 2.4e-6) tmp = Float64(t_0 + Float64(Float64(r * r) * Float64(-0.25 * Float64(w * w)))); else tmp = Float64(Float64(3.0 - Float64(Float64(r * w) * Float64(w * Float64(r * 0.375)))) - 4.5); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 2.0 / (r * r); tmp = 0.0; if (r <= 2.5e-143) tmp = t_0; elseif (r <= 2.4e-6) tmp = t_0 + ((r * r) * (-0.25 * (w * w))); else tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5; end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[r, 2.5e-143], t$95$0, If[LessEqual[r, 2.4e-6], N[(t$95$0 + N[(N[(r * r), $MachinePrecision] * N[(-0.25 * N[(w * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(3.0 - N[(N[(r * w), $MachinePrecision] * N[(w * N[(r * 0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;r \leq 2.5 \cdot 10^{-143}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;r \leq 2.4 \cdot 10^{-6}:\\
\;\;\;\;t\_0 + \left(r \cdot r\right) \cdot \left(-0.25 \cdot \left(w \cdot w\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(3 - \left(r \cdot w\right) \cdot \left(w \cdot \left(r \cdot 0.375\right)\right)\right) - 4.5\\
\end{array}
\end{array}
if r < 2.5000000000000001e-143Initial program 90.0%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified98.8%
Taylor expanded in r around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6454.5%
Simplified54.5%
if 2.5000000000000001e-143 < r < 2.3999999999999999e-6Initial program 99.8%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified99.8%
Taylor expanded in v around inf
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6496.6%
Simplified96.6%
Taylor expanded in r around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6496.6%
Simplified96.6%
if 2.3999999999999999e-6 < r Initial program 91.9%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr95.0%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in v around 0
*-commutativeN/A
*-lowering-*.f6490.6%
Simplified90.6%
Taylor expanded in r around inf
Simplified90.6%
Final simplification66.8%
(FPCore (v w r) :precision binary64 (+ (/ 2.0 (* r r)) (+ (* (+ 0.375 (* v -0.25)) (/ (* (* r w) (* r w)) (+ v -1.0))) -1.5)))
double code(double v, double w, double r) {
return (2.0 / (r * r)) + (((0.375 + (v * -0.25)) * (((r * w) * (r * w)) / (v + -1.0))) + -1.5);
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
code = (2.0d0 / (r * r)) + (((0.375d0 + (v * (-0.25d0))) * (((r * w) * (r * w)) / (v + (-1.0d0)))) + (-1.5d0))
end function
public static double code(double v, double w, double r) {
return (2.0 / (r * r)) + (((0.375 + (v * -0.25)) * (((r * w) * (r * w)) / (v + -1.0))) + -1.5);
}
def code(v, w, r): return (2.0 / (r * r)) + (((0.375 + (v * -0.25)) * (((r * w) * (r * w)) / (v + -1.0))) + -1.5)
function code(v, w, r) return Float64(Float64(2.0 / Float64(r * r)) + Float64(Float64(Float64(0.375 + Float64(v * -0.25)) * Float64(Float64(Float64(r * w) * Float64(r * w)) / Float64(v + -1.0))) + -1.5)) end
function tmp = code(v, w, r) tmp = (2.0 / (r * r)) + (((0.375 + (v * -0.25)) * (((r * w) * (r * w)) / (v + -1.0))) + -1.5); end
code[v_, w_, r_] := N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(0.375 + N[(v * -0.25), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(r * w), $MachinePrecision] * N[(r * w), $MachinePrecision]), $MachinePrecision] / N[(v + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.5), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{r \cdot r} + \left(\left(0.375 + v \cdot -0.25\right) \cdot \frac{\left(r \cdot w\right) \cdot \left(r \cdot w\right)}{v + -1} + -1.5\right)
\end{array}
Initial program 91.4%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified99.2%
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.8%
Applied egg-rr99.8%
(FPCore (v w r) :precision binary64 (+ (/ 2.0 (* r r)) (+ -1.5 (* (+ 0.375 (* v -0.25)) (/ (* r (* w (* r w))) (+ v -1.0))))))
double code(double v, double w, double r) {
return (2.0 / (r * r)) + (-1.5 + ((0.375 + (v * -0.25)) * ((r * (w * (r * w))) / (v + -1.0))));
}
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) + ((0.375d0 + (v * (-0.25d0))) * ((r * (w * (r * w))) / (v + (-1.0d0)))))
end function
public static double code(double v, double w, double r) {
return (2.0 / (r * r)) + (-1.5 + ((0.375 + (v * -0.25)) * ((r * (w * (r * w))) / (v + -1.0))));
}
def code(v, w, r): return (2.0 / (r * r)) + (-1.5 + ((0.375 + (v * -0.25)) * ((r * (w * (r * w))) / (v + -1.0))))
function code(v, w, r) return Float64(Float64(2.0 / Float64(r * r)) + Float64(-1.5 + Float64(Float64(0.375 + Float64(v * -0.25)) * Float64(Float64(r * Float64(w * Float64(r * w))) / Float64(v + -1.0))))) end
function tmp = code(v, w, r) tmp = (2.0 / (r * r)) + (-1.5 + ((0.375 + (v * -0.25)) * ((r * (w * (r * w))) / (v + -1.0)))); end
code[v_, w_, r_] := N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + N[(-1.5 + N[(N[(0.375 + N[(v * -0.25), $MachinePrecision]), $MachinePrecision] * N[(N[(r * N[(w * N[(r * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(v + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{r \cdot r} + \left(-1.5 + \left(0.375 + v \cdot -0.25\right) \cdot \frac{r \cdot \left(w \cdot \left(r \cdot w\right)\right)}{v + -1}\right)
\end{array}
Initial program 91.4%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified99.2%
Final simplification99.2%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<= r 8.5e-34)
t_0
(if (<= r 4e+90)
(* (* r r) (* (* w w) -0.375))
(if (<= r 1.35e+106) (+ t_0 -1.5) (* (* r w) (* (* r w) -0.25)))))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (r <= 8.5e-34) {
tmp = t_0;
} else if (r <= 4e+90) {
tmp = (r * r) * ((w * w) * -0.375);
} else if (r <= 1.35e+106) {
tmp = t_0 + -1.5;
} else {
tmp = (r * w) * ((r * w) * -0.25);
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: t_0
real(8) :: tmp
t_0 = 2.0d0 / (r * r)
if (r <= 8.5d-34) then
tmp = t_0
else if (r <= 4d+90) then
tmp = (r * r) * ((w * w) * (-0.375d0))
else if (r <= 1.35d+106) then
tmp = t_0 + (-1.5d0)
else
tmp = (r * w) * ((r * w) * (-0.25d0))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (r <= 8.5e-34) {
tmp = t_0;
} else if (r <= 4e+90) {
tmp = (r * r) * ((w * w) * -0.375);
} else if (r <= 1.35e+106) {
tmp = t_0 + -1.5;
} else {
tmp = (r * w) * ((r * w) * -0.25);
}
return tmp;
}
def code(v, w, r): t_0 = 2.0 / (r * r) tmp = 0 if r <= 8.5e-34: tmp = t_0 elif r <= 4e+90: tmp = (r * r) * ((w * w) * -0.375) elif r <= 1.35e+106: tmp = t_0 + -1.5 else: tmp = (r * w) * ((r * w) * -0.25) return tmp
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (r <= 8.5e-34) tmp = t_0; elseif (r <= 4e+90) tmp = Float64(Float64(r * r) * Float64(Float64(w * w) * -0.375)); elseif (r <= 1.35e+106) tmp = Float64(t_0 + -1.5); else tmp = Float64(Float64(r * w) * Float64(Float64(r * w) * -0.25)); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 2.0 / (r * r); tmp = 0.0; if (r <= 8.5e-34) tmp = t_0; elseif (r <= 4e+90) tmp = (r * r) * ((w * w) * -0.375); elseif (r <= 1.35e+106) tmp = t_0 + -1.5; else tmp = (r * w) * ((r * w) * -0.25); end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[r, 8.5e-34], t$95$0, If[LessEqual[r, 4e+90], N[(N[(r * r), $MachinePrecision] * N[(N[(w * w), $MachinePrecision] * -0.375), $MachinePrecision]), $MachinePrecision], If[LessEqual[r, 1.35e+106], N[(t$95$0 + -1.5), $MachinePrecision], N[(N[(r * w), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] * -0.25), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;r \leq 8.5 \cdot 10^{-34}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;r \leq 4 \cdot 10^{+90}:\\
\;\;\;\;\left(r \cdot r\right) \cdot \left(\left(w \cdot w\right) \cdot -0.375\right)\\
\mathbf{elif}\;r \leq 1.35 \cdot 10^{+106}:\\
\;\;\;\;t\_0 + -1.5\\
\mathbf{else}:\\
\;\;\;\;\left(r \cdot w\right) \cdot \left(\left(r \cdot w\right) \cdot -0.25\right)\\
\end{array}
\end{array}
if r < 8.5000000000000001e-34Initial program 91.0%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified98.9%
Taylor expanded in r around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6457.9%
Simplified57.9%
if 8.5000000000000001e-34 < r < 3.99999999999999987e90Initial program 96.0%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr96.1%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in v around 0
*-commutativeN/A
*-lowering-*.f6492.6%
Simplified92.6%
Taylor expanded in r around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.3%
Simplified69.3%
if 3.99999999999999987e90 < r < 1.35000000000000003e106Initial program 100.0%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified100.0%
Taylor expanded in r around 0
Simplified72.3%
if 1.35000000000000003e106 < r Initial program 88.3%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified100.0%
Taylor expanded in v around inf
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6476.2%
Simplified76.2%
Taylor expanded in r around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6466.1%
Simplified66.1%
associate-*r*N/A
swap-sqrN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6474.2%
Applied egg-rr74.2%
Final simplification61.5%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<= r 9.2e-34)
t_0
(if (<= r 4e+90)
(* (* r r) (* (* w w) -0.375))
(if (<= r 1.35e+106) (+ t_0 -1.5) (* r (* w (* w (* r -0.25)))))))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (r <= 9.2e-34) {
tmp = t_0;
} else if (r <= 4e+90) {
tmp = (r * r) * ((w * w) * -0.375);
} else if (r <= 1.35e+106) {
tmp = t_0 + -1.5;
} else {
tmp = r * (w * (w * (r * -0.25)));
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: t_0
real(8) :: tmp
t_0 = 2.0d0 / (r * r)
if (r <= 9.2d-34) then
tmp = t_0
else if (r <= 4d+90) then
tmp = (r * r) * ((w * w) * (-0.375d0))
else if (r <= 1.35d+106) then
tmp = t_0 + (-1.5d0)
else
tmp = r * (w * (w * (r * (-0.25d0))))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (r <= 9.2e-34) {
tmp = t_0;
} else if (r <= 4e+90) {
tmp = (r * r) * ((w * w) * -0.375);
} else if (r <= 1.35e+106) {
tmp = t_0 + -1.5;
} else {
tmp = r * (w * (w * (r * -0.25)));
}
return tmp;
}
def code(v, w, r): t_0 = 2.0 / (r * r) tmp = 0 if r <= 9.2e-34: tmp = t_0 elif r <= 4e+90: tmp = (r * r) * ((w * w) * -0.375) elif r <= 1.35e+106: tmp = t_0 + -1.5 else: tmp = r * (w * (w * (r * -0.25))) return tmp
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (r <= 9.2e-34) tmp = t_0; elseif (r <= 4e+90) tmp = Float64(Float64(r * r) * Float64(Float64(w * w) * -0.375)); elseif (r <= 1.35e+106) tmp = Float64(t_0 + -1.5); else tmp = Float64(r * Float64(w * Float64(w * Float64(r * -0.25)))); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 2.0 / (r * r); tmp = 0.0; if (r <= 9.2e-34) tmp = t_0; elseif (r <= 4e+90) tmp = (r * r) * ((w * w) * -0.375); elseif (r <= 1.35e+106) tmp = t_0 + -1.5; else tmp = r * (w * (w * (r * -0.25))); end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[r, 9.2e-34], t$95$0, If[LessEqual[r, 4e+90], N[(N[(r * r), $MachinePrecision] * N[(N[(w * w), $MachinePrecision] * -0.375), $MachinePrecision]), $MachinePrecision], If[LessEqual[r, 1.35e+106], N[(t$95$0 + -1.5), $MachinePrecision], N[(r * N[(w * N[(w * N[(r * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;r \leq 9.2 \cdot 10^{-34}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;r \leq 4 \cdot 10^{+90}:\\
\;\;\;\;\left(r \cdot r\right) \cdot \left(\left(w \cdot w\right) \cdot -0.375\right)\\
\mathbf{elif}\;r \leq 1.35 \cdot 10^{+106}:\\
\;\;\;\;t\_0 + -1.5\\
\mathbf{else}:\\
\;\;\;\;r \cdot \left(w \cdot \left(w \cdot \left(r \cdot -0.25\right)\right)\right)\\
\end{array}
\end{array}
if r < 9.20000000000000045e-34Initial program 91.0%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified98.9%
Taylor expanded in r around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6457.9%
Simplified57.9%
if 9.20000000000000045e-34 < r < 3.99999999999999987e90Initial program 96.0%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr96.1%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in v around 0
*-commutativeN/A
*-lowering-*.f6492.6%
Simplified92.6%
Taylor expanded in r around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.3%
Simplified69.3%
if 3.99999999999999987e90 < r < 1.35000000000000003e106Initial program 100.0%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified100.0%
Taylor expanded in r around 0
Simplified72.3%
if 1.35000000000000003e106 < r Initial program 88.3%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified100.0%
Taylor expanded in v around inf
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6476.2%
Simplified76.2%
Taylor expanded in r around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6466.1%
Simplified66.1%
*-commutativeN/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6466.6%
Applied egg-rr66.6%
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6474.1%
Applied egg-rr74.1%
Final simplification61.4%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<= r 9e-34)
t_0
(if (<= r 3.8e+90)
(* (* r r) (* (* w w) -0.375))
(if (<= r 1.35e+106) (+ t_0 -1.5) (* -0.25 (* r (* r (* w w)))))))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (r <= 9e-34) {
tmp = t_0;
} else if (r <= 3.8e+90) {
tmp = (r * r) * ((w * w) * -0.375);
} else if (r <= 1.35e+106) {
tmp = t_0 + -1.5;
} else {
tmp = -0.25 * (r * (r * (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) :: t_0
real(8) :: tmp
t_0 = 2.0d0 / (r * r)
if (r <= 9d-34) then
tmp = t_0
else if (r <= 3.8d+90) then
tmp = (r * r) * ((w * w) * (-0.375d0))
else if (r <= 1.35d+106) then
tmp = t_0 + (-1.5d0)
else
tmp = (-0.25d0) * (r * (r * (w * w)))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (r <= 9e-34) {
tmp = t_0;
} else if (r <= 3.8e+90) {
tmp = (r * r) * ((w * w) * -0.375);
} else if (r <= 1.35e+106) {
tmp = t_0 + -1.5;
} else {
tmp = -0.25 * (r * (r * (w * w)));
}
return tmp;
}
def code(v, w, r): t_0 = 2.0 / (r * r) tmp = 0 if r <= 9e-34: tmp = t_0 elif r <= 3.8e+90: tmp = (r * r) * ((w * w) * -0.375) elif r <= 1.35e+106: tmp = t_0 + -1.5 else: tmp = -0.25 * (r * (r * (w * w))) return tmp
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (r <= 9e-34) tmp = t_0; elseif (r <= 3.8e+90) tmp = Float64(Float64(r * r) * Float64(Float64(w * w) * -0.375)); elseif (r <= 1.35e+106) tmp = Float64(t_0 + -1.5); else tmp = Float64(-0.25 * Float64(r * Float64(r * Float64(w * w)))); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 2.0 / (r * r); tmp = 0.0; if (r <= 9e-34) tmp = t_0; elseif (r <= 3.8e+90) tmp = (r * r) * ((w * w) * -0.375); elseif (r <= 1.35e+106) tmp = t_0 + -1.5; else tmp = -0.25 * (r * (r * (w * w))); end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[r, 9e-34], t$95$0, If[LessEqual[r, 3.8e+90], N[(N[(r * r), $MachinePrecision] * N[(N[(w * w), $MachinePrecision] * -0.375), $MachinePrecision]), $MachinePrecision], If[LessEqual[r, 1.35e+106], N[(t$95$0 + -1.5), $MachinePrecision], N[(-0.25 * N[(r * N[(r * N[(w * w), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;r \leq 9 \cdot 10^{-34}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;r \leq 3.8 \cdot 10^{+90}:\\
\;\;\;\;\left(r \cdot r\right) \cdot \left(\left(w \cdot w\right) \cdot -0.375\right)\\
\mathbf{elif}\;r \leq 1.35 \cdot 10^{+106}:\\
\;\;\;\;t\_0 + -1.5\\
\mathbf{else}:\\
\;\;\;\;-0.25 \cdot \left(r \cdot \left(r \cdot \left(w \cdot w\right)\right)\right)\\
\end{array}
\end{array}
if r < 9.00000000000000085e-34Initial program 91.0%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified98.9%
Taylor expanded in r around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6457.9%
Simplified57.9%
if 9.00000000000000085e-34 < r < 3.8000000000000001e90Initial program 96.0%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr96.1%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in v around 0
*-commutativeN/A
*-lowering-*.f6492.6%
Simplified92.6%
Taylor expanded in r around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.3%
Simplified69.3%
if 3.8000000000000001e90 < r < 1.35000000000000003e106Initial program 100.0%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified100.0%
Taylor expanded in r around 0
Simplified72.3%
if 1.35000000000000003e106 < r Initial program 88.3%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified100.0%
Taylor expanded in v around inf
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6476.2%
Simplified76.2%
Taylor expanded in r around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6466.1%
Simplified66.1%
associate-*r*N/A
swap-sqrN/A
*-lowering-*.f64N/A
swap-sqrN/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6473.1%
Applied egg-rr73.1%
Final simplification61.3%
(FPCore (v w r)
:precision binary64
(let* ((t_0 (/ 2.0 (* r r))))
(if (<= v 7.885e-29)
(- (- (+ 3.0 t_0) (* (* r w) (* w (* r 0.375)))) 4.5)
(+ t_0 (+ -1.5 (* (* r w) (* (* r w) -0.25)))))))
double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (v <= 7.885e-29) {
tmp = ((3.0 + t_0) - ((r * w) * (w * (r * 0.375)))) - 4.5;
} else {
tmp = t_0 + (-1.5 + ((r * w) * ((r * w) * -0.25)));
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: t_0
real(8) :: tmp
t_0 = 2.0d0 / (r * r)
if (v <= 7.885d-29) then
tmp = ((3.0d0 + t_0) - ((r * w) * (w * (r * 0.375d0)))) - 4.5d0
else
tmp = t_0 + ((-1.5d0) + ((r * w) * ((r * w) * (-0.25d0))))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double t_0 = 2.0 / (r * r);
double tmp;
if (v <= 7.885e-29) {
tmp = ((3.0 + t_0) - ((r * w) * (w * (r * 0.375)))) - 4.5;
} else {
tmp = t_0 + (-1.5 + ((r * w) * ((r * w) * -0.25)));
}
return tmp;
}
def code(v, w, r): t_0 = 2.0 / (r * r) tmp = 0 if v <= 7.885e-29: tmp = ((3.0 + t_0) - ((r * w) * (w * (r * 0.375)))) - 4.5 else: tmp = t_0 + (-1.5 + ((r * w) * ((r * w) * -0.25))) return tmp
function code(v, w, r) t_0 = Float64(2.0 / Float64(r * r)) tmp = 0.0 if (v <= 7.885e-29) tmp = Float64(Float64(Float64(3.0 + t_0) - Float64(Float64(r * w) * Float64(w * Float64(r * 0.375)))) - 4.5); else tmp = Float64(t_0 + Float64(-1.5 + Float64(Float64(r * w) * Float64(Float64(r * w) * -0.25)))); end return tmp end
function tmp_2 = code(v, w, r) t_0 = 2.0 / (r * r); tmp = 0.0; if (v <= 7.885e-29) tmp = ((3.0 + t_0) - ((r * w) * (w * (r * 0.375)))) - 4.5; else tmp = t_0 + (-1.5 + ((r * w) * ((r * w) * -0.25))); end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[v, 7.885e-29], N[(N[(N[(3.0 + t$95$0), $MachinePrecision] - N[(N[(r * w), $MachinePrecision] * N[(w * N[(r * 0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision], N[(t$95$0 + N[(-1.5 + N[(N[(r * w), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{r \cdot r}\\
\mathbf{if}\;v \leq 7.885 \cdot 10^{-29}:\\
\;\;\;\;\left(\left(3 + t\_0\right) - \left(r \cdot w\right) \cdot \left(w \cdot \left(r \cdot 0.375\right)\right)\right) - 4.5\\
\mathbf{else}:\\
\;\;\;\;t\_0 + \left(-1.5 + \left(r \cdot w\right) \cdot \left(\left(r \cdot w\right) \cdot -0.25\right)\right)\\
\end{array}
\end{array}
if v < 7.88500000000000052e-29Initial program 93.6%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr93.5%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in v around 0
*-commutativeN/A
*-lowering-*.f6498.0%
Simplified98.0%
if 7.88500000000000052e-29 < v Initial program 82.3%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified99.9%
Taylor expanded in v around inf
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6481.1%
Simplified81.1%
associate-*r*N/A
swap-sqrN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6498.9%
Applied egg-rr98.9%
Final simplification98.2%
(FPCore (v w r) :precision binary64 (let* ((t_0 (* (* r w) (* w (* r 0.375))))) (if (<= r 0.0054) (- (- (/ 2.0 (* r r)) t_0) 4.5) (- (- 3.0 t_0) 4.5))))
double code(double v, double w, double r) {
double t_0 = (r * w) * (w * (r * 0.375));
double tmp;
if (r <= 0.0054) {
tmp = ((2.0 / (r * r)) - t_0) - 4.5;
} else {
tmp = (3.0 - t_0) - 4.5;
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: t_0
real(8) :: tmp
t_0 = (r * w) * (w * (r * 0.375d0))
if (r <= 0.0054d0) then
tmp = ((2.0d0 / (r * r)) - t_0) - 4.5d0
else
tmp = (3.0d0 - t_0) - 4.5d0
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double t_0 = (r * w) * (w * (r * 0.375));
double tmp;
if (r <= 0.0054) {
tmp = ((2.0 / (r * r)) - t_0) - 4.5;
} else {
tmp = (3.0 - t_0) - 4.5;
}
return tmp;
}
def code(v, w, r): t_0 = (r * w) * (w * (r * 0.375)) tmp = 0 if r <= 0.0054: tmp = ((2.0 / (r * r)) - t_0) - 4.5 else: tmp = (3.0 - t_0) - 4.5 return tmp
function code(v, w, r) t_0 = Float64(Float64(r * w) * Float64(w * Float64(r * 0.375))) tmp = 0.0 if (r <= 0.0054) tmp = Float64(Float64(Float64(2.0 / Float64(r * r)) - t_0) - 4.5); else tmp = Float64(Float64(3.0 - t_0) - 4.5); end return tmp end
function tmp_2 = code(v, w, r) t_0 = (r * w) * (w * (r * 0.375)); tmp = 0.0; if (r <= 0.0054) tmp = ((2.0 / (r * r)) - t_0) - 4.5; else tmp = (3.0 - t_0) - 4.5; end tmp_2 = tmp; end
code[v_, w_, r_] := Block[{t$95$0 = N[(N[(r * w), $MachinePrecision] * N[(w * N[(r * 0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[r, 0.0054], N[(N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] - t$95$0), $MachinePrecision] - 4.5), $MachinePrecision], N[(N[(3.0 - t$95$0), $MachinePrecision] - 4.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(r \cdot w\right) \cdot \left(w \cdot \left(r \cdot 0.375\right)\right)\\
\mathbf{if}\;r \leq 0.0054:\\
\;\;\;\;\left(\frac{2}{r \cdot r} - t\_0\right) - 4.5\\
\mathbf{else}:\\
\;\;\;\;\left(3 - t\_0\right) - 4.5\\
\end{array}
\end{array}
if r < 0.0054000000000000003Initial program 91.3%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr91.3%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in v around 0
*-commutativeN/A
*-lowering-*.f6497.3%
Simplified97.3%
Taylor expanded in r around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6489.5%
Simplified89.5%
if 0.0054000000000000003 < r Initial program 91.7%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr94.9%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in v around 0
*-commutativeN/A
*-lowering-*.f6490.9%
Simplified90.9%
Taylor expanded in r around inf
Simplified90.9%
(FPCore (v w r) :precision binary64 (if (<= r 1.6e-5) (+ (/ 2.0 (* r r)) (+ -1.5 (* (* r w) (* (* r w) -0.25)))) (- (- 3.0 (* (* r w) (* w (* r 0.375)))) 4.5)))
double code(double v, double w, double r) {
double tmp;
if (r <= 1.6e-5) {
tmp = (2.0 / (r * r)) + (-1.5 + ((r * w) * ((r * w) * -0.25)));
} else {
tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5;
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: tmp
if (r <= 1.6d-5) then
tmp = (2.0d0 / (r * r)) + ((-1.5d0) + ((r * w) * ((r * w) * (-0.25d0))))
else
tmp = (3.0d0 - ((r * w) * (w * (r * 0.375d0)))) - 4.5d0
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 1.6e-5) {
tmp = (2.0 / (r * r)) + (-1.5 + ((r * w) * ((r * w) * -0.25)));
} else {
tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5;
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 1.6e-5: tmp = (2.0 / (r * r)) + (-1.5 + ((r * w) * ((r * w) * -0.25))) else: tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5 return tmp
function code(v, w, r) tmp = 0.0 if (r <= 1.6e-5) tmp = Float64(Float64(2.0 / Float64(r * r)) + Float64(-1.5 + Float64(Float64(r * w) * Float64(Float64(r * w) * -0.25)))); else tmp = Float64(Float64(3.0 - Float64(Float64(r * w) * Float64(w * Float64(r * 0.375)))) - 4.5); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 1.6e-5) tmp = (2.0 / (r * r)) + (-1.5 + ((r * w) * ((r * w) * -0.25))); else tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5; end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 1.6e-5], N[(N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision] + N[(-1.5 + N[(N[(r * w), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(3.0 - N[(N[(r * w), $MachinePrecision] * N[(w * N[(r * 0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 1.6 \cdot 10^{-5}:\\
\;\;\;\;\frac{2}{r \cdot r} + \left(-1.5 + \left(r \cdot w\right) \cdot \left(\left(r \cdot w\right) \cdot -0.25\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(3 - \left(r \cdot w\right) \cdot \left(w \cdot \left(r \cdot 0.375\right)\right)\right) - 4.5\\
\end{array}
\end{array}
if r < 1.59999999999999993e-5Initial program 91.3%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified99.0%
Taylor expanded in v around inf
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6483.7%
Simplified83.7%
associate-*r*N/A
swap-sqrN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6493.3%
Applied egg-rr93.3%
if 1.59999999999999993e-5 < r Initial program 91.9%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr95.0%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in v around 0
*-commutativeN/A
*-lowering-*.f6490.6%
Simplified90.6%
Taylor expanded in r around inf
Simplified90.6%
Final simplification92.7%
(FPCore (v w r)
:precision binary64
(if (<= r 7.3e-34)
(/ 2.0 (* r r))
(if (<= r 1.35e+154)
(+ -1.5 (* (* r r) (* (* w w) -0.375)))
(* (* r w) (* (* r w) -0.25)))))
double code(double v, double w, double r) {
double tmp;
if (r <= 7.3e-34) {
tmp = 2.0 / (r * r);
} else if (r <= 1.35e+154) {
tmp = -1.5 + ((r * r) * ((w * w) * -0.375));
} else {
tmp = (r * w) * ((r * w) * -0.25);
}
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.3d-34) then
tmp = 2.0d0 / (r * r)
else if (r <= 1.35d+154) then
tmp = (-1.5d0) + ((r * r) * ((w * w) * (-0.375d0)))
else
tmp = (r * w) * ((r * w) * (-0.25d0))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 7.3e-34) {
tmp = 2.0 / (r * r);
} else if (r <= 1.35e+154) {
tmp = -1.5 + ((r * r) * ((w * w) * -0.375));
} else {
tmp = (r * w) * ((r * w) * -0.25);
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 7.3e-34: tmp = 2.0 / (r * r) elif r <= 1.35e+154: tmp = -1.5 + ((r * r) * ((w * w) * -0.375)) else: tmp = (r * w) * ((r * w) * -0.25) return tmp
function code(v, w, r) tmp = 0.0 if (r <= 7.3e-34) tmp = Float64(2.0 / Float64(r * r)); elseif (r <= 1.35e+154) tmp = Float64(-1.5 + Float64(Float64(r * r) * Float64(Float64(w * w) * -0.375))); else tmp = Float64(Float64(r * w) * Float64(Float64(r * w) * -0.25)); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 7.3e-34) tmp = 2.0 / (r * r); elseif (r <= 1.35e+154) tmp = -1.5 + ((r * r) * ((w * w) * -0.375)); else tmp = (r * w) * ((r * w) * -0.25); end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 7.3e-34], N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision], If[LessEqual[r, 1.35e+154], N[(-1.5 + N[(N[(r * r), $MachinePrecision] * N[(N[(w * w), $MachinePrecision] * -0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(r * w), $MachinePrecision] * N[(N[(r * w), $MachinePrecision] * -0.25), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 7.3 \cdot 10^{-34}:\\
\;\;\;\;\frac{2}{r \cdot r}\\
\mathbf{elif}\;r \leq 1.35 \cdot 10^{+154}:\\
\;\;\;\;-1.5 + \left(r \cdot r\right) \cdot \left(\left(w \cdot w\right) \cdot -0.375\right)\\
\mathbf{else}:\\
\;\;\;\;\left(r \cdot w\right) \cdot \left(\left(r \cdot w\right) \cdot -0.25\right)\\
\end{array}
\end{array}
if r < 7.29999999999999996e-34Initial program 91.0%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified98.9%
Taylor expanded in r around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6457.9%
Simplified57.9%
if 7.29999999999999996e-34 < r < 1.35000000000000003e154Initial program 95.4%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr97.7%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in v around 0
*-commutativeN/A
*-lowering-*.f6491.1%
Simplified91.1%
Taylor expanded in r around inf
mul-1-negN/A
distribute-rgt-inN/A
distribute-neg-inN/A
associate-*r*N/A
*-commutativeN/A
distribute-lft-neg-inN/A
metadata-evalN/A
unsub-negN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified86.6%
if 1.35000000000000003e154 < r Initial program 86.8%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified100.0%
Taylor expanded in v around inf
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.7%
Simplified67.7%
Taylor expanded in r around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.7%
Simplified67.7%
associate-*r*N/A
swap-sqrN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6479.7%
Applied egg-rr79.7%
Final simplification64.5%
(FPCore (v w r) :precision binary64 (if (<= r 9.2e-34) (/ 2.0 (* r r)) (- (- 3.0 (* (* r w) (* w (* r 0.375)))) 4.5)))
double code(double v, double w, double r) {
double tmp;
if (r <= 9.2e-34) {
tmp = 2.0 / (r * r);
} else {
tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5;
}
return tmp;
}
real(8) function code(v, w, r)
real(8), intent (in) :: v
real(8), intent (in) :: w
real(8), intent (in) :: r
real(8) :: tmp
if (r <= 9.2d-34) then
tmp = 2.0d0 / (r * r)
else
tmp = (3.0d0 - ((r * w) * (w * (r * 0.375d0)))) - 4.5d0
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 9.2e-34) {
tmp = 2.0 / (r * r);
} else {
tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5;
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 9.2e-34: tmp = 2.0 / (r * r) else: tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5 return tmp
function code(v, w, r) tmp = 0.0 if (r <= 9.2e-34) tmp = Float64(2.0 / Float64(r * r)); else tmp = Float64(Float64(3.0 - Float64(Float64(r * w) * Float64(w * Float64(r * 0.375)))) - 4.5); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 9.2e-34) tmp = 2.0 / (r * r); else tmp = (3.0 - ((r * w) * (w * (r * 0.375)))) - 4.5; end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 9.2e-34], N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision], N[(N[(3.0 - N[(N[(r * w), $MachinePrecision] * N[(w * N[(r * 0.375), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 4.5), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 9.2 \cdot 10^{-34}:\\
\;\;\;\;\frac{2}{r \cdot r}\\
\mathbf{else}:\\
\;\;\;\;\left(3 - \left(r \cdot w\right) \cdot \left(w \cdot \left(r \cdot 0.375\right)\right)\right) - 4.5\\
\end{array}
\end{array}
if r < 9.20000000000000045e-34Initial program 91.0%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified98.9%
Taylor expanded in r around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6457.9%
Simplified57.9%
if 9.20000000000000045e-34 < r Initial program 92.6%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr95.5%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in v around 0
*-commutativeN/A
*-lowering-*.f6491.5%
Simplified91.5%
Taylor expanded in r around inf
Simplified88.4%
(FPCore (v w r) :precision binary64 (if (<= r 7.3e-34) (/ 2.0 (* r r)) (* (* r r) (* (* w w) -0.375))))
double code(double v, double w, double r) {
double tmp;
if (r <= 7.3e-34) {
tmp = 2.0 / (r * r);
} else {
tmp = (r * r) * ((w * 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 <= 7.3d-34) then
tmp = 2.0d0 / (r * r)
else
tmp = (r * r) * ((w * w) * (-0.375d0))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 7.3e-34) {
tmp = 2.0 / (r * r);
} else {
tmp = (r * r) * ((w * w) * -0.375);
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 7.3e-34: tmp = 2.0 / (r * r) else: tmp = (r * r) * ((w * w) * -0.375) return tmp
function code(v, w, r) tmp = 0.0 if (r <= 7.3e-34) tmp = Float64(2.0 / Float64(r * r)); else tmp = Float64(Float64(r * r) * Float64(Float64(w * w) * -0.375)); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 7.3e-34) tmp = 2.0 / (r * r); else tmp = (r * r) * ((w * w) * -0.375); end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 7.3e-34], N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision], N[(N[(r * r), $MachinePrecision] * N[(N[(w * w), $MachinePrecision] * -0.375), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 7.3 \cdot 10^{-34}:\\
\;\;\;\;\frac{2}{r \cdot r}\\
\mathbf{else}:\\
\;\;\;\;\left(r \cdot r\right) \cdot \left(\left(w \cdot w\right) \cdot -0.375\right)\\
\end{array}
\end{array}
if r < 7.29999999999999996e-34Initial program 91.0%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified98.9%
Taylor expanded in r around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6457.9%
Simplified57.9%
if 7.29999999999999996e-34 < r Initial program 92.6%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr95.5%
associate-*l*N/A
associate-*r*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
un-div-invN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
associate-*l*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in v around 0
*-commutativeN/A
*-lowering-*.f6491.5%
Simplified91.5%
Taylor expanded in r around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6463.0%
Simplified63.0%
Final simplification59.2%
(FPCore (v w r) :precision binary64 (if (<= r 9.2e-34) (/ 2.0 (* r r)) (* w (* (* r r) (* w -0.25)))))
double code(double v, double w, double r) {
double tmp;
if (r <= 9.2e-34) {
tmp = 2.0 / (r * r);
} else {
tmp = w * ((r * r) * (w * -0.25));
}
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 <= 9.2d-34) then
tmp = 2.0d0 / (r * r)
else
tmp = w * ((r * r) * (w * (-0.25d0)))
end if
code = tmp
end function
public static double code(double v, double w, double r) {
double tmp;
if (r <= 9.2e-34) {
tmp = 2.0 / (r * r);
} else {
tmp = w * ((r * r) * (w * -0.25));
}
return tmp;
}
def code(v, w, r): tmp = 0 if r <= 9.2e-34: tmp = 2.0 / (r * r) else: tmp = w * ((r * r) * (w * -0.25)) return tmp
function code(v, w, r) tmp = 0.0 if (r <= 9.2e-34) tmp = Float64(2.0 / Float64(r * r)); else tmp = Float64(w * Float64(Float64(r * r) * Float64(w * -0.25))); end return tmp end
function tmp_2 = code(v, w, r) tmp = 0.0; if (r <= 9.2e-34) tmp = 2.0 / (r * r); else tmp = w * ((r * r) * (w * -0.25)); end tmp_2 = tmp; end
code[v_, w_, r_] := If[LessEqual[r, 9.2e-34], N[(2.0 / N[(r * r), $MachinePrecision]), $MachinePrecision], N[(w * N[(N[(r * r), $MachinePrecision] * N[(w * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;r \leq 9.2 \cdot 10^{-34}:\\
\;\;\;\;\frac{2}{r \cdot r}\\
\mathbf{else}:\\
\;\;\;\;w \cdot \left(\left(r \cdot r\right) \cdot \left(w \cdot -0.25\right)\right)\\
\end{array}
\end{array}
if r < 9.20000000000000045e-34Initial program 91.0%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified98.9%
Taylor expanded in r around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6457.9%
Simplified57.9%
if 9.20000000000000045e-34 < r Initial program 92.6%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified99.9%
Taylor expanded in v around inf
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6479.3%
Simplified79.3%
Taylor expanded in r around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6457.7%
Simplified57.7%
*-commutativeN/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6457.9%
Applied egg-rr57.9%
Final simplification57.9%
(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 91.4%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified99.2%
Taylor expanded in r around 0
Simplified57.3%
(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 91.4%
associate--l-N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate--l+N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
Simplified99.2%
Taylor expanded in r around 0
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6444.6%
Simplified44.6%
herbie shell --seed 2024152
(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))