restore Clenshaw for AbstractQuasiVector

Author: dlfivefifty

src/ClassicalOrthogonalPolynomials.jl

@@ -39,8 +39,9 @@ import ContinuumArrays: Basis, Weight, basis_axes, @simplify, Identity, Abstract
     AffineQuasiVector, AffineMap, AbstractWeightLayout, AbstractWeightedBasisLayout, WeightedBasisLayout, WeightedBasisLayouts, demap, AbstractBasisLayout, BasisLayout,
     checkpoints, weight, unweighted, MappedBasisLayouts, sum_layout, invmap, plan_ldiv, layout_broadcasted, MappedBasisLayout, SubBasisLayout, broadcastbasis_layout,
     plan_grid_transform, plan_transform, MAX_PLOT_POINTS, MulPlan, grammatrix, AdjointBasisLayout, grammatrix_layout, plan_transform_layout, _cumsum
-import FastTransforms: Λ, forwardrecurrence, forwardrecurrence!, _forwardrecurrence!, clenshaw, clenshaw!,
-                        _forwardrecurrence_next, _clenshaw_next, check_clenshaw_recurrences, ChebyshevGrid, chebyshevpoints, Plan, ScaledPlan, th_cheb2leg
+import FastTransforms: Λ, ChebyshevGrid, chebyshevpoints, Plan, ScaledPlan, th_cheb2leg
+import RecurrenceRelatioships: forwardrecurrence, forwardrecurrence!, _forwardrecurrence!, clenshaw, clenshaw!,
+                        _forwardrecurrence_next, _clenshaw_next, check_clenshaw_recurrences
 
 import FastGaussQuadrature: jacobimoment
 

src/clenshaw.jl

@@ -1,4 +1,6 @@
 
+
+
 # Assume 1 normalization
 _p0(A) = one(eltype(A))
 
@@ -123,4 +125,10 @@ function _broadcasted_layout_broadcasted_mul(::Tuple{AbstractWeightLayout,Polyno
     Q = OrthogonalPolynomial(w)
     a = (w .* Q) * (Q \ v)
     a .* P
+end
+
+# constructor for Clenshaw
+function Clenshaw(a::AbstractQuasiVector, X::AbstractQuasiMatrix)
+    P,c = arguments(a)
+    Clenshaw(paddeddata(c), recurrencecoefficients(P)..., jacobimatrix(X), _p0(P))
 end