Move shared code here, export spaces

Author: dlfivefifty

REQUIRE

@@ -1,7 +1,6 @@
 julia 0.7
 FastGaussQuadrature 0.3.2
 FastTransforms 0.4.2
-RecipesBase 0.5
 DualNumbers 0.4
 ToeplitzMatrices 0.4
 Calculus 0.1.15

src/ApproxFunBase.jl

@@ -96,9 +96,6 @@ export pad!, pad, chop!, sample,
        complexroots, roots, svfft, isvfft,
        reverseorientation, jumplocations
 
-##Testing
-export bisectioninv
-
 export .., Interval, ChebyshevInterval, leftendpoint, rightendpoint, endpoints, cache
 
 

src/Domains/Domains.jl

@@ -39,3 +39,29 @@ end
 # multivariate domainxs
 
 include("multivariate.jl")
+
+function affine_setdiff(d::Domain, ptsin::UnionDomain)    
+    pts=Number.(elements(ptsin))
+    isempty(pts) && return d
+    tol=sqrt(eps(arclength(d)))
+    da=leftendpoint(d)
+    isapprox(da,pts[1];atol=tol) && popfirst!(pts)
+    isempty(pts) && return d
+    db=rightendpoint(d)
+    isapprox(db,pts[end];atol=tol) && pop!(pts)
+
+    sort!(pts)
+    leftendpoint(d) > rightendpoint(d) && reverse!(pts)
+    filter!(p->p ∈ d,pts)
+
+    isempty(pts) && return d
+    length(pts) == 1 && return d \ pts[1]
+
+    ret = Array{Domain}(undef, length(pts)+1)
+    ret[1] = Domain(leftendpoint(d) .. pts[1])
+    for k = 2:length(pts)
+        ret[k] = Domain(pts[k-1]..pts[k])
+    end
+    ret[end] = Domain(pts[end] .. rightendpoint(d))
+    UnionDomain(ret)
+end

src/Spaces/DiracSpace.jl

@@ -1,4 +1,5 @@
-export DiracDelta, KroneckerDelta
+export DiracDelta, KroneckerDelta, DiracSpace, PointSpace
+
 for TYP in (:DiracSpace,:PointSpace)
     @eval begin
         struct $TYP{T,D,R} <: Space{D,R}

src/specialfunctions.jl

@@ -537,6 +537,8 @@ end
 extrema(f::Fun{PiecewiseSpace{SV,DD,RR},T}) where {SV,DD<:UnionDomain,RR<:Real,T<:Real} =
     mapreduce(extrema,(x,y)->extrema([x...;y...]),components(f))
 
+function complexroots end
+
 function roots(f::Fun{P}) where P<:PiecewiseSpace
     rts=mapreduce(roots,vcat,components(f))
     k=1