
(FPCore (x tau) :precision binary32 (let* ((t_1 (* (* x PI) tau))) (* (/ (sin t_1) t_1) (/ (sin (* x PI)) (* x PI)))))
float code(float x, float tau) {
float t_1 = (x * ((float) M_PI)) * tau;
return (sinf(t_1) / t_1) * (sinf((x * ((float) M_PI))) / (x * ((float) M_PI)));
}
function code(x, tau) t_1 = Float32(Float32(x * Float32(pi)) * tau) return Float32(Float32(sin(t_1) / t_1) * Float32(sin(Float32(x * Float32(pi))) / Float32(x * Float32(pi)))) end
function tmp = code(x, tau) t_1 = (x * single(pi)) * tau; tmp = (sin(t_1) / t_1) * (sin((x * single(pi))) / (x * single(pi))); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(x \cdot \pi\right) \cdot tau\\
\frac{\sin t\_1}{t\_1} \cdot \frac{\sin \left(x \cdot \pi\right)}{x \cdot \pi}
\end{array}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 24 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x tau) :precision binary32 (let* ((t_1 (* (* x PI) tau))) (* (/ (sin t_1) t_1) (/ (sin (* x PI)) (* x PI)))))
float code(float x, float tau) {
float t_1 = (x * ((float) M_PI)) * tau;
return (sinf(t_1) / t_1) * (sinf((x * ((float) M_PI))) / (x * ((float) M_PI)));
}
function code(x, tau) t_1 = Float32(Float32(x * Float32(pi)) * tau) return Float32(Float32(sin(t_1) / t_1) * Float32(sin(Float32(x * Float32(pi))) / Float32(x * Float32(pi)))) end
function tmp = code(x, tau) t_1 = (x * single(pi)) * tau; tmp = (sin(t_1) / t_1) * (sin((x * single(pi))) / (x * single(pi))); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(x \cdot \pi\right) \cdot tau\\
\frac{\sin t\_1}{t\_1} \cdot \frac{\sin \left(x \cdot \pi\right)}{x \cdot \pi}
\end{array}
\end{array}
(FPCore (x tau) :precision binary32 (let* ((t_1 (* x (* tau PI)))) (* (/ (sin t_1) t_1) (/ (sin (* x PI)) (* x PI)))))
float code(float x, float tau) {
float t_1 = x * (tau * ((float) M_PI));
return (sinf(t_1) / t_1) * (sinf((x * ((float) M_PI))) / (x * ((float) M_PI)));
}
function code(x, tau) t_1 = Float32(x * Float32(tau * Float32(pi))) return Float32(Float32(sin(t_1) / t_1) * Float32(sin(Float32(x * Float32(pi))) / Float32(x * Float32(pi)))) end
function tmp = code(x, tau) t_1 = x * (tau * single(pi)); tmp = (sin(t_1) / t_1) * (sin((x * single(pi))) / (x * single(pi))); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(tau \cdot \pi\right)\\
\frac{\sin t\_1}{t\_1} \cdot \frac{\sin \left(x \cdot \pi\right)}{x \cdot \pi}
\end{array}
\end{array}
Initial program 97.7%
Taylor expanded in x around inf
lower-/.f32N/A
lower-sin.f32N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3297.8
Applied rewrites97.8%
(FPCore (x tau) :precision binary32 (let* ((t_1 (* PI (* x tau)))) (/ (* (sin (* x PI)) (sin t_1)) (* t_1 (* x PI)))))
float code(float x, float tau) {
float t_1 = ((float) M_PI) * (x * tau);
return (sinf((x * ((float) M_PI))) * sinf(t_1)) / (t_1 * (x * ((float) M_PI)));
}
function code(x, tau) t_1 = Float32(Float32(pi) * Float32(x * tau)) return Float32(Float32(sin(Float32(x * Float32(pi))) * sin(t_1)) / Float32(t_1 * Float32(x * Float32(pi)))) end
function tmp = code(x, tau) t_1 = single(pi) * (x * tau); tmp = (sin((x * single(pi))) * sin(t_1)) / (t_1 * (x * single(pi))); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \pi \cdot \left(x \cdot tau\right)\\
\frac{\sin \left(x \cdot \pi\right) \cdot \sin t\_1}{t\_1 \cdot \left(x \cdot \pi\right)}
\end{array}
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
associate-*l/N/A
lift-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
associate-*l*N/A
associate-/r*N/A
*-commutativeN/A
lower-/.f32N/A
Applied rewrites97.6%
lift-/.f32N/A
lift-/.f32N/A
associate-/r*N/A
lift-*.f32N/A
lift-/.f32N/A
associate-*r/N/A
lift-*.f32N/A
associate-/r*N/A
Applied rewrites97.6%
Final simplification97.6%
(FPCore (x tau) :precision binary32 (/ (* (sin (* x PI)) (sin (* PI (* x tau)))) (* PI (* (* x tau) (* x PI)))))
float code(float x, float tau) {
return (sinf((x * ((float) M_PI))) * sinf((((float) M_PI) * (x * tau)))) / (((float) M_PI) * ((x * tau) * (x * ((float) M_PI))));
}
function code(x, tau) return Float32(Float32(sin(Float32(x * Float32(pi))) * sin(Float32(Float32(pi) * Float32(x * tau)))) / Float32(Float32(pi) * Float32(Float32(x * tau) * Float32(x * Float32(pi))))) end
function tmp = code(x, tau) tmp = (sin((x * single(pi))) * sin((single(pi) * (x * tau)))) / (single(pi) * ((x * tau) * (x * single(pi)))); end
\begin{array}{l}
\\
\frac{\sin \left(x \cdot \pi\right) \cdot \sin \left(\pi \cdot \left(x \cdot tau\right)\right)}{\pi \cdot \left(\left(x \cdot tau\right) \cdot \left(x \cdot \pi\right)\right)}
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
associate-*l/N/A
lift-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
associate-*l*N/A
associate-/r*N/A
*-commutativeN/A
lower-/.f32N/A
Applied rewrites97.6%
lift-/.f32N/A
lift-/.f32N/A
associate-/r*N/A
lift-*.f32N/A
lift-/.f32N/A
associate-*r/N/A
lift-*.f32N/A
associate-/r*N/A
Applied rewrites97.6%
lift-*.f32N/A
*-commutativeN/A
lift-neg.f32N/A
lift-*.f32N/A
distribute-lft-neg-inN/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
associate-*l*N/A
lower-*.f32N/A
lower-neg.f32N/A
lower-*.f3297.4
lift-*.f32N/A
*-commutativeN/A
lift-*.f3297.4
Applied rewrites97.4%
Final simplification97.4%
(FPCore (x tau) :precision binary32 (let* ((t_1 (* PI (* x tau)))) (* (sin (* x PI)) (/ (sin t_1) (* x (* PI t_1))))))
float code(float x, float tau) {
float t_1 = ((float) M_PI) * (x * tau);
return sinf((x * ((float) M_PI))) * (sinf(t_1) / (x * (((float) M_PI) * t_1)));
}
function code(x, tau) t_1 = Float32(Float32(pi) * Float32(x * tau)) return Float32(sin(Float32(x * Float32(pi))) * Float32(sin(t_1) / Float32(x * Float32(Float32(pi) * t_1)))) end
function tmp = code(x, tau) t_1 = single(pi) * (x * tau); tmp = sin((x * single(pi))) * (sin(t_1) / (x * (single(pi) * t_1))); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \pi \cdot \left(x \cdot tau\right)\\
\sin \left(x \cdot \pi\right) \cdot \frac{\sin t\_1}{x \cdot \left(\pi \cdot t\_1\right)}
\end{array}
\end{array}
Initial program 97.7%
lift-*.f32N/A
*-commutativeN/A
lift-/.f32N/A
div-invN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
Applied rewrites97.3%
Final simplification97.3%
(FPCore (x tau) :precision binary32 (let* ((t_1 (* tau (* x PI)))) (* (/ (sin (* x PI)) (* PI (* x t_1))) (sin t_1))))
float code(float x, float tau) {
float t_1 = tau * (x * ((float) M_PI));
return (sinf((x * ((float) M_PI))) / (((float) M_PI) * (x * t_1))) * sinf(t_1);
}
function code(x, tau) t_1 = Float32(tau * Float32(x * Float32(pi))) return Float32(Float32(sin(Float32(x * Float32(pi))) / Float32(Float32(pi) * Float32(x * t_1))) * sin(t_1)) end
function tmp = code(x, tau) t_1 = tau * (x * single(pi)); tmp = (sin((x * single(pi))) / (single(pi) * (x * t_1))) * sin(t_1); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := tau \cdot \left(x \cdot \pi\right)\\
\frac{\sin \left(x \cdot \pi\right)}{\pi \cdot \left(x \cdot t\_1\right)} \cdot \sin t\_1
\end{array}
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
associate-*l/N/A
lift-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
associate-*l*N/A
associate-/r*N/A
*-commutativeN/A
lower-/.f32N/A
Applied rewrites97.6%
lift-/.f32N/A
lift-/.f32N/A
associate-/r*N/A
lift-*.f32N/A
lift-/.f32N/A
associate-*r/N/A
lift-*.f32N/A
associate-/r*N/A
Applied rewrites97.6%
Applied rewrites97.3%
Final simplification97.3%
(FPCore (x tau) :precision binary32 (let* ((t_1 (* PI (* x tau)))) (* (sin t_1) (/ (sin (* x PI)) (* x (* PI t_1))))))
float code(float x, float tau) {
float t_1 = ((float) M_PI) * (x * tau);
return sinf(t_1) * (sinf((x * ((float) M_PI))) / (x * (((float) M_PI) * t_1)));
}
function code(x, tau) t_1 = Float32(Float32(pi) * Float32(x * tau)) return Float32(sin(t_1) * Float32(sin(Float32(x * Float32(pi))) / Float32(x * Float32(Float32(pi) * t_1)))) end
function tmp = code(x, tau) t_1 = single(pi) * (x * tau); tmp = sin(t_1) * (sin((x * single(pi))) / (x * (single(pi) * t_1))); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \pi \cdot \left(x \cdot tau\right)\\
\sin t\_1 \cdot \frac{\sin \left(x \cdot \pi\right)}{x \cdot \left(\pi \cdot t\_1\right)}
\end{array}
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
div-invN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
Applied rewrites97.3%
Final simplification97.3%
(FPCore (x tau) :precision binary32 (* (sin (* x (* tau PI))) (/ (sin (* x PI)) (* tau (* x (* x (* PI PI)))))))
float code(float x, float tau) {
return sinf((x * (tau * ((float) M_PI)))) * (sinf((x * ((float) M_PI))) / (tau * (x * (x * (((float) M_PI) * ((float) M_PI))))));
}
function code(x, tau) return Float32(sin(Float32(x * Float32(tau * Float32(pi)))) * Float32(sin(Float32(x * Float32(pi))) / Float32(tau * Float32(x * Float32(x * Float32(Float32(pi) * Float32(pi))))))) end
function tmp = code(x, tau) tmp = sin((x * (tau * single(pi)))) * (sin((x * single(pi))) / (tau * (x * (x * (single(pi) * single(pi)))))); end
\begin{array}{l}
\\
\sin \left(x \cdot \left(tau \cdot \pi\right)\right) \cdot \frac{\sin \left(x \cdot \pi\right)}{tau \cdot \left(x \cdot \left(x \cdot \left(\pi \cdot \pi\right)\right)\right)}
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
associate-/r*N/A
lift-/.f32N/A
frac-timesN/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
lift-*.f32N/A
Applied rewrites97.2%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f32N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f3275.6
Applied rewrites75.6%
Taylor expanded in tau around 0
Applied rewrites61.1%
Taylor expanded in x around inf
associate-/l*N/A
lower-*.f32N/A
lower-sin.f32N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-/.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
unpow2N/A
associate-*l*N/A
lower-*.f32N/A
lower-*.f32N/A
Applied rewrites97.1%
(FPCore (x tau) :precision binary32 (* (sin (* x (* tau PI))) (/ (sin (* x PI)) (* PI (* PI (* tau (* x x)))))))
float code(float x, float tau) {
return sinf((x * (tau * ((float) M_PI)))) * (sinf((x * ((float) M_PI))) / (((float) M_PI) * (((float) M_PI) * (tau * (x * x)))));
}
function code(x, tau) return Float32(sin(Float32(x * Float32(tau * Float32(pi)))) * Float32(sin(Float32(x * Float32(pi))) / Float32(Float32(pi) * Float32(Float32(pi) * Float32(tau * Float32(x * x)))))) end
function tmp = code(x, tau) tmp = sin((x * (tau * single(pi)))) * (sin((x * single(pi))) / (single(pi) * (single(pi) * (tau * (x * x))))); end
\begin{array}{l}
\\
\sin \left(x \cdot \left(tau \cdot \pi\right)\right) \cdot \frac{\sin \left(x \cdot \pi\right)}{\pi \cdot \left(\pi \cdot \left(tau \cdot \left(x \cdot x\right)\right)\right)}
\end{array}
Initial program 97.7%
Taylor expanded in x around inf
associate-/l*N/A
lower-*.f32N/A
lower-sin.f32N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-/r*N/A
associate-/l/N/A
lower-/.f32N/A
Applied rewrites96.9%
Final simplification96.9%
(FPCore (x tau)
:precision binary32
(*
(sin (* PI (* x tau)))
(/
(fma
(* x x)
(fma
x
(* x (/ (* 0.008333333333333333 (* PI (* PI PI))) tau))
(/ (* PI -0.16666666666666666) tau))
(/ 1.0 (* tau PI)))
x)))
float code(float x, float tau) {
return sinf((((float) M_PI) * (x * tau))) * (fmaf((x * x), fmaf(x, (x * ((0.008333333333333333f * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))) / tau)), ((((float) M_PI) * -0.16666666666666666f) / tau)), (1.0f / (tau * ((float) M_PI)))) / x);
}
function code(x, tau) return Float32(sin(Float32(Float32(pi) * Float32(x * tau))) * Float32(fma(Float32(x * x), fma(x, Float32(x * Float32(Float32(Float32(0.008333333333333333) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))) / tau)), Float32(Float32(Float32(pi) * Float32(-0.16666666666666666)) / tau)), Float32(Float32(1.0) / Float32(tau * Float32(pi)))) / x)) end
\begin{array}{l}
\\
\sin \left(\pi \cdot \left(x \cdot tau\right)\right) \cdot \frac{\mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x, x \cdot \frac{0.008333333333333333 \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)}{tau}, \frac{\pi \cdot -0.16666666666666666}{tau}\right), \frac{1}{tau \cdot \pi}\right)}{x}
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
div-invN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
Applied rewrites97.3%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
lower-/.f32N/A
Applied rewrites83.4%
Taylor expanded in x around 0
lower-/.f32N/A
Applied rewrites90.3%
Final simplification90.3%
(FPCore (x tau)
:precision binary32
(*
(sin (* PI (* x tau)))
(/
(fma
(/ (* x x) tau)
(fma
0.008333333333333333
(* x (* x (* PI (* PI PI))))
(* PI -0.16666666666666666))
(/ 1.0 (* tau PI)))
x)))
float code(float x, float tau) {
return sinf((((float) M_PI) * (x * tau))) * (fmaf(((x * x) / tau), fmaf(0.008333333333333333f, (x * (x * (((float) M_PI) * (((float) M_PI) * ((float) M_PI))))), (((float) M_PI) * -0.16666666666666666f)), (1.0f / (tau * ((float) M_PI)))) / x);
}
function code(x, tau) return Float32(sin(Float32(Float32(pi) * Float32(x * tau))) * Float32(fma(Float32(Float32(x * x) / tau), fma(Float32(0.008333333333333333), Float32(x * Float32(x * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))))), Float32(Float32(pi) * Float32(-0.16666666666666666))), Float32(Float32(1.0) / Float32(tau * Float32(pi)))) / x)) end
\begin{array}{l}
\\
\sin \left(\pi \cdot \left(x \cdot tau\right)\right) \cdot \frac{\mathsf{fma}\left(\frac{x \cdot x}{tau}, \mathsf{fma}\left(0.008333333333333333, x \cdot \left(x \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right), \pi \cdot -0.16666666666666666\right), \frac{1}{tau \cdot \pi}\right)}{x}
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
div-invN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
Applied rewrites97.3%
Taylor expanded in x around 0
lower-/.f32N/A
Applied rewrites90.3%
Final simplification90.3%
(FPCore (x tau) :precision binary32 (let* ((t_1 (* tau (* x PI)))) (* (/ (sin t_1) t_1) (fma (* PI PI) (* (* x x) -0.16666666666666666) 1.0))))
float code(float x, float tau) {
float t_1 = tau * (x * ((float) M_PI));
return (sinf(t_1) / t_1) * fmaf((((float) M_PI) * ((float) M_PI)), ((x * x) * -0.16666666666666666f), 1.0f);
}
function code(x, tau) t_1 = Float32(tau * Float32(x * Float32(pi))) return Float32(Float32(sin(t_1) / t_1) * fma(Float32(Float32(pi) * Float32(pi)), Float32(Float32(x * x) * Float32(-0.16666666666666666)), Float32(1.0))) end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := tau \cdot \left(x \cdot \pi\right)\\
\frac{\sin t\_1}{t\_1} \cdot \mathsf{fma}\left(\pi \cdot \pi, \left(x \cdot x\right) \cdot -0.16666666666666666, 1\right)
\end{array}
\end{array}
Initial program 97.7%
Taylor expanded in x around 0
*-commutativeN/A
associate-*r*N/A
+-commutativeN/A
*-commutativeN/A
associate-*l*N/A
lower-fma.f32N/A
unpow2N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
*-commutativeN/A
lower-*.f32N/A
unpow2N/A
lower-*.f3284.0
Applied rewrites84.0%
Final simplification84.0%
(FPCore (x tau) :precision binary32 (let* ((t_1 (* x (* tau PI)))) (* (/ (sin t_1) t_1) (fma (* PI PI) (* (* x x) -0.16666666666666666) 1.0))))
float code(float x, float tau) {
float t_1 = x * (tau * ((float) M_PI));
return (sinf(t_1) / t_1) * fmaf((((float) M_PI) * ((float) M_PI)), ((x * x) * -0.16666666666666666f), 1.0f);
}
function code(x, tau) t_1 = Float32(x * Float32(tau * Float32(pi))) return Float32(Float32(sin(t_1) / t_1) * fma(Float32(Float32(pi) * Float32(pi)), Float32(Float32(x * x) * Float32(-0.16666666666666666)), Float32(1.0))) end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(tau \cdot \pi\right)\\
\frac{\sin t\_1}{t\_1} \cdot \mathsf{fma}\left(\pi \cdot \pi, \left(x \cdot x\right) \cdot -0.16666666666666666, 1\right)
\end{array}
\end{array}
Initial program 97.7%
Taylor expanded in x around inf
lower-/.f32N/A
lower-sin.f32N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3297.8
Applied rewrites97.8%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
lower-fma.f32N/A
unpow2N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
*-commutativeN/A
lower-*.f32N/A
unpow2N/A
lower-*.f3284.0
Applied rewrites84.0%
(FPCore (x tau) :precision binary32 (* (sin (* x PI)) (fma -0.16666666666666666 (* x (* tau (* tau PI))) (/ 1.0 (* x PI)))))
float code(float x, float tau) {
return sinf((x * ((float) M_PI))) * fmaf(-0.16666666666666666f, (x * (tau * (tau * ((float) M_PI)))), (1.0f / (x * ((float) M_PI))));
}
function code(x, tau) return Float32(sin(Float32(x * Float32(pi))) * fma(Float32(-0.16666666666666666), Float32(x * Float32(tau * Float32(tau * Float32(pi)))), Float32(Float32(1.0) / Float32(x * Float32(pi))))) end
\begin{array}{l}
\\
\sin \left(x \cdot \pi\right) \cdot \mathsf{fma}\left(-0.16666666666666666, x \cdot \left(tau \cdot \left(tau \cdot \pi\right)\right), \frac{1}{x \cdot \pi}\right)
\end{array}
Initial program 97.7%
lift-*.f32N/A
*-commutativeN/A
lift-/.f32N/A
div-invN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
Applied rewrites97.3%
Taylor expanded in tau around 0
lower-fma.f32N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
unpow2N/A
associate-*l*N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-/.f32N/A
lower-*.f32N/A
lower-PI.f3276.6
Applied rewrites76.6%
Final simplification76.6%
(FPCore (x tau)
:precision binary32
(*
(/ (+ (/ 1.0 (* tau PI)) (* (* x x) (/ (* PI -0.16666666666666666) tau))) x)
(*
x
(fma
(* PI (* x (* x (* PI PI))))
(* -0.16666666666666666 (* tau (* tau tau)))
(* tau PI)))))
float code(float x, float tau) {
return (((1.0f / (tau * ((float) M_PI))) + ((x * x) * ((((float) M_PI) * -0.16666666666666666f) / tau))) / x) * (x * fmaf((((float) M_PI) * (x * (x * (((float) M_PI) * ((float) M_PI))))), (-0.16666666666666666f * (tau * (tau * tau))), (tau * ((float) M_PI))));
}
function code(x, tau) return Float32(Float32(Float32(Float32(Float32(1.0) / Float32(tau * Float32(pi))) + Float32(Float32(x * x) * Float32(Float32(Float32(pi) * Float32(-0.16666666666666666)) / tau))) / x) * Float32(x * fma(Float32(Float32(pi) * Float32(x * Float32(x * Float32(Float32(pi) * Float32(pi))))), Float32(Float32(-0.16666666666666666) * Float32(tau * Float32(tau * tau))), Float32(tau * Float32(pi))))) end
\begin{array}{l}
\\
\frac{\frac{1}{tau \cdot \pi} + \left(x \cdot x\right) \cdot \frac{\pi \cdot -0.16666666666666666}{tau}}{x} \cdot \left(x \cdot \mathsf{fma}\left(\pi \cdot \left(x \cdot \left(x \cdot \left(\pi \cdot \pi\right)\right)\right), -0.16666666666666666 \cdot \left(tau \cdot \left(tau \cdot tau\right)\right), tau \cdot \pi\right)\right)
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
div-invN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
Applied rewrites97.3%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
lower-/.f32N/A
Applied rewrites83.4%
Applied rewrites83.4%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
lower-*.f32N/A
+-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
Applied rewrites75.9%
Final simplification75.9%
(FPCore (x tau)
:precision binary32
(*
(/ (fma (* x x) (/ (* PI -0.16666666666666666) tau) (/ 1.0 (* tau PI))) x)
(*
x
(fma
(* x x)
(* (* PI (* PI PI)) (* tau (* -0.16666666666666666 (* tau tau))))
(* tau PI)))))
float code(float x, float tau) {
return (fmaf((x * x), ((((float) M_PI) * -0.16666666666666666f) / tau), (1.0f / (tau * ((float) M_PI)))) / x) * (x * fmaf((x * x), ((((float) M_PI) * (((float) M_PI) * ((float) M_PI))) * (tau * (-0.16666666666666666f * (tau * tau)))), (tau * ((float) M_PI))));
}
function code(x, tau) return Float32(Float32(fma(Float32(x * x), Float32(Float32(Float32(pi) * Float32(-0.16666666666666666)) / tau), Float32(Float32(1.0) / Float32(tau * Float32(pi)))) / x) * Float32(x * fma(Float32(x * x), Float32(Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))) * Float32(tau * Float32(Float32(-0.16666666666666666) * Float32(tau * tau)))), Float32(tau * Float32(pi))))) end
\begin{array}{l}
\\
\frac{\mathsf{fma}\left(x \cdot x, \frac{\pi \cdot -0.16666666666666666}{tau}, \frac{1}{tau \cdot \pi}\right)}{x} \cdot \left(x \cdot \mathsf{fma}\left(x \cdot x, \left(\pi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(tau \cdot \left(-0.16666666666666666 \cdot \left(tau \cdot tau\right)\right)\right), tau \cdot \pi\right)\right)
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
div-invN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
Applied rewrites97.3%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
lower-/.f32N/A
Applied rewrites83.4%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
lower-*.f32N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites75.9%
Final simplification75.9%
(FPCore (x tau)
:precision binary32
(*
(/ (fma (* x x) (/ (* PI -0.16666666666666666) tau) (/ 1.0 (* tau PI))) x)
(*
tau
(fma
-0.16666666666666666
(* (* PI (* PI PI)) (* x (* tau (* tau (* x x)))))
(* x PI)))))
float code(float x, float tau) {
return (fmaf((x * x), ((((float) M_PI) * -0.16666666666666666f) / tau), (1.0f / (tau * ((float) M_PI)))) / x) * (tau * fmaf(-0.16666666666666666f, ((((float) M_PI) * (((float) M_PI) * ((float) M_PI))) * (x * (tau * (tau * (x * x))))), (x * ((float) M_PI))));
}
function code(x, tau) return Float32(Float32(fma(Float32(x * x), Float32(Float32(Float32(pi) * Float32(-0.16666666666666666)) / tau), Float32(Float32(1.0) / Float32(tau * Float32(pi)))) / x) * Float32(tau * fma(Float32(-0.16666666666666666), Float32(Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))) * Float32(x * Float32(tau * Float32(tau * Float32(x * x))))), Float32(x * Float32(pi))))) end
\begin{array}{l}
\\
\frac{\mathsf{fma}\left(x \cdot x, \frac{\pi \cdot -0.16666666666666666}{tau}, \frac{1}{tau \cdot \pi}\right)}{x} \cdot \left(tau \cdot \mathsf{fma}\left(-0.16666666666666666, \left(\pi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(x \cdot \left(tau \cdot \left(tau \cdot \left(x \cdot x\right)\right)\right)\right), x \cdot \pi\right)\right)
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
div-invN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
Applied rewrites97.3%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
lower-/.f32N/A
Applied rewrites83.4%
Taylor expanded in tau around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites75.7%
Final simplification75.7%
(FPCore (x tau) :precision binary32 (fma (* (* x PI) (* PI (fma tau (* tau -0.16666666666666666) -0.16666666666666666))) x 1.0))
float code(float x, float tau) {
return fmaf(((x * ((float) M_PI)) * (((float) M_PI) * fmaf(tau, (tau * -0.16666666666666666f), -0.16666666666666666f))), x, 1.0f);
}
function code(x, tau) return fma(Float32(Float32(x * Float32(pi)) * Float32(Float32(pi) * fma(tau, Float32(tau * Float32(-0.16666666666666666)), Float32(-0.16666666666666666)))), x, Float32(1.0)) end
\begin{array}{l}
\\
\mathsf{fma}\left(\left(x \cdot \pi\right) \cdot \left(\pi \cdot \mathsf{fma}\left(tau, tau \cdot -0.16666666666666666, -0.16666666666666666\right)\right), x, 1\right)
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
associate-/r*N/A
lift-/.f32N/A
frac-timesN/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
lift-*.f32N/A
Applied rewrites97.2%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f32N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f3275.6
Applied rewrites75.6%
Applied rewrites75.6%
Final simplification75.6%
(FPCore (x tau) :precision binary32 (fma (* x x) (* (* PI PI) (fma -0.16666666666666666 (* tau tau) -0.16666666666666666)) 1.0))
float code(float x, float tau) {
return fmaf((x * x), ((((float) M_PI) * ((float) M_PI)) * fmaf(-0.16666666666666666f, (tau * tau), -0.16666666666666666f)), 1.0f);
}
function code(x, tau) return fma(Float32(x * x), Float32(Float32(Float32(pi) * Float32(pi)) * fma(Float32(-0.16666666666666666), Float32(tau * tau), Float32(-0.16666666666666666))), Float32(1.0)) end
\begin{array}{l}
\\
\mathsf{fma}\left(x \cdot x, \left(\pi \cdot \pi\right) \cdot \mathsf{fma}\left(-0.16666666666666666, tau \cdot tau, -0.16666666666666666\right), 1\right)
\end{array}
Initial program 97.7%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f32N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f3275.6
Applied rewrites75.6%
(FPCore (x tau) :precision binary32 (fma (* x x) (* tau (* -0.16666666666666666 (* tau (* PI PI)))) 1.0))
float code(float x, float tau) {
return fmaf((x * x), (tau * (-0.16666666666666666f * (tau * (((float) M_PI) * ((float) M_PI))))), 1.0f);
}
function code(x, tau) return fma(Float32(x * x), Float32(tau * Float32(Float32(-0.16666666666666666) * Float32(tau * Float32(Float32(pi) * Float32(pi))))), Float32(1.0)) end
\begin{array}{l}
\\
\mathsf{fma}\left(x \cdot x, tau \cdot \left(-0.16666666666666666 \cdot \left(tau \cdot \left(\pi \cdot \pi\right)\right)\right), 1\right)
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
associate-/r*N/A
lift-/.f32N/A
frac-timesN/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
lift-*.f32N/A
Applied rewrites97.2%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f32N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f3275.6
Applied rewrites75.6%
Taylor expanded in tau around inf
Applied rewrites66.5%
(FPCore (x tau) :precision binary32 (fma (* x (* x (* PI PI))) -0.16666666666666666 1.0))
float code(float x, float tau) {
return fmaf((x * (x * (((float) M_PI) * ((float) M_PI)))), -0.16666666666666666f, 1.0f);
}
function code(x, tau) return fma(Float32(x * Float32(x * Float32(Float32(pi) * Float32(pi)))), Float32(-0.16666666666666666), Float32(1.0)) end
\begin{array}{l}
\\
\mathsf{fma}\left(x \cdot \left(x \cdot \left(\pi \cdot \pi\right)\right), -0.16666666666666666, 1\right)
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
associate-/r*N/A
lift-/.f32N/A
frac-timesN/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
lift-*.f32N/A
Applied rewrites97.2%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f32N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f3275.6
Applied rewrites75.6%
Taylor expanded in tau around 0
Applied rewrites61.1%
Applied rewrites61.1%
(FPCore (x tau) :precision binary32 (fma x (* (* PI PI) (* x -0.16666666666666666)) 1.0))
float code(float x, float tau) {
return fmaf(x, ((((float) M_PI) * ((float) M_PI)) * (x * -0.16666666666666666f)), 1.0f);
}
function code(x, tau) return fma(x, Float32(Float32(Float32(pi) * Float32(pi)) * Float32(x * Float32(-0.16666666666666666))), Float32(1.0)) end
\begin{array}{l}
\\
\mathsf{fma}\left(x, \left(\pi \cdot \pi\right) \cdot \left(x \cdot -0.16666666666666666\right), 1\right)
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
associate-/r*N/A
lift-/.f32N/A
frac-timesN/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
lift-*.f32N/A
Applied rewrites97.2%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f32N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f3275.6
Applied rewrites75.6%
Taylor expanded in tau around 0
Applied rewrites61.1%
Applied rewrites61.1%
(FPCore (x tau) :precision binary32 (fma x (* (* x PI) (* PI -0.16666666666666666)) 1.0))
float code(float x, float tau) {
return fmaf(x, ((x * ((float) M_PI)) * (((float) M_PI) * -0.16666666666666666f)), 1.0f);
}
function code(x, tau) return fma(x, Float32(Float32(x * Float32(pi)) * Float32(Float32(pi) * Float32(-0.16666666666666666))), Float32(1.0)) end
\begin{array}{l}
\\
\mathsf{fma}\left(x, \left(x \cdot \pi\right) \cdot \left(\pi \cdot -0.16666666666666666\right), 1\right)
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
associate-/r*N/A
lift-/.f32N/A
frac-timesN/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
lift-*.f32N/A
Applied rewrites97.2%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f32N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f3275.6
Applied rewrites75.6%
Taylor expanded in tau around 0
Applied rewrites61.1%
Applied rewrites61.1%
Final simplification61.1%
(FPCore (x tau) :precision binary32 (fma x (* x (* (* PI PI) -0.16666666666666666)) 1.0))
float code(float x, float tau) {
return fmaf(x, (x * ((((float) M_PI) * ((float) M_PI)) * -0.16666666666666666f)), 1.0f);
}
function code(x, tau) return fma(x, Float32(x * Float32(Float32(Float32(pi) * Float32(pi)) * Float32(-0.16666666666666666))), Float32(1.0)) end
\begin{array}{l}
\\
\mathsf{fma}\left(x, x \cdot \left(\left(\pi \cdot \pi\right) \cdot -0.16666666666666666\right), 1\right)
\end{array}
Initial program 97.7%
lift-*.f32N/A
lift-/.f32N/A
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
associate-/r*N/A
lift-/.f32N/A
frac-timesN/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
lift-*.f32N/A
Applied rewrites97.2%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f32N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
+-commutativeN/A
lower-fma.f32N/A
unpow2N/A
lower-*.f3275.6
Applied rewrites75.6%
Taylor expanded in tau around 0
Applied rewrites61.1%
Final simplification61.1%
(FPCore (x tau) :precision binary32 1.0)
float code(float x, float tau) {
return 1.0f;
}
real(4) function code(x, tau)
real(4), intent (in) :: x
real(4), intent (in) :: tau
code = 1.0e0
end function
function code(x, tau) return Float32(1.0) end
function tmp = code(x, tau) tmp = single(1.0); end
\begin{array}{l}
\\
1
\end{array}
Initial program 97.7%
Taylor expanded in x around 0
Applied rewrites60.1%
herbie shell --seed 2024222
(FPCore (x tau)
:name "Lanczos kernel"
:precision binary32
:pre (and (and (<= 1e-5 x) (<= x 1.0)) (and (<= 1.0 tau) (<= tau 5.0)))
(* (/ (sin (* (* x PI) tau)) (* (* x PI) tau)) (/ (sin (* x PI)) (* x PI))))