Move out singular integrals (#129)

* start working on fast Stieltjes

* Update test_stieltjes.jl

* faststietljes

* move out SIEs

* increase cov

* Update test_normalized.jl

* Update test_normalized.jl

* Update Project.toml

* add UltrasphericalWeight(::LegendreWeight)

Author: dlfivefifty

Project.toml

@@ -1,7 +1,7 @@
 name = "ClassicalOrthogonalPolynomials"
 uuid = "b30e2e7b-c4ee-47da-9d5f-2c5c27239acd"
 authors = ["Sheehan Olver <[email protected]>"]
-version = "0.7.6"
+version = "0.8"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
@@ -25,21 +25,21 @@ QuasiArrays = "c4ea9172-b204-11e9-377d-29865faadc5c"
 SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
 
 [compat]
-ArrayLayouts = "0.8, 1"
+ArrayLayouts = "1"
 BandedMatrices = "0.17.17"
 BlockArrays = "0.16.9"
-BlockBandedMatrices = "0.11.6, 0.12"
+BlockBandedMatrices = "0.12"
 ContinuumArrays = "0.12.4"
 DomainSets = "0.5.6, 0.6"
 FFTW = "1.1"
 FastGaussQuadrature = "0.4.3, 0.5"
 FastTransforms = "0.15.2"
-FillArrays = "0.13, 1"
+FillArrays = "1"
 HypergeometricFunctions = "0.3.4"
 InfiniteArrays = "0.12.11"
 InfiniteLinearAlgebra = "0.6.16"
 IntervalSets = "0.5, 0.6, 0.7"
-LazyArrays = "0.22.18, 1"
+LazyArrays = "1.0.1"
 LazyBandedMatrices = "0.8.5"
 QuasiArrays = "0.9.6"
 SpecialFunctions = "1.0, 2"

src/ClassicalOrthogonalPolynomials.jl

@@ -35,7 +35,7 @@ import QuasiArrays: cardinality, checkindex, QuasiAdjoint, QuasiTranspose, Inclu
 import InfiniteArrays: OneToInf, InfAxes, Infinity, AbstractInfUnitRange, InfiniteCardinal, InfRanges
 import InfiniteLinearAlgebra: chop!, chop, pad, choplength, compatible_resize!
 import ContinuumArrays: Basis, Weight, basis, @simplify, Identity, AbstractAffineQuasiVector, ProjectionFactorization,
-    inbounds_getindex, grid, plotgrid, _plotgrid, _grid, transform_ldiv, TransformFactorization, QInfAxes, broadcastbasis, ExpansionLayout, basismap,
+    grid, plotgrid, _plotgrid, _grid, transform_ldiv, TransformFactorization, QInfAxes, broadcastbasis, ExpansionLayout, basismap,
     AffineQuasiVector, AffineMap, WeightLayout, AbstractWeightedBasisLayout, WeightedBasisLayout, WeightedBasisLayouts, demap, AbstractBasisLayout, BasisLayout,
     checkpoints, weight, unweighted, MappedBasisLayouts, __sum, invmap, plan_ldiv, layout_broadcasted, MappedBasisLayout, SubBasisLayout, _broadcastbasis,
     plan_transform, plan_grid_transform, MAX_PLOT_POINTS, MulPlan
@@ -322,7 +322,6 @@ include("classical/chebyshev.jl")
 include("classical/ultraspherical.jl")
 include("classical/laguerre.jl")
 include("classical/fourier.jl")
-include("stieltjes.jl")
 include("roots.jl")
 
 end # module

src/classical/ultraspherical.jl

@@ -11,6 +11,7 @@ end
 
 UltrasphericalWeight{T}(λ) where T = UltrasphericalWeight{T,typeof(λ)}(λ)
 UltrasphericalWeight(λ) = UltrasphericalWeight{float(typeof(λ)),typeof(λ)}(λ)
+UltrasphericalWeight(::LegendreWeight{T}) where T = UltrasphericalWeight(one(T)/2)
 
 ==(a::UltrasphericalWeight, b::UltrasphericalWeight) = a.λ == b.λ