Merge pull request #2 from JuliaApproximation/dl/domainsets002_2

DomainSets 0.0.2 support

Author: dlfivefifty

Project.toml

@@ -1,6 +1,6 @@
 name = "ApproxFunBase"
 uuid = "fbd15aa5-315a-5a7d-a8a4-24992e37be05"
-version = "0.0.4"
+version = "0.0.5"
 
 [deps]
 AbstractFFTs = "621f4979-c628-5d54-868e-fcf4e3e8185c"
@@ -35,7 +35,7 @@ BlockArrays = "0.8.0"
 BlockBandedMatrices = "0.4.0"
 Calculus = "≥ 0.1.15"
 DSP = "≥ 0.5.1"
-DomainSets = "0.0.1"
+DomainSets = "0.0.2"
 DualNumbers = "0.6.2"
 FFTW = "≥ 0.2.4"
 FastGaussQuadrature = "≥ 0.3.2"

src/ApproxFunBase.jl

@@ -8,7 +8,7 @@ import StaticArrays, Calculus
 
 import DomainSets: Domain, indomain, UnionDomain, ProductDomain, FullSpace, Point, elements, DifferenceDomain,
             Interval, ChebyshevInterval, boundary, ∂, rightendpoint, leftendpoint,
-            dimension
+            dimension, WrappedDomain
 
 
 

src/Domains/Domains.jl

@@ -40,8 +40,13 @@ end
 
 include("multivariate.jl")
 
-function affine_setdiff(d::Domain, ptsin::UnionDomain)    
-    pts=Number.(elements(ptsin))
+affine_setdiff(d::Domain, ptsin::UnionDomain) = 
+    _affine_setdiff(d, Number.(elements(ptsin)))
+
+affine_setdiff(d::Domain, ptsin::WrappedDomain{<:AbstractVector}) = 
+    _affine_setdiff(d, ptsin.domain)
+
+function _affine_setdiff(d::Domain, pts)    
     isempty(pts) && return d
     tol=sqrt(eps(arclength(d)))
     da=leftendpoint(d)