comment out matrix transforms, add chebyshev point tests

Author: dlfivefifty

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=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
+#
 
 
 ## Chebyshev U

test/chebyshevtests.jl

@@ -1,6 +1,19 @@
 using FastTransforms, Compat, Compat.Test
 
 @testset "Chebyshev transform"  begin
+    @testset "Chebyshev points" begin
+        @test chebyshevpoints(10; kind=1) == chebyshevpoints(Float64, 10; kind=1)
+        @test chebyshevpoints(10; kind=2) == chebyshevpoints(Float64, 10; kind=2)
+        for T in (Float32, Float64, ComplexF32, ComplexF64)
+            @test chebyshevpoints(T, 0, kind=1) == chebyshevpoints(T, 0, kind=2) == T[]
+            @test chebyshevpoints(T, 1, kind=1) == chebyshevpoints(T, 1, kind=2) == T[]
+
+            n = 20
+            @inferred(chebyshevpoints(T, n, kind=1))
+            @inferred(chebyshevpoints(T, n, kind=1))
+        end
+    end
+
     @testset "Chebyshev first kind points <-> first kind coefficients" begin
         for T in (Float32, Float64, ComplexF32, ComplexF64)
             n = 20