Merge pull request #44 from simonsobs/patch_lmax

adjust lmax to not be crazy high
simonsobs · Feb 26, 2022 · 0d414bc · 0d414bc
2 parents 68360ce + c709ff9
commit 0d414bc
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/transforms.jl b/src/transforms.jl
@@ -76,7 +76,7 @@ function create_sht_band(m::Enmap)
 end
 
 """Generate an estimate of the Nyquist frequency for a map"""
-getlmax(wcs) = 3 * fullringsize(wcs)
+getlmax(wcs) = fullringsize(wcs) ÷ 2
 
 """perform forward SHT of an intensity map"""
 function map2alm(input_map::Enmap{T,2}; lmax=nothing, mmax=lmax) where T

diff --git a/test/test_transforms.jl b/test/test_transforms.jl
@@ -18,7 +18,7 @@ end
     ref_alms = data[:,1] + 1im .* data[:,2]
     @test (alms.alm ≈ ref_alms)
 
-    @test length(map2alm(m).alm) ==5995
+    @test length(map2alm(m).alm) == 190
 
     alms = map2alm(m[6:end-2, 5:end-3]; lmax=18)
     data = readdlm("data/simple_analytic_sht_sliced.txt")