enable use of fftw threads

Adds fftw_init_threads & fftw_plan_with_nthreads functions.

base/signal_fftw.jl

@@ -1,8 +1,5 @@
 ## FFT: Implement fft by calling fftw.
 
-_jl_libfftw = dlopen("libfftw3")
-_jl_libfftwf = dlopen("libfftw3f")
-
 ## Direction of FFT
 
 const _jl_FFTW_FORWARD = int32(-1)
@@ -35,6 +32,20 @@ const _jl_FFTW_RODFT11 = int32(10)
 
 ## Julia wrappers around FFTW functions
 
+# Threads
+
+function fftw_init_threads()
+    stat = ccall(dlsym(_jl_libfftw,:fftw_init_threads), Int32, ())
+    @assert stat != 0
+    stat = ccall(dlsym(_jl_libfftwf,:fftwf_init_threads), Int32, ())
+    @assert stat != 0
+end
+
+function fftw_plan_with_nthreads(nthreads::Int)
+    ccall(dlsym(_jl_libfftw,:fftw_plan_with_nthreads), Void, (Int32,), nthreads)
+    ccall(dlsym(_jl_libfftwf,:fftwf_plan_with_nthreads), Void, (Int32,), nthreads)
+end
+
 # Execute
 
 _jl_fftw_execute(precision::Union(Type{Float64}, Type{Complex128}), plan) =

base/start_image.jl

@@ -17,8 +17,8 @@ _jl_librandom = dlopen("librandom"); _jl_librandom_init();
 # Optional libraries
 const _jl_libblas = dlopen(_jl_libblas_name)
 const _jl_liblapack = _jl_libblas
-_jl_libfftw = dlopen("libfftw3")
-_jl_libfftwf = dlopen("libfftw3f")
+const _jl_libfftw = dlopen("libfftw3_threads")
+const _jl_libfftwf = dlopen("libfftw3f_threads")
 
 ##_jl_libglpk = dlopen("libglpk")
 ##_jl_libglpk = dlopen("libglpk_wrapper")

deps/Makefile

@@ -519,7 +519,7 @@ fftw-$(FFTW_VER)-single/configure: fftw-$(FFTW_VER).tar.gz
 	touch $@
 fftw-$(FFTW_VER)-single/config.status: fftw-$(FFTW_VER)-single/configure
 	cd fftw-$(FFTW_VER)-single && \
-	./configure --prefix=$(abspath $(USR)) $(FFTW_CONFIG) --enable-shared --disable-fortran --enable-sse --enable-single --enable-threads && \
+	./configure --prefix=$(abspath $(USR)) $(FFTW_CONFIG) --enable-sse --enable-single && \
 	$(MAKE) clean
 	touch $@
 $(FFTW_SINGLE_OBJ_TARGET): fftw-$(FFTW_VER)-single/config.status