
(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 15 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 (* (+ m -1.0) (/ (+ v (* m (+ m -1.0))) v)))
double code(double m, double v) {
return (m + -1.0) * ((v + (m * (m + -1.0))) / v);
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
code = (m + (-1.0d0)) * ((v + (m * (m + (-1.0d0)))) / v)
end function
public static double code(double m, double v) {
return (m + -1.0) * ((v + (m * (m + -1.0))) / v);
}
def code(m, v): return (m + -1.0) * ((v + (m * (m + -1.0))) / v)
function code(m, v) return Float64(Float64(m + -1.0) * Float64(Float64(v + Float64(m * Float64(m + -1.0))) / v)) end
function tmp = code(m, v) tmp = (m + -1.0) * ((v + (m * (m + -1.0))) / v); end
code[m_, v_] := N[(N[(m + -1.0), $MachinePrecision] * N[(N[(v + N[(m * N[(m + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(m + -1\right) \cdot \frac{v + m \cdot \left(m + -1\right)}{v}
\end{array}
Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in v around 0 99.9%
neg-mul-199.9%
+-commutative99.9%
unsub-neg99.9%
Simplified99.9%
Final simplification99.9%
(FPCore (m v) :precision binary64 (if (<= m 1.25e-18) (+ -1.0 (/ m v)) (* (+ m -1.0) (* m (/ (+ m -1.0) v)))))
double code(double m, double v) {
double tmp;
if (m <= 1.25e-18) {
tmp = -1.0 + (m / v);
} else {
tmp = (m + -1.0) * (m * ((m + -1.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.25d-18) then
tmp = (-1.0d0) + (m / v)
else
tmp = (m + (-1.0d0)) * (m * ((m + (-1.0d0)) / v))
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.25e-18) {
tmp = -1.0 + (m / v);
} else {
tmp = (m + -1.0) * (m * ((m + -1.0) / v));
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.25e-18: tmp = -1.0 + (m / v) else: tmp = (m + -1.0) * (m * ((m + -1.0) / v)) return tmp
function code(m, v) tmp = 0.0 if (m <= 1.25e-18) tmp = Float64(-1.0 + Float64(m / v)); else tmp = Float64(Float64(m + -1.0) * Float64(m * Float64(Float64(m + -1.0) / v))); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.25e-18) tmp = -1.0 + (m / v); else tmp = (m + -1.0) * (m * ((m + -1.0) / v)); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.25e-18], N[(-1.0 + N[(m / v), $MachinePrecision]), $MachinePrecision], N[(N[(m + -1.0), $MachinePrecision] * N[(m * N[(N[(m + -1.0), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.25 \cdot 10^{-18}:\\
\;\;\;\;-1 + \frac{m}{v}\\
\mathbf{else}:\\
\;\;\;\;\left(m + -1\right) \cdot \left(m \cdot \frac{m + -1}{v}\right)\\
\end{array}
\end{array}
if m < 1.25000000000000009e-18Initial program 100.0%
*-commutative100.0%
sub-neg100.0%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in m around 0 99.8%
+-commutative40.5%
distribute-lft-in40.5%
div-inv40.7%
*-rgt-identity40.7%
Applied egg-rr100.0%
Taylor expanded in v around 0 100.0%
if 1.25000000000000009e-18 < 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%
neg-mul-199.9%
+-commutative99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in v around 0 99.9%
associate-*r/99.9%
Simplified99.9%
Final simplification99.9%
(FPCore (m v) :precision binary64 (if (<= m 1.25e-18) (+ -1.0 (/ m v)) (* (- 1.0 m) (* (- 1.0 m) (/ m v)))))
double code(double m, double v) {
double tmp;
if (m <= 1.25e-18) {
tmp = -1.0 + (m / v);
} else {
tmp = (1.0 - m) * ((1.0 - m) * (m / 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.25d-18) then
tmp = (-1.0d0) + (m / v)
else
tmp = (1.0d0 - m) * ((1.0d0 - m) * (m / v))
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.25e-18) {
tmp = -1.0 + (m / v);
} else {
tmp = (1.0 - m) * ((1.0 - m) * (m / v));
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.25e-18: tmp = -1.0 + (m / v) else: tmp = (1.0 - m) * ((1.0 - m) * (m / v)) return tmp
function code(m, v) tmp = 0.0 if (m <= 1.25e-18) tmp = Float64(-1.0 + Float64(m / v)); else tmp = Float64(Float64(1.0 - m) * Float64(Float64(1.0 - m) * Float64(m / v))); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.25e-18) tmp = -1.0 + (m / v); else tmp = (1.0 - m) * ((1.0 - m) * (m / v)); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.25e-18], N[(-1.0 + N[(m / v), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 - m), $MachinePrecision] * N[(N[(1.0 - m), $MachinePrecision] * N[(m / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.25 \cdot 10^{-18}:\\
\;\;\;\;-1 + \frac{m}{v}\\
\mathbf{else}:\\
\;\;\;\;\left(1 - m\right) \cdot \left(\left(1 - m\right) \cdot \frac{m}{v}\right)\\
\end{array}
\end{array}
if m < 1.25000000000000009e-18Initial program 100.0%
*-commutative100.0%
sub-neg100.0%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in m around 0 99.8%
+-commutative40.5%
distribute-lft-in40.5%
div-inv40.7%
*-rgt-identity40.7%
Applied egg-rr100.0%
Taylor expanded in v around 0 100.0%
if 1.25000000000000009e-18 < 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%
neg-mul-199.9%
+-commutative99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in v around 0 99.9%
associate-*r/99.9%
Simplified99.9%
Taylor expanded in m around 0 99.9%
+-commutative99.9%
neg-mul-199.9%
sub-neg99.9%
div-sub99.9%
associate-*r/99.9%
associate-*l/99.9%
Simplified99.9%
Final simplification99.9%
(FPCore (m v) :precision binary64 (if (<= m 1.0) (* (- 1.0 m) (+ -1.0 (/ m v))) (* (/ m v) (+ m -1.0))))
double code(double m, double v) {
double tmp;
if (m <= 1.0) {
tmp = (1.0 - m) * (-1.0 + (m / v));
} else {
tmp = (m / v) * (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 (m <= 1.0d0) then
tmp = (1.0d0 - m) * ((-1.0d0) + (m / v))
else
tmp = (m / v) * (m + (-1.0d0))
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.0) {
tmp = (1.0 - m) * (-1.0 + (m / v));
} else {
tmp = (m / v) * (m + -1.0);
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.0: tmp = (1.0 - m) * (-1.0 + (m / v)) else: tmp = (m / v) * (m + -1.0) return tmp
function code(m, v) tmp = 0.0 if (m <= 1.0) tmp = Float64(Float64(1.0 - m) * Float64(-1.0 + Float64(m / v))); else tmp = Float64(Float64(m / v) * Float64(m + -1.0)); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.0) tmp = (1.0 - m) * (-1.0 + (m / v)); else tmp = (m / v) * (m + -1.0); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.0], N[(N[(1.0 - m), $MachinePrecision] * N[(-1.0 + N[(m / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(m / v), $MachinePrecision] * N[(m + -1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1:\\
\;\;\;\;\left(1 - m\right) \cdot \left(-1 + \frac{m}{v}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{m}{v} \cdot \left(m + -1\right)\\
\end{array}
\end{array}
if m < 1Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in m around 0 97.0%
if 1 < 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%
neg-mul-199.9%
+-commutative99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in v around 0 99.9%
associate-*r/99.9%
Simplified99.9%
Taylor expanded in m around 0 0.1%
associate-*r*0.1%
*-commutative0.1%
div-inv0.1%
frac-2neg0.1%
distribute-rgt-neg-in0.1%
add-sqr-sqrt0.0%
sqrt-unprod84.2%
sqr-neg84.2%
sqrt-unprod84.2%
add-sqr-sqrt84.2%
Applied egg-rr84.2%
associate-/l*84.2%
neg-sub084.2%
associate--r-84.2%
metadata-eval84.2%
+-commutative84.2%
associate-*r/84.2%
associate-*l/84.2%
Simplified84.2%
Final simplification90.4%
(FPCore (m v) :precision binary64 (if (<= m 1.0) (/ (* (- 1.0 m) (- m v)) v) (* (/ m v) (+ m -1.0))))
double code(double m, double v) {
double tmp;
if (m <= 1.0) {
tmp = ((1.0 - m) * (m - v)) / v;
} else {
tmp = (m / v) * (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 (m <= 1.0d0) then
tmp = ((1.0d0 - m) * (m - v)) / v
else
tmp = (m / v) * (m + (-1.0d0))
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.0) {
tmp = ((1.0 - m) * (m - v)) / v;
} else {
tmp = (m / v) * (m + -1.0);
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.0: tmp = ((1.0 - m) * (m - v)) / v else: tmp = (m / v) * (m + -1.0) return tmp
function code(m, v) tmp = 0.0 if (m <= 1.0) tmp = Float64(Float64(Float64(1.0 - m) * Float64(m - v)) / v); else tmp = Float64(Float64(m / v) * Float64(m + -1.0)); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.0) tmp = ((1.0 - m) * (m - v)) / v; else tmp = (m / v) * (m + -1.0); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.0], N[(N[(N[(1.0 - m), $MachinePrecision] * N[(m - v), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision], N[(N[(m / v), $MachinePrecision] * N[(m + -1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1:\\
\;\;\;\;\frac{\left(1 - m\right) \cdot \left(m - v\right)}{v}\\
\mathbf{else}:\\
\;\;\;\;\frac{m}{v} \cdot \left(m + -1\right)\\
\end{array}
\end{array}
if m < 1Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in m around 0 97.0%
Taylor expanded in v around 0 97.0%
+-commutative97.0%
associate-*r*97.0%
neg-mul-197.0%
distribute-rgt-out97.0%
unsub-neg97.0%
Simplified97.0%
if 1 < 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%
neg-mul-199.9%
+-commutative99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in v around 0 99.9%
associate-*r/99.9%
Simplified99.9%
Taylor expanded in m around 0 0.1%
associate-*r*0.1%
*-commutative0.1%
div-inv0.1%
frac-2neg0.1%
distribute-rgt-neg-in0.1%
add-sqr-sqrt0.0%
sqrt-unprod84.2%
sqr-neg84.2%
sqrt-unprod84.2%
add-sqr-sqrt84.2%
Applied egg-rr84.2%
associate-/l*84.2%
neg-sub084.2%
associate--r-84.2%
metadata-eval84.2%
+-commutative84.2%
associate-*r/84.2%
associate-*l/84.2%
Simplified84.2%
Final simplification90.4%
(FPCore (m v) :precision binary64 (if (<= m 1.0) (/ (* (- 1.0 m) (- m v)) v) (* (* m (/ m v)) (+ m -1.0))))
double code(double m, double v) {
double tmp;
if (m <= 1.0) {
tmp = ((1.0 - m) * (m - v)) / v;
} else {
tmp = (m * (m / v)) * (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 (m <= 1.0d0) then
tmp = ((1.0d0 - m) * (m - v)) / v
else
tmp = (m * (m / v)) * (m + (-1.0d0))
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.0) {
tmp = ((1.0 - m) * (m - v)) / v;
} else {
tmp = (m * (m / v)) * (m + -1.0);
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.0: tmp = ((1.0 - m) * (m - v)) / v else: tmp = (m * (m / v)) * (m + -1.0) return tmp
function code(m, v) tmp = 0.0 if (m <= 1.0) tmp = Float64(Float64(Float64(1.0 - m) * Float64(m - v)) / v); else tmp = Float64(Float64(m * Float64(m / v)) * Float64(m + -1.0)); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.0) tmp = ((1.0 - m) * (m - v)) / v; else tmp = (m * (m / v)) * (m + -1.0); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.0], N[(N[(N[(1.0 - m), $MachinePrecision] * N[(m - v), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision], N[(N[(m * N[(m / v), $MachinePrecision]), $MachinePrecision] * N[(m + -1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1:\\
\;\;\;\;\frac{\left(1 - m\right) \cdot \left(m - v\right)}{v}\\
\mathbf{else}:\\
\;\;\;\;\left(m \cdot \frac{m}{v}\right) \cdot \left(m + -1\right)\\
\end{array}
\end{array}
if m < 1Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in m around 0 97.0%
Taylor expanded in v around 0 97.0%
+-commutative97.0%
associate-*r*97.0%
neg-mul-197.0%
distribute-rgt-out97.0%
unsub-neg97.0%
Simplified97.0%
if 1 < 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%
neg-mul-199.9%
+-commutative99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in v around 0 99.9%
associate-*r/99.9%
Simplified99.9%
Taylor expanded in m around inf 98.4%
neg-mul-198.4%
distribute-neg-frac298.4%
Simplified98.4%
Final simplification97.7%
(FPCore (m v) :precision binary64 (if (<= m 1.0) (/ (* (- 1.0 m) (- m v)) v) (/ (* m (* m (+ m -1.0))) v)))
double code(double m, double v) {
double tmp;
if (m <= 1.0) {
tmp = ((1.0 - m) * (m - v)) / v;
} else {
tmp = (m * (m * (m + -1.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.0d0) then
tmp = ((1.0d0 - m) * (m - v)) / v
else
tmp = (m * (m * (m + (-1.0d0)))) / v
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.0) {
tmp = ((1.0 - m) * (m - v)) / v;
} else {
tmp = (m * (m * (m + -1.0))) / v;
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.0: tmp = ((1.0 - m) * (m - v)) / v else: tmp = (m * (m * (m + -1.0))) / v return tmp
function code(m, v) tmp = 0.0 if (m <= 1.0) tmp = Float64(Float64(Float64(1.0 - m) * Float64(m - v)) / v); else tmp = Float64(Float64(m * Float64(m * Float64(m + -1.0))) / v); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.0) tmp = ((1.0 - m) * (m - v)) / v; else tmp = (m * (m * (m + -1.0))) / v; end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.0], N[(N[(N[(1.0 - m), $MachinePrecision] * N[(m - v), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision], N[(N[(m * N[(m * N[(m + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1:\\
\;\;\;\;\frac{\left(1 - m\right) \cdot \left(m - v\right)}{v}\\
\mathbf{else}:\\
\;\;\;\;\frac{m \cdot \left(m \cdot \left(m + -1\right)\right)}{v}\\
\end{array}
\end{array}
if m < 1Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in m around 0 97.0%
Taylor expanded in v around 0 97.0%
+-commutative97.0%
associate-*r*97.0%
neg-mul-197.0%
distribute-rgt-out97.0%
unsub-neg97.0%
Simplified97.0%
if 1 < 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%
neg-mul-199.9%
+-commutative99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in v around 0 99.9%
associate-*r/99.9%
Simplified99.9%
Taylor expanded in m around inf 98.4%
neg-mul-198.4%
distribute-neg-frac298.4%
Simplified98.4%
associate-*r*98.4%
*-commutative98.4%
distribute-frac-neg298.4%
distribute-frac-neg98.4%
associate-*r/98.4%
Applied egg-rr98.4%
Final simplification97.8%
(FPCore (m v) :precision binary64 (* (- 1.0 m) (- -1.0 (* m (/ (+ m -1.0) v)))))
double code(double m, double v) {
return (1.0 - m) * (-1.0 - (m * ((m + -1.0) / v)));
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
code = (1.0d0 - m) * ((-1.0d0) - (m * ((m + (-1.0d0)) / v)))
end function
public static double code(double m, double v) {
return (1.0 - m) * (-1.0 - (m * ((m + -1.0) / v)));
}
def code(m, v): return (1.0 - m) * (-1.0 - (m * ((m + -1.0) / v)))
function code(m, v) return Float64(Float64(1.0 - m) * Float64(-1.0 - Float64(m * Float64(Float64(m + -1.0) / v)))) end
function tmp = code(m, v) tmp = (1.0 - m) * (-1.0 - (m * ((m + -1.0) / v))); end
code[m_, v_] := N[(N[(1.0 - m), $MachinePrecision] * N[(-1.0 - N[(m * N[(N[(m + -1.0), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 - m\right) \cdot \left(-1 - m \cdot \frac{m + -1}{v}\right)
\end{array}
Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
Final simplification99.8%
(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.8%
metadata-eval99.8%
Simplified99.8%
clear-num99.8%
un-div-inv99.9%
Applied egg-rr99.9%
Final simplification99.9%
(FPCore (m v) :precision binary64 (if (<= m 1.0) (+ -1.0 (/ m v)) (* (/ m v) (+ m -1.0))))
double code(double m, double v) {
double tmp;
if (m <= 1.0) {
tmp = -1.0 + (m / v);
} else {
tmp = (m / v) * (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 (m <= 1.0d0) then
tmp = (-1.0d0) + (m / v)
else
tmp = (m / v) * (m + (-1.0d0))
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.0) {
tmp = -1.0 + (m / v);
} else {
tmp = (m / v) * (m + -1.0);
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.0: tmp = -1.0 + (m / v) else: tmp = (m / v) * (m + -1.0) return tmp
function code(m, v) tmp = 0.0 if (m <= 1.0) tmp = Float64(-1.0 + Float64(m / v)); else tmp = Float64(Float64(m / v) * Float64(m + -1.0)); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.0) tmp = -1.0 + (m / v); else tmp = (m / v) * (m + -1.0); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.0], N[(-1.0 + N[(m / v), $MachinePrecision]), $MachinePrecision], N[(N[(m / v), $MachinePrecision] * N[(m + -1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1:\\
\;\;\;\;-1 + \frac{m}{v}\\
\mathbf{else}:\\
\;\;\;\;\frac{m}{v} \cdot \left(m + -1\right)\\
\end{array}
\end{array}
if m < 1Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in m around 0 96.7%
+-commutative42.8%
distribute-lft-in42.8%
div-inv43.0%
*-rgt-identity43.0%
Applied egg-rr96.9%
Taylor expanded in v around 0 96.9%
if 1 < 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%
neg-mul-199.9%
+-commutative99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in v around 0 99.9%
associate-*r/99.9%
Simplified99.9%
Taylor expanded in m around 0 0.1%
associate-*r*0.1%
*-commutative0.1%
div-inv0.1%
frac-2neg0.1%
distribute-rgt-neg-in0.1%
add-sqr-sqrt0.0%
sqrt-unprod84.2%
sqr-neg84.2%
sqrt-unprod84.2%
add-sqr-sqrt84.2%
Applied egg-rr84.2%
associate-/l*84.2%
neg-sub084.2%
associate--r-84.2%
metadata-eval84.2%
+-commutative84.2%
associate-*r/84.2%
associate-*l/84.2%
Simplified84.2%
Final simplification90.3%
(FPCore (m v) :precision binary64 (if (<= m 2.05e-153) -1.0 (+ m (/ m v))))
double code(double m, double v) {
double tmp;
if (m <= 2.05e-153) {
tmp = -1.0;
} else {
tmp = m + (m / v);
}
return tmp;
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
real(8) :: tmp
if (m <= 2.05d-153) then
tmp = -1.0d0
else
tmp = m + (m / v)
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 2.05e-153) {
tmp = -1.0;
} else {
tmp = m + (m / v);
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 2.05e-153: tmp = -1.0 else: tmp = m + (m / v) return tmp
function code(m, v) tmp = 0.0 if (m <= 2.05e-153) tmp = -1.0; else tmp = Float64(m + Float64(m / v)); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 2.05e-153) tmp = -1.0; else tmp = m + (m / v); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 2.05e-153], -1.0, N[(m + N[(m / v), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 2.05 \cdot 10^{-153}:\\
\;\;\;\;-1\\
\mathbf{else}:\\
\;\;\;\;m + \frac{m}{v}\\
\end{array}
\end{array}
if m < 2.05e-153Initial program 100.0%
*-commutative100.0%
sub-neg100.0%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in m around 0 80.6%
if 2.05e-153 < m Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in m around 0 70.0%
Taylor expanded in m around inf 61.2%
+-commutative61.2%
distribute-lft-in61.2%
div-inv61.3%
*-rgt-identity61.3%
Applied egg-rr61.3%
Final simplification65.9%
(FPCore (m v) :precision binary64 (if (<= m 1.2e-46) -1.0 m))
double code(double m, double v) {
double tmp;
if (m <= 1.2e-46) {
tmp = -1.0;
} else {
tmp = 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.2d-46) then
tmp = -1.0d0
else
tmp = m
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.2e-46) {
tmp = -1.0;
} else {
tmp = m;
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.2e-46: tmp = -1.0 else: tmp = m return tmp
function code(m, v) tmp = 0.0 if (m <= 1.2e-46) tmp = -1.0; else tmp = m; end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.2e-46) tmp = -1.0; else tmp = m; end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.2e-46], -1.0, m]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.2 \cdot 10^{-46}:\\
\;\;\;\;-1\\
\mathbf{else}:\\
\;\;\;\;m\\
\end{array}
\end{array}
if m < 1.20000000000000007e-46Initial program 100.0%
*-commutative100.0%
sub-neg100.0%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in m around 0 63.1%
if 1.20000000000000007e-46 < m Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in m around 0 62.3%
Taylor expanded in m around inf 61.9%
Taylor expanded in v around inf 5.5%
Final simplification28.3%
(FPCore (m v) :precision binary64 (+ -1.0 (/ m v)))
double code(double m, double v) {
return -1.0 + (m / v);
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
code = (-1.0d0) + (m / v)
end function
public static double code(double m, double v) {
return -1.0 + (m / v);
}
def code(m, v): return -1.0 + (m / v)
function code(m, v) return Float64(-1.0 + Float64(m / v)) end
function tmp = code(m, v) tmp = -1.0 + (m / v); end
code[m_, v_] := N[(-1.0 + N[(m / v), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-1 + \frac{m}{v}
\end{array}
Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in m around 0 77.1%
+-commutative51.2%
distribute-lft-in51.2%
div-inv51.3%
*-rgt-identity51.3%
Applied egg-rr77.2%
Taylor expanded in v around 0 77.2%
Final simplification77.2%
(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.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in v around inf 28.0%
neg-mul-128.0%
neg-sub028.0%
associate--r-28.0%
metadata-eval28.0%
Simplified28.0%
Final simplification28.0%
(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.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in m around 0 25.3%
Final simplification25.3%
herbie shell --seed 2024067
(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)))