reorg

Author: dlfivefifty

src/Disk/.DS_Store

@@ -0,0 +1,30 @@
+using ApproxFun
+f = (x,y) -> x*y+cos(y-0.1)+sin(x)+1; ff = Fun((r,θ) -> f(r*cos(θ),r*sin(θ)), (-1..1) × PeriodicSegment());
+ApproxFunBase.coefficientmatrix(ff)
+
+struct ChebyshevDisk{V,T} <: Space{Disk{V},T}
+    domain::Disk{V}
+end
+
+ChebyshevDisk(d::Disk{V}) where V = ChebyshevDisk{V,eltype(V)}(d)
+ChebyshevDisk() = ChebyshevDisk(Disk())
+
+rectspace(_) = Chebyshev() * Laurent()
+
+
+function points(S::ChebyshevDisk, N)
+    pts = points(rectspace(S), N)
+    fromcanonical.(Ref(S.domain), polar.(pts))
+end
+
+plan_transform(S::DuffyTriangle, n::AbstractVector) = 
+    TransformPlan(S, plan_transform(rectspace(S),n), Val{false})
+plan_itransform(S::DuffyTriangle, n::AbstractVector) = 
+    ITransformPlan(S, plan_itransform(rectspace(S),n), Val{false})
+
+
+*(P::TransformPlan{<:Any,<:DuffyTriangle}, v::AbstractArray) = checkerboard(P.plan*v)
+*(P::ITransformPlan{<:Any,<:DuffyTriangle}, v::AbstractArray) = P.plan*icheckerboard(v)
+
+evaluate(cfs::AbstractVector, S::DuffyTriangle, x) = evaluate(cfs, rectspace(S), ipolar(tocanonical(S.domain,x)))
+

src/Disk/Disk.jl

@@ -1,9 +1,7 @@
-__precompile__()
-
 module MultivariateOrthogonalPolynomials
 using Base, RecipesBase, ApproxFun, BandedMatrices, BlockArrays,
     FastTransforms, FastGaussQuadrature, StaticArrays, FillArrays,
-    LinearAlgebra, Libdl, SpecialFunctions, LazyArrays
+    LinearAlgebra, Libdl, SpecialFunctions, LazyArrays, InfiniteArrays
 
 # package code goes here
 import Base: values,getindex,setindex!,*, +, -, ==,<,<=,>,
@@ -15,16 +13,16 @@ import BandedMatrices: inbands_getindex, inbands_setindex!
 import BlockArrays: blocksizes, BlockSizes, getblock, global2blockindex
 
 # ApproxFun general import
-import ApproxFun: BandedMatrix, order, blocksize,
+import ApproxFunBase: BandedMatrix, order, blocksize,
                   linesum,complexlength, BandedBlockBandedMatrix,
-                  real, eps, isapproxinteger, ∞, FiniteRange, DFunction,
+                  real, eps, isapproxinteger, FiniteRange, DFunction,
                   TransformPlan, ITransformPlan
 
 # Domains import
-import ApproxFun: fromcanonical, tocanonical, domainscompatible
+import ApproxFunBase: fromcanonical, tocanonical, domainscompatible
 
 # Operator import
-import ApproxFun:    bandwidths,SpaceOperator, ConversionWrapper, DerivativeWrapper,
+import ApproxFunBase:    bandwidths,SpaceOperator, ConversionWrapper, DerivativeWrapper,
                   rangespace, domainspace, InterlaceOperator,
                   promotedomainspace,  CalculusOperator, interlace, Multiplication,
                    choosedomainspace, SubOperator, ZeroOperator,
@@ -34,7 +32,7 @@ import ApproxFun:    bandwidths,SpaceOperator, ConversionWrapper, DerivativeWrap
 
 
 # Spaces import
-import ApproxFun: PolynomialSpace, ConstantSpace, NoSpace, prectype,
+import ApproxFunBase:   ConstantSpace, NoSpace, prectype,
                     SumSpace,PiecewiseSpace, ArraySpace, @containsconstants,
                     UnsetSpace, canonicalspace, canonicaldomain, domain, evaluate,
                     AnyDomain, plan_transform,plan_itransform,
@@ -51,26 +49,25 @@ import ApproxFun: Domain2d, DirectSumSpace, AbstractProductSpace, factor,
                     Tensorizer, tensorizer, block, blockstart, blockstop, blocklengths,
                     domaintensorizer, rangetensorizer, blockrange, Block, BlockRange1
 
-
-# Jacobi import
-import ApproxFun: jacobip, JacobiSD
-
 # Singularities
-import ApproxFun: WeightSpace, weight
+import ApproxFunBase: WeightSpace, weight
 
 # Vec is for two points
-import ApproxFun: Vec
+import ApproxFunBase: Vec
+
+# Jacobi import
+import ApproxFunOrthogonalPolynomials: jacobip, JacobiSD, PolynomialSpace
 
 
 
-include("c_tri2cheb.jl")
+include("c_transforms.jl")
 
-include("Triangle.jl")
-include("DirichletTriangle.jl")
+include("Triangle/Triangle.jl")
+include("Disk/Disk.jl")
 
 include("clenshaw.jl")
 
-# include("SphericalHarmonics.jl")
+# include("Sphere/SphericalHarmonics.jl")
 
 include("show.jl")
 include("plot.jl")

src/DiskSpace.jl renamed to src/Disk/DiskSpace.jl

@@ -1381,3 +1381,6 @@ function *(A_in::Operator, f::Fun{<:JacobiTriangle})
     v = pad(coefficients(f),size(M,2))
     Fun(rangespace(A),M*v)
 end
+
+
+include("DirichletTriangle.jl")