
(FPCore (x y z a) :precision binary64 (+ x (- (tan (+ y z)) (tan a))))
double code(double x, double y, double z, double a) {
return x + (tan((y + z)) - tan(a));
}
real(8) function code(x, y, z, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: a
code = x + (tan((y + z)) - tan(a))
end function
public static double code(double x, double y, double z, double a) {
return x + (Math.tan((y + z)) - Math.tan(a));
}
def code(x, y, z, a): return x + (math.tan((y + z)) - math.tan(a))
function code(x, y, z, a) return Float64(x + Float64(tan(Float64(y + z)) - tan(a))) end
function tmp = code(x, y, z, a) tmp = x + (tan((y + z)) - tan(a)); end
code[x_, y_, z_, a_] := N[(x + N[(N[Tan[N[(y + z), $MachinePrecision]], $MachinePrecision] - N[Tan[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(\tan \left(y + z\right) - \tan a\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 8 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z a) :precision binary64 (+ x (- (tan (+ y z)) (tan a))))
double code(double x, double y, double z, double a) {
return x + (tan((y + z)) - tan(a));
}
real(8) function code(x, y, z, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: a
code = x + (tan((y + z)) - tan(a))
end function
public static double code(double x, double y, double z, double a) {
return x + (Math.tan((y + z)) - Math.tan(a));
}
def code(x, y, z, a): return x + (math.tan((y + z)) - math.tan(a))
function code(x, y, z, a) return Float64(x + Float64(tan(Float64(y + z)) - tan(a))) end
function tmp = code(x, y, z, a) tmp = x + (tan((y + z)) - tan(a)); end
code[x_, y_, z_, a_] := N[(x + N[(N[Tan[N[(y + z), $MachinePrecision]], $MachinePrecision] - N[Tan[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(\tan \left(y + z\right) - \tan a\right)
\end{array}
(FPCore (x y z a) :precision binary64 (+ (/ (fma (+ (tan y) (tan z)) (cos a) (* (- (* (tan y) (tan z)) 1.0) (sin a))) (* (fma (- (tan z)) (tan y) 1.0) (cos a))) x))
double code(double x, double y, double z, double a) {
return (fma((tan(y) + tan(z)), cos(a), (((tan(y) * tan(z)) - 1.0) * sin(a))) / (fma(-tan(z), tan(y), 1.0) * cos(a))) + x;
}
function code(x, y, z, a) return Float64(Float64(fma(Float64(tan(y) + tan(z)), cos(a), Float64(Float64(Float64(tan(y) * tan(z)) - 1.0) * sin(a))) / Float64(fma(Float64(-tan(z)), tan(y), 1.0) * cos(a))) + x) end
code[x_, y_, z_, a_] := N[(N[(N[(N[(N[Tan[y], $MachinePrecision] + N[Tan[z], $MachinePrecision]), $MachinePrecision] * N[Cos[a], $MachinePrecision] + N[(N[(N[(N[Tan[y], $MachinePrecision] * N[Tan[z], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision] * N[Sin[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[((-N[Tan[z], $MachinePrecision]) * N[Tan[y], $MachinePrecision] + 1.0), $MachinePrecision] * N[Cos[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}
\\
\frac{\mathsf{fma}\left(\tan y + \tan z, \cos a, \left(\tan y \cdot \tan z - 1\right) \cdot \sin a\right)}{\mathsf{fma}\left(-\tan z, \tan y, 1\right) \cdot \cos a} + x
\end{array}
Initial program 78.9%
lift--.f64N/A
sub-negN/A
lift-tan.f64N/A
lift-+.f64N/A
tan-sumN/A
lift-tan.f64N/A
tan-quotN/A
distribute-neg-fracN/A
frac-addN/A
lower-/.f64N/A
Applied rewrites99.8%
lift-fma.f64N/A
+-commutativeN/A
lift-neg.f64N/A
cancel-sign-sub-invN/A
lift-tan.f64N/A
lift-tan.f64N/A
lower--.f64N/A
lift-tan.f64N/A
lift-tan.f64N/A
*-commutativeN/A
lower-*.f6499.8
Applied rewrites99.8%
Final simplification99.8%
(FPCore (x y z a)
:precision binary64
(let* ((t_0 (+ (tan y) (tan z))) (t_1 (- x (- (tan a) (/ t_0 1.0)))))
(if (<= (tan a) -0.005)
t_1
(if (<= (tan a) 2e-9)
(- (/ t_0 (fma (- (tan z)) (tan y) 1.0)) (- x))
t_1))))
double code(double x, double y, double z, double a) {
double t_0 = tan(y) + tan(z);
double t_1 = x - (tan(a) - (t_0 / 1.0));
double tmp;
if (tan(a) <= -0.005) {
tmp = t_1;
} else if (tan(a) <= 2e-9) {
tmp = (t_0 / fma(-tan(z), tan(y), 1.0)) - -x;
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, a) t_0 = Float64(tan(y) + tan(z)) t_1 = Float64(x - Float64(tan(a) - Float64(t_0 / 1.0))) tmp = 0.0 if (tan(a) <= -0.005) tmp = t_1; elseif (tan(a) <= 2e-9) tmp = Float64(Float64(t_0 / fma(Float64(-tan(z)), tan(y), 1.0)) - Float64(-x)); else tmp = t_1; end return tmp end
code[x_, y_, z_, a_] := Block[{t$95$0 = N[(N[Tan[y], $MachinePrecision] + N[Tan[z], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(x - N[(N[Tan[a], $MachinePrecision] - N[(t$95$0 / 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[Tan[a], $MachinePrecision], -0.005], t$95$1, If[LessEqual[N[Tan[a], $MachinePrecision], 2e-9], N[(N[(t$95$0 / N[((-N[Tan[z], $MachinePrecision]) * N[Tan[y], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] - (-x)), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \tan y + \tan z\\
t_1 := x - \left(\tan a - \frac{t\_0}{1}\right)\\
\mathbf{if}\;\tan a \leq -0.005:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\tan a \leq 2 \cdot 10^{-9}:\\
\;\;\;\;\frac{t\_0}{\mathsf{fma}\left(-\tan z, \tan y, 1\right)} - \left(-x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (tan.f64 a) < -0.0050000000000000001 or 2.00000000000000012e-9 < (tan.f64 a) Initial program 82.0%
lift-tan.f64N/A
lift-+.f64N/A
tan-sumN/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-tan.f64N/A
lower-tan.f64N/A
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-fma.f64N/A
lower-neg.f64N/A
lower-tan.f64N/A
lower-tan.f6499.7
Applied rewrites99.7%
Taylor expanded in z around 0
Applied rewrites82.5%
if -0.0050000000000000001 < (tan.f64 a) < 2.00000000000000012e-9Initial program 75.1%
lift-+.f64N/A
+-commutativeN/A
lift--.f64N/A
associate-+l-N/A
lower--.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower--.f6475.1
Applied rewrites75.1%
Taylor expanded in a around 0
mul-1-negN/A
lower-neg.f6475.1
Applied rewrites75.1%
lift-tan.f64N/A
lift-+.f64N/A
tan-sumN/A
lift-tan.f64N/A
lift-tan.f64N/A
lift-+.f64N/A
lift-tan.f64N/A
lift-tan.f64N/A
cancel-sign-sub-invN/A
lift-neg.f64N/A
+-commutativeN/A
lift-fma.f64N/A
lower-/.f6499.1
lift-+.f64N/A
+-commutativeN/A
lower-+.f6499.1
Applied rewrites99.1%
Final simplification90.0%
(FPCore (x y z a) :precision binary64 (+ (- (/ (+ (tan y) (tan z)) (- 1.0 (* (tan y) (tan z)))) (tan a)) x))
double code(double x, double y, double z, double a) {
return (((tan(y) + tan(z)) / (1.0 - (tan(y) * tan(z)))) - tan(a)) + x;
}
real(8) function code(x, y, z, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: a
code = (((tan(y) + tan(z)) / (1.0d0 - (tan(y) * tan(z)))) - tan(a)) + x
end function
public static double code(double x, double y, double z, double a) {
return (((Math.tan(y) + Math.tan(z)) / (1.0 - (Math.tan(y) * Math.tan(z)))) - Math.tan(a)) + x;
}
def code(x, y, z, a): return (((math.tan(y) + math.tan(z)) / (1.0 - (math.tan(y) * math.tan(z)))) - math.tan(a)) + x
function code(x, y, z, a) return Float64(Float64(Float64(Float64(tan(y) + tan(z)) / Float64(1.0 - Float64(tan(y) * tan(z)))) - tan(a)) + x) end
function tmp = code(x, y, z, a) tmp = (((tan(y) + tan(z)) / (1.0 - (tan(y) * tan(z)))) - tan(a)) + x; end
code[x_, y_, z_, a_] := N[(N[(N[(N[(N[Tan[y], $MachinePrecision] + N[Tan[z], $MachinePrecision]), $MachinePrecision] / N[(1.0 - N[(N[Tan[y], $MachinePrecision] * N[Tan[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[Tan[a], $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}
\\
\left(\frac{\tan y + \tan z}{1 - \tan y \cdot \tan z} - \tan a\right) + x
\end{array}
Initial program 78.9%
lift-tan.f64N/A
lift-+.f64N/A
tan-sumN/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-tan.f64N/A
lower-tan.f64N/A
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-fma.f64N/A
lower-neg.f64N/A
lower-tan.f64N/A
lower-tan.f6499.8
Applied rewrites99.8%
lift-fma.f64N/A
+-commutativeN/A
lift-neg.f64N/A
cancel-sign-sub-invN/A
lift-tan.f64N/A
lift-tan.f64N/A
lower--.f64N/A
lift-tan.f64N/A
lift-tan.f64N/A
*-commutativeN/A
lower-*.f6499.8
Applied rewrites99.8%
Final simplification99.8%
(FPCore (x y z a) :precision binary64 (- x (- (tan a) (/ (+ (tan y) (tan z)) 1.0))))
double code(double x, double y, double z, double a) {
return x - (tan(a) - ((tan(y) + tan(z)) / 1.0));
}
real(8) function code(x, y, z, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: a
code = x - (tan(a) - ((tan(y) + tan(z)) / 1.0d0))
end function
public static double code(double x, double y, double z, double a) {
return x - (Math.tan(a) - ((Math.tan(y) + Math.tan(z)) / 1.0));
}
def code(x, y, z, a): return x - (math.tan(a) - ((math.tan(y) + math.tan(z)) / 1.0))
function code(x, y, z, a) return Float64(x - Float64(tan(a) - Float64(Float64(tan(y) + tan(z)) / 1.0))) end
function tmp = code(x, y, z, a) tmp = x - (tan(a) - ((tan(y) + tan(z)) / 1.0)); end
code[x_, y_, z_, a_] := N[(x - N[(N[Tan[a], $MachinePrecision] - N[(N[(N[Tan[y], $MachinePrecision] + N[Tan[z], $MachinePrecision]), $MachinePrecision] / 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x - \left(\tan a - \frac{\tan y + \tan z}{1}\right)
\end{array}
Initial program 78.9%
lift-tan.f64N/A
lift-+.f64N/A
tan-sumN/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-tan.f64N/A
lower-tan.f64N/A
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-fma.f64N/A
lower-neg.f64N/A
lower-tan.f64N/A
lower-tan.f6499.8
Applied rewrites99.8%
Taylor expanded in z around 0
Applied rewrites79.6%
Final simplification79.6%
(FPCore (x y z a) :precision binary64 (+ (- (tan (+ y z)) (/ (sin a) (cos a))) x))
double code(double x, double y, double z, double a) {
return (tan((y + z)) - (sin(a) / cos(a))) + x;
}
real(8) function code(x, y, z, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: a
code = (tan((y + z)) - (sin(a) / cos(a))) + x
end function
public static double code(double x, double y, double z, double a) {
return (Math.tan((y + z)) - (Math.sin(a) / Math.cos(a))) + x;
}
def code(x, y, z, a): return (math.tan((y + z)) - (math.sin(a) / math.cos(a))) + x
function code(x, y, z, a) return Float64(Float64(tan(Float64(y + z)) - Float64(sin(a) / cos(a))) + x) end
function tmp = code(x, y, z, a) tmp = (tan((y + z)) - (sin(a) / cos(a))) + x; end
code[x_, y_, z_, a_] := N[(N[(N[Tan[N[(y + z), $MachinePrecision]], $MachinePrecision] - N[(N[Sin[a], $MachinePrecision] / N[Cos[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}
\\
\left(\tan \left(y + z\right) - \frac{\sin a}{\cos a}\right) + x
\end{array}
Initial program 78.9%
lift-tan.f64N/A
tan-quotN/A
lower-/.f64N/A
lower-sin.f64N/A
lower-cos.f6478.9
Applied rewrites78.9%
Final simplification78.9%
(FPCore (x y z a) :precision binary64 (+ (- (tan (+ y z)) (tan a)) x))
double code(double x, double y, double z, double a) {
return (tan((y + z)) - tan(a)) + x;
}
real(8) function code(x, y, z, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: a
code = (tan((y + z)) - tan(a)) + x
end function
public static double code(double x, double y, double z, double a) {
return (Math.tan((y + z)) - Math.tan(a)) + x;
}
def code(x, y, z, a): return (math.tan((y + z)) - math.tan(a)) + x
function code(x, y, z, a) return Float64(Float64(tan(Float64(y + z)) - tan(a)) + x) end
function tmp = code(x, y, z, a) tmp = (tan((y + z)) - tan(a)) + x; end
code[x_, y_, z_, a_] := N[(N[(N[Tan[N[(y + z), $MachinePrecision]], $MachinePrecision] - N[Tan[a], $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}
\\
\left(\tan \left(y + z\right) - \tan a\right) + x
\end{array}
Initial program 78.9%
Final simplification78.9%
(FPCore (x y z a) :precision binary64 (- (tan (+ y z)) (- x)))
double code(double x, double y, double z, double a) {
return tan((y + z)) - -x;
}
real(8) function code(x, y, z, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: a
code = tan((y + z)) - -x
end function
public static double code(double x, double y, double z, double a) {
return Math.tan((y + z)) - -x;
}
def code(x, y, z, a): return math.tan((y + z)) - -x
function code(x, y, z, a) return Float64(tan(Float64(y + z)) - Float64(-x)) end
function tmp = code(x, y, z, a) tmp = tan((y + z)) - -x; end
code[x_, y_, z_, a_] := N[(N[Tan[N[(y + z), $MachinePrecision]], $MachinePrecision] - (-x)), $MachinePrecision]
\begin{array}{l}
\\
\tan \left(y + z\right) - \left(-x\right)
\end{array}
Initial program 78.9%
lift-+.f64N/A
+-commutativeN/A
lift--.f64N/A
associate-+l-N/A
lower--.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower--.f6478.8
Applied rewrites78.8%
Taylor expanded in a around 0
mul-1-negN/A
lower-neg.f6447.2
Applied rewrites47.2%
Final simplification47.2%
(FPCore (x y z a) :precision binary64 (* -1.0 (- x)))
double code(double x, double y, double z, double a) {
return -1.0 * -x;
}
real(8) function code(x, y, z, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: a
code = (-1.0d0) * -x
end function
public static double code(double x, double y, double z, double a) {
return -1.0 * -x;
}
def code(x, y, z, a): return -1.0 * -x
function code(x, y, z, a) return Float64(-1.0 * Float64(-x)) end
function tmp = code(x, y, z, a) tmp = -1.0 * -x; end
code[x_, y_, z_, a_] := N[(-1.0 * (-x)), $MachinePrecision]
\begin{array}{l}
\\
-1 \cdot \left(-x\right)
\end{array}
Initial program 78.9%
lift-+.f64N/A
+-commutativeN/A
lift--.f64N/A
associate-+l-N/A
lower--.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower--.f6478.8
Applied rewrites78.8%
Taylor expanded in x around -inf
associate-*r*N/A
mul-1-negN/A
lower-*.f64N/A
lower-neg.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
Applied rewrites78.7%
Taylor expanded in x around inf
Applied rewrites30.8%
Final simplification30.8%
herbie shell --seed 2024255
(FPCore (x y z a)
:name "tan-example (used to crash)"
:precision binary64
:pre (and (and (and (or (== x 0.0) (and (<= 0.5884142 x) (<= x 505.5909))) (or (and (<= -1.796658e+308 y) (<= y -9.425585e-310)) (and (<= 1.284938e-309 y) (<= y 1.751224e+308)))) (or (and (<= -1.776707e+308 z) (<= z -8.599796e-310)) (and (<= 3.293145e-311 z) (<= z 1.725154e+308)))) (or (and (<= -1.796658e+308 a) (<= a -9.425585e-310)) (and (<= 1.284938e-309 a) (<= a 1.751224e+308))))
(+ x (- (tan (+ y z)) (tan a))))