| Alternative 1 | |
|---|---|
| Error | 1.3 |
| Cost | 2441 |
(FPCore (u v t1) :precision binary64 (/ (* (- t1) v) (* (+ t1 u) (+ t1 u))))
(FPCore (u v t1)
:precision binary64
(let* ((t_1 (/ v (+ t1 u))))
(if (<= t1 -4.5e-250)
(/ t_1 (- -1.0 (/ u t1)))
(if (<= t1 9.5e-230)
(/ (* v (/ t1 u)) (- u))
(* t_1 (/ (- t1) (+ t1 u)))))))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 = v / (t1 + u);
double tmp;
if (t1 <= -4.5e-250) {
tmp = t_1 / (-1.0 - (u / t1));
} else if (t1 <= 9.5e-230) {
tmp = (v * (t1 / u)) / -u;
} else {
tmp = t_1 * (-t1 / (t1 + u));
}
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 = v / (t1 + u)
if (t1 <= (-4.5d-250)) then
tmp = t_1 / ((-1.0d0) - (u / t1))
else if (t1 <= 9.5d-230) then
tmp = (v * (t1 / u)) / -u
else
tmp = t_1 * (-t1 / (t1 + u))
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 = v / (t1 + u);
double tmp;
if (t1 <= -4.5e-250) {
tmp = t_1 / (-1.0 - (u / t1));
} else if (t1 <= 9.5e-230) {
tmp = (v * (t1 / u)) / -u;
} else {
tmp = t_1 * (-t1 / (t1 + u));
}
return tmp;
}
def code(u, v, t1): return (-t1 * v) / ((t1 + u) * (t1 + u))
def code(u, v, t1): t_1 = v / (t1 + u) tmp = 0 if t1 <= -4.5e-250: tmp = t_1 / (-1.0 - (u / t1)) elif t1 <= 9.5e-230: tmp = (v * (t1 / u)) / -u else: tmp = t_1 * (-t1 / (t1 + u)) 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(v / Float64(t1 + u)) tmp = 0.0 if (t1 <= -4.5e-250) tmp = Float64(t_1 / Float64(-1.0 - Float64(u / t1))); elseif (t1 <= 9.5e-230) tmp = Float64(Float64(v * Float64(t1 / u)) / Float64(-u)); else tmp = Float64(t_1 * Float64(Float64(-t1) / Float64(t1 + u))); 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 = v / (t1 + u); tmp = 0.0; if (t1 <= -4.5e-250) tmp = t_1 / (-1.0 - (u / t1)); elseif (t1 <= 9.5e-230) tmp = (v * (t1 / u)) / -u; else tmp = t_1 * (-t1 / (t1 + u)); 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[(v / N[(t1 + u), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t1, -4.5e-250], N[(t$95$1 / N[(-1.0 - N[(u / t1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t1, 9.5e-230], N[(N[(v * N[(t1 / u), $MachinePrecision]), $MachinePrecision] / (-u)), $MachinePrecision], N[(t$95$1 * N[((-t1) / N[(t1 + u), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\frac{\left(-t1\right) \cdot v}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\begin{array}{l}
t_1 := \frac{v}{t1 + u}\\
\mathbf{if}\;t1 \leq -4.5 \cdot 10^{-250}:\\
\;\;\;\;\frac{t_1}{-1 - \frac{u}{t1}}\\
\mathbf{elif}\;t1 \leq 9.5 \cdot 10^{-230}:\\
\;\;\;\;\frac{v \cdot \frac{t1}{u}}{-u}\\
\mathbf{else}:\\
\;\;\;\;t_1 \cdot \frac{-t1}{t1 + u}\\
\end{array}
Results
if t1 < -4.49999999999999993e-250Initial program 18.6
Simplified0.7
[Start]18.6 | \[ \frac{\left(-t1\right) \cdot v}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\] |
|---|---|
*-commutative [=>]18.6 | \[ \frac{\color{blue}{v \cdot \left(-t1\right)}}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\] |
associate-/l* [=>]16.2 | \[ \color{blue}{\frac{v}{\frac{\left(t1 + u\right) \cdot \left(t1 + u\right)}{-t1}}}
\] |
associate-*r/ [<=]3.5 | \[ \frac{v}{\color{blue}{\left(t1 + u\right) \cdot \frac{t1 + u}{-t1}}}
\] |
associate-/r* [=>]0.7 | \[ \color{blue}{\frac{\frac{v}{t1 + u}}{\frac{t1 + u}{-t1}}}
\] |
neg-mul-1 [=>]0.7 | \[ \frac{\frac{v}{t1 + u}}{\frac{t1 + u}{\color{blue}{-1 \cdot t1}}}
\] |
associate-/l/ [<=]0.7 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\frac{\frac{t1 + u}{t1}}{-1}}}
\] |
metadata-eval [<=]0.7 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{0 - 1}}}
\] |
mul0-lft [<=]6.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{0 \cdot \frac{t1 + u}{t1}} - 1}}
\] |
associate-*r/ [=>]0.7 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{0 \cdot \left(t1 + u\right)}{t1}} - 1}}
\] |
mul0-lft [=>]0.7 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{0}}{t1} - 1}}
\] |
*-inverses [<=]0.7 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{0}{t1} - \color{blue}{\frac{t1}{t1}}}}
\] |
div-sub [<=]0.7 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{0 - t1}{t1}}}}
\] |
neg-sub0 [<=]0.7 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{-t1}}{t1}}}
\] |
neg-mul-1 [=>]0.7 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{-1 \cdot t1}}{t1}}}
\] |
*-commutative [=>]0.7 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{t1 \cdot -1}}{t1}}}
\] |
associate-/l* [=>]0.7 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{t1}{\frac{t1}{-1}}}}}
\] |
associate-/l* [<=]0.7 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\frac{\frac{t1 + u}{t1} \cdot \frac{t1}{-1}}{t1}}}
\] |
*-commutative [=>]0.7 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1}{-1} \cdot \frac{t1 + u}{t1}}}{t1}}
\] |
times-frac [<=]14.6 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1 \cdot \left(t1 + u\right)}{-1 \cdot t1}}}{t1}}
\] |
neg-mul-1 [<=]14.6 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 \cdot \left(t1 + u\right)}{\color{blue}{-t1}}}{t1}}
\] |
associate-/l* [=>]0.8 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1}{\frac{-t1}{t1 + u}}}}{t1}}
\] |
if -4.49999999999999993e-250 < t1 < 9.5000000000000004e-230Initial program 16.5
Simplified6.1
[Start]16.5 | \[ \frac{\left(-t1\right) \cdot v}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\] |
|---|---|
*-commutative [=>]16.5 | \[ \frac{\color{blue}{v \cdot \left(-t1\right)}}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\] |
associate-/l* [=>]15.8 | \[ \color{blue}{\frac{v}{\frac{\left(t1 + u\right) \cdot \left(t1 + u\right)}{-t1}}}
\] |
associate-*r/ [<=]3.7 | \[ \frac{v}{\color{blue}{\left(t1 + u\right) \cdot \frac{t1 + u}{-t1}}}
\] |
associate-/r* [=>]6.1 | \[ \color{blue}{\frac{\frac{v}{t1 + u}}{\frac{t1 + u}{-t1}}}
\] |
neg-mul-1 [=>]6.1 | \[ \frac{\frac{v}{t1 + u}}{\frac{t1 + u}{\color{blue}{-1 \cdot t1}}}
\] |
associate-/l/ [<=]6.1 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\frac{\frac{t1 + u}{t1}}{-1}}}
\] |
metadata-eval [<=]6.1 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{0 - 1}}}
\] |
mul0-lft [<=]35.0 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{0 \cdot \frac{t1 + u}{t1}} - 1}}
\] |
associate-*r/ [=>]6.1 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{0 \cdot \left(t1 + u\right)}{t1}} - 1}}
\] |
mul0-lft [=>]6.1 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{0}}{t1} - 1}}
\] |
*-inverses [<=]6.1 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{0}{t1} - \color{blue}{\frac{t1}{t1}}}}
\] |
div-sub [<=]6.1 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{0 - t1}{t1}}}}
\] |
neg-sub0 [<=]6.1 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{-t1}}{t1}}}
\] |
neg-mul-1 [=>]6.1 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{-1 \cdot t1}}{t1}}}
\] |
*-commutative [=>]6.1 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{t1 \cdot -1}}{t1}}}
\] |
associate-/l* [=>]6.1 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{t1}{\frac{t1}{-1}}}}}
\] |
associate-/l* [<=]6.1 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\frac{\frac{t1 + u}{t1} \cdot \frac{t1}{-1}}{t1}}}
\] |
*-commutative [=>]6.1 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1}{-1} \cdot \frac{t1 + u}{t1}}}{t1}}
\] |
times-frac [<=]27.2 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1 \cdot \left(t1 + u\right)}{-1 \cdot t1}}}{t1}}
\] |
neg-mul-1 [<=]27.2 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 \cdot \left(t1 + u\right)}{\color{blue}{-t1}}}{t1}}
\] |
associate-/l* [=>]6.1 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1}{\frac{-t1}{t1 + u}}}}{t1}}
\] |
Taylor expanded in t1 around 0 16.5
Simplified16.5
[Start]16.5 | \[ -1 \cdot \frac{t1 \cdot v}{{u}^{2}}
\] |
|---|---|
associate-*r/ [=>]16.5 | \[ \color{blue}{\frac{-1 \cdot \left(t1 \cdot v\right)}{{u}^{2}}}
\] |
associate-*r* [=>]16.5 | \[ \frac{\color{blue}{\left(-1 \cdot t1\right) \cdot v}}{{u}^{2}}
\] |
neg-mul-1 [<=]16.5 | \[ \frac{\color{blue}{\left(-t1\right)} \cdot v}{{u}^{2}}
\] |
unpow2 [=>]16.5 | \[ \frac{\left(-t1\right) \cdot v}{\color{blue}{u \cdot u}}
\] |
Taylor expanded in t1 around 0 16.5
Simplified12.5
[Start]16.5 | \[ -1 \cdot \frac{t1 \cdot v}{{u}^{2}}
\] |
|---|---|
unpow2 [=>]16.5 | \[ -1 \cdot \frac{t1 \cdot v}{\color{blue}{u \cdot u}}
\] |
associate-*r/ [=>]16.5 | \[ \color{blue}{\frac{-1 \cdot \left(t1 \cdot v\right)}{u \cdot u}}
\] |
associate-*l/ [<=]16.8 | \[ \color{blue}{\frac{-1}{u \cdot u} \cdot \left(t1 \cdot v\right)}
\] |
associate-*r* [=>]16.3 | \[ \color{blue}{\left(\frac{-1}{u \cdot u} \cdot t1\right) \cdot v}
\] |
associate-/r/ [<=]15.8 | \[ \color{blue}{\frac{-1}{\frac{u \cdot u}{t1}}} \cdot v
\] |
associate-*l/ [=>]15.8 | \[ \color{blue}{\frac{-1 \cdot v}{\frac{u \cdot u}{t1}}}
\] |
neg-mul-1 [<=]15.8 | \[ \frac{\color{blue}{-v}}{\frac{u \cdot u}{t1}}
\] |
associate-/r/ [=>]16.6 | \[ \color{blue}{\frac{-v}{u \cdot u} \cdot t1}
\] |
*-commutative [=>]16.6 | \[ \color{blue}{t1 \cdot \frac{-v}{u \cdot u}}
\] |
associate-/r* [=>]12.5 | \[ t1 \cdot \color{blue}{\frac{\frac{-v}{u}}{u}}
\] |
Applied egg-rr12.3
Applied egg-rr9.5
if 9.5000000000000004e-230 < t1 Initial program 18.6
Simplified0.8
[Start]18.6 | \[ \frac{\left(-t1\right) \cdot v}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\] |
|---|---|
times-frac [=>]0.8 | \[ \color{blue}{\frac{-t1}{t1 + u} \cdot \frac{v}{t1 + u}}
\] |
Final simplification1.6
| Alternative 1 | |
|---|---|
| Error | 1.3 |
| Cost | 2441 |
| Alternative 2 | |
|---|---|
| Error | 14.5 |
| Cost | 1040 |
| Alternative 3 | |
|---|---|
| Error | 1.7 |
| Cost | 969 |
| Alternative 4 | |
|---|---|
| Error | 14.8 |
| Cost | 841 |
| Alternative 5 | |
|---|---|
| Error | 14.7 |
| Cost | 777 |
| Alternative 6 | |
|---|---|
| Error | 13.8 |
| Cost | 777 |
| Alternative 7 | |
|---|---|
| Error | 14.7 |
| Cost | 777 |
| Alternative 8 | |
|---|---|
| Error | 20.8 |
| Cost | 713 |
| Alternative 9 | |
|---|---|
| Error | 21.2 |
| Cost | 712 |
| Alternative 10 | |
|---|---|
| Error | 20.6 |
| Cost | 712 |
| Alternative 11 | |
|---|---|
| Error | 27.8 |
| Cost | 520 |
| Alternative 12 | |
|---|---|
| Error | 50.5 |
| Cost | 456 |
| Alternative 13 | |
|---|---|
| Error | 54.8 |
| Cost | 192 |
herbie shell --seed 2023032
(FPCore (u v t1)
:name "Rosa's DopplerBench"
:precision binary64
(/ (* (- t1) v) (* (+ t1 u) (+ t1 u))))