| Alternative 1 | |
|---|---|
| Error | 16.1 |
| Cost | 1305 |
(FPCore (u v t1) :precision binary64 (/ (* (- t1) v) (* (+ t1 u) (+ t1 u))))
(FPCore (u v t1) :precision binary64 (/ (/ v (+ 1.0 (/ u t1))) (- (- 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) {
return (v / (1.0 + (u / t1))) / (-u - t1);
}
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
code = (v / (1.0d0 + (u / t1))) / (-u - t1)
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) {
return (v / (1.0 + (u / t1))) / (-u - t1);
}
def code(u, v, t1): return (-t1 * v) / ((t1 + u) * (t1 + u))
def code(u, v, t1): return (v / (1.0 + (u / t1))) / (-u - t1)
function code(u, v, t1) return Float64(Float64(Float64(-t1) * v) / Float64(Float64(t1 + u) * Float64(t1 + u))) end
function code(u, v, t1) return Float64(Float64(v / Float64(1.0 + Float64(u / t1))) / Float64(Float64(-u) - t1)) end
function tmp = code(u, v, t1) tmp = (-t1 * v) / ((t1 + u) * (t1 + u)); end
function tmp = code(u, v, t1) tmp = (v / (1.0 + (u / t1))) / (-u - t1); 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_] := N[(N[(v / N[(1.0 + N[(u / t1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[((-u) - t1), $MachinePrecision]), $MachinePrecision]
\frac{\left(-t1\right) \cdot v}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\frac{\frac{v}{1 + \frac{u}{t1}}}{\left(-u\right) - t1}
Results
Initial program 18.3
Simplified1.5
[Start]18.3 | \[ \frac{\left(-t1\right) \cdot v}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\] |
|---|---|
*-commutative [=>]18.3 | \[ \frac{\color{blue}{v \cdot \left(-t1\right)}}{\left(t1 + u\right) \cdot \left(t1 + u\right)}
\] |
associate-/l* [=>]15.9 | \[ \color{blue}{\frac{v}{\frac{\left(t1 + u\right) \cdot \left(t1 + u\right)}{-t1}}}
\] |
associate-*r/ [<=]3.3 | \[ \frac{v}{\color{blue}{\left(t1 + u\right) \cdot \frac{t1 + u}{-t1}}}
\] |
associate-/r* [=>]1.5 | \[ \color{blue}{\frac{\frac{v}{t1 + u}}{\frac{t1 + u}{-t1}}}
\] |
neg-mul-1 [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{t1 + u}{\color{blue}{-1 \cdot t1}}}
\] |
associate-/l/ [<=]1.5 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\frac{\frac{t1 + u}{t1}}{-1}}}
\] |
metadata-eval [<=]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{0 - 1}}}
\] |
mul0-lft [<=]8.9 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{0 \cdot \frac{t1 + u}{t1}} - 1}}
\] |
associate-*r/ [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{0 \cdot \left(t1 + u\right)}{t1}} - 1}}
\] |
mul0-lft [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{0}}{t1} - 1}}
\] |
*-inverses [<=]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{0}{t1} - \color{blue}{\frac{t1}{t1}}}}
\] |
div-sub [<=]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{0 - t1}{t1}}}}
\] |
neg-sub0 [<=]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{-t1}}{t1}}}
\] |
neg-mul-1 [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{-1 \cdot t1}}{t1}}}
\] |
*-commutative [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\frac{\color{blue}{t1 \cdot -1}}{t1}}}
\] |
associate-/l* [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 + u}{t1}}{\color{blue}{\frac{t1}{\frac{t1}{-1}}}}}
\] |
associate-/l* [<=]1.5 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\frac{\frac{t1 + u}{t1} \cdot \frac{t1}{-1}}{t1}}}
\] |
*-commutative [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1}{-1} \cdot \frac{t1 + u}{t1}}}{t1}}
\] |
times-frac [<=]16.3 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1 \cdot \left(t1 + u\right)}{-1 \cdot t1}}}{t1}}
\] |
neg-mul-1 [<=]16.3 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1 \cdot \left(t1 + u\right)}{\color{blue}{-t1}}}{t1}}
\] |
associate-/l* [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1}{\frac{-t1}{t1 + u}}}}{t1}}
\] |
associate-/r/ [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{t1}{-t1} \cdot \left(t1 + u\right)}}{t1}}
\] |
neg-mul-1 [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1}{\color{blue}{-1 \cdot t1}} \cdot \left(t1 + u\right)}{t1}}
\] |
*-commutative [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{t1}{\color{blue}{t1 \cdot -1}} \cdot \left(t1 + u\right)}{t1}}
\] |
associate-/r* [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\frac{\frac{t1}{t1}}{-1}} \cdot \left(t1 + u\right)}{t1}}
\] |
*-inverses [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\frac{\color{blue}{1}}{-1} \cdot \left(t1 + u\right)}{t1}}
\] |
metadata-eval [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{-1} \cdot \left(t1 + u\right)}{t1}}
\] |
neg-mul-1 [<=]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{-\left(t1 + u\right)}}{t1}}
\] |
distribute-neg-in [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\left(-t1\right) + \left(-u\right)}}{t1}}
\] |
sub-neg [<=]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{\left(-t1\right) - u}}{t1}}
\] |
div-sub [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\frac{-t1}{t1} - \frac{u}{t1}}}
\] |
neg-sub0 [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\frac{\color{blue}{0 - t1}}{t1} - \frac{u}{t1}}
\] |
div-sub [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{\left(\frac{0}{t1} - \frac{t1}{t1}\right)} - \frac{u}{t1}}
\] |
mul0-lft [<=]1.5 | \[ \frac{\frac{v}{t1 + u}}{\left(\frac{\color{blue}{0 \cdot \left(t1 + u\right)}}{t1} - \frac{t1}{t1}\right) - \frac{u}{t1}}
\] |
associate-*r/ [<=]8.9 | \[ \frac{\frac{v}{t1 + u}}{\left(\color{blue}{0 \cdot \frac{t1 + u}{t1}} - \frac{t1}{t1}\right) - \frac{u}{t1}}
\] |
mul0-lft [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\left(\color{blue}{0} - \frac{t1}{t1}\right) - \frac{u}{t1}}
\] |
*-inverses [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\left(0 - \color{blue}{1}\right) - \frac{u}{t1}}
\] |
metadata-eval [=>]1.5 | \[ \frac{\frac{v}{t1 + u}}{\color{blue}{-1} - \frac{u}{t1}}
\] |
Applied egg-rr3.2
Simplified1.4
[Start]3.2 | \[ \frac{\frac{-1}{\frac{u}{t1} + 1}}{\left(-u\right) - t1} \cdot \left(-v\right)
\] |
|---|---|
associate-*l/ [=>]1.4 | \[ \color{blue}{\frac{\frac{-1}{\frac{u}{t1} + 1} \cdot \left(-v\right)}{\left(-u\right) - t1}}
\] |
*-commutative [<=]1.4 | \[ \frac{\color{blue}{\left(-v\right) \cdot \frac{-1}{\frac{u}{t1} + 1}}}{\left(-u\right) - t1}
\] |
distribute-lft-neg-out [=>]1.4 | \[ \frac{\color{blue}{-v \cdot \frac{-1}{\frac{u}{t1} + 1}}}{\left(-u\right) - t1}
\] |
/-rgt-identity [<=]1.4 | \[ \frac{-\color{blue}{\frac{v}{1}} \cdot \frac{-1}{\frac{u}{t1} + 1}}{\left(-u\right) - t1}
\] |
associate-/r/ [<=]1.4 | \[ \frac{-\color{blue}{\frac{v}{\frac{1}{\frac{-1}{\frac{u}{t1} + 1}}}}}{\left(-u\right) - t1}
\] |
associate-/r/ [=>]1.4 | \[ \frac{-\frac{v}{\color{blue}{\frac{1}{-1} \cdot \left(\frac{u}{t1} + 1\right)}}}{\left(-u\right) - t1}
\] |
metadata-eval [=>]1.4 | \[ \frac{-\frac{v}{\color{blue}{-1} \cdot \left(\frac{u}{t1} + 1\right)}}{\left(-u\right) - t1}
\] |
mul-1-neg [=>]1.4 | \[ \frac{-\frac{v}{\color{blue}{-\left(\frac{u}{t1} + 1\right)}}}{\left(-u\right) - t1}
\] |
+-commutative [<=]1.4 | \[ \frac{-\frac{v}{-\color{blue}{\left(1 + \frac{u}{t1}\right)}}}{\left(-u\right) - t1}
\] |
distribute-neg-in [=>]1.4 | \[ \frac{-\frac{v}{\color{blue}{\left(-1\right) + \left(-\frac{u}{t1}\right)}}}{\left(-u\right) - t1}
\] |
metadata-eval [=>]1.4 | \[ \frac{-\frac{v}{\color{blue}{-1} + \left(-\frac{u}{t1}\right)}}{\left(-u\right) - t1}
\] |
sub-neg [<=]1.4 | \[ \frac{-\frac{v}{\color{blue}{-1 - \frac{u}{t1}}}}{\left(-u\right) - t1}
\] |
distribute-neg-frac [=>]1.4 | \[ \frac{\color{blue}{\frac{-v}{-1 - \frac{u}{t1}}}}{\left(-u\right) - t1}
\] |
Taylor expanded in v around 0 1.4
Final simplification1.4
| Alternative 1 | |
|---|---|
| Error | 16.1 |
| Cost | 1305 |
| Alternative 2 | |
|---|---|
| Error | 15.5 |
| Cost | 1304 |
| Alternative 3 | |
|---|---|
| Error | 15.5 |
| Cost | 1304 |
| Alternative 4 | |
|---|---|
| Error | 15.5 |
| Cost | 1040 |
| Alternative 5 | |
|---|---|
| Error | 15.3 |
| Cost | 777 |
| Alternative 6 | |
|---|---|
| Error | 13.3 |
| Cost | 777 |
| Alternative 7 | |
|---|---|
| Error | 1.4 |
| Cost | 768 |
| Alternative 8 | |
|---|---|
| Error | 22.1 |
| Cost | 713 |
| Alternative 9 | |
|---|---|
| Error | 22.1 |
| Cost | 713 |
| Alternative 10 | |
|---|---|
| Error | 21.9 |
| Cost | 713 |
| Alternative 11 | |
|---|---|
| Error | 1.5 |
| Cost | 704 |
| Alternative 12 | |
|---|---|
| Error | 27.8 |
| Cost | 585 |
| Alternative 13 | |
|---|---|
| Error | 27.9 |
| Cost | 521 |
| Alternative 14 | |
|---|---|
| Error | 27.8 |
| Cost | 520 |
| Alternative 15 | |
|---|---|
| Error | 25.0 |
| Cost | 384 |
| Alternative 16 | |
|---|---|
| Error | 53.5 |
| Cost | 192 |
herbie shell --seed 2022354
(FPCore (u v t1)
:name "Rosa's DopplerBench"
:precision binary64
(/ (* (- t1) v) (* (+ t1 u) (+ t1 u))))