FastTransform 0.7 compatibility

Author: dlfivefifty

src/FastTransforms.jl

@@ -92,8 +92,8 @@ include("hermite.jl")
 #cheb2leg(x...)=th_cheb2leg(x...)
 leg2chebu(x...) = th_leg2chebu(x...)
 ultra2ultra(x...) = th_ultra2ultra(x...)
-jac2jac(x...) = th_jac2jac(x...)
-plan_jac2jac(x...) = th_jac2jacplan(x...)
+jac2jac(x...) = jjt(x...)
+plan_jac2jac(x...) = plan_jjt(x...)
 plan_ultra2ultra(x...) = th_ultra2ultraplan(x...)
 
 include("hierarchical.jl")
@@ -106,4 +106,20 @@ include("gaunt.jl")
 include("precompile.jl")
 _precompile_()
 
+
+##
+# FastTransform v0.7 compatibility
+##
+
+export plan_cheb2jac, plan_jac2cheb, cheb2jac, jac2cheb
+
+plan_cheb2jac(v...) = plan_icjt(v...)
+plan_jac2cheb(v...) = plan_cjt(v...)
+cheb2jac(v...) = icjt(v...)
+jac2cheb(v...) = cjt(v...)
+
+plan_cheb2leg(::Type{T}, n::Int) where T = th_cheb2legplan(T, n)
+plan_leg2cheb(::Type{T}, n::Int) where T = th_leg2chebplan(T, n)
+
+
 end # module

src/chebyshevtransform.jl

@@ -119,51 +119,51 @@ end
 
 
 # Matrix inputs
-#
-#
-# function chebyshevtransform!(X::AbstractMatrix{T}; kind::Integer=1) where T<:fftwNumber
-#     if kind == 1
-#         if size(X) == (1,1)
-#             X
-#         else
-#             X=r2r!(X,REDFT10)
-#             X[:,1]/=2;X[1,:]/=2;
-#             lmul!(1/(size(X,1)*size(X,2)),X)
-#         end
-#     elseif kind == 2
-#         if size(X) == (1,1)
-#             X
-#         else
-#             X=r2r!(X,REDFT00)
-#             lmul!(1/((size(X,1)-1)*(size(X,2)-1)),X)
-#             X[:,1]/=2;X[:,end]/=2
-#             X[1,:]/=2;X[end,:]/=2
-#             X
-#         end
-#     end
-# end
-#
-# function ichebyshevtransform!(X::AbstractMatrix{T}; kind::Integer=1) where T<:fftwNumber
-#     if kind == 1
-#         if size(X) == (1,1)
-#             X
-#         else
-#             X[1,:]*=2;X[:,1]*=2
-#             X = r2r(X,REDFT01)
-#             lmul!(0.25, X)
-#         end
-#     elseif kind == 2
-#         if size(X) == (1,1)
-#             X
-#         else
-#             X[1,:]*=2;X[end,:]*=2;X[:,1]*=2;X[:,end]*=2
-#             X=chebyshevtransform!(X;kind=kind)
-#             X[1,:]*=2;X[end,:]*=2;X[:,1]*=2;X[:,end]*=2
-#             lmul!((size(X,1)-1)*(size(X,2)-1)/4,X)
-#         end
-#     end
-# end
-#
+
+
+function chebyshevtransform!(X::AbstractMatrix{T}; kind::Integer=1) where T<:fftwNumber
+    if kind == 1
+        if size(X) == (1,1)
+            X
+        else
+            X= FFTW.r2r!(X,REDFT10)
+            X[:,1]/=2;X[1,:]/=2;
+            lmul!(1/(size(X,1)*size(X,2)),X)
+        end
+    elseif kind == 2
+        if size(X) == (1,1)
+            X
+        else
+            X= FFTW.r2r!(X,REDFT00)
+            lmul!(1/((size(X,1)-1)*(size(X,2)-1)),X)
+            X[:,1]/=2;X[:,end]/=2
+            X[1,:]/=2;X[end,:]/=2
+            X
+        end
+    end
+end
+
+function ichebyshevtransform!(X::AbstractMatrix{T}; kind::Integer=1) where T<:fftwNumber
+    if kind == 1
+        if size(X) == (1,1)
+            X
+        else
+            X[1,:]*=2;X[:,1]*=2
+            X =  FFTW.r2r(X,REDFT01)
+            lmul!(0.25, X)
+        end
+    elseif kind == 2
+        if size(X) == (1,1)
+            X
+        else
+            X[1,:]*=2;X[end,:]*=2;X[:,1]*=2;X[:,end]*=2
+            X=chebyshevtransform!(X;kind=kind)
+            X[1,:]*=2;X[end,:]*=2;X[:,1]*=2;X[:,end]*=2
+            lmul!((size(X,1)-1)*(size(X,2)-1)/4,X)
+        end
+    end
+end
+
 
 
 ## Chebyshev U