make ChebyshevToLegendrePlan consistent

th_cheb2legplan now creates the ChebyshevToLegendrePlan.

Author: MikaelSlevinsky

docs/api/FastTransforms.md

@@ -13,7 +13,7 @@ See also [`icjt`](#method__icjt.1) and [`jjt`](#method__jjt.1).
 
 
 *source:*
-[FastTransforms/src/cjt.jl:121](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/cjt.jl#L121)
+[FastTransforms/src/cjt.jl:121](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/cjt.jl#L121)
 
 ---
 
@@ -23,7 +23,7 @@ Calculates the Gaunt coefficients in 64-bit floating-point arithmetic.
 
 
 *source:*
-[FastTransforms/src/gaunt.jl:24](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/gaunt.jl#L24)
+[FastTransforms/src/gaunt.jl:24](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/gaunt.jl#L24)
 
 ---
 
@@ -43,7 +43,7 @@ This is a Julia implementation of the stable recurrence described in:
 
 
 *source:*
-[FastTransforms/src/gaunt.jl:14](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/gaunt.jl#L14)
+[FastTransforms/src/gaunt.jl:14](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/gaunt.jl#L14)
 
 ---
 
@@ -56,7 +56,7 @@ See also [`cjt`](#method__cjt.1) and [`jjt`](#method__jjt.1).
 
 
 *source:*
-[FastTransforms/src/cjt.jl:129](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/cjt.jl#L129)
+[FastTransforms/src/cjt.jl:129](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/cjt.jl#L129)
 
 ---
 
@@ -69,7 +69,7 @@ See also [`cjt`](#method__cjt.1) and [`icjt`](#method__icjt.1).
 
 
 *source:*
-[FastTransforms/src/cjt.jl:137](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/cjt.jl#L137)
+[FastTransforms/src/cjt.jl:137](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/cjt.jl#L137)
 
 ---
 
@@ -88,7 +88,7 @@ Optionally:
 
 
 *source:*
-[FastTransforms/src/cjt.jl:158](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/cjt.jl#L158)
+[FastTransforms/src/cjt.jl:158](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/cjt.jl#L158)
 
 ---
 
@@ -107,7 +107,7 @@ Optionally:
 
 
 *source:*
-[FastTransforms/src/cjt.jl:172](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/cjt.jl#L172)
+[FastTransforms/src/cjt.jl:172](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/cjt.jl#L172)
 
 ## Internal
 
@@ -122,7 +122,7 @@ Modified Chebyshev moments of the first kind with respect to the Jacobi weight:
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:399](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L399)
+[FastTransforms/src/specialfunctions.jl:399](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L399)
 
 ---
 
@@ -135,7 +135,7 @@ Modified Chebyshev moments of the second kind with respect to the Jacobi weight:
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:417](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L417)
+[FastTransforms/src/specialfunctions.jl:417](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L417)
 
 ---
 
@@ -145,7 +145,7 @@ Compute weights of the Clenshaw—Curtis quadrature rule with a Jacobi weight.
 
 
 *source:*
-[FastTransforms/src/clenshawcurtis.jl:12](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/clenshawcurtis.jl#L12)
+[FastTransforms/src/clenshawcurtis.jl:12](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/clenshawcurtis.jl#L12)
 
 ---
 
@@ -155,7 +155,7 @@ Compute nodes and weights of the Clenshaw—Curtis quadrature rule with a Jacobi
 
 
 *source:*
-[FastTransforms/src/clenshawcurtis.jl:6](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/clenshawcurtis.jl#L6)
+[FastTransforms/src/clenshawcurtis.jl:6](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/clenshawcurtis.jl#L6)
 
 ---
 
@@ -165,7 +165,7 @@ Compute Jacobi expansion coefficients in Pₙ^(α-1,β) given Jacobi expansion c
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:467](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L467)
+[FastTransforms/src/specialfunctions.jl:467](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L467)
 
 ---
 
@@ -175,7 +175,7 @@ Compute Jacobi expansion coefficients in Pₙ^(α-1,α-1) given Jacobi expansion
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:489](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L489)
+[FastTransforms/src/specialfunctions.jl:489](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L489)
 
 ---
 
@@ -185,7 +185,7 @@ Compute Jacobi expansion coefficients in Pₙ^(α,β-1) given Jacobi expansion c
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:478](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L478)
+[FastTransforms/src/specialfunctions.jl:478](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L478)
 
 ---
 
@@ -195,7 +195,7 @@ Compute nodes and weights of Fejer's first quadrature rule with a Jacobi weight.
 
 
 *source:*
-[FastTransforms/src/fejer.jl:7](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/fejer.jl#L7)
+[FastTransforms/src/fejer.jl:7](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/fejer.jl#L7)
 
 ---
 
@@ -205,7 +205,7 @@ Compute nodes and weights of Fejer's second quadrature rule with a Jacobi weight
 
 
 *source:*
-[FastTransforms/src/fejer.jl:12](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/fejer.jl#L12)
+[FastTransforms/src/fejer.jl:12](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/fejer.jl#L12)
 
 ---
 
@@ -215,7 +215,7 @@ Compute weights of Fejer's first quadrature rule with a Jacobi weight.
 
 
 *source:*
-[FastTransforms/src/fejer.jl:21](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/fejer.jl#L21)
+[FastTransforms/src/fejer.jl:21](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/fejer.jl#L21)
 
 ---
 
@@ -225,7 +225,7 @@ Compute weights of Fejer's second quadrature rule with a Jacobi weight.
 
 
 *source:*
-[FastTransforms/src/fejer.jl:26](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/fejer.jl#L26)
+[FastTransforms/src/fejer.jl:26](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/fejer.jl#L26)
 
 ---
 
@@ -235,7 +235,7 @@ Compute a typed 0.5.
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:12](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L12)
+[FastTransforms/src/specialfunctions.jl:12](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L12)
 
 ---
 
@@ -245,7 +245,7 @@ Compute Jacobi expansion coefficients in Pₙ^(α+1,β) given Jacobi expansion c
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:432](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L432)
+[FastTransforms/src/specialfunctions.jl:432](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L432)
 
 ---
 
@@ -255,7 +255,7 @@ Compute Jacobi expansion coefficients in Pₙ^(α+1,α+1) given Jacobi expansion
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:454](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L454)
+[FastTransforms/src/specialfunctions.jl:454](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L454)
 
 ---
 
@@ -265,7 +265,7 @@ Compute Jacobi expansion coefficients in Pₙ^(α,β+1) given Jacobi expansion c
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:443](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L443)
+[FastTransforms/src/specialfunctions.jl:443](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L443)
 
 ---
 
@@ -275,7 +275,7 @@ Pochhammer symbol (x)_n = Γ(x+n)/Γ(x) for the rising factorial.
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:32](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L32)
+[FastTransforms/src/specialfunctions.jl:32](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L32)
 
 ---
 
@@ -285,7 +285,7 @@ Stirling series for Γ(z).
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:63](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L63)
+[FastTransforms/src/specialfunctions.jl:63](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L63)
 
 ---
 
@@ -295,7 +295,7 @@ Compute a typed 2.
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:20](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L20)
+[FastTransforms/src/specialfunctions.jl:20](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L20)
 
 ---
 
@@ -307,7 +307,7 @@ For 64-bit floating-point arithmetic, the Lambda function uses the asymptotic se
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:147](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L147)
+[FastTransforms/src/specialfunctions.jl:147](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L147)
 
 ---
 
@@ -317,7 +317,7 @@ The Lambda function Λ(z) = Γ(z+½)/Γ(z+1) for the ratio of gamma functions.
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:141](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L141)
+[FastTransforms/src/specialfunctions.jl:141](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L141)
 
 ---
 
@@ -327,7 +327,7 @@ The Lambda function Λ(z,λ₁,λ₂) = Γ(z+λ₁)/Γ(z+λ₂) for the ratio of
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:160](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L160)
+[FastTransforms/src/specialfunctions.jl:160](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L160)
 
 ---
 
@@ -337,7 +337,7 @@ The Kronecker δ function.
 
 
 *source:*
-[FastTransforms/src/specialfunctions.jl:26](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/specialfunctions.jl#L26)
+[FastTransforms/src/specialfunctions.jl:26](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/specialfunctions.jl#L26)
 
 ---
 
@@ -351,5 +351,5 @@ where the Hankel matrix `H` is non-negative definite. This allows a Cholesky dec
 
 
 *source:*
-[FastTransforms/src/toeplitzhankel.jl:8](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/b6e4fdf35a2160f0d9948e6565fd69c2177f4d82/src/toeplitzhankel.jl#L8)
+[FastTransforms/src/toeplitzhankel.jl:8](https://github.com/MikaelSlevinsky/FastTransforms.jl/tree/cf090e3655e959c23556d1af681c3378ab9c3081/src/toeplitzhankel.jl#L8)
 

src/toeplitzhankel.jl

@@ -156,31 +156,26 @@ function jac2jacTH{S}(::Type{S},n,α,β,γ,δ)
     T,H,DL,DR
 end
 
-th_leg2chebplan{S}(::Type{S},n)=ToeplitzHankelPlan(leg2chebTH(S,n)...,ones(S,n))
-th_cheb2legplan{S}(::Type{S},n)=ToeplitzHankelPlan(cheb2legTH(S,n)...)
-th_leg2chebuplan{S}(::Type{S},n)=ToeplitzHankelPlan(leg2chebuTH(S,n)...,1:n,ones(S,n))
-th_ultra2ultraplan{S}(::Type{S},n,λ₁,λ₂)=ToeplitzHankelPlan(ultra2ultraTH(S,n,λ₁,λ₂)...)
-th_jac2jacplan{S}(::Type{S},n,α,β,γ,δ)=ToeplitzHankelPlan(jac2jacTH(S,n,α,β,γ,δ)...)
-
-
 immutable ChebyshevToLegendrePlan{TH}
     toeplitzhankel::TH
 end
 
-ChebyshevToLegendrePlan{S}(::Type{S},n) = ChebyshevToLegendrePlan(th_cheb2legplan(S,n))
-
-
 function *(P::ChebyshevToLegendrePlan,v::AbstractVector)
     w = zero(v)
     S,n = eltype(v),length(v)
     w[1:2:end] = -one(S)./(one(S):two(S):n)./(-one(S):two(S):n-two(S))
     [dot(w,v);P.toeplitzhankel*view(v,2:n)]
 end
 
+th_leg2chebplan{S}(::Type{S},n)=ToeplitzHankelPlan(leg2chebTH(S,n)...,ones(S,n))
+th_cheb2legplan{S}(::Type{S},n)=ChebyshevToLegendrePlan(ToeplitzHankelPlan(cheb2legTH(S,n)...))
+th_leg2chebuplan{S}(::Type{S},n)=ToeplitzHankelPlan(leg2chebuTH(S,n)...,1:n,ones(S,n))
+th_ultra2ultraplan{S}(::Type{S},n,λ₁,λ₂)=ToeplitzHankelPlan(ultra2ultraTH(S,n,λ₁,λ₂)...)
+th_jac2jacplan{S}(::Type{S},n,α,β,γ,δ)=ToeplitzHankelPlan(jac2jacTH(S,n,α,β,γ,δ)...)
 
-th_leg2cheb(v)=th_leg2chebplan(eltype(v),length(v))*v
-th_cheb2leg(v) = ChebyshevToLegendrePlan(eltype(v),length(v))*v
 
+th_leg2cheb(v)=th_leg2chebplan(eltype(v),length(v))*v
+th_cheb2leg(v) = th_cheb2legplan(eltype(v),length(v))*v
 th_leg2chebu(v)=th_leg2chebuplan(eltype(v),length(v))*v
 th_ultra2ultra(v,λ₁,λ₂)=th_ultra2ultraplan(eltype(v),length(v),λ₁,λ₂)*v
 th_jac2jac(v,α,β,γ,δ)=th_jac2jacplan(eltype(v),length(v),α,β,γ,δ)*v