tests pass locally!

Author: MikaelSlevinsky

.gitignore

@@ -1,3 +1,3 @@
 docs/build/
 docs/site/
-FastTransforms
+deps/FastTransforms

Project.toml

@@ -1,27 +1,21 @@
 name = "FastTransforms"
 uuid = "057dd010-8810-581a-b7be-e3fc3b93f78c"
-version = "0.5"
+version = "0.5.0"
 
 [deps]
 AbstractFFTs = "621f4979-c628-5d54-868e-fcf4e3e8185c"
-Compat = "34da2185-b29b-5c13-b0c7-acf172513d20"
 DSP = "717857b8-e6f2-59f4-9121-6e50c889abd2"
 FFTW = "7a1cc6ca-52ef-59f5-83cd-3a7055c09341"
-HierarchicalMatrices = "7c893195-952b-5c83-bb6e-be12f22eed51"
+Libdl = "8f399da3-3557-5675-b5ff-fb832c97cbdb"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
-LowRankApprox = "898213cb-b102-5a47-900c-97e73b919f73"
-ProgressMeter = "92933f4c-e287-5a05-a399-4b506db050ca"
 SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
+Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 ToeplitzMatrices = "c751599d-da0a-543b-9d20-d0a503d91d24"
 
 [compat]
 AbstractFFTs = "≥ 0.3.1"
-Compat = "≥ 0.17.0"
 DSP = "≥ 0.4.0"
 FFTW = "≥ 0.0.4"
-HierarchicalMatrices = "≥ 0.1.1"
-LowRankApprox = "≥ 0.1.1"
-ProgressMeter = "≥ 0.3.4"
 SpecialFunctions = "≥ 0.3.4"
 ToeplitzMatrices = "≥ 0.4.0"
-julia = "0.7, 1"
+julia = "≥ 1.0"

examples/chebyshev.jl

@@ -7,7 +7,7 @@ using FastTransforms
 
 # first kind points -> first kind polynomials
 n = 20
-p_1 = [sinpi((n-2k-1)/2n) for k=0:n-1]
+p_1 = chebyshevpoints(Float64, n; kind=1)
 f = exp.(p_1)
 f̌ = chebyshevtransform(f; kind=1)
 
@@ -18,7 +18,7 @@ f̃(0.1) ≈ exp(0.1)
 ichebyshevtransform(f̌; kind=1) ≈ exp.(p_1)
 
 # second kind points -> first kind polynomials
-p_2 = [cospi(k/(n-1)) for k=0:n-1]
+p_2 = chebyshevpoints(Float64, n; kind=2)
 f = exp.(p_2)
 f̌ = chebyshevtransform(f; kind=2)
 
@@ -31,7 +31,7 @@ ichebyshevtransform(f̌; kind=2) ≈ exp.(p_2)
 
 # first kind points -> second kind polynomials
 n = 20
-p_1 = [sinpi((n-2k-1)/2n) for k=0:n-1]
+p_1 = chebyshevpoints(Float64, n; kind=1)
 f = exp.(p_1)
 f̌ = chebyshevutransform(f; kind=1)
 f̃ = x -> [sin((k+1)*acos(x))/sin(acos(x)) for k=0:n-1]' * f̌
@@ -42,7 +42,7 @@ ichebyshevutransform(f̌; kind=1) ≈ exp.(p_1)
 
 
 # second kind points -> second kind polynomials
-p_2 = [cospi(k/(n-1)) for k=1:n-2]
+p_2 = chebyshevpoints(Float64, n; kind=2)[2:n-1]
 f = exp.(p_2)
 f̌ = chebyshevutransform(f; kind=2)
 f̃ = x -> [sin((k+1)*acos(x))/sin(acos(x)) for k=0:n-3]' * f̌

examples/sphere.jl

@@ -51,35 +51,34 @@ P4 = x -> (35*x^4-30*x^2+3)/8
 # On the tensor product grid, our function samples are:
 F = [(P4(z(θ,φ)⋅y) - P4(x⋅y))/(z(θ,φ)⋅y - x⋅y) for θ in θ, φ in φ]
 
-P = plan_sph2fourier(F);
-PA = FastTransforms.plan_analysis(F);
+P = plan_sph2fourier(Float64, N)
+PA = plan_sph_analysis(Float64, N, M)
 
 # Its spherical harmonic coefficients demonstrate that it is degree-3:
-V = zero(F);
-A_mul_B!(V, PA, F);
+V = PA*F
 U3 = P\V
 
 # Similarly, on the tensor product grid, the Legendre polynomial P₄(z⋅y) is:
 F = [P4(z(θ,φ)⋅y) for θ in θ, φ in φ]
 
 # Its spherical harmonic coefficients demonstrate that it is exact-degree-4:
-A_mul_B!(V, PA, F);
+V = PA*F
 U4 = P\V
 
-nrm1 = vecnorm(U4);
+nrm1 = norm(U4);
 
 # Finally, the Legendre polynomial P₄(z⋅x) is aligned with the grid:
 F = [P4(z(θ,φ)⋅x) for θ in θ, φ in φ]
 
 # It only has one nonnegligible spherical harmonic coefficient.
 # Can you spot it?
-A_mul_B!(V, PA, F);
+V = PA*F
 U4 = P\V
 
 # That nonnegligible coefficient should be approximately √(2π/(4+1/2)),
 # since the convention in this library is to orthonormalize.
 
-nrm2 = vecnorm(U4);
+nrm2 = norm(U4);
 
 # Note that the integrals of both functions P₄(z⋅y) and P₄(z⋅x) and their
 # L²(𝕊²) norms are the same because of rotational invariance. The integral of