Refactor into multiple packages (#658)

* move out files

* move out files

* Update ApproxFun.jl

* move tests back here

* using works

* Fractional operator types

* work on tests

* tests pass

* ApproxFunSingularities

* move nullspace test

* Project.toml, NEWS

* export sample

* `@time` to prevent timeout

* More `@time`

* Update FractionalTest.jl

* Skip fractionaltests for now

Author: dlfivefifty

.DS_Store

@@ -3,6 +3,24 @@
 
 ### notes on release changes, ongoing development, and future planned work
 
+#### 0.11
+ - Split into ApproxFunBase.jl, ApproxFunFourier.jl, ApproxFunOrthogonalPolynomials.jl, and ApproxFunSingularities.jl
+
+#### 0.10.4
+ - Update for latest BandedMatrices.jl
+ - Support for symmetric discretizations
+ - Bug fixes
+
+#### 0.10.3
+ - Bug fixes
+ - Update for latest BlockBandedMatrices.jl
+
+#### 0.10.2
+ - Bug fixes for special functions on large intervals
+
+#### 0.10.1 
+ - Update for latest BlockBandedMatrices.jl
+ 
 #### 0.10.0
  - Use DomainSets.jl for representing domains
  - Use InfiniteArrays.jl for representing block sizes

NEWS.md

@@ -0,0 +1,62 @@
+name = "ApproxFun"
+uuid = "28f2ccd6-bb30-5033-b560-165f7b14dc2f"
+version = "0.11"
+
+[deps]
+AbstractFFTs = "621f4979-c628-5d54-868e-fcf4e3e8185c"
+ApproxFunBase = "fbd15aa5-315a-5a7d-a8a4-24992e37be05"
+ApproxFunFourier = "59844689-9c9d-51bf-9583-5b794ec66d30"
+ApproxFunOrthogonalPolynomials = "b70543e2-c0d9-56b8-a290-0d4d6d4de211"
+ApproxFunSingularities = "f8fcb915-6b99-5be2-b79a-d6dbef8e6e7e"
+BandedMatrices = "aae01518-5342-5314-be14-df237901396f"
+BlockArrays = "8e7c35d0-a365-5155-bbbb-fb81a777f24e"
+BlockBandedMatrices = "ffab5731-97b5-5995-9138-79e8c1846df0"
+Calculus = "49dc2e85-a5d0-5ad3-a950-438e2897f1b9"
+DSP = "717857b8-e6f2-59f4-9121-6e50c889abd2"
+DomainSets = "5b8099bc-c8ec-5219-889f-1d9e522a28bf"
+DualNumbers = "fa6b7ba4-c1ee-5f82-b5fc-ecf0adba8f74"
+FFTW = "7a1cc6ca-52ef-59f5-83cd-3a7055c09341"
+FastGaussQuadrature = "442a2c76-b920-505d-bb47-c5924d526838"
+FastTransforms = "057dd010-8810-581a-b7be-e3fc3b93f78c"
+FillArrays = "1a297f60-69ca-5386-bcde-b61e274b549b"
+InfiniteArrays = "4858937d-0d70-526a-a4dd-2d5cb5dd786c"
+IntervalSets = "8197267c-284f-5f27-9208-e0e47529a953"
+LazyArrays = "5078a376-72f3-5289-bfd5-ec5146d43c02"
+LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
+LowRankApprox = "898213cb-b102-5a47-900c-97e73b919f73"
+Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
+RecipesBase = "3cdcf5f2-1ef4-517c-9805-6587b60abb01"
+Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
+SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
+SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
+StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
+Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
+ToeplitzMatrices = "c751599d-da0a-543b-9d20-d0a503d91d24"
+
+[compat]
+AbstractFFTs = "≥ 0.3.1"
+ApproxFunBase = "≥ 0.0.3"
+ApproxFunFourier = "≥ 0.0.1"
+ApproxFunSingularities = "≥ 0.0.1"
+ApproxFunOrthogonalPolynomails = "≥ 0.0.2"
+BandedMatrices = "≥ 0.9.0"
+BlockArrays = "≥ 0.8.0"
+BlockBandedMatrices = "≥ 0.4.0"
+Calculus = "≥ 0.1.15"
+DSP = "≥ 0.5.1"
+DomainSets = "≥ 0.0.1"
+DualNumbers = "≥ 0.4.0"
+FFTW = "≥ 0.2.4"
+FastGaussQuadrature = "≥ 0.3.2"
+FastTransforms = "≥ 0.4.2"
+FillArrays = "≥ 0.5.0"
+InfiniteArrays = "≥ 0.0.3"
+IntervalSets = "≥ 0.3.1"
+LazyArrays = "≥ 0.8.0"
+LowRankApprox = "≥ 0.1.4"
+RecipesBase = "≥ 0.5.0"
+Reexport = "≥ 0.2.0"
+SpecialFunctions = "≥ 0.7.0"
+StaticArrays = "≥ 0.8.3"
+ToeplitzMatrices = "≥ 0.4.0"
+julia = "≥ 0.7.0"

Project.toml

@@ -1,12 +1,60 @@
 __precompile__()
 
 module ApproxFun
-    using Base, RecipesBase, FastGaussQuadrature, FastTransforms, DualNumbers,
-            BlockArrays, BandedMatrices, BlockBandedMatrices, DomainSets, IntervalSets,
-            SpecialFunctions, AbstractFFTs, FFTW, SpecialFunctions, DSP,
-            LinearAlgebra, LowRankApprox, SparseArrays, FillArrays, InfiniteArrays #, Arpack
-    import StaticArrays, ToeplitzMatrices, Calculus
-
+    using Base, Reexport, BlockArrays, BandedMatrices, BlockBandedMatrices, DomainSets, IntervalSets,
+            SpecialFunctions, AbstractFFTs, FFTW, SpecialFunctions, DSP, DualNumbers, FastTransforms,
+            LinearAlgebra, SparseArrays, LowRankApprox, FillArrays, InfiniteArrays, RecipesBase,
+            FFTW, AbstractFFTs #, Arpack
+
+import Calculus
+
+@reexport using ApproxFunBase    
+@reexport using ApproxFunFourier
+@reexport using ApproxFunOrthogonalPolynomials
+@reexport using ApproxFunSingularities
+
+import ApproxFunBase: normalize!, flipsign, FiniteRange, Fun, MatrixFun, UnsetSpace, VFun, RowVector,
+                    UnivariateSpace, AmbiguousSpace, SumSpace, SubSpace, WeightSpace, NoSpace, Space,
+                    HeavisideSpace, PointSpace,
+                    IntervalOrSegment, RaggedMatrix, AlmostBandedMatrix,
+                    AnyDomain, ZeroSpace, ArraySpace, TrivialInterlacer, BlockInterlacer, 
+                    AbstractTransformPlan, TransformPlan, ITransformPlan,
+                    ConcreteConversion, ConcreteMultiplication, ConcreteDerivative, ConcreteIntegral, CalculusOperator,
+                    ConcreteVolterra, Volterra, VolterraWrapper,
+                    MultiplicationWrapper, ConversionWrapper, DerivativeWrapper, Evaluation, EvaluationWrapper,
+                    Conversion, defaultConversion, defaultcoefficients, default_Fun, Multiplication, Derivative, Integral, bandwidths, 
+                    ConcreteEvaluation, ConcreteDefiniteLineIntegral, ConcreteDefiniteIntegral, ConcreteIntegral,
+                    DefiniteLineIntegral, DefiniteIntegral, ConcreteDefiniteIntegral, ConcreteDefiniteLineIntegral, IntegralWrapper,
+                    ReverseOrientation, ReverseOrientationWrapper, ReverseWrapper, Reverse, NegateEven, 
+                    Dirichlet, ConcreteDirichlet, DirichletWrapper,
+                    TridiagonalOperator, SubOperator, Space, @containsconstants, spacescompatible,
+                    hasfasttransform, canonicalspace, domain, setdomain, prectype, domainscompatible, 
+                    plan_transform, plan_itransform, plan_transform!, plan_itransform!, transform, itransform, hasfasttransform, 
+                    CanonicalTransformPlan, ICanonicalTransformPlan,
+                    Integral, 
+                    domainspace, rangespace, boundary, 
+                    union_rule, conversion_rule, maxspace_rule, conversion_type, maxspace, hasconversion, points, 
+                    rdirichlet, ldirichlet, lneumann, rneumann, ivp, bvp, 
+                    linesum, differentiate, integrate, linebilinearform, bilinearform, 
+                    UnsetNumber, coefficienttimes, subspace_coefficients, sumspacecoefficients, specialfunctionnormalizationpoint,
+                    Segment, IntervalOrSegmentDomain, PiecewiseSegment, isambiguous, Vec, eps, isperiodic,
+                    arclength, complexlength,
+                    invfromcanonicalD, fromcanonical, tocanonical, fromcanonicalD, tocanonicalD, canonicaldomain, setcanonicaldomain, mappoint,
+                    reverseorientation, checkpoints, evaluate, mul_coefficients, coefficients, coefficientmatrix, isconvertible,
+                    clenshaw, ClenshawPlan, sineshaw,
+                    toeplitz_getindex, toeplitz_axpy!, sym_toeplitz_axpy!, hankel_axpy!, ToeplitzOperator, SymToeplitzOperator, hankel_getindex, 
+                    SpaceOperator, ZeroOperator, InterlaceOperator,
+                    interlace!, reverseeven!, negateeven!, cfstype, pad!, alternatesign!, mobius,
+                    extremal_args, hesseneigvals, chebyshev_clenshaw, recA, recB, recC, roots,splitatroots,
+                    chebmult_getindex, intpow, alternatingsum,
+                    domaintype, diagindshift, rangetype, weight, isapproxinteger, default_Dirichlet, scal!, dotu,
+                    components, promoterangespace, promotedomainspace, choosedomainspace,
+                    block, blockstart, blockstop, blocklengths, isblockbanded, pointscompatible,
+                    AbstractProductSpace, MultivariateFun, BivariateSpace, 
+                    @wrapperstructure, @wrapperspaces, @wrapper, @calculus_operator, resizedata!, slnorm,
+                    sample
+
+import ApproxFunOrthogonalPolynomials: order
 
 import DomainSets: Domain, indomain, UnionDomain, ProductDomain, FullSpace, Point, elements, DifferenceDomain,
             Interval, ChebyshevInterval, boundary, ∂, rightendpoint, leftendpoint,
@@ -16,7 +64,6 @@ import AbstractFFTs: Plan, fft, ifft
 import FFTW: plan_r2r!, fftwNumber, REDFT10, REDFT01, REDFT00, RODFT00, R2HC, HC2R,
                 r2r!, r2r,  plan_fft, plan_ifft, plan_ifft!, plan_fft!
 
-
 import Base: values, convert, getindex, setindex!, *, +, -, ==, <, <=, >, |, !, !=, eltype, iterate,
                 >=, /, ^, \, ∪, transpose, size, reindex, tail, broadcast, broadcast!, copyto!, copy, to_index, (:),
                 similar, map, vcat, hcat, hvcat, show, summary, stride, sum, cumsum, sign, imag, conj, inv,
@@ -32,7 +79,7 @@ import Base: values, convert, getindex, setindex!, *, +, -, ==, <, <=, >, |, !,
 import Base.Broadcast: BroadcastStyle, Broadcasted, AbstractArrayStyle, broadcastable,
                         DefaultArrayStyle, broadcasted
 
-import Statistics: mean
+
 
 import LinearAlgebra: BlasInt, BlasFloat, norm, ldiv!, mul!, det, eigvals, dot, cross,
                         qr, qr!, rank, isdiag, istril, istriu, issymmetric, ishermitian,
@@ -64,7 +111,6 @@ import BlockArrays: nblocks, blocksize, global2blockindex, globalrange, BlockSiz
 
 import BandedMatrices: bandrange, bandshift,
                         inbands_getindex, inbands_setindex!, bandwidth, AbstractBandedMatrix,
-                        flipsign,
                         colstart, colstop, colrange, rowstart, rowstop, rowrange,
                         bandwidths, _BandedMatrix, BandedMatrix
 
@@ -82,11 +128,24 @@ import FillArrays: AbstractFill, getindex_value
 import LazyArrays: cache
 import InfiniteArrays: Infinity, InfRanges, AbstractInfUnitRange, OneToInf
 
+"""
+`Curve` Represents a domain defined by the image of a Fun.  Example
+usage would be
+
+```julia
+x=Fun(1..2)
+Curve(exp(im*x))  # represents an arc
+```
+"""
+const Curve{S,T} = Union{IntervalCurve{S,T},PeriodicCurve{S,T}}
+Curve(f::Fun{<:Space{<:PeriodicDomain}}) = PeriodicCurve(f)
+
+#TODO: Make type stable
+Curve(f::Fun{<:Space{<:ChebyshevInterval}}) = IntervalCurve(f) 
 
+export Curve
 
 
-# convenience for 1-d block ranges
-const BlockRange1 = BlockRange{1,Tuple{UnitRange{Int}}}
 
 import Base: view
 
@@ -95,42 +154,16 @@ import StaticArrays: StaticArray, SVector
 
 import IntervalSets: (..), endpoints
 
-const Vec{d,T} = SVector{d,T}
-
-export pad!, pad, chop!, sample,
-       complexroots, roots, svfft, isvfft,
-       reverseorientation, jumplocations
 
 ##Testing
 export bisectioninv
 
-export .., Interval, ChebyshevInterval, leftendpoint, rightendpoint, endpoints
-
-
-
-include("LinearAlgebra/LinearAlgebra.jl")
-
-
-include("Fun/Fun.jl")
-
-
-include("Domains/Domains.jl")
-include("Multivariate/Multivariate.jl")
-include("Operators/Operator.jl")
-
-include("Spaces/Spaces.jl")
-
-
-
 
 ## Further extra features
 
-include("PDE/PDE.jl")
-include("Caching/caching.jl")
 include("Extras/Extras.jl")
 include("Plot/Plot.jl")
 include("docs.jl")
-include("testing.jl")
 
 #include("precompile.jl")
 #_precompile_()