Move SymmetricEigen to ApproxFunBase (#248)

* Move symeigen to AFB

* bump AFB to v0.8.23

* bump version to 0.6.28

Author: jishnub

Project.toml

@@ -20,7 +20,7 @@ StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 
 [compat]
-ApproxFunBase = "0.8.22"
+ApproxFunBase = "0.8.23"
 ApproxFunBaseTest = "0.1"
 Aqua = "0.5"
 BandedMatrices = "0.16, 0.17"

src/ApproxFunOrthogonalPolynomials.jl

@@ -44,7 +44,8 @@ import ApproxFunBase: Fun, SubSpace, WeightSpace, NoSpace, HeavisideSpace,
                     ℓ⁰, recα, recβ, recγ, ℵ₀, ∞, RectDomain,
                     assert_integer, supportsinplacetransform, ContinuousSpace, SpecialEvalPtType,
                     isleftendpoint, isrightendpoint, evaluation_point, boundaryfn, ldiffbc, rdiffbc,
-                    LeftEndPoint, RightEndPoint, normalizedspace, promotedomainspace, eigs
+                    LeftEndPoint, RightEndPoint, normalizedspace, promotedomainspace,
+                    bandmatrices_eigen, SymmetricEigensystem, SkewSymmetricEigensystem
 
 import DomainSets: Domain, indomain, UnionDomain, FullSpace, Point,
             Interval, ChebyshevInterval, boundary, rightendpoint, leftendpoint,
@@ -84,6 +85,8 @@ using StaticArrays: SVector
 
 import LinearAlgebra: isdiag, eigvals, eigen
 
+export bandmatrices_eigen, SymmetricEigensystem, SkewSymmetricEigensystem
+
 points(d::IntervalOrSegmentDomain{T},n::Integer) where {T} =
     _maybefromcanonical(d, chebyshevpoints(float(real(eltype(T))), n))  # eltype to handle point
 _maybefromcanonical(d, pts) = fromcanonical.(Ref(d), pts)
@@ -140,7 +143,6 @@ include("Spaces/Spaces.jl")
 include("roots.jl")
 include("specialfunctions.jl")
 include("fastops.jl")
-include("symeigen.jl")
 include("show.jl")
 
 end