Move SparseArrays and Statistics to extensions

Author: jishnub

Project.toml

@@ -1,6 +1,6 @@
 name = "ApproxFunBase"
 uuid = "fbd15aa5-315a-5a7d-a8a4-24992e37be05"
-version = "0.8.43"
+version = "0.8.44"
 
 [deps]
 AbstractFFTs = "621f4979-c628-5d54-868e-fcf4e3e8185c"
@@ -26,6 +26,14 @@ StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
+[weakdeps]
+SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
+Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
+
+[extensions]
+ApproxFunBaseSparseArraysExt = "SparseArrays"
+ApproxFunBaseStatisticsExt = "Statistics"
+
 [compat]
 AbstractFFTs = "0.5, 1"
 Aqua = "0.6"
@@ -63,6 +71,7 @@ Infinities = "e1ba4f0e-776d-440f-acd9-e1d2e9742647"
 OddEvenIntegers = "8d37c425-f37a-4ca2-9b9d-a61bc06559d2"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 Static = "aedffcd0-7271-4cad-89d0-dc628f76c6d3"
+Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 
 [targets]
-test = ["Aqua", "Random", "FastTransforms", "FastGaussQuadrature", "HalfIntegers", "Infinities", "OddEvenIntegers", "Static"]
+test = ["Aqua", "Random", "FastTransforms", "FastGaussQuadrature", "HalfIntegers", "Infinities", "OddEvenIntegers", "Static", "Statistics"]

ext/ApproxFunBaseSparseArraysExt.jl

@@ -0,0 +1,36 @@
+module ApproxFunBaseSparseArraysExt
+
+using SparseArrays
+import SparseArrays: blockdiag
+using ApproxFunBase
+using ApproxFunBase: promote_eltypeof
+
+##TODO: unify with other blockdiag
+function blockdiag(d1::AbstractVector{T}, d2::AbstractVector{T}) where T<:Operator
+    if isempty(d1) && isempty(d2)
+        error("Empty blockdiag")
+    end
+    if isempty(d1)
+        TT=promote_eltypeof(d2)
+    elseif isempty(d2)
+        TT=promote_eltypeof(d1)
+    else
+        TT=promote_type(promote_eltypeof(d1),
+                        promote_eltypeof(d2))
+    end
+
+      D=zeros(Operator{TT},length(d1)+length(d2),2)
+      D[1:length(d1),1]=d1
+      D[length(d1)+1:end,2]=d2
+      D
+end
+
+function blockdiag(a::Operator, b::Operator)
+	blockdiag(Operator{promote_type(eltype(a),eltype(b))}[a],
+			Operator{promote_type(eltype(a),eltype(b))}[b])
+end
+
+blockdiag(A::PlusOperator) = mapreduce(blockdiag, +, A.ops)
+blockdiag(A::TimesOperator) = mapreduce(blockdiag, .*, A.ops)
+
+end

ext/ApproxFunBaseStatisticsExt.jl

@@ -0,0 +1,9 @@
+module ApproxFunBaseStatisticsExt
+
+import Statistics: mean
+using ApproxFunBase
+using ApproxFunBase: IntervalOrSegment
+
+mean(d::IntervalOrSegment) = ApproxFunBase.mean(d)
+
+end

src/ApproxFunBase.jl

@@ -1,7 +1,7 @@
 module ApproxFunBase
 using BlockArrays, BandedMatrices, BlockBandedMatrices, DomainSets,
               IntervalSets, SpecialFunctions, AbstractFFTs, FFTW,
-              SpecialFunctions, DSP, DualNumbers, LinearAlgebra, SparseArrays,
+              SpecialFunctions, DSP, DualNumbers, LinearAlgebra,
               LowRankApprox, FillArrays, InfiniteArrays, InfiniteLinearAlgebra
 
 import Calculus
@@ -39,8 +39,6 @@ import Base: values, convert, getindex, setindex!, *, +, -, ==, <, <=, >, |, !,
 import Base.Broadcast: BroadcastStyle, Broadcasted, AbstractArrayStyle,
               broadcastable, DefaultArrayStyle, broadcasted
 
-import Statistics: mean
-
 import Combinatorics: multiexponents
 
 import LinearAlgebra: BlasInt, BlasFloat, norm, ldiv!, mul!, det, cross,
@@ -49,10 +47,6 @@ import LinearAlgebra: BlasInt, BlasFloat, norm, ldiv!, mul!, det, cross,
               nullspace, Hermitian, Symmetric, adjoint, transpose, char_uplo,
               axpy!, eigvals
 
-import SparseArrays: blockdiag
-
-# import Arpack: eigs
-
 # we need to import all special functions to use Calculus.symbolic_derivatives_1arg
 # we can't do importall Base as we replace some Base definitions
 import SpecialFunctions: airy, besselh,
@@ -155,4 +149,9 @@ include("testing.jl")
 include("specialfunctions.jl")
 include("show.jl")
 
+if !isdefined(Base, :get_extension)
+    include("../ext/ApproxFunBaseSparseArraysExt.jl")
+    include("../ext/ApproxFunBaseStatisticsExt.jl")
+end
+
 end #module

src/Domains/Segment.jl

@@ -68,6 +68,7 @@ issubset(a::Segment,b::Segment) = leftendpoint(a)∈b && rightendpoint(a)∈b
 arclength(d::AbstractInterval) = width(d)
 arclength(d::Segment) = norm(complexlength(d))
 complexlength(d::IntervalOrSegment) = rightendpoint(d)-leftendpoint(d)
+# ApproxFunBase.mean != Statistics.mean, as the latter is defined in the Statistics extension
 mean(d::IntervalOrSegment) = (rightendpoint(d)+leftendpoint(d))/2
 angle(d::IntervalOrSegment) = angle(complexlength(d))
 sign(d::IntervalOrSegment) = sign(complexlength(d))

src/Operators/Operator.jl

@@ -329,13 +329,6 @@ istril(A::Operator) = bandwidth(A, 2) <= 0
 
 include("SubOperator.jl")
 
-
-#
-# sparse(B::Operator,n::Integer)=sparse(BandedMatrix(B,n))
-# sparse(B::Operator,n::AbstractRange,m::AbstractRange)=sparse(BandedMatrix(B,n,m))
-# sparse(B::Operator,n::Colon,m::AbstractRange)=sparse(BandedMatrix(B,n,m))
-# sparse(B::Operator,n::AbstractRange,m::Colon)=sparse(BandedMatrix(B,n,m))
-
 ## getindex
 
 

src/Operators/systems.jl

@@ -33,29 +33,6 @@ function diagm_container(size, kv::Pair{<:Integer,<:AbstractVector{<:Operator}}.
     zeros(Operator{T}, n, n)
 end
 
-##TODO: unify with other blockdiag
-function blockdiag(d1::AbstractVector{T},d2::AbstractVector{T}) where T<:Operator
-    if isempty(d1)&&isempty(d2)
-        error("Empty blockdiag")
-    end
-    if isempty(d1)
-        TT=promote_eltypeof(d2)
-    elseif isempty(d2)
-        TT=promote_eltypeof(d1)
-    else
-        TT=promote_type(promote_eltypeof(d1),
-                        promote_eltypeof(d2))
-    end
-
-      D=zeros(Operator{TT},length(d1)+length(d2),2)
-      D[1:length(d1),1]=d1
-      D[length(d1)+1:end,2]=d2
-      D
-end
-
-blockdiag(a::Operator,b::Operator) = blockdiag(Operator{promote_type(eltype(a),eltype(b))}[a],
-                                                Operator{promote_type(eltype(a),eltype(b))}[b])
-
 ## broadcase
 
 broadcast(::typeof(*),A::AbstractArray{N},D::Operator) where {N<:Number} =

src/Spaces/ProductSpaceOperators.jl

@@ -55,9 +55,6 @@ end
 continuity(d::UnionDomain,k) = continuity(Space(d),k)
 continuity(d) = continuity(d,0)
 
-blockdiag(A::PlusOperator) = mapreduce(blockdiag, +, A.ops)
-blockdiag(A::TimesOperator) = mapreduce(blockdiag, .*, A.ops)
-
 # TODO: general wrappers
 
 Evaluation(S::SumSpace,x,order) =