Merge pull request #40 from slimgroup/fix-dft

fix iffshift in centered apply_ifft

Author: mloubout

Diff for: Project.toml

@@ -1,7 +1,7 @@
 name = "JOLI"
 uuid = "bb331ad6-a1cf-11e9-23da-9bcb53c69f6f"
 authors = ["Henryk Modzelewski <[email protected]>"]
-version = "0.8.0"
+version = "0.8.1"
 
 [deps]
 Distributed = "8ba89e20-285c-5b6f-9357-94700520ee1b"

Diff for: src/joLinearFunctionConstructors/joDFT.jl

@@ -45,7 +45,7 @@ module joDFT_etc
         mp=prod(ms)
         rv=reshape(v,msc)
         pf=plan_ifft(rv,dims)
-        if centered rv=fftshift(rv,dims) end
+        if centered rv=ifftshift(rv,dims) end
         rv=(pf*rv)*sqrt(mp)
         rv=reshape(rv,(prod(ms),nvc))
         rv=jo_convert(rdt,rv,false)
@@ -92,7 +92,7 @@ module joDFT_etc
         mp=prod(ms)
         rv=reshape(v,msc)
         pf=plan_ifft(rv,dims)
-        if centered rv=fftshift(rv,dims) end
+        if centered rv=ifftshift(rv,dims) end
         rv=(pf*rv)*sqrt(mp)
         rv=reshape(rv,(prod(ms),nvc))
         rv=jo_convert(rdt,rv,false)

Diff for: test/test_joDFT.jl

@@ -42,6 +42,11 @@ for t=1:T # start test loop
         @test norm((adjoint(A2)*A2)*vmv2-vmv2)<joTol
     end
 
+    ### test https://github.com/slimgroup/JOLI.jl/issues/39
+    A1=joDFT(m;DDT=ComplexF64,centered=true)
+    verbose && println("$tsname ($m[,$m]) - centered")
+    @test A1 * v1 == vec(A1 * reshape(v1, :, 1))
+
     A1=joDFT(m;planned=false,DDT=ComplexF64)
     A2=joDFT(m,m;planned=false,DDT=ComplexF64)
     verbose && println("$tsname ($m[,$m]) - not planned")