Work on Cone transform

Author: dlfivefifty

src/Cone/Cone.jl

@@ -182,8 +182,11 @@ points(::Cone, n) = conemap.(points(rectspace(DuffyCone()), n))
 checkpoints(::Cone) = conemap.(checkpoints(rectspace(DuffyCone())))
 
 
-plan_transform(S::DuffyCone, n::AbstractVector) = 
-    TransformPlan(S, plan_transform(rectspace(S),n), Val{false})
+function plan_transform(S::DuffyCone, n::AbstractVector) 
+    rs = rectspace(S)
+    P = TransformPlan(sp,((plan_transform(sp.spaces[1],T,N),N), (plan_transform(sp.spaces[2],T,M),M)), Val{false})
+    TransformPlan(S, P, Val{false})
+end
 plan_itransform(S::DuffyCone, n::AbstractVector) = 
     ITransformPlan(S, plan_itransform(rectspace(S),n), Val{false})
 

test/test_cone.jl

@@ -1,4 +1,4 @@
-using ApproxFun, MultivariateOrthogonalPolynomials, Test
+using ApproxFun, MultivariateOrthogonalPolynomials, StaticArrays, Test
 import MultivariateOrthogonalPolynomials: rectspace, totensor, duffy2legendreconic!, legendre2duffyconic!, c_plan_rottriangle, plan_transform
 
 @testset "Conic" begin
@@ -70,6 +70,8 @@ end
     @testset "rectspace" begin
         rs = rectspace(DuffyCone())
         @test points(rs,10) isa Vector{SVector{3,Float64}}
+        @test_broken @inferred(checkpoints(rs))
+        @test checkpoints(rs) isa Vector{SVector{3,Float64}}
     end
 
     @testset "DuffyCone" begin