
(FPCore (m v) :precision binary64 (* (- (/ (* m (- 1.0 m)) v) 1.0) m))
double code(double m, double v) {
return (((m * (1.0 - m)) / v) - 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) * m
end function
public static double code(double m, double v) {
return (((m * (1.0 - m)) / v) - 1.0) * m;
}
def code(m, v): return (((m * (1.0 - m)) / v) - 1.0) * m
function code(m, v) return Float64(Float64(Float64(Float64(m * Float64(1.0 - m)) / v) - 1.0) * m) end
function tmp = code(m, v) tmp = (((m * (1.0 - m)) / v) - 1.0) * m; end
code[m_, v_] := N[(N[(N[(N[(m * N[(1.0 - m), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision] - 1.0), $MachinePrecision] * m), $MachinePrecision]
\begin{array}{l}
\\
\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot m
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 12 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (m v) :precision binary64 (* (- (/ (* m (- 1.0 m)) v) 1.0) m))
double code(double m, double v) {
return (((m * (1.0 - m)) / v) - 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) * m
end function
public static double code(double m, double v) {
return (((m * (1.0 - m)) / v) - 1.0) * m;
}
def code(m, v): return (((m * (1.0 - m)) / v) - 1.0) * m
function code(m, v) return Float64(Float64(Float64(Float64(m * Float64(1.0 - m)) / v) - 1.0) * m) end
function tmp = code(m, v) tmp = (((m * (1.0 - m)) / v) - 1.0) * m; end
code[m_, v_] := N[(N[(N[(N[(m * N[(1.0 - m), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision] - 1.0), $MachinePrecision] * m), $MachinePrecision]
\begin{array}{l}
\\
\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot m
\end{array}
(FPCore (m v) :precision binary64 (* m (+ (/ (* m (- 1.0 m)) v) -1.0)))
double code(double m, double v) {
return m * (((m * (1.0 - m)) / v) + -1.0);
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
code = m * (((m * (1.0d0 - m)) / v) + (-1.0d0))
end function
public static double code(double m, double v) {
return m * (((m * (1.0 - m)) / v) + -1.0);
}
def code(m, v): return m * (((m * (1.0 - m)) / v) + -1.0)
function code(m, v) return Float64(m * Float64(Float64(Float64(m * Float64(1.0 - m)) / v) + -1.0)) end
function tmp = code(m, v) tmp = m * (((m * (1.0 - m)) / v) + -1.0); end
code[m_, v_] := N[(m * N[(N[(N[(m * N[(1.0 - m), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
m \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} + -1\right)
\end{array}
Initial program 99.8%
Final simplification99.8%
(FPCore (m v)
:precision binary64
(let* ((t_0 (* m (/ m v))))
(if (<= m 6.4e-224)
(- m)
(if (<= m 3.8e-157)
t_0
(if (<= m 6e-153)
(- m)
(if (<= m 1.9e-127)
t_0
(if (<= m 8.6e-116)
(- m)
(if (<= m 4.7e-14) (/ m (/ v m)) (* m (/ m (- v)))))))))))
double code(double m, double v) {
double t_0 = m * (m / v);
double tmp;
if (m <= 6.4e-224) {
tmp = -m;
} else if (m <= 3.8e-157) {
tmp = t_0;
} else if (m <= 6e-153) {
tmp = -m;
} else if (m <= 1.9e-127) {
tmp = t_0;
} else if (m <= 8.6e-116) {
tmp = -m;
} else if (m <= 4.7e-14) {
tmp = m / (v / m);
} 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) :: t_0
real(8) :: tmp
t_0 = m * (m / v)
if (m <= 6.4d-224) then
tmp = -m
else if (m <= 3.8d-157) then
tmp = t_0
else if (m <= 6d-153) then
tmp = -m
else if (m <= 1.9d-127) then
tmp = t_0
else if (m <= 8.6d-116) then
tmp = -m
else if (m <= 4.7d-14) then
tmp = m / (v / m)
else
tmp = m * (m / -v)
end if
code = tmp
end function
public static double code(double m, double v) {
double t_0 = m * (m / v);
double tmp;
if (m <= 6.4e-224) {
tmp = -m;
} else if (m <= 3.8e-157) {
tmp = t_0;
} else if (m <= 6e-153) {
tmp = -m;
} else if (m <= 1.9e-127) {
tmp = t_0;
} else if (m <= 8.6e-116) {
tmp = -m;
} else if (m <= 4.7e-14) {
tmp = m / (v / m);
} else {
tmp = m * (m / -v);
}
return tmp;
}
def code(m, v): t_0 = m * (m / v) tmp = 0 if m <= 6.4e-224: tmp = -m elif m <= 3.8e-157: tmp = t_0 elif m <= 6e-153: tmp = -m elif m <= 1.9e-127: tmp = t_0 elif m <= 8.6e-116: tmp = -m elif m <= 4.7e-14: tmp = m / (v / m) else: tmp = m * (m / -v) return tmp
function code(m, v) t_0 = Float64(m * Float64(m / v)) tmp = 0.0 if (m <= 6.4e-224) tmp = Float64(-m); elseif (m <= 3.8e-157) tmp = t_0; elseif (m <= 6e-153) tmp = Float64(-m); elseif (m <= 1.9e-127) tmp = t_0; elseif (m <= 8.6e-116) tmp = Float64(-m); elseif (m <= 4.7e-14) tmp = Float64(m / Float64(v / m)); else tmp = Float64(m * Float64(m / Float64(-v))); end return tmp end
function tmp_2 = code(m, v) t_0 = m * (m / v); tmp = 0.0; if (m <= 6.4e-224) tmp = -m; elseif (m <= 3.8e-157) tmp = t_0; elseif (m <= 6e-153) tmp = -m; elseif (m <= 1.9e-127) tmp = t_0; elseif (m <= 8.6e-116) tmp = -m; elseif (m <= 4.7e-14) tmp = m / (v / m); else tmp = m * (m / -v); end tmp_2 = tmp; end
code[m_, v_] := Block[{t$95$0 = N[(m * N[(m / v), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[m, 6.4e-224], (-m), If[LessEqual[m, 3.8e-157], t$95$0, If[LessEqual[m, 6e-153], (-m), If[LessEqual[m, 1.9e-127], t$95$0, If[LessEqual[m, 8.6e-116], (-m), If[LessEqual[m, 4.7e-14], N[(m / N[(v / m), $MachinePrecision]), $MachinePrecision], N[(m * N[(m / (-v)), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := m \cdot \frac{m}{v}\\
\mathbf{if}\;m \leq 6.4 \cdot 10^{-224}:\\
\;\;\;\;-m\\
\mathbf{elif}\;m \leq 3.8 \cdot 10^{-157}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;m \leq 6 \cdot 10^{-153}:\\
\;\;\;\;-m\\
\mathbf{elif}\;m \leq 1.9 \cdot 10^{-127}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;m \leq 8.6 \cdot 10^{-116}:\\
\;\;\;\;-m\\
\mathbf{elif}\;m \leq 4.7 \cdot 10^{-14}:\\
\;\;\;\;\frac{m}{\frac{v}{m}}\\
\mathbf{else}:\\
\;\;\;\;m \cdot \frac{m}{-v}\\
\end{array}
\end{array}
if m < 6.3999999999999997e-224 or 3.8000000000000002e-157 < m < 6e-153 or 1.90000000000000001e-127 < m < 8.5999999999999994e-116Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in m around 0 81.7%
neg-mul-181.7%
Simplified81.7%
if 6.3999999999999997e-224 < m < 3.8000000000000002e-157 or 6e-153 < m < 1.90000000000000001e-127Initial program 99.8%
*-commutative99.8%
sub-neg99.8%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
clear-num99.8%
un-div-inv99.8%
Applied egg-rr99.8%
Taylor expanded in v around 0 61.6%
Taylor expanded in m around 0 61.6%
if 8.5999999999999994e-116 < m < 4.7000000000000002e-14Initial program 99.6%
*-commutative99.6%
sub-neg99.6%
associate-/l*99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in v around 0 89.1%
associate-/l*89.1%
Simplified89.1%
unpow289.1%
associate-*r*89.0%
clear-num89.0%
div-inv89.1%
clear-num89.0%
un-div-inv89.2%
Applied egg-rr89.2%
Taylor expanded in m around 0 89.3%
if 4.7000000000000002e-14 < m Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
clear-num99.8%
un-div-inv99.8%
Applied egg-rr99.8%
Taylor expanded in v around 0 99.9%
associate-*r/99.9%
*-commutative99.9%
Applied egg-rr99.9%
Taylor expanded in m around 0 3.0%
associate-*l/3.0%
*-un-lft-identity3.0%
frac-2neg3.0%
remove-double-div3.0%
metadata-eval3.0%
frac-2neg3.0%
add-sqr-sqrt0.0%
sqrt-unprod73.0%
frac-times73.0%
metadata-eval73.0%
metadata-eval73.0%
frac-times73.0%
sqrt-unprod73.0%
add-sqr-sqrt73.0%
clear-num73.0%
/-rgt-identity73.0%
Applied egg-rr73.0%
(FPCore (m v) :precision binary64 (if (<= m 1.65e-110) (* m (+ -1.0 (* m (/ 1.0 v)))) (/ m (/ v (* m (- 1.0 m))))))
double code(double m, double v) {
double tmp;
if (m <= 1.65e-110) {
tmp = m * (-1.0 + (m * (1.0 / v)));
} else {
tmp = m / (v / (m * (1.0 - 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.65d-110) then
tmp = m * ((-1.0d0) + (m * (1.0d0 / v)))
else
tmp = m / (v / (m * (1.0d0 - m)))
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.65e-110) {
tmp = m * (-1.0 + (m * (1.0 / v)));
} else {
tmp = m / (v / (m * (1.0 - m)));
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.65e-110: tmp = m * (-1.0 + (m * (1.0 / v))) else: tmp = m / (v / (m * (1.0 - m))) return tmp
function code(m, v) tmp = 0.0 if (m <= 1.65e-110) tmp = Float64(m * Float64(-1.0 + Float64(m * Float64(1.0 / v)))); else tmp = Float64(m / Float64(v / Float64(m * Float64(1.0 - m)))); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.65e-110) tmp = m * (-1.0 + (m * (1.0 / v))); else tmp = m / (v / (m * (1.0 - m))); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.65e-110], N[(m * N[(-1.0 + N[(m * N[(1.0 / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(m / N[(v / N[(m * N[(1.0 - m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.65 \cdot 10^{-110}:\\
\;\;\;\;m \cdot \left(-1 + m \cdot \frac{1}{v}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{m}{\frac{v}{m \cdot \left(1 - m\right)}}\\
\end{array}
\end{array}
if m < 1.65e-110Initial program 99.8%
*-commutative99.8%
sub-neg99.8%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in m around 0 99.9%
if 1.65e-110 < m Initial program 99.8%
*-commutative99.8%
sub-neg99.8%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in v around 0 98.1%
associate-/l*98.1%
Simplified98.1%
unpow298.1%
associate-*r*98.1%
clear-num98.0%
div-inv98.1%
clear-num98.1%
un-div-inv98.2%
Applied egg-rr98.2%
Taylor expanded in v around 0 98.2%
Final simplification98.7%
(FPCore (m v) :precision binary64 (if (<= m 1.65e-110) (* m (+ -1.0 (* m (/ 1.0 v)))) (* m (/ (* m (- 1.0 m)) v))))
double code(double m, double v) {
double tmp;
if (m <= 1.65e-110) {
tmp = m * (-1.0 + (m * (1.0 / v)));
} else {
tmp = m * ((m * (1.0 - 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.65d-110) then
tmp = m * ((-1.0d0) + (m * (1.0d0 / v)))
else
tmp = m * ((m * (1.0d0 - m)) / v)
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.65e-110) {
tmp = m * (-1.0 + (m * (1.0 / v)));
} else {
tmp = m * ((m * (1.0 - m)) / v);
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.65e-110: tmp = m * (-1.0 + (m * (1.0 / v))) else: tmp = m * ((m * (1.0 - m)) / v) return tmp
function code(m, v) tmp = 0.0 if (m <= 1.65e-110) tmp = Float64(m * Float64(-1.0 + Float64(m * Float64(1.0 / v)))); else tmp = Float64(m * Float64(Float64(m * Float64(1.0 - m)) / v)); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.65e-110) tmp = m * (-1.0 + (m * (1.0 / v))); else tmp = m * ((m * (1.0 - m)) / v); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.65e-110], N[(m * N[(-1.0 + N[(m * N[(1.0 / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(m * N[(N[(m * N[(1.0 - m), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.65 \cdot 10^{-110}:\\
\;\;\;\;m \cdot \left(-1 + m \cdot \frac{1}{v}\right)\\
\mathbf{else}:\\
\;\;\;\;m \cdot \frac{m \cdot \left(1 - m\right)}{v}\\
\end{array}
\end{array}
if m < 1.65e-110Initial program 99.8%
*-commutative99.8%
sub-neg99.8%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in m around 0 99.9%
if 1.65e-110 < m Initial program 99.8%
*-commutative99.8%
sub-neg99.8%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
clear-num99.7%
un-div-inv99.8%
Applied egg-rr99.8%
Taylor expanded in v around 0 98.1%
Final simplification98.7%
(FPCore (m v) :precision binary64 (if (<= m 1.65e-110) (* m (+ -1.0 (* m (/ 1.0 v)))) (* m (* m (/ (- 1.0 m) v)))))
double code(double m, double v) {
double tmp;
if (m <= 1.65e-110) {
tmp = m * (-1.0 + (m * (1.0 / v)));
} else {
tmp = m * (m * ((1.0 - 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.65d-110) then
tmp = m * ((-1.0d0) + (m * (1.0d0 / v)))
else
tmp = m * (m * ((1.0d0 - m) / v))
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.65e-110) {
tmp = m * (-1.0 + (m * (1.0 / v)));
} else {
tmp = m * (m * ((1.0 - m) / v));
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.65e-110: tmp = m * (-1.0 + (m * (1.0 / v))) else: tmp = m * (m * ((1.0 - m) / v)) return tmp
function code(m, v) tmp = 0.0 if (m <= 1.65e-110) tmp = Float64(m * Float64(-1.0 + Float64(m * Float64(1.0 / v)))); else tmp = Float64(m * Float64(m * Float64(Float64(1.0 - m) / v))); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.65e-110) tmp = m * (-1.0 + (m * (1.0 / v))); else tmp = m * (m * ((1.0 - m) / v)); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.65e-110], N[(m * N[(-1.0 + N[(m * N[(1.0 / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(m * N[(m * N[(N[(1.0 - m), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.65 \cdot 10^{-110}:\\
\;\;\;\;m \cdot \left(-1 + m \cdot \frac{1}{v}\right)\\
\mathbf{else}:\\
\;\;\;\;m \cdot \left(m \cdot \frac{1 - m}{v}\right)\\
\end{array}
\end{array}
if m < 1.65e-110Initial program 99.8%
*-commutative99.8%
sub-neg99.8%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in m around 0 99.9%
if 1.65e-110 < m Initial program 99.8%
*-commutative99.8%
sub-neg99.8%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
clear-num99.7%
un-div-inv99.8%
Applied egg-rr99.8%
Taylor expanded in v around 0 98.1%
associate-*r/98.1%
Simplified98.1%
Final simplification98.7%
(FPCore (m v) :precision binary64 (if (<= m 1.65e-110) (* m (+ -1.0 (/ m v))) (* m (* m (/ (- 1.0 m) v)))))
double code(double m, double v) {
double tmp;
if (m <= 1.65e-110) {
tmp = m * (-1.0 + (m / v));
} else {
tmp = m * (m * ((1.0 - 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.65d-110) then
tmp = m * ((-1.0d0) + (m / v))
else
tmp = m * (m * ((1.0d0 - m) / v))
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 1.65e-110) {
tmp = m * (-1.0 + (m / v));
} else {
tmp = m * (m * ((1.0 - m) / v));
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 1.65e-110: tmp = m * (-1.0 + (m / v)) else: tmp = m * (m * ((1.0 - m) / v)) return tmp
function code(m, v) tmp = 0.0 if (m <= 1.65e-110) tmp = Float64(m * Float64(-1.0 + Float64(m / v))); else tmp = Float64(m * Float64(m * Float64(Float64(1.0 - m) / v))); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 1.65e-110) tmp = m * (-1.0 + (m / v)); else tmp = m * (m * ((1.0 - m) / v)); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 1.65e-110], N[(m * N[(-1.0 + N[(m / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(m * N[(m * N[(N[(1.0 - m), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.65 \cdot 10^{-110}:\\
\;\;\;\;m \cdot \left(-1 + \frac{m}{v}\right)\\
\mathbf{else}:\\
\;\;\;\;m \cdot \left(m \cdot \frac{1 - m}{v}\right)\\
\end{array}
\end{array}
if m < 1.65e-110Initial program 99.8%
Taylor expanded in m around 0 99.8%
if 1.65e-110 < m Initial program 99.8%
*-commutative99.8%
sub-neg99.8%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
clear-num99.7%
un-div-inv99.8%
Applied egg-rr99.8%
Taylor expanded in v around 0 98.1%
associate-*r/98.1%
Simplified98.1%
Final simplification98.7%
(FPCore (m v) :precision binary64 (if (<= m 3.8e-10) (* m (+ -1.0 (/ m v))) (/ m (/ (/ v (- m)) m))))
double code(double m, double v) {
double tmp;
if (m <= 3.8e-10) {
tmp = m * (-1.0 + (m / v));
} else {
tmp = m / ((v / -m) / m);
}
return tmp;
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
real(8) :: tmp
if (m <= 3.8d-10) then
tmp = m * ((-1.0d0) + (m / v))
else
tmp = m / ((v / -m) / m)
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 3.8e-10) {
tmp = m * (-1.0 + (m / v));
} else {
tmp = m / ((v / -m) / m);
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 3.8e-10: tmp = m * (-1.0 + (m / v)) else: tmp = m / ((v / -m) / m) return tmp
function code(m, v) tmp = 0.0 if (m <= 3.8e-10) tmp = Float64(m * Float64(-1.0 + Float64(m / v))); else tmp = Float64(m / Float64(Float64(v / Float64(-m)) / m)); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 3.8e-10) tmp = m * (-1.0 + (m / v)); else tmp = m / ((v / -m) / m); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 3.8e-10], N[(m * N[(-1.0 + N[(m / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(m / N[(N[(v / (-m)), $MachinePrecision] / m), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 3.8 \cdot 10^{-10}:\\
\;\;\;\;m \cdot \left(-1 + \frac{m}{v}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{m}{\frac{\frac{v}{-m}}{m}}\\
\end{array}
\end{array}
if m < 3.7999999999999998e-10Initial program 99.8%
Taylor expanded in m around 0 99.3%
if 3.7999999999999998e-10 < 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%
associate-/l*99.9%
Simplified99.9%
unpow299.9%
associate-*r*99.9%
clear-num99.8%
div-inv99.9%
clear-num99.9%
un-div-inv99.9%
Applied egg-rr99.9%
Taylor expanded in m around inf 96.8%
associate-*r/96.8%
mul-1-neg96.8%
Simplified96.8%
Final simplification98.0%
(FPCore (m v) :precision binary64 (if (<= m 5e-11) (* m (+ -1.0 (/ m v))) (* m (/ m (- v)))))
double code(double m, double v) {
double tmp;
if (m <= 5e-11) {
tmp = m * (-1.0 + (m / v));
} 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 <= 5d-11) then
tmp = m * ((-1.0d0) + (m / v))
else
tmp = m * (m / -v)
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (m <= 5e-11) {
tmp = m * (-1.0 + (m / v));
} else {
tmp = m * (m / -v);
}
return tmp;
}
def code(m, v): tmp = 0 if m <= 5e-11: tmp = m * (-1.0 + (m / v)) else: tmp = m * (m / -v) return tmp
function code(m, v) tmp = 0.0 if (m <= 5e-11) tmp = Float64(m * Float64(-1.0 + Float64(m / v))); else tmp = Float64(m * Float64(m / Float64(-v))); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (m <= 5e-11) tmp = m * (-1.0 + (m / v)); else tmp = m * (m / -v); end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[m, 5e-11], N[(m * N[(-1.0 + N[(m / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(m * N[(m / (-v)), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 5 \cdot 10^{-11}:\\
\;\;\;\;m \cdot \left(-1 + \frac{m}{v}\right)\\
\mathbf{else}:\\
\;\;\;\;m \cdot \frac{m}{-v}\\
\end{array}
\end{array}
if m < 5.00000000000000018e-11Initial program 99.8%
Taylor expanded in m around 0 99.5%
if 5.00000000000000018e-11 < m Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.9%
metadata-eval99.9%
Simplified99.9%
clear-num99.8%
un-div-inv99.8%
Applied egg-rr99.8%
Taylor expanded in v around 0 99.9%
associate-*r/99.9%
*-commutative99.9%
Applied egg-rr99.9%
Taylor expanded in m around 0 1.8%
associate-*l/1.8%
*-un-lft-identity1.8%
frac-2neg1.8%
remove-double-div1.8%
metadata-eval1.8%
frac-2neg1.8%
add-sqr-sqrt0.0%
sqrt-unprod74.1%
frac-times74.1%
metadata-eval74.1%
metadata-eval74.1%
frac-times74.1%
sqrt-unprod74.1%
add-sqr-sqrt74.1%
clear-num74.1%
/-rgt-identity74.1%
Applied egg-rr74.1%
Final simplification86.4%
(FPCore (m v) :precision binary64 (* m (+ (/ m (/ v (- 1.0 m))) -1.0)))
double code(double m, double v) {
return 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 = m * ((m / (v / (1.0d0 - m))) + (-1.0d0))
end function
public static double code(double m, double v) {
return m * ((m / (v / (1.0 - m))) + -1.0);
}
def code(m, v): return m * ((m / (v / (1.0 - m))) + -1.0)
function code(m, v) return Float64(m * Float64(Float64(m / Float64(v / Float64(1.0 - m))) + -1.0)) end
function tmp = code(m, v) tmp = m * ((m / (v / (1.0 - m))) + -1.0); end
code[m_, v_] := N[(m * N[(N[(m / N[(v / N[(1.0 - m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
m \cdot \left(\frac{m}{\frac{v}{1 - m}} + -1\right)
\end{array}
Initial program 99.8%
*-commutative99.8%
sub-neg99.8%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
clear-num99.8%
un-div-inv99.8%
Applied egg-rr99.8%
(FPCore (m v) :precision binary64 (* m (+ -1.0 (* m (/ (- 1.0 m) v)))))
double code(double m, double v) {
return m * (-1.0 + (m * ((1.0 - m) / v)));
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
code = m * ((-1.0d0) + (m * ((1.0d0 - m) / v)))
end function
public static double code(double m, double v) {
return m * (-1.0 + (m * ((1.0 - m) / v)));
}
def code(m, v): return m * (-1.0 + (m * ((1.0 - m) / v)))
function code(m, v) return Float64(m * Float64(-1.0 + Float64(m * Float64(Float64(1.0 - m) / v)))) end
function tmp = code(m, v) tmp = m * (-1.0 + (m * ((1.0 - m) / v))); end
code[m_, v_] := N[(m * N[(-1.0 + N[(m * N[(N[(1.0 - m), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
m \cdot \left(-1 + m \cdot \frac{1 - m}{v}\right)
\end{array}
Initial program 99.8%
*-commutative99.8%
sub-neg99.8%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
Final simplification99.8%
(FPCore (m v) :precision binary64 (if (<= v 4.2e-208) (* m (/ m v)) (- m)))
double code(double m, double v) {
double tmp;
if (v <= 4.2e-208) {
tmp = m * (m / v);
} 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 (v <= 4.2d-208) then
tmp = m * (m / v)
else
tmp = -m
end if
code = tmp
end function
public static double code(double m, double v) {
double tmp;
if (v <= 4.2e-208) {
tmp = m * (m / v);
} else {
tmp = -m;
}
return tmp;
}
def code(m, v): tmp = 0 if v <= 4.2e-208: tmp = m * (m / v) else: tmp = -m return tmp
function code(m, v) tmp = 0.0 if (v <= 4.2e-208) tmp = Float64(m * Float64(m / v)); else tmp = Float64(-m); end return tmp end
function tmp_2 = code(m, v) tmp = 0.0; if (v <= 4.2e-208) tmp = m * (m / v); else tmp = -m; end tmp_2 = tmp; end
code[m_, v_] := If[LessEqual[v, 4.2e-208], N[(m * N[(m / v), $MachinePrecision]), $MachinePrecision], (-m)]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq 4.2 \cdot 10^{-208}:\\
\;\;\;\;m \cdot \frac{m}{v}\\
\mathbf{else}:\\
\;\;\;\;-m\\
\end{array}
\end{array}
if v < 4.20000000000000024e-208Initial program 99.8%
*-commutative99.8%
sub-neg99.8%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
clear-num99.8%
un-div-inv99.8%
Applied egg-rr99.8%
Taylor expanded in v around 0 97.5%
Taylor expanded in m around 0 46.6%
if 4.20000000000000024e-208 < v Initial program 99.9%
*-commutative99.9%
sub-neg99.9%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in m around 0 37.2%
neg-mul-137.2%
Simplified37.2%
(FPCore (m v) :precision binary64 (- m))
double code(double m, double v) {
return -m;
}
real(8) function code(m, v)
real(8), intent (in) :: m
real(8), intent (in) :: v
code = -m
end function
public static double code(double m, double v) {
return -m;
}
def code(m, v): return -m
function code(m, v) return Float64(-m) end
function tmp = code(m, v) tmp = -m; end
code[m_, v_] := (-m)
\begin{array}{l}
\\
-m
\end{array}
Initial program 99.8%
*-commutative99.8%
sub-neg99.8%
associate-/l*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in m around 0 25.2%
neg-mul-125.2%
Simplified25.2%
herbie shell --seed 2024107
(FPCore (m v)
:name "a 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) m))