| Alternative 1 | |
|---|---|
| Error | 14.1 |
| Cost | 1172 |
(FPCore (u v t1) :precision binary64 (/ (* (- t1) v) (* (+ t1 u) (+ t1 u))))
(FPCore (u v t1)
:precision binary64
(let* ((t_1 (- -1.0 (/ u t1))))
(if (or (<= t1 -3.9e-104) (not (<= t1 1.04e-200)))
(/ (/ v (+ u t1)) t_1)
(/ v (* t_1 (+ u t1))))))double code(double u, double v, double t1) {
return (-t1 * v) / ((t1 + u) * (t1 + u));
}
double code(double u, double v, double t1) {
double t_1 = -1.0 - (u / t1);
double tmp;
if ((t1 <= -3.9e-104) || !(t1 <= 1.04e-200)) {
tmp = (v / (u + t1)) / t_1;
} else {
tmp = v / (t_1 * (u + t1));
}
return tmp;
}
real(8) function code(u, v, t1)
real(8), intent (in) :: u
real(8), intent (in) :: v
real(8), intent (in) :: t1
code = (-t1 * v) / ((t1 + u) * (t1 + u))
end function
real(8) function code(u, v, t1)
real(8), intent (in) :: u
real(8), intent (in) :: v
real(8), intent (in) :: t1
real(8) :: t_1
real(8) :: tmp
t_1 = (-1.0d0) - (u / t1)
if ((t1 <= (-3.9d-104)) .or. (.not. (t1 <= 1.04d-200))) then
tmp = (v / (u + t1)) / t_1
else
tmp = v / (t_1 * (u + t1))
end if
code = tmp
end function
public static double code(double u, double v, double t1) {
return (-t1 * v) / ((t1 + u) * (t1 + u));
}
public static double code(double u, double v, double t1) {
double t_1 = -1.0 - (u / t1);
double tmp;
if ((t1 <= -3.9e-104) || !(t1 <= 1.04e-200)) {
tmp = (v / (u + t1)) / t_1;
} else {
tmp = v / (t_1 * (u + t1));
}
return tmp;
}
def code(u, v, t1): return (-t1 * v) / ((t1 + u) * (t1 + u))
def code(u, v, t1): t_1 = -1.0 - (u / t1) tmp = 0 if (t1 <= -3.9e-104) or not (t1 <= 1.04e-200): tmp = (v / (u + t1)) / t_1 else: tmp = v / (t_1 * (u + t1)) return tmp
function code(u, v, t1) return Float64(Float64(Float64(-t1) * v) / Float64(Float64(t1 + u) * Float64(t1 + u))) end
function code(u, v, t1) t_1 = Float64(-1.0 - Float64(u / t1)) tmp = 0.0 if ((t1 <= -3.9e-104) || !(t1 <= 1.04e-200)) tmp = Float64(Float64(v / Float64(u + t1)) / t_1); else tmp = Float64(v / Float64(t_1 * Float64(u + t1))); end return tmp end
function tmp = code(u, v, t1) tmp = (-t1 * v) / ((t1 + u) * (t1 + u)); end
function tmp_2 = code(u, v, t1) t_1 = -1.0 - (u / t1); tmp = 0.0; if ((t1 <= -3.9e-104) || ~((t1 <= 1.04e-200))) tmp = (v / (u + t1)) / t_1; else tmp = v / (t_1 * (u + t1)); end tmp_2 = tmp; end
code[u_, v_, t1_] := N[(N[((-t1) * v), $MachinePrecision] / N[(N[(t1 + u), $MachinePrecision] * N[(t1 + u), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[u_, v_, t1_] := Block[{t$95$1 = N[(-1.0 - N[(u / t1), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t1, -3.9e-104], N[Not[LessEqual[t1, 1.04e-200]], $MachinePrecision]], N[(N[(v / N[(u + t1), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision], N[(v / N[(t$95$1 * N[(u + t1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\frac{\left(-t1\right) \cdot v}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\begin{array}{l}
t_1 := -1 - \frac{u}{t1}\\
\mathbf{if}\;t1 \leq -3.9 \cdot 10^{-104} \lor \neg \left(t1 \leq 1.04 \cdot 10^{-200}\right):\\
\;\;\;\;\frac{\frac{v}{u + t1}}{t_1}\\
\mathbf{else}:\\
\;\;\;\;\frac{v}{t_1 \cdot \left(u + t1\right)}\\
\end{array}
Results
if t1 < -3.9000000000000002e-104 or 1.03999999999999993e-200 < t1 Initial program 19.3
Simplified0.4
[Start]19.3 | \[ \frac{\left(-t1\right) \cdot v}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\] |
|---|---|
*-commutative [=>]19.3 | \[ \frac{\color{blue}{v \cdot \left(-t1\right)}}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\] |
associate-/l* [=>]16.4 | \[ \color{blue}{\frac{v}{\frac{\left(t1 + u\right) \cdot \left(t1 + u\right)}{-t1}}}
\] |
associate-*r/ [<=]3.0 | \[ \frac{v}{\color{blue}{\left(t1 + u\right) \cdot \frac{t1 + u}{-t1}}}
\] |
associate-/r* [=>]0.4 | \[ \color{blue}{\frac{\frac{v}{t1 + u}}{\frac{t1 + u}{-t1}}}
\] |
neg-mul-1 [=>]0.4 | \[ \frac{\frac{v}{t1 + u}}{\frac{t1 + u}{\color{blue}{-1 \cdot t1}}}
\] |
associate-/l/ [<=]0.4 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\frac{\frac{t1 + u}{t1}}{-1}}}
\] |
metadata-eval [<=]0.4 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{0 - 1}}}
\] |
mul0-lft [<=]3.3 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{0 \cdot \frac{t1 + u}{t1}} - 1}}
\] |
associate-*r/ [=>]0.4 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{0 \cdot \left(t1 + u\right)}{t1}} - 1}}
\] |
mul0-lft [=>]0.4 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{0}}{t1} - 1}}
\] |
*-inverses [<=]0.4 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{0}{t1} - \color{blue}{\frac{t1}{t1}}}}
\] |
div-sub [<=]0.4 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{0 - t1}{t1}}}}
\] |
neg-sub0 [<=]0.4 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{-t1}}{t1}}}
\] |
neg-mul-1 [=>]0.4 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{-1 \cdot t1}}{t1}}}
\] |
*-commutative [=>]0.4 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{t1 \cdot -1}}{t1}}}
\] |
associate-/l* [=>]0.4 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{t1}{\frac{t1}{-1}}}}}
\] |
associate-/l* [<=]0.4 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\frac{\frac{t1 + u}{t1} \cdot \frac{t1}{-1}}{t1}}}
\] |
*-commutative [=>]0.4 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1}{-1} \cdot \frac{t1 + u}{t1}}}{t1}}
\] |
times-frac [<=]14.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1 \cdot \left(t1 + u\right)}{-1 \cdot t1}}}{t1}}
\] |
neg-mul-1 [<=]14.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 \cdot \left(t1 + u\right)}{\color{blue}{-t1}}}{t1}}
\] |
associate-/l* [=>]0.4 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1}{\frac{-t1}{t1 + u}}}}{t1}}
\] |
if -3.9000000000000002e-104 < t1 < 1.03999999999999993e-200Initial program 15.0
Simplified5.0
[Start]15.0 | \[ \frac{\left(-t1\right) \cdot v}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\] |
|---|---|
*-commutative [=>]15.0 | \[ \frac{\color{blue}{v \cdot \left(-t1\right)}}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\] |
associate-/l* [=>]14.6 | \[ \color{blue}{\frac{v}{\frac{\left(t1 + u\right) \cdot \left(t1 + u\right)}{-t1}}}
\] |
associate-*r/ [<=]4.2 | \[ \frac{v}{\color{blue}{\left(t1 + u\right) \cdot \frac{t1 + u}{-t1}}}
\] |
associate-/r* [=>]5.0 | \[ \color{blue}{\frac{\frac{v}{t1 + u}}{\frac{t1 + u}{-t1}}}
\] |
neg-mul-1 [=>]5.0 | \[ \frac{\frac{v}{t1 + u}}{\frac{t1 + u}{\color{blue}{-1 \cdot t1}}}
\] |
associate-/l/ [<=]5.0 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\frac{\frac{t1 + u}{t1}}{-1}}}
\] |
metadata-eval [<=]5.0 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{0 - 1}}}
\] |
mul0-lft [<=]28.3 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{0 \cdot \frac{t1 + u}{t1}} - 1}}
\] |
associate-*r/ [=>]5.0 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{0 \cdot \left(t1 + u\right)}{t1}} - 1}}
\] |
mul0-lft [=>]5.0 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{0}}{t1} - 1}}
\] |
*-inverses [<=]5.0 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{0}{t1} - \color{blue}{\frac{t1}{t1}}}}
\] |
div-sub [<=]5.0 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{0 - t1}{t1}}}}
\] |
neg-sub0 [<=]5.0 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{-t1}}{t1}}}
\] |
neg-mul-1 [=>]5.0 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{-1 \cdot t1}}{t1}}}
\] |
*-commutative [=>]5.0 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{t1 \cdot -1}}{t1}}}
\] |
associate-/l* [=>]5.0 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{t1}{\frac{t1}{-1}}}}}
\] |
associate-/l* [<=]5.0 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\frac{\frac{t1 + u}{t1} \cdot \frac{t1}{-1}}{t1}}}
\] |
*-commutative [=>]5.0 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1}{-1} \cdot \frac{t1 + u}{t1}}}{t1}}
\] |
times-frac [<=]20.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1 \cdot \left(t1 + u\right)}{-1 \cdot t1}}}{t1}}
\] |
neg-mul-1 [<=]20.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 \cdot \left(t1 + u\right)}{\color{blue}{-t1}}}{t1}}
\] |
associate-/l* [=>]5.1 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1}{\frac{-t1}{t1 + u}}}}{t1}}
\] |
Taylor expanded in v around 0 4.2
Simplified4.2
[Start]4.2 | \[ -1 \cdot \frac{v}{\left(t1 + u\right) \cdot \left(1 + \frac{u}{t1}\right)}
\] |
|---|---|
associate-*r/ [=>]4.2 | \[ \color{blue}{\frac{-1 \cdot v}{\left(t1 + u\right) \cdot \left(1 + \frac{u}{t1}\right)}}
\] |
+-commutative [=>]4.2 | \[ \frac{-1 \cdot v}{\left(t1 + u\right) \cdot \color{blue}{\left(\frac{u}{t1} + 1\right)}}
\] |
*-commutative [<=]4.2 | \[ \frac{-1 \cdot v}{\color{blue}{\left(\frac{u}{t1} + 1\right) \cdot \left(t1 + u\right)}}
\] |
times-frac [=>]5.1 | \[ \color{blue}{\frac{-1}{\frac{u}{t1} + 1} \cdot \frac{v}{t1 + u}}
\] |
/-rgt-identity [<=]5.1 | \[ \color{blue}{\frac{\frac{-1}{\frac{u}{t1} + 1}}{1}} \cdot \frac{v}{t1 + u}
\] |
*-commutative [=>]5.1 | \[ \color{blue}{\frac{v}{t1 + u} \cdot \frac{\frac{-1}{\frac{u}{t1} + 1}}{1}}
\] |
associate-*r/ [=>]5.1 | \[ \color{blue}{\frac{\frac{v}{t1 + u} \cdot \frac{-1}{\frac{u}{t1} + 1}}{1}}
\] |
associate-/l* [=>]5.1 | \[ \color{blue}{\frac{\frac{v}{t1 + u}}{\frac{1}{\frac{-1}{\frac{u}{t1} + 1}}}}
\] |
associate-/r/ [=>]5.0 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\frac{1}{-1} \cdot \left(\frac{u}{t1} + 1\right)}}
\] |
metadata-eval [=>]5.0 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{-1} \cdot \left(\frac{u}{t1} + 1\right)}
\] |
mul-1-neg [=>]5.0 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{-\left(\frac{u}{t1} + 1\right)}}
\] |
+-commutative [<=]5.0 | \[ \frac{\frac{v}{t1 + u}}{-\color{blue}{\left(1 + \frac{u}{t1}\right)}}
\] |
distribute-neg-in [=>]5.0 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\left(-1\right) + \left(-\frac{u}{t1}\right)}}
\] |
metadata-eval [=>]5.0 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{-1} + \left(-\frac{u}{t1}\right)}
\] |
sub-neg [<=]5.0 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{-1 - \frac{u}{t1}}}
\] |
associate-/r* [<=]4.2 | \[ \color{blue}{\frac{v}{\left(t1 + u\right) \cdot \left(-1 - \frac{u}{t1}\right)}}
\] |
Final simplification1.3
| Alternative 1 | |
|---|---|
| Error | 14.1 |
| Cost | 1172 |
| Alternative 2 | |
|---|---|
| Error | 13.9 |
| Cost | 1105 |
| Alternative 3 | |
|---|---|
| Error | 13.7 |
| Cost | 1105 |
| Alternative 4 | |
|---|---|
| Error | 13.8 |
| Cost | 1105 |
| Alternative 5 | |
|---|---|
| Error | 13.7 |
| Cost | 1104 |
| Alternative 6 | |
|---|---|
| Error | 13.8 |
| Cost | 1042 |
| Alternative 7 | |
|---|---|
| Error | 13.9 |
| Cost | 1041 |
| Alternative 8 | |
|---|---|
| Error | 20.6 |
| Cost | 713 |
| Alternative 9 | |
|---|---|
| Error | 20.6 |
| Cost | 712 |
| Alternative 10 | |
|---|---|
| Error | 3.3 |
| Cost | 704 |
| Alternative 11 | |
|---|---|
| Error | 1.4 |
| Cost | 704 |
| Alternative 12 | |
|---|---|
| Error | 27.7 |
| Cost | 521 |
| Alternative 13 | |
|---|---|
| Error | 24.8 |
| Cost | 384 |
| Alternative 14 | |
|---|---|
| Error | 30.4 |
| Cost | 256 |
herbie shell --seed 2023031
(FPCore (u v t1)
:name "Rosa's DopplerBench"
:precision binary64
(/ (* (- t1) v) (* (+ t1 u) (+ t1 u))))