Skip to content

perf: pre-compute FFT plan #16

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
wants to merge 2 commits into
base: master
Choose a base branch
from
Open

perf: pre-compute FFT plan #16

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
82a75ca
perf: pre-compute FFT plan
johnnychen94 Oct 17, 2022
58354b6
ImageSmooth v0.1.2
johnnychen94 Oct 17, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion Project.toml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
name = "ImageSmooth"
uuid = "a14a6c0f-317a-4379-b4a0-4ba4e1e9b5ed"
authors = ["Johnny Chen <[email protected]>", "JKay <[email protected]>"]
version = "0.1.1"
version = "0.1.2"

[deps]
FFTW = "7a1cc6ca-52ef-59f5-83cd-3a7055c09341"
Expand Down
12 changes: 7 additions & 5 deletions src/algorithms/l0_smooth.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -128,6 +128,9 @@ function (f::L0Smooth)(out::GenericGrayImage,
t¹ = trues(N, M)
ℱ𝑆 = similar(ℱ𝐼)

# precompute the FFT plan so that we get fast FFT inside the iteration
F = plan_fft!(ℱ𝐼, (1, 2))
IF = plan_ifft!(ℱ𝐼, (1, 2))
while 𝛽 < 𝛽max
# Computing (ℎ, 𝑣) via solving equation (9) in [1]
# We get the solution (12) in [1] through following process
Expand All @@ -152,14 +155,13 @@ function (f::L0Smooth)(out::GenericGrayImage,
# 𝛥₁ᵀ() and 𝛥₂ᵀ() indicate the transposition of forward difference along horizontal axis and vertical axis
fdiff!(𝛥₁ᵀℎ, 𝛥₁𝑆, dims = 2, rev=true, boundary=:periodic)
fdiff!(𝛥₂ᵀ𝑣, 𝛥₂𝑆, dims = 1, rev=true, boundary=:periodic)
@. 𝛥₁ᵀℎ = -𝛥₁ᵀℎ
@. 𝛥₂ᵀ𝑣 = -𝛥₂ᵀ𝑣

# Computing S via equation (8) in [1]
@. Normin = complex(𝛥₁ᵀℎ + 𝛥₂ᵀ𝑣)
fft!(Normin, (1, 2))
@. Normin = complex(-𝛥₁ᵀℎ - 𝛥₂ᵀ𝑣)
F * Normin # fft!(Normin, (1, 2))

@. ℱ𝑆 = (ℱ𝐼 + 𝛽 * Normin) / (1 + 𝛽 * Denormin)
ifft!(ℱ𝑆, (1, 2))
IF * ℱ𝑆 # ifft!(ℱ𝑆, (1, 2))
@. 𝑆 = real(ℱ𝑆)

𝛽 = 𝛽 * 𝜅
Expand Down