| Alternative 1 | |
|---|---|
| Error | 9.0 |
| Cost | 4172 |
(FPCore (x y z t a) :precision binary64 (+ x (* (- y z) (/ (- t x) (- a z)))))
(FPCore (x y z t a)
:precision binary64
(let* ((t_1
(+
(* (- y z) (- (/ t (- z a))))
(* x (+ (/ y (- z a)) (- 1.0 (/ z (- z a)))))))
(t_2 (+ x (* (- y z) (/ (- t x) (- a z))))))
(if (<= t_2 -4e-188)
t_1
(if (<= t_2 0.0) (+ (* (+ -1.0 (/ y z)) (- t)) (* y (/ x z))) t_1))))double code(double x, double y, double z, double t, double a) {
return x + ((y - z) * ((t - x) / (a - z)));
}
double code(double x, double y, double z, double t, double a) {
double t_1 = ((y - z) * -(t / (z - a))) + (x * ((y / (z - a)) + (1.0 - (z / (z - a)))));
double t_2 = x + ((y - z) * ((t - x) / (a - z)));
double tmp;
if (t_2 <= -4e-188) {
tmp = t_1;
} else if (t_2 <= 0.0) {
tmp = ((-1.0 + (y / z)) * -t) + (y * (x / z));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = x + ((y - z) * ((t - x) / (a - z)))
end function
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = ((y - z) * -(t / (z - a))) + (x * ((y / (z - a)) + (1.0d0 - (z / (z - a)))))
t_2 = x + ((y - z) * ((t - x) / (a - z)))
if (t_2 <= (-4d-188)) then
tmp = t_1
else if (t_2 <= 0.0d0) then
tmp = (((-1.0d0) + (y / z)) * -t) + (y * (x / z))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
return x + ((y - z) * ((t - x) / (a - z)));
}
public static double code(double x, double y, double z, double t, double a) {
double t_1 = ((y - z) * -(t / (z - a))) + (x * ((y / (z - a)) + (1.0 - (z / (z - a)))));
double t_2 = x + ((y - z) * ((t - x) / (a - z)));
double tmp;
if (t_2 <= -4e-188) {
tmp = t_1;
} else if (t_2 <= 0.0) {
tmp = ((-1.0 + (y / z)) * -t) + (y * (x / z));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): return x + ((y - z) * ((t - x) / (a - z)))
def code(x, y, z, t, a): t_1 = ((y - z) * -(t / (z - a))) + (x * ((y / (z - a)) + (1.0 - (z / (z - a))))) t_2 = x + ((y - z) * ((t - x) / (a - z))) tmp = 0 if t_2 <= -4e-188: tmp = t_1 elif t_2 <= 0.0: tmp = ((-1.0 + (y / z)) * -t) + (y * (x / z)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) return Float64(x + Float64(Float64(y - z) * Float64(Float64(t - x) / Float64(a - z)))) end
function code(x, y, z, t, a) t_1 = Float64(Float64(Float64(y - z) * Float64(-Float64(t / Float64(z - a)))) + Float64(x * Float64(Float64(y / Float64(z - a)) + Float64(1.0 - Float64(z / Float64(z - a)))))) t_2 = Float64(x + Float64(Float64(y - z) * Float64(Float64(t - x) / Float64(a - z)))) tmp = 0.0 if (t_2 <= -4e-188) tmp = t_1; elseif (t_2 <= 0.0) tmp = Float64(Float64(Float64(-1.0 + Float64(y / z)) * Float64(-t)) + Float64(y * Float64(x / z))); else tmp = t_1; end return tmp end
function tmp = code(x, y, z, t, a) tmp = x + ((y - z) * ((t - x) / (a - z))); end
function tmp_2 = code(x, y, z, t, a) t_1 = ((y - z) * -(t / (z - a))) + (x * ((y / (z - a)) + (1.0 - (z / (z - a))))); t_2 = x + ((y - z) * ((t - x) / (a - z))); tmp = 0.0; if (t_2 <= -4e-188) tmp = t_1; elseif (t_2 <= 0.0) tmp = ((-1.0 + (y / z)) * -t) + (y * (x / z)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(N[(y - z), $MachinePrecision] * (-N[(t / N[(z - a), $MachinePrecision]), $MachinePrecision])), $MachinePrecision] + N[(x * N[(N[(y / N[(z - a), $MachinePrecision]), $MachinePrecision] + N[(1.0 - N[(z / N[(z - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -4e-188], t$95$1, If[LessEqual[t$95$2, 0.0], N[(N[(N[(-1.0 + N[(y / z), $MachinePrecision]), $MachinePrecision] * (-t)), $MachinePrecision] + N[(y * N[(x / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
x + \left(y - z\right) \cdot \frac{t - x}{a - z}
\begin{array}{l}
t_1 := \left(y - z\right) \cdot \left(-\frac{t}{z - a}\right) + x \cdot \left(\frac{y}{z - a} + \left(1 - \frac{z}{z - a}\right)\right)\\
t_2 := x + \left(y - z\right) \cdot \frac{t - x}{a - z}\\
\mathbf{if}\;t_2 \leq -4 \cdot 10^{-188}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t_2 \leq 0:\\
\;\;\;\;\left(-1 + \frac{y}{z}\right) \cdot \left(-t\right) + y \cdot \frac{x}{z}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
Results
if (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < -3.9999999999999998e-188 or 0.0 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) Initial program 6.9
Simplified6.9
[Start]6.9 | \[ x + \left(y - z\right) \cdot \frac{t - x}{a - z}
\] |
|---|---|
rational.json-simplify-50 [=>]6.9 | \[ x + \left(y - z\right) \cdot \color{blue}{\frac{-\left(t - x\right)}{z - a}}
\] |
rational.json-simplify-12 [=>]6.9 | \[ x + \left(y - z\right) \cdot \frac{\color{blue}{0 - \left(t - x\right)}}{z - a}
\] |
rational.json-simplify-45 [=>]6.9 | \[ x + \left(y - z\right) \cdot \frac{\color{blue}{x - \left(t - 0\right)}}{z - a}
\] |
rational.json-simplify-5 [=>]6.9 | \[ x + \left(y - z\right) \cdot \frac{x - \color{blue}{t}}{z - a}
\] |
Taylor expanded in x around 0 14.2
Simplified5.1
[Start]14.2 | \[ -1 \cdot \frac{t \cdot \left(y - z\right)}{z - a} + \left(\left(\frac{y}{z - a} + 1\right) - \frac{z}{z - a}\right) \cdot x
\] |
|---|---|
rational.json-simplify-49 [=>]5.9 | \[ -1 \cdot \color{blue}{\left(\left(y - z\right) \cdot \frac{t}{z - a}\right)} + \left(\left(\frac{y}{z - a} + 1\right) - \frac{z}{z - a}\right) \cdot x
\] |
rational.json-simplify-43 [=>]5.9 | \[ \color{blue}{\left(y - z\right) \cdot \left(\frac{t}{z - a} \cdot -1\right)} + \left(\left(\frac{y}{z - a} + 1\right) - \frac{z}{z - a}\right) \cdot x
\] |
rational.json-simplify-9 [=>]5.9 | \[ \left(y - z\right) \cdot \color{blue}{\left(-\frac{t}{z - a}\right)} + \left(\left(\frac{y}{z - a} + 1\right) - \frac{z}{z - a}\right) \cdot x
\] |
rational.json-simplify-2 [=>]5.9 | \[ \left(y - z\right) \cdot \left(-\frac{t}{z - a}\right) + \color{blue}{x \cdot \left(\left(\frac{y}{z - a} + 1\right) - \frac{z}{z - a}\right)}
\] |
rational.json-simplify-1 [=>]5.9 | \[ \left(y - z\right) \cdot \left(-\frac{t}{z - a}\right) + x \cdot \left(\color{blue}{\left(1 + \frac{y}{z - a}\right)} - \frac{z}{z - a}\right)
\] |
rational.json-simplify-48 [=>]5.1 | \[ \left(y - z\right) \cdot \left(-\frac{t}{z - a}\right) + x \cdot \color{blue}{\left(\frac{y}{z - a} + \left(1 - \frac{z}{z - a}\right)\right)}
\] |
if -3.9999999999999998e-188 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 0.0Initial program 55.7
Simplified55.8
[Start]55.7 | \[ x + \left(y - z\right) \cdot \frac{t - x}{a - z}
\] |
|---|---|
rational.json-simplify-50 [=>]55.7 | \[ x + \left(y - z\right) \cdot \color{blue}{\frac{-\left(t - x\right)}{z - a}}
\] |
rational.json-simplify-8 [=>]55.7 | \[ x + \left(y - z\right) \cdot \frac{\color{blue}{\left(t - x\right) \cdot -1}}{z - a}
\] |
rational.json-simplify-49 [=>]55.7 | \[ x + \left(y - z\right) \cdot \color{blue}{\left(-1 \cdot \frac{t - x}{z - a}\right)}
\] |
rational.json-simplify-43 [<=]55.7 | \[ x + \color{blue}{\frac{t - x}{z - a} \cdot \left(\left(y - z\right) \cdot -1\right)}
\] |
rational.json-simplify-5 [<=]55.7 | \[ x + \frac{t - x}{\color{blue}{\left(z - a\right) - 0}} \cdot \left(\left(y - z\right) \cdot -1\right)
\] |
rational.json-simplify-50 [=>]55.7 | \[ x + \color{blue}{\frac{-\left(t - x\right)}{0 - \left(z - a\right)}} \cdot \left(\left(y - z\right) \cdot -1\right)
\] |
rational.json-simplify-12 [<=]55.7 | \[ x + \frac{-\left(t - x\right)}{\color{blue}{-\left(z - a\right)}} \cdot \left(\left(y - z\right) \cdot -1\right)
\] |
rational.json-simplify-10 [=>]55.7 | \[ x + \frac{-\left(t - x\right)}{\color{blue}{\frac{z - a}{-1}}} \cdot \left(\left(y - z\right) \cdot -1\right)
\] |
rational.json-simplify-61 [=>]55.7 | \[ x + \color{blue}{\frac{-1}{\frac{z - a}{-\left(t - x\right)}}} \cdot \left(\left(y - z\right) \cdot -1\right)
\] |
rational.json-simplify-9 [=>]55.7 | \[ x + \frac{-1}{\frac{z - a}{-\left(t - x\right)}} \cdot \color{blue}{\left(-\left(y - z\right)\right)}
\] |
rational.json-simplify-10 [=>]55.7 | \[ x + \frac{-1}{\frac{z - a}{-\left(t - x\right)}} \cdot \color{blue}{\frac{y - z}{-1}}
\] |
rational.json-simplify-55 [=>]55.8 | \[ x + \color{blue}{\frac{\frac{y - z}{-1}}{\frac{\frac{z - a}{-\left(t - x\right)}}{-1}}}
\] |
rational.json-simplify-10 [<=]55.8 | \[ x + \frac{\color{blue}{-\left(y - z\right)}}{\frac{\frac{z - a}{-\left(t - x\right)}}{-1}}
\] |
rational.json-simplify-12 [=>]55.8 | \[ x + \frac{\color{blue}{0 - \left(y - z\right)}}{\frac{\frac{z - a}{-\left(t - x\right)}}{-1}}
\] |
rational.json-simplify-45 [=>]55.8 | \[ x + \frac{\color{blue}{z - \left(y - 0\right)}}{\frac{\frac{z - a}{-\left(t - x\right)}}{-1}}
\] |
rational.json-simplify-5 [=>]55.8 | \[ x + \frac{z - \color{blue}{y}}{\frac{\frac{z - a}{-\left(t - x\right)}}{-1}}
\] |
rational.json-simplify-46 [<=]55.8 | \[ x + \frac{z - y}{\color{blue}{\frac{z - a}{\left(-\left(t - x\right)\right) \cdot -1}}}
\] |
rational.json-simplify-2 [=>]55.8 | \[ x + \frac{z - y}{\frac{z - a}{\color{blue}{-1 \cdot \left(-\left(t - x\right)\right)}}}
\] |
rational.json-simplify-46 [=>]55.8 | \[ x + \frac{z - y}{\color{blue}{\frac{\frac{z - a}{-1}}{-\left(t - x\right)}}}
\] |
rational.json-simplify-10 [<=]55.8 | \[ x + \frac{z - y}{\frac{\color{blue}{-\left(z - a\right)}}{-\left(t - x\right)}}
\] |
rational.json-simplify-12 [=>]55.8 | \[ x + \frac{z - y}{\frac{\color{blue}{0 - \left(z - a\right)}}{-\left(t - x\right)}}
\] |
rational.json-simplify-45 [=>]55.8 | \[ x + \frac{z - y}{\frac{\color{blue}{a - \left(z - 0\right)}}{-\left(t - x\right)}}
\] |
rational.json-simplify-5 [=>]55.8 | \[ x + \frac{z - y}{\frac{a - \color{blue}{z}}{-\left(t - x\right)}}
\] |
rational.json-simplify-12 [=>]55.8 | \[ x + \frac{z - y}{\frac{a - z}{\color{blue}{0 - \left(t - x\right)}}}
\] |
Taylor expanded in z around inf 45.5
Simplified42.4
[Start]45.5 | \[ \left(-1 \cdot \left(x - t\right) + \left(\frac{y \cdot \left(x - t\right)}{z} + x\right)\right) - \frac{a \cdot \left(x - t\right)}{z}
\] |
|---|---|
rational.json-simplify-1 [=>]45.5 | \[ \left(-1 \cdot \left(x - t\right) + \color{blue}{\left(x + \frac{y \cdot \left(x - t\right)}{z}\right)}\right) - \frac{a \cdot \left(x - t\right)}{z}
\] |
rational.json-simplify-41 [=>]45.5 | \[ \color{blue}{\left(x + \left(\frac{y \cdot \left(x - t\right)}{z} + -1 \cdot \left(x - t\right)\right)\right)} - \frac{a \cdot \left(x - t\right)}{z}
\] |
rational.json-simplify-49 [=>]45.3 | \[ \left(x + \left(\color{blue}{\left(x - t\right) \cdot \frac{y}{z}} + -1 \cdot \left(x - t\right)\right)\right) - \frac{a \cdot \left(x - t\right)}{z}
\] |
rational.json-simplify-51 [=>]45.4 | \[ \left(x + \color{blue}{\left(x - t\right) \cdot \left(-1 + \frac{y}{z}\right)}\right) - \frac{a \cdot \left(x - t\right)}{z}
\] |
rational.json-simplify-49 [=>]42.4 | \[ \left(x + \left(x - t\right) \cdot \left(-1 + \frac{y}{z}\right)\right) - \color{blue}{\left(x - t\right) \cdot \frac{a}{z}}
\] |
Taylor expanded in a around 0 57.0
Taylor expanded in x around 0 26.2
Simplified22.6
[Start]26.2 | \[ \frac{y \cdot x}{z} + -1 \cdot \left(\left(\frac{y}{z} - 1\right) \cdot t\right)
\] |
|---|---|
rational.json-simplify-1 [=>]26.2 | \[ \color{blue}{-1 \cdot \left(\left(\frac{y}{z} - 1\right) \cdot t\right) + \frac{y \cdot x}{z}}
\] |
rational.json-simplify-15 [<=]26.2 | \[ -1 \cdot \left(\color{blue}{\left(\frac{y}{z} + -1\right)} \cdot t\right) + \frac{y \cdot x}{z}
\] |
rational.json-simplify-43 [=>]26.2 | \[ \color{blue}{\left(\frac{y}{z} + -1\right) \cdot \left(t \cdot -1\right)} + \frac{y \cdot x}{z}
\] |
rational.json-simplify-1 [=>]26.2 | \[ \color{blue}{\left(-1 + \frac{y}{z}\right)} \cdot \left(t \cdot -1\right) + \frac{y \cdot x}{z}
\] |
rational.json-simplify-9 [=>]26.2 | \[ \left(-1 + \frac{y}{z}\right) \cdot \color{blue}{\left(-t\right)} + \frac{y \cdot x}{z}
\] |
rational.json-simplify-2 [=>]26.2 | \[ \left(-1 + \frac{y}{z}\right) \cdot \left(-t\right) + \frac{\color{blue}{x \cdot y}}{z}
\] |
rational.json-simplify-49 [=>]22.6 | \[ \left(-1 + \frac{y}{z}\right) \cdot \left(-t\right) + \color{blue}{y \cdot \frac{x}{z}}
\] |
Final simplification8.0
| Alternative 1 | |
|---|---|
| Error | 9.0 |
| Cost | 4172 |
| Alternative 2 | |
|---|---|
| Error | 9.9 |
| Cost | 3660 |
| Alternative 3 | |
|---|---|
| Error | 10.7 |
| Cost | 3532 |
| Alternative 4 | |
|---|---|
| Error | 10.0 |
| Cost | 3532 |
| Alternative 5 | |
|---|---|
| Error | 10.0 |
| Cost | 3532 |
| Alternative 6 | |
|---|---|
| Error | 19.3 |
| Cost | 1364 |
| Alternative 7 | |
|---|---|
| Error | 25.8 |
| Cost | 1300 |
| Alternative 8 | |
|---|---|
| Error | 36.4 |
| Cost | 1240 |
| Alternative 9 | |
|---|---|
| Error | 30.8 |
| Cost | 1240 |
| Alternative 10 | |
|---|---|
| Error | 30.8 |
| Cost | 1240 |
| Alternative 11 | |
|---|---|
| Error | 30.7 |
| Cost | 1240 |
| Alternative 12 | |
|---|---|
| Error | 29.5 |
| Cost | 1240 |
| Alternative 13 | |
|---|---|
| Error | 29.2 |
| Cost | 1240 |
| Alternative 14 | |
|---|---|
| Error | 25.9 |
| Cost | 1236 |
| Alternative 15 | |
|---|---|
| Error | 25.9 |
| Cost | 1236 |
| Alternative 16 | |
|---|---|
| Error | 25.8 |
| Cost | 1236 |
| Alternative 17 | |
|---|---|
| Error | 18.5 |
| Cost | 1100 |
| Alternative 18 | |
|---|---|
| Error | 26.5 |
| Cost | 972 |
| Alternative 19 | |
|---|---|
| Error | 26.8 |
| Cost | 840 |
| Alternative 20 | |
|---|---|
| Error | 33.9 |
| Cost | 712 |
| Alternative 21 | |
|---|---|
| Error | 28.7 |
| Cost | 712 |
| Alternative 22 | |
|---|---|
| Error | 28.7 |
| Cost | 712 |
| Alternative 23 | |
|---|---|
| Error | 36.5 |
| Cost | 328 |
| Alternative 24 | |
|---|---|
| Error | 45.8 |
| Cost | 64 |
herbie shell --seed 2023074
(FPCore (x y z t a)
:name "Numeric.Signal:interpolate from hsignal-0.2.7.1"
:precision binary64
(+ x (* (- y z) (/ (- t x) (- a z)))))