updates NFCT and NFFCt accuracy

Author: mischmi96

Project.toml

@@ -1,7 +1,7 @@
 name = "GroupedTransforms"
 uuid = "76248ec8-aa36-4347-850f-ece2d673bc9f"
 authors = ["Michael Schmischke <[email protected]> and Felix Bartel <[email protected]>"]
-version = "1.1.1"
+version = "1.1.2"
 
 [deps]
 Aqua = "4c88cf16-eb10-579e-8560-4a9242c79595"

src/NFCTtools.jl

@@ -118,7 +118,7 @@ function get_transform(bandwidths::Vector{Int}, X::Array{Float64})::Int64
 
     mask = nfct_mask(bandwidths)
 
-    plan = NFCT(Tuple(bandwidths), M)
+    plan = NFCT(Tuple(bandwidths), M, Tuple(2 * collect(bandwidths)), 5)
     plan.x = X
     factor_sqrt = sqrt(2)^d
 

src/NFFTtools.jl

@@ -118,7 +118,7 @@ function get_transform(bandwidths::Vector{Int}, X::Array{Float64})::Int64
 
     mask = nfft_mask(bandwidths)
     N = Tuple(bandwidths)
-    plan = NFFT(N, M, Tuple(2 * collect(N)), 6)
+    plan = NFFT(N, M, Tuple(2 * collect(N)), 5)
     plan.x = X
 
     function trafo(fhat::Vector{ComplexF64})::Vector{ComplexF64}

test/nfct_U.jl

@@ -51,7 +51,7 @@ f_direct = F_direct * vec(fhat)
 # compare results ####################################################
 
 error = norm(f - f_direct)
-@test error < 1e-8
+@test error < 1e-5
 
 # generate random function values ####################################
 
@@ -68,4 +68,4 @@ fhat_direct = F_direct' * y
 # compare results ####################################################
 
 error = norm(vec(fhat) - fhat_direct)
-@test error < 1e-8
+@test error < 1e-5

test/nfct_ds.jl

@@ -49,7 +49,7 @@ f_direct = F_direct * vec(fhat)
 # compare results ####################################################
 
 error = norm(f - f_direct)
-@test error < 1e-8
+@test error < 1e-5
 
 # generate random function values ####################################
 
@@ -66,4 +66,4 @@ fhat_direct = F_direct' * y
 # compare results ####################################################
 
 error = norm(vec(fhat) - fhat_direct)
-@test error < 1e-8
+@test error < 1e-5