Fix spin-weighted with |m| = 1

Close #139 and #136

Author: MikaelSlevinsky

Project.toml

@@ -1,6 +1,6 @@
 name = "FastTransforms"
 uuid = "057dd010-8810-581a-b7be-e3fc3b93f78c"
-version = "0.12.1"
+version = "0.12.2"
 
 [deps]
 AbstractFFTs = "621f4979-c628-5d54-868e-fcf4e3e8185c"
@@ -25,7 +25,7 @@ BinaryProvider = "0.5"
 DSP = "0.6"
 FFTW = "1"
 FastGaussQuadrature = "0.4"
-FastTransforms_jll = "0.5.0"
+FastTransforms_jll = "0.5.1"
 FillArrays = "0.9, 0.10, 0.11"
 Reexport = "0.2, 1.0"
 SpecialFunctions = "0.10, 1"

deps/build.jl

@@ -42,4 +42,6 @@ if ft_build_from_source == "true"
         print_error()
     end
     println("FastTransforms built from source.")
+else
+    println("FastTransforms using precompiled binaries.")
 end

src/libfasttransforms.jl

@@ -10,13 +10,14 @@ else
     using FastTransforms_jll
 end
 
-function ft_fftw_plan_with_nthreads(n::Integer)
-    ccall((:ft_fftw_plan_with_nthreads, libfasttransforms), Cvoid, (Cint, ), n)
-end
+ft_set_num_threads(n::Integer) = ccall((:ft_set_num_threads, libfasttransforms), Cvoid, (Cint, ), n)
+ft_fftw_plan_with_nthreads(n::Integer) = ccall((:ft_fftw_plan_with_nthreads, libfasttransforms), Cvoid, (Cint, ), n)
 
 function __init__()
+    n = ceil(Int, Sys.CPU_THREADS/2)
+    ft_set_num_threads(n)
     ccall((:ft_fftw_init_threads, libfasttransforms), Cint, ())
-    ft_fftw_plan_with_nthreads(ceil(Int, Sys.CPU_THREADS/2))
+    ft_fftw_plan_with_nthreads(n)
 end
 
 
@@ -49,8 +50,6 @@ function renew!(x::Array{BigFloat})
     return x
 end
 
-set_num_threads(n::Integer) = ccall((:ft_set_num_threads, libfasttransforms), Cvoid, (Cint, ), n)
-
 function horner!(c::Vector{Float64}, x::Vector{Float64}, f::Vector{Float64})
     @assert length(x) == length(f)
     ccall((:ft_horner, libfasttransforms), Cvoid, (Cint, Ptr{Float64}, Cint, Cint, Ptr{Float64}, Ptr{Float64}), length(c), c, 1, length(x), x, f)