Revision 73a9e57fff6f0f4457a1a350a761b3c9bc3f5c3f authored by Tobias on 11 February 2022, 13:52:40 UTC, committed by Tobias on 11 February 2022, 13:52:40 UTC
1 parent a8bba44
constructors.jl
@testset "Constructors" begin
@test_throws ArgumentError NFFTPlan(zeros(1,4), (2,2))
@test_throws ArgumentError NFFTPlan(zeros(2,4), (2,3))
p = NFFTPlan(zeros(2,4), (2,2))
pCopy = copy(p)
for n in fieldnames(typeof(p))
println(n)
if n ∉ [:tmpVec, :forwardFFT, :backwardFFT, :blocks, :nodesInBlock]
@test getfield(p,n) == getfield(pCopy,n)
end
end
@show p
## test range error
x = [-0.6 0.9; 0.5 -0.5]
@test_throws ArgumentError NFFTPlan(x, (2,2))
## test nodes!(p, tr)
Nx = 32
trj1 = rand(2, 1000) .- 0.5
trj2 = rand(2, 1000) .- 0.5
p1 = NFFT.NFFTPlan(trj1, (Nx, Nx))
p2 = NFFT.NFFTPlan(trj2, (Nx, Nx))
NFFT.nodes!(p2, trj1)
@test p1.N == p2.N
@test p1.NOut == p2.NOut
@test p1.M == p2.M
@test all(p1.x .== p2.x)
@test p1.n == p2.n
@test p1.dims == p2.dims
for n in fieldnames(typeof(p1.params))
@test getfield(p1.params,n) == getfield(p2.params,n)
end
@test p1.windowLUT == p2.windowLUT
@test p1.windowHatInvLUT == p2.windowHatInvLUT
@test p1.B == p2.B
for n in fieldnames(typeof(p1.forwardFFT))
println(n)
if n != :pinv
@test getfield(p1.forwardFFT,n) == getfield(p1.forwardFFT,n)
@test getfield(p1.backwardFFT,n) == getfield(p1.backwardFFT,n)
end
end
end
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...