zymon/RoomAcoustics.jl
1
0
Fork 0
Code Issues Pull requests Releases Activity

rm random stuff

Browse source
This commit is contained in:
zymon 2024-06-11 14:57:43 +02:00
parent f0c2574ff4
commit f0fa05284e
1 changed files with 0 additions and 28 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

28
random/sinc.jl
View file

@ -1,28 +0,0 @@
# @INFO: there is no SIMD version for sinc (sinpi etc.)
# We need a sinc kernel that can be SIMD vectorized
# Below is a sinple comparision of accuracy between a few approaches
using GLMakie
pix(x) = π * x
sinc2(x) = x |> pix |> (x -> iszero(x) ? 1.0 : sin(x)/x)
sinc3(x) = x |> pix |> (x -> x + eps(x)) |> (x -> sin(x)/x)
fs = 480e3
Δt = 1 / fs
t = -0.5:Δt:0.5
y_sinc = t .|> sinc;
y_sinc2 = t .|> sinc2;
y_sinc3 = t .|> sinc3;
@show maximum(abs, y_sinc .- y_sinc2)
@show maximum(abs, y_sinc .- y_sinc3)
fig = Figure()
ax = Axis(fig[1, 1])
lines!(ax, t, abs.(y_sinc - y_sinc2) .|> log10; label="if sinc")
lines!(ax, t, abs.(y_sinc - y_sinc3) .|> log10; label="eps sinc")
ylims!(ax, -16, -12)
xlims!(ax, (t |> extrema)...)
axislegend(ax)