Explicitly use th_cheb2leg to support large dimensions (#126)

Author: dlfivefifty

src/ClassicalOrthogonalPolynomials.jl

@@ -40,7 +40,7 @@ import ContinuumArrays: Basis, Weight, basis, @simplify, Identity, AbstractAffin
     checkpoints, weight, unweighted, MappedBasisLayouts, __sum, invmap, plan_ldiv, layout_broadcasted, MappedBasisLayout, SubBasisLayout, _broadcastbasis,
     plan_transform, plan_grid_transform, MAX_PLOT_POINTS, MulPlan
 import FastTransforms: Λ, forwardrecurrence, forwardrecurrence!, _forwardrecurrence!, clenshaw, clenshaw!,
-                        _forwardrecurrence_next, _clenshaw_next, check_clenshaw_recurrences, ChebyshevGrid, chebyshevpoints, Plan, ScaledPlan
+                        _forwardrecurrence_next, _clenshaw_next, check_clenshaw_recurrences, ChebyshevGrid, chebyshevpoints, Plan, ScaledPlan, th_cheb2leg
 
 import FastGaussQuadrature: jacobimoment
 

src/classical/legendre.jl

@@ -82,7 +82,7 @@ ldiv(P::Legendre{V}, f::AbstractQuasiFill{T,1}) where {T,V} = _op_ldiv(P, f)
 function transform_ldiv(::Legendre{V}, f::Union{AbstractQuasiVector,AbstractQuasiMatrix}) where V
     T = ChebyshevT{V}()
     dat = T \ f
-    pad(cheb2leg(paddeddata(dat)), axes(dat)...)
+    pad(th_cheb2leg(paddeddata(dat)), axes(dat)...)
 end
 
 struct LegendreTransformPlan{T, CHEB2LEG, DCT} <: Plan{T}