Use cld instead of ceil & division

Sagnac · Jul 2, 2023 · 1f6adcd · 1f6adcd
1 parent f5f6185
commit 1f6adcd
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/WavefrontError.jl b/src/WavefrontError.jl
@@ -75,7 +75,7 @@ end
 
 # synthesis function
 function Λ(ΔW, a, v, n_max; finesse::Int)
-    finesse = finesse ∈ 1:100 ? finesse : ceil(Int, 100 / √ length(a))
+    finesse::Int = finesse ∈ 1:100 ? finesse : cld(100, √ length(a))
     ρ, θ = polar(n_max, n_max; finesse)
     # construct the estimated wavefront error
     ΔWp = ΔW.(ρ', θ)

diff --git a/src/Zernike.jl b/src/Zernike.jl
@@ -79,7 +79,7 @@ function polar(m::Int, n::Int; finesse::Int = 100)
 end
 
 # radial order
-get_n(j::Int)::Int = ceil((-3 + √(9 + 8j)) / 2)
+get_n(j::Int)::Int = cld(-3 + √(9 + 8j), 2)
 # azimuthal frequency
 get_m(j::Int, n::Int) = 2j - (n + 2)n
 # ISO / ANSI / OSA standard single mode-ordering index