
(FPCore (m v) :precision binary64 (* (- (/ (* m (- 1.0 m)) v) 1.0) (- 1.0 m)))
double code(double m, double v) {
return (((m * (1.0 - m)) / v) - 1.0) * (1.0 - m);
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
code = (((m * (1.0d0 - m)) / v) - 1.0d0) * (1.0d0 - m)
end function
public static double code(double m, double v) {
return (((m * (1.0 - m)) / v) - 1.0) * (1.0 - m);
}
def code(m, v): return (((m * (1.0 - m)) / v) - 1.0) * (1.0 - m)
function code(m, v) return Float64(Float64(Float64(Float64(m * Float64(1.0 - m)) / v) - 1.0) * Float64(1.0 - m)) end
function tmp = code(m, v) tmp = (((m * (1.0 - m)) / v) - 1.0) * (1.0 - m); end
code[m_, v_] := N[(N[(N[(N[(m * N[(1.0 - m), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision] - 1.0), $MachinePrecision] * N[(1.0 - m), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \left(1 - m\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 9 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (m v) :precision binary64 (* (- (/ (* m (- 1.0 m)) v) 1.0) (- 1.0 m)))
double code(double m, double v) {
return (((m * (1.0 - m)) / v) - 1.0) * (1.0 - m);
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
code = (((m * (1.0d0 - m)) / v) - 1.0d0) * (1.0d0 - m)
end function
public static double code(double m, double v) {
return (((m * (1.0 - m)) / v) - 1.0) * (1.0 - m);
}
def code(m, v): return (((m * (1.0 - m)) / v) - 1.0) * (1.0 - m)
function code(m, v) return Float64(Float64(Float64(Float64(m * Float64(1.0 - m)) / v) - 1.0) * Float64(1.0 - m)) end
function tmp = code(m, v) tmp = (((m * (1.0 - m)) / v) - 1.0) * (1.0 - m); end
code[m_, v_] := N[(N[(N[(N[(m * N[(1.0 - m), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision] - 1.0), $MachinePrecision] * N[(1.0 - m), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \left(1 - m\right)
\end{array}
(FPCore (m v) :precision binary64 (* (- 1.0 m) (+ (/ m (/ v (- 1.0 m))) -1.0)))
double code(double m, double v) {
return (1.0 - m) * ((m / (v / (1.0 - m))) + -1.0);
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
code = (1.0d0 - m) * ((m / (v / (1.0d0 - m))) + (-1.0d0))
end function
public static double code(double m, double v) {
return (1.0 - m) * ((m / (v / (1.0 - m))) + -1.0);
}
def code(m, v): return (1.0 - m) * ((m / (v / (1.0 - m))) + -1.0)
function code(m, v) return Float64(Float64(1.0 - m) * Float64(Float64(m / Float64(v / Float64(1.0 - m))) + -1.0)) end
function tmp = code(m, v) tmp = (1.0 - m) * ((m / (v / (1.0 - m))) + -1.0); end
code[m_, v_] := N[(N[(1.0 - m), $MachinePrecision] * N[(N[(m / N[(v / N[(1.0 - m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 - m\right) \cdot \left(\frac{m}{\frac{v}{1 - m}} + -1\right)
\end{array}
Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Final simplification99.9%
(FPCore (m v) :precision binary64 (if (<= m 1.7e-15) (+ -1.0 (+ m (/ m v))) (/ (* m (+ 1.0 (* m (+ m -2.0)))) v)))
double code(double m, double v) {
double tmp;
if (m <= 1.7e-15) {
tmp = -1.0 + (m + (m / v));
} else {
tmp = (m * (1.0 + (m * (m + -2.0)))) / v;
}
return tmp;
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
real(8) :: tmp
if (m <= 1.7d-15) then
tmp = (-1.0d0) + (m + (m / v))
else
tmp = (m * (1.0d0 + (m * (m + (-2.0d0))))) / v
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.7e-15) {
tmp = -1.0 + (m + (m / v));
} else {
tmp = (m * (1.0 + (m * (m + -2.0)))) / v;
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.7e-15: tmp = -1.0 + (m + (m / v)) else: tmp = (m * (1.0 + (m * (m + -2.0)))) / v return tmp
function code(m, v) tmp = 0.0 if (m <= 1.7e-15) tmp = Float64(-1.0 + Float64(m + Float64(m / v))); else tmp = Float64(Float64(m * Float64(1.0 + Float64(m * Float64(m + -2.0)))) / v); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.7e-15) tmp = -1.0 + (m + (m / v)); else tmp = (m * (1.0 + (m * (m + -2.0)))) / v; end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.7e-15], N[(-1.0 + N[(m + N[(m / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(m * N[(1.0 + N[(m * N[(m + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.7 \cdot 10^{-15}:\\
\;\;\;\;-1 + \left(m + \frac{m}{v}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{m \cdot \left(1 + m \cdot \left(m + -2\right)\right)}{v}\\
\end{array}
\end{array}
if m < 1.7e-15Initial program 100.0%
*-commutative100.0%
sub-neg100.0%
associate-/l*100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in m around 0 99.8%
+-commutative99.8%
distribute-lft-in99.8%
un-div-inv100.0%
*-rgt-identity100.0%
Applied egg-rr100.0%
if 1.7e-15 < m Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in v around 0 99.9%
Taylor expanded in m around 0 99.9%
+-commutative99.9%
unpow299.9%
distribute-rgt-out99.9%
Simplified99.9%
Final simplification99.9%
(FPCore (m v) :precision binary64 (if (<= m 2.4) (+ -1.0 (+ m (/ m v))) (* m (/ m (/ v (+ m -2.0))))))
double code(double m, double v) {
double tmp;
if (m <= 2.4) {
tmp = -1.0 + (m + (m / v));
} else {
tmp = m * (m / (v / (m + -2.0)));
}
return tmp;
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
real(8) :: tmp
if (m <= 2.4d0) then
tmp = (-1.0d0) + (m + (m / v))
else
tmp = m * (m / (v / (m + (-2.0d0))))
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 2.4) {
tmp = -1.0 + (m + (m / v));
} else {
tmp = m * (m / (v / (m + -2.0)));
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 2.4: tmp = -1.0 + (m + (m / v)) else: tmp = m * (m / (v / (m + -2.0))) return tmp
function code(m, v) tmp = 0.0 if (m <= 2.4) tmp = Float64(-1.0 + Float64(m + Float64(m / v))); else tmp = Float64(m * Float64(m / Float64(v / Float64(m + -2.0)))); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 2.4) tmp = -1.0 + (m + (m / v)); else tmp = m * (m / (v / (m + -2.0))); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 2.4], N[(-1.0 + N[(m + N[(m / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(m * N[(m / N[(v / N[(m + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 2.4:\\
\;\;\;\;-1 + \left(m + \frac{m}{v}\right)\\
\mathbf{else}:\\
\;\;\;\;m \cdot \frac{m}{\frac{v}{m + -2}}\\
\end{array}
\end{array}
if m < 2.39999999999999991Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in m around 0 97.8%
+-commutative97.8%
distribute-lft-in97.8%
un-div-inv97.9%
*-rgt-identity97.9%
Applied egg-rr97.9%
if 2.39999999999999991 < m Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in v around 0 99.9%
Taylor expanded in m around inf 99.4%
+-commutative99.4%
unpow299.4%
distribute-rgt-out99.4%
Simplified99.4%
div-inv99.4%
associate-*l*99.4%
Applied egg-rr99.4%
un-div-inv99.4%
associate-/l*99.5%
Applied egg-rr99.5%
Final simplification98.7%
(FPCore (m v) :precision binary64 (if (<= m 1.65) (* (- 1.0 m) (+ -1.0 (/ m v))) (* m (/ m (/ v (+ m -2.0))))))
double code(double m, double v) {
double tmp;
if (m <= 1.65) {
tmp = (1.0 - m) * (-1.0 + (m / v));
} else {
tmp = m * (m / (v / (m + -2.0)));
}
return tmp;
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
real(8) :: tmp
if (m <= 1.65d0) then
tmp = (1.0d0 - m) * ((-1.0d0) + (m / v))
else
tmp = m * (m / (v / (m + (-2.0d0))))
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.65) {
tmp = (1.0 - m) * (-1.0 + (m / v));
} else {
tmp = m * (m / (v / (m + -2.0)));
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.65: tmp = (1.0 - m) * (-1.0 + (m / v)) else: tmp = m * (m / (v / (m + -2.0))) return tmp
function code(m, v) tmp = 0.0 if (m <= 1.65) tmp = Float64(Float64(1.0 - m) * Float64(-1.0 + Float64(m / v))); else tmp = Float64(m * Float64(m / Float64(v / Float64(m + -2.0)))); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.65) tmp = (1.0 - m) * (-1.0 + (m / v)); else tmp = m * (m / (v / (m + -2.0))); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.65], N[(N[(1.0 - m), $MachinePrecision] * N[(-1.0 + N[(m / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(m * N[(m / N[(v / N[(m + -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.65:\\
\;\;\;\;\left(1 - m\right) \cdot \left(-1 + \frac{m}{v}\right)\\
\mathbf{else}:\\
\;\;\;\;m \cdot \frac{m}{\frac{v}{m + -2}}\\
\end{array}
\end{array}
if m < 1.6499999999999999Initial program 99.9%
Taylor expanded in m around 0 98.0%
if 1.6499999999999999 < m Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in v around 0 99.9%
Taylor expanded in m around inf 99.4%
+-commutative99.4%
unpow299.4%
distribute-rgt-out99.4%
Simplified99.4%
div-inv99.4%
associate-*l*99.4%
Applied egg-rr99.4%
un-div-inv99.4%
associate-/l*99.5%
Applied egg-rr99.5%
Final simplification98.8%
(FPCore (m v) :precision binary64 (if (<= m 1.65) (* (- 1.0 m) (+ -1.0 (/ m v))) (/ (* m (+ m -2.0)) (/ v m))))
double code(double m, double v) {
double tmp;
if (m <= 1.65) {
tmp = (1.0 - m) * (-1.0 + (m / v));
} else {
tmp = (m * (m + -2.0)) / (v / m);
}
return tmp;
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
real(8) :: tmp
if (m <= 1.65d0) then
tmp = (1.0d0 - m) * ((-1.0d0) + (m / v))
else
tmp = (m * (m + (-2.0d0))) / (v / m)
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.65) {
tmp = (1.0 - m) * (-1.0 + (m / v));
} else {
tmp = (m * (m + -2.0)) / (v / m);
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.65: tmp = (1.0 - m) * (-1.0 + (m / v)) else: tmp = (m * (m + -2.0)) / (v / m) return tmp
function code(m, v) tmp = 0.0 if (m <= 1.65) tmp = Float64(Float64(1.0 - m) * Float64(-1.0 + Float64(m / v))); else tmp = Float64(Float64(m * Float64(m + -2.0)) / Float64(v / m)); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.65) tmp = (1.0 - m) * (-1.0 + (m / v)); else tmp = (m * (m + -2.0)) / (v / m); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.65], N[(N[(1.0 - m), $MachinePrecision] * N[(-1.0 + N[(m / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(m * N[(m + -2.0), $MachinePrecision]), $MachinePrecision] / N[(v / m), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.65:\\
\;\;\;\;\left(1 - m\right) \cdot \left(-1 + \frac{m}{v}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{m \cdot \left(m + -2\right)}{\frac{v}{m}}\\
\end{array}
\end{array}
if m < 1.6499999999999999Initial program 99.9%
Taylor expanded in m around 0 98.0%
if 1.6499999999999999 < m Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in v around 0 99.9%
Taylor expanded in m around inf 99.4%
+-commutative99.4%
unpow299.4%
distribute-rgt-out99.4%
Simplified99.4%
associate-/l*99.5%
associate-/r/99.4%
Applied egg-rr99.4%
*-commutative99.4%
clear-num99.5%
un-div-inv99.5%
Applied egg-rr99.5%
Final simplification98.8%
(FPCore (m v) :precision binary64 (+ -1.0 (+ m (/ m v))))
double code(double m, double v) {
return -1.0 + (m + (m / v));
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
code = (-1.0d0) + (m + (m / v))
end function
public static double code(double m, double v) {
return -1.0 + (m + (m / v));
}
def code(m, v): return -1.0 + (m + (m / v))
function code(m, v) return Float64(-1.0 + Float64(m + Float64(m / v))) end
function tmp = code(m, v) tmp = -1.0 + (m + (m / v)); end
code[m_, v_] := N[(-1.0 + N[(m + N[(m / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-1 + \left(m + \frac{m}{v}\right)
\end{array}
Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in m around 0 74.7%
+-commutative74.7%
distribute-lft-in74.7%
un-div-inv74.7%
*-rgt-identity74.7%
Applied egg-rr74.7%
Final simplification74.7%
(FPCore (m v) :precision binary64 (if (<= v 8e-144) (/ m v) (+ m -1.0)))
double code(double m, double v) {
double tmp;
if (v <= 8e-144) {
tmp = m / v;
} else {
tmp = m + -1.0;
}
return tmp;
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
real(8) :: tmp
if (v <= 8d-144) then
tmp = m / v
else
tmp = m + (-1.0d0)
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (v <= 8e-144) {
tmp = m / v;
} else {
tmp = m + -1.0;
}
return tmp;
}
def code(m, v): tmp = 0 if v <= 8e-144: tmp = m / v else: tmp = m + -1.0 return tmp
function code(m, v) tmp = 0.0 if (v <= 8e-144) tmp = Float64(m / v); else tmp = Float64(m + -1.0); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (v <= 8e-144) tmp = m / v; else tmp = m + -1.0; end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[v, 8e-144], N[(m / v), $MachinePrecision], N[(m + -1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq 8 \cdot 10^{-144}:\\
\;\;\;\;\frac{m}{v}\\
\mathbf{else}:\\
\;\;\;\;m + -1\\
\end{array}
\end{array}
if v < 7.9999999999999996e-144Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in v around 0 89.7%
Taylor expanded in m around 0 89.7%
+-commutative89.7%
unpow289.7%
distribute-rgt-out89.7%
Simplified89.7%
Taylor expanded in m around 0 75.3%
if 7.9999999999999996e-144 < v Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in v around inf 44.7%
neg-mul-144.7%
neg-sub044.7%
associate--r-44.7%
metadata-eval44.7%
Simplified44.7%
Final simplification61.9%
(FPCore (m v) :precision binary64 (+ m -1.0))
double code(double m, double v) {
return m + -1.0;
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
code = m + (-1.0d0)
end function
public static double code(double m, double v) {
return m + -1.0;
}
def code(m, v): return m + -1.0
function code(m, v) return Float64(m + -1.0) end
function tmp = code(m, v) tmp = m + -1.0; end
code[m_, v_] := N[(m + -1.0), $MachinePrecision]
\begin{array}{l}
\\
m + -1
\end{array}
Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in v around inf 26.1%
neg-mul-126.1%
neg-sub026.1%
associate--r-26.1%
metadata-eval26.1%
Simplified26.1%
Final simplification26.1%
(FPCore (m v) :precision binary64 -1.0)
double code(double m, double v) {
return -1.0;
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
code = -1.0d0
end function
public static double code(double m, double v) {
return -1.0;
}
def code(m, v): return -1.0
function code(m, v) return -1.0 end
function tmp = code(m, v) tmp = -1.0; end
code[m_, v_] := -1.0
\begin{array}{l}
\\
-1
\end{array}
Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in m around 0 23.7%
Final simplification23.7%
herbie shell --seed 2023309
(FPCore (m v)
:name "b parameter of renormalized beta distribution"
:precision binary64
:pre (and (and (< 0.0 m) (< 0.0 v)) (< v 0.25))
(* (- (/ (* m (- 1.0 m)) v) 1.0) (- 1.0 m)))