fix fftw threads

Project.toml

@@ -9,6 +9,7 @@ Hwloc = "0e44f5e4-bd66-52a0-8798-143a42290a1d"
 JLD2 = "033835bb-8acc-5ee8-8aae-3f567f8a3819"
 PaddedViews = "5432bcbf-9aad-5242-b902-cca2824c8663"
 QuadGK = "1fd47b50-473d-5c70-9696-f719f8f3bcdc"
+Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
 SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
 Tullio = "bc48ee85-29a4-5162-ae0b-a64e1601d4bc"
 UnPack = "3a884ed6-31ef-47d7-9d2a-63182c4928ed"

src/QuantumFluidSpectra.jl

@@ -1,8 +1,9 @@
 module QuantumFluidSpectra
 
+using Reexport
 using Tullio
 using Hwloc
-using FFTW
+@reexport using FFTW
 using SpecialFunctions
 using PaddedViews
 using UnPack

test/runtests.jl

@@ -3,7 +3,7 @@ using Test
 
 @testset "2D Analysis" begin 
     # Velocity and Helmholtz tests
-    n = 100
+    n = 256
     L = (1,1)
     N = (n,n)
     X,K,dX,dK = xk_arrays(L,N)
@@ -24,7 +24,7 @@ using Test
 
     ## Orthogonality
     vidotvc = Vi[1].*Vc[1] .+ Vi[2].*Vc[2]
-    @test maximum(abs.(vidotvc)) < 1e-10
+    @test maximum(abs.(vidotvc)) < 1e-8
 
     ## Projective
     @test Vi[1] .+ Vc[1] ≈ vx
@@ -36,7 +36,7 @@ using Test
 
     Natoms = sum(abs2.(ψ))*dx^2
     AC = auto_correlate(ψ,X,K)
-    Nr,Nim = AC[101,101] .|> (real,imag) 
+    Nr,Nim = AC[n+1,n+1] .|> (real,imag) 
 
     @test Natoms ≈ Nr
     @test Nim ≈ 0.0 
@@ -49,7 +49,7 @@ end
 
 @testset "3D Analysis" begin 
     # Velocity and Helmholtz tests
-    n = 20
+    n = 64
     L = (1,1,1)
     N = (n,n,n)
     X,K,dX,dK = xk_arrays(L,N)
@@ -84,7 +84,7 @@ end
 
     Natoms = sum(abs2.(ψ))*dx^3
     AC = auto_correlate(ψ,X,K)
-    Nr,Nim = AC[21,21,21] .|> (real,imag) 
+    Nr,Nim = AC[n+1,n+1,n+1] .|> (real,imag) 
 
     @test Natoms ≈ Nr
     @test Nim ≈ 0.0