
(FPCore (x y) :precision binary64 (+ x (/ (fabs (- y x)) 2.0)))
double code(double x, double y) {
return x + (fabs((y - x)) / 2.0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x + (abs((y - x)) / 2.0d0)
end function
public static double code(double x, double y) {
return x + (Math.abs((y - x)) / 2.0);
}
def code(x, y): return x + (math.fabs((y - x)) / 2.0)
function code(x, y) return Float64(x + Float64(abs(Float64(y - x)) / 2.0)) end
function tmp = code(x, y) tmp = x + (abs((y - x)) / 2.0); end
code[x_, y_] := N[(x + N[(N[Abs[N[(y - x), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{\left|y - x\right|}{2}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 9 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (+ x (/ (fabs (- y x)) 2.0)))
double code(double x, double y) {
return x + (fabs((y - x)) / 2.0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x + (abs((y - x)) / 2.0d0)
end function
public static double code(double x, double y) {
return x + (Math.abs((y - x)) / 2.0);
}
def code(x, y): return x + (math.fabs((y - x)) / 2.0)
function code(x, y) return Float64(x + Float64(abs(Float64(y - x)) / 2.0)) end
function tmp = code(x, y) tmp = x + (abs((y - x)) / 2.0); end
code[x_, y_] := N[(x + N[(N[Abs[N[(y - x), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{\left|y - x\right|}{2}
\end{array}
(FPCore (x y) :precision binary64 (fma (fabs (- y x)) 0.5 x))
double code(double x, double y) {
return fma(fabs((y - x)), 0.5, x);
}
function code(x, y) return fma(abs(Float64(y - x)), 0.5, x) end
code[x_, y_] := N[(N[Abs[N[(y - x), $MachinePrecision]], $MachinePrecision] * 0.5 + x), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\left|y - x\right|, 0.5, x\right)
\end{array}
Initial program 99.9%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
div-invN/A
lower-fma.f64N/A
lift-fabs.f64N/A
neg-fabsN/A
lower-fabs.f64N/A
lift--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
sub-negN/A
lower--.f64N/A
metadata-eval99.9
Applied rewrites99.9%
Final simplification99.9%
(FPCore (x y) :precision binary64 (if (<= (+ (/ (fabs (- y x)) 2.0) x) -5e-216) (* 0.5 x) (fma (- y x) -0.5 x)))
double code(double x, double y) {
double tmp;
if (((fabs((y - x)) / 2.0) + x) <= -5e-216) {
tmp = 0.5 * x;
} else {
tmp = fma((y - x), -0.5, x);
}
return tmp;
}
function code(x, y) tmp = 0.0 if (Float64(Float64(abs(Float64(y - x)) / 2.0) + x) <= -5e-216) tmp = Float64(0.5 * x); else tmp = fma(Float64(y - x), -0.5, x); end return tmp end
code[x_, y_] := If[LessEqual[N[(N[(N[Abs[N[(y - x), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision] + x), $MachinePrecision], -5e-216], N[(0.5 * x), $MachinePrecision], N[(N[(y - x), $MachinePrecision] * -0.5 + x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{\left|y - x\right|}{2} + x \leq -5 \cdot 10^{-216}:\\
\;\;\;\;0.5 \cdot x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(y - x, -0.5, x\right)\\
\end{array}
\end{array}
if (+.f64 x (/.f64 (fabs.f64 (-.f64 y x)) #s(literal 2 binary64))) < -5.00000000000000021e-216Initial program 100.0%
Taylor expanded in x around 0
*-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
lower-fabs.f64N/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
sub-negN/A
lower--.f641.3
Applied rewrites1.3%
Applied rewrites98.8%
Taylor expanded in y around 0
Applied rewrites98.9%
if -5.00000000000000021e-216 < (+.f64 x (/.f64 (fabs.f64 (-.f64 y x)) #s(literal 2 binary64))) Initial program 99.9%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
div-invN/A
lower-fma.f64N/A
lift-fabs.f64N/A
neg-fabsN/A
lower-fabs.f64N/A
lift--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
sub-negN/A
lower--.f64N/A
metadata-eval99.9
Applied rewrites99.9%
lift-fma.f64N/A
Applied rewrites69.7%
Final simplification75.7%
(FPCore (x y) :precision binary64 (if (<= (+ (/ (fabs (- y x)) 2.0) x) -5e-216) (* 0.5 x) (* -0.5 y)))
double code(double x, double y) {
double tmp;
if (((fabs((y - x)) / 2.0) + x) <= -5e-216) {
tmp = 0.5 * x;
} else {
tmp = -0.5 * y;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (((abs((y - x)) / 2.0d0) + x) <= (-5d-216)) then
tmp = 0.5d0 * x
else
tmp = (-0.5d0) * y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (((Math.abs((y - x)) / 2.0) + x) <= -5e-216) {
tmp = 0.5 * x;
} else {
tmp = -0.5 * y;
}
return tmp;
}
def code(x, y): tmp = 0 if ((math.fabs((y - x)) / 2.0) + x) <= -5e-216: tmp = 0.5 * x else: tmp = -0.5 * y return tmp
function code(x, y) tmp = 0.0 if (Float64(Float64(abs(Float64(y - x)) / 2.0) + x) <= -5e-216) tmp = Float64(0.5 * x); else tmp = Float64(-0.5 * y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (((abs((y - x)) / 2.0) + x) <= -5e-216) tmp = 0.5 * x; else tmp = -0.5 * y; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[N[(N[(N[Abs[N[(y - x), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision] + x), $MachinePrecision], -5e-216], N[(0.5 * x), $MachinePrecision], N[(-0.5 * y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{\left|y - x\right|}{2} + x \leq -5 \cdot 10^{-216}:\\
\;\;\;\;0.5 \cdot x\\
\mathbf{else}:\\
\;\;\;\;-0.5 \cdot y\\
\end{array}
\end{array}
if (+.f64 x (/.f64 (fabs.f64 (-.f64 y x)) #s(literal 2 binary64))) < -5.00000000000000021e-216Initial program 100.0%
Taylor expanded in x around 0
*-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
lower-fabs.f64N/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
sub-negN/A
lower--.f641.3
Applied rewrites1.3%
Applied rewrites98.8%
Taylor expanded in y around 0
Applied rewrites98.9%
if -5.00000000000000021e-216 < (+.f64 x (/.f64 (fabs.f64 (-.f64 y x)) #s(literal 2 binary64))) Initial program 99.9%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
div-invN/A
lower-fma.f64N/A
lift-fabs.f64N/A
neg-fabsN/A
lower-fabs.f64N/A
lift--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
sub-negN/A
lower--.f64N/A
metadata-eval99.9
Applied rewrites99.9%
lift-fma.f64N/A
Applied rewrites69.7%
Taylor expanded in y around inf
lower-*.f6435.0
Applied rewrites35.0%
Final simplification48.2%
(FPCore (x y) :precision binary64 (if (<= x -1.5e-165) (* -0.5 (- y x)) (if (<= x 60000.0) (fma (fabs (- y)) 0.5 x) (fma 1.5 x (* -0.5 y)))))
double code(double x, double y) {
double tmp;
if (x <= -1.5e-165) {
tmp = -0.5 * (y - x);
} else if (x <= 60000.0) {
tmp = fma(fabs(-y), 0.5, x);
} else {
tmp = fma(1.5, x, (-0.5 * y));
}
return tmp;
}
function code(x, y) tmp = 0.0 if (x <= -1.5e-165) tmp = Float64(-0.5 * Float64(y - x)); elseif (x <= 60000.0) tmp = fma(abs(Float64(-y)), 0.5, x); else tmp = fma(1.5, x, Float64(-0.5 * y)); end return tmp end
code[x_, y_] := If[LessEqual[x, -1.5e-165], N[(-0.5 * N[(y - x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 60000.0], N[(N[Abs[(-y)], $MachinePrecision] * 0.5 + x), $MachinePrecision], N[(1.5 * x + N[(-0.5 * y), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.5 \cdot 10^{-165}:\\
\;\;\;\;-0.5 \cdot \left(y - x\right)\\
\mathbf{elif}\;x \leq 60000:\\
\;\;\;\;\mathsf{fma}\left(\left|-y\right|, 0.5, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1.5, x, -0.5 \cdot y\right)\\
\end{array}
\end{array}
if x < -1.49999999999999989e-165Initial program 100.0%
Taylor expanded in x around 0
*-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
lower-fabs.f64N/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
sub-negN/A
lower--.f6434.9
Applied rewrites34.9%
Applied rewrites84.8%
if -1.49999999999999989e-165 < x < 6e4Initial program 99.9%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
div-invN/A
lower-fma.f64N/A
lift-fabs.f64N/A
neg-fabsN/A
lower-fabs.f64N/A
lift--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
sub-negN/A
lower--.f64N/A
metadata-eval99.9
Applied rewrites99.9%
Taylor expanded in y around inf
mul-1-negN/A
lower-neg.f6483.8
Applied rewrites83.8%
if 6e4 < x Initial program 99.8%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
div-invN/A
lower-fma.f64N/A
lift-fabs.f64N/A
neg-fabsN/A
lower-fabs.f64N/A
lift--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
sub-negN/A
lower--.f64N/A
metadata-eval99.8
Applied rewrites99.8%
lift-fma.f64N/A
Applied rewrites85.2%
Taylor expanded in y around 0
+-commutativeN/A
associate-+r+N/A
distribute-rgt1-inN/A
metadata-evalN/A
lower-fma.f64N/A
lower-*.f6485.4
Applied rewrites85.4%
Final simplification84.6%
(FPCore (x y) :precision binary64 (if (<= x -1.5e-165) (* -0.5 (- y x)) (if (<= x 60000.0) (fma (fabs (- y)) 0.5 x) (fma (- y x) -0.5 x))))
double code(double x, double y) {
double tmp;
if (x <= -1.5e-165) {
tmp = -0.5 * (y - x);
} else if (x <= 60000.0) {
tmp = fma(fabs(-y), 0.5, x);
} else {
tmp = fma((y - x), -0.5, x);
}
return tmp;
}
function code(x, y) tmp = 0.0 if (x <= -1.5e-165) tmp = Float64(-0.5 * Float64(y - x)); elseif (x <= 60000.0) tmp = fma(abs(Float64(-y)), 0.5, x); else tmp = fma(Float64(y - x), -0.5, x); end return tmp end
code[x_, y_] := If[LessEqual[x, -1.5e-165], N[(-0.5 * N[(y - x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 60000.0], N[(N[Abs[(-y)], $MachinePrecision] * 0.5 + x), $MachinePrecision], N[(N[(y - x), $MachinePrecision] * -0.5 + x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.5 \cdot 10^{-165}:\\
\;\;\;\;-0.5 \cdot \left(y - x\right)\\
\mathbf{elif}\;x \leq 60000:\\
\;\;\;\;\mathsf{fma}\left(\left|-y\right|, 0.5, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(y - x, -0.5, x\right)\\
\end{array}
\end{array}
if x < -1.49999999999999989e-165Initial program 100.0%
Taylor expanded in x around 0
*-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
lower-fabs.f64N/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
sub-negN/A
lower--.f6434.9
Applied rewrites34.9%
Applied rewrites84.8%
if -1.49999999999999989e-165 < x < 6e4Initial program 99.9%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
div-invN/A
lower-fma.f64N/A
lift-fabs.f64N/A
neg-fabsN/A
lower-fabs.f64N/A
lift--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
sub-negN/A
lower--.f64N/A
metadata-eval99.9
Applied rewrites99.9%
Taylor expanded in y around inf
mul-1-negN/A
lower-neg.f6483.8
Applied rewrites83.8%
if 6e4 < x Initial program 99.8%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
div-invN/A
lower-fma.f64N/A
lift-fabs.f64N/A
neg-fabsN/A
lower-fabs.f64N/A
lift--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
sub-negN/A
lower--.f64N/A
metadata-eval99.8
Applied rewrites99.8%
lift-fma.f64N/A
Applied rewrites85.2%
Final simplification84.5%
(FPCore (x y) :precision binary64 (if (<= x -1.5e-165) (* -0.5 (- y x)) (if (<= x 60000.0) (* (fabs (- y x)) 0.5) (fma (- y x) -0.5 x))))
double code(double x, double y) {
double tmp;
if (x <= -1.5e-165) {
tmp = -0.5 * (y - x);
} else if (x <= 60000.0) {
tmp = fabs((y - x)) * 0.5;
} else {
tmp = fma((y - x), -0.5, x);
}
return tmp;
}
function code(x, y) tmp = 0.0 if (x <= -1.5e-165) tmp = Float64(-0.5 * Float64(y - x)); elseif (x <= 60000.0) tmp = Float64(abs(Float64(y - x)) * 0.5); else tmp = fma(Float64(y - x), -0.5, x); end return tmp end
code[x_, y_] := If[LessEqual[x, -1.5e-165], N[(-0.5 * N[(y - x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 60000.0], N[(N[Abs[N[(y - x), $MachinePrecision]], $MachinePrecision] * 0.5), $MachinePrecision], N[(N[(y - x), $MachinePrecision] * -0.5 + x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.5 \cdot 10^{-165}:\\
\;\;\;\;-0.5 \cdot \left(y - x\right)\\
\mathbf{elif}\;x \leq 60000:\\
\;\;\;\;\left|y - x\right| \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(y - x, -0.5, x\right)\\
\end{array}
\end{array}
if x < -1.49999999999999989e-165Initial program 100.0%
Taylor expanded in x around 0
*-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
lower-fabs.f64N/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
sub-negN/A
lower--.f6434.9
Applied rewrites34.9%
Applied rewrites84.8%
if -1.49999999999999989e-165 < x < 6e4Initial program 99.9%
Taylor expanded in x around 0
*-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
lower-fabs.f64N/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
sub-negN/A
lower--.f6483.0
Applied rewrites83.0%
if 6e4 < x Initial program 99.8%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
div-invN/A
lower-fma.f64N/A
lift-fabs.f64N/A
neg-fabsN/A
lower-fabs.f64N/A
lift--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
sub-negN/A
lower--.f64N/A
metadata-eval99.8
Applied rewrites99.8%
lift-fma.f64N/A
Applied rewrites85.2%
Final simplification84.2%
(FPCore (x y) :precision binary64 (if (<= x -6.2e-186) (* 0.5 x) (if (<= x 4.6e-20) (* -0.5 y) (* 1.5 x))))
double code(double x, double y) {
double tmp;
if (x <= -6.2e-186) {
tmp = 0.5 * x;
} else if (x <= 4.6e-20) {
tmp = -0.5 * y;
} else {
tmp = 1.5 * x;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= (-6.2d-186)) then
tmp = 0.5d0 * x
else if (x <= 4.6d-20) then
tmp = (-0.5d0) * y
else
tmp = 1.5d0 * x
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= -6.2e-186) {
tmp = 0.5 * x;
} else if (x <= 4.6e-20) {
tmp = -0.5 * y;
} else {
tmp = 1.5 * x;
}
return tmp;
}
def code(x, y): tmp = 0 if x <= -6.2e-186: tmp = 0.5 * x elif x <= 4.6e-20: tmp = -0.5 * y else: tmp = 1.5 * x return tmp
function code(x, y) tmp = 0.0 if (x <= -6.2e-186) tmp = Float64(0.5 * x); elseif (x <= 4.6e-20) tmp = Float64(-0.5 * y); else tmp = Float64(1.5 * x); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= -6.2e-186) tmp = 0.5 * x; elseif (x <= 4.6e-20) tmp = -0.5 * y; else tmp = 1.5 * x; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, -6.2e-186], N[(0.5 * x), $MachinePrecision], If[LessEqual[x, 4.6e-20], N[(-0.5 * y), $MachinePrecision], N[(1.5 * x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -6.2 \cdot 10^{-186}:\\
\;\;\;\;0.5 \cdot x\\
\mathbf{elif}\;x \leq 4.6 \cdot 10^{-20}:\\
\;\;\;\;-0.5 \cdot y\\
\mathbf{else}:\\
\;\;\;\;1.5 \cdot x\\
\end{array}
\end{array}
if x < -6.20000000000000018e-186Initial program 100.0%
Taylor expanded in x around 0
*-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
lower-fabs.f64N/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
sub-negN/A
lower--.f6439.3
Applied rewrites39.3%
Applied rewrites80.1%
Taylor expanded in y around 0
Applied rewrites61.2%
if -6.20000000000000018e-186 < x < 4.5999999999999998e-20Initial program 99.9%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
div-invN/A
lower-fma.f64N/A
lift-fabs.f64N/A
neg-fabsN/A
lower-fabs.f64N/A
lift--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
sub-negN/A
lower--.f64N/A
metadata-eval99.9
Applied rewrites99.9%
lift-fma.f64N/A
Applied rewrites66.0%
Taylor expanded in y around inf
lower-*.f6449.2
Applied rewrites49.2%
if 4.5999999999999998e-20 < x Initial program 99.8%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
div-invN/A
lower-fma.f64N/A
lift-fabs.f64N/A
neg-fabsN/A
lower-fabs.f64N/A
lift--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
sub-negN/A
lower--.f64N/A
metadata-eval99.8
Applied rewrites99.8%
lift-fma.f64N/A
Applied rewrites80.2%
Taylor expanded in y around 0
distribute-rgt1-inN/A
metadata-evalN/A
lower-*.f6471.4
Applied rewrites71.4%
(FPCore (x y) :precision binary64 (if (<= x 4.6e-20) (* -0.5 (- y x)) (* 1.5 x)))
double code(double x, double y) {
double tmp;
if (x <= 4.6e-20) {
tmp = -0.5 * (y - x);
} else {
tmp = 1.5 * x;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= 4.6d-20) then
tmp = (-0.5d0) * (y - x)
else
tmp = 1.5d0 * x
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 4.6e-20) {
tmp = -0.5 * (y - x);
} else {
tmp = 1.5 * x;
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 4.6e-20: tmp = -0.5 * (y - x) else: tmp = 1.5 * x return tmp
function code(x, y) tmp = 0.0 if (x <= 4.6e-20) tmp = Float64(-0.5 * Float64(y - x)); else tmp = Float64(1.5 * x); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 4.6e-20) tmp = -0.5 * (y - x); else tmp = 1.5 * x; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 4.6e-20], N[(-0.5 * N[(y - x), $MachinePrecision]), $MachinePrecision], N[(1.5 * x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 4.6 \cdot 10^{-20}:\\
\;\;\;\;-0.5 \cdot \left(y - x\right)\\
\mathbf{else}:\\
\;\;\;\;1.5 \cdot x\\
\end{array}
\end{array}
if x < 4.5999999999999998e-20Initial program 100.0%
Taylor expanded in x around 0
*-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
lower-fabs.f64N/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
sub-negN/A
lower--.f6462.0
Applied rewrites62.0%
Applied rewrites66.0%
if 4.5999999999999998e-20 < x Initial program 99.8%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
div-invN/A
lower-fma.f64N/A
lift-fabs.f64N/A
neg-fabsN/A
lower-fabs.f64N/A
lift--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
sub-negN/A
lower--.f64N/A
metadata-eval99.8
Applied rewrites99.8%
lift-fma.f64N/A
Applied rewrites80.2%
Taylor expanded in y around 0
distribute-rgt1-inN/A
metadata-evalN/A
lower-*.f6471.4
Applied rewrites71.4%
Final simplification67.7%
(FPCore (x y) :precision binary64 (* 0.5 x))
double code(double x, double y) {
return 0.5 * x;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 0.5d0 * x
end function
public static double code(double x, double y) {
return 0.5 * x;
}
def code(x, y): return 0.5 * x
function code(x, y) return Float64(0.5 * x) end
function tmp = code(x, y) tmp = 0.5 * x; end
code[x_, y_] := N[(0.5 * x), $MachinePrecision]
\begin{array}{l}
\\
0.5 \cdot x
\end{array}
Initial program 99.9%
Taylor expanded in x around 0
*-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
lower-fabs.f64N/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
sub-negN/A
lower--.f6455.6
Applied rewrites55.6%
Applied rewrites52.2%
Taylor expanded in y around 0
Applied rewrites27.4%
herbie shell --seed 2024276
(FPCore (x y)
:name "Graphics.Rendering.Chart.Plot.AreaSpots:renderSpotLegend from Chart-1.5.3"
:precision binary64
(+ x (/ (fabs (- y x)) 2.0)))