CompatHelper: bump compat for "HypergeometricFunctions" to "0.3" (#137)

* CompatHelper: bump compat for "HypergeometricFunctions" to "0.3"

* test alternate without dual numbers?

* revert change

* Ok epsilon, you win

Co-authored-by: github-actions[bot] <41898282+github-actions[bot]@users.noreply.github.com>
Co-authored-by: Sheehan Olver <solver@mac.com>
Co-authored-by: Mikael Slevinsky <Richard.Slevinsky@umanitoba.ca>

Project.toml

@@ -26,7 +26,7 @@ ApproxFunSingularities = "0.1.6"
 BandedMatrices = "0.14.2, 0.15"
 DomainSets = "0.1, 0.2, 0.3"
 DualNumbers = "0.6"
-HypergeometricFunctions = "0.1, 0.2"
+HypergeometricFunctions = "0.1, 0.2, 0.3"
 LowRankApprox = "0.4"
 SpecialFunctions = "0.7, 0.8, 0.9, 0.10"
 julia = "1.3"

src/stieltjesmoment.jl

@@ -30,11 +30,12 @@ function logjacobimoment(α::Real,β::Real,n::Int,z)
     if n == 0
         2normalization(0,α,β)*(log(z-1)-dualpart(_₂F₁(dual(zero(α)+eps(α+β),one(β)),α+1,α+β+2,x)))
         # For testing purposes only, should be equivalent to above within radius of convergence
-        #2normalization(0,α,β)*(log(z-1)-(α+1)/(α+β+2)*x.*_₃F₂(α+2,α+β+3,x))
+        # 2normalization(0,α,β)*(log(z-1)-(α+1)/(α+β+2)*x*_₃F₂(α+2,α+β+3,x))
     else
         -2normalization(n,α,β)/n*(-x)^n*_₂F₁(n,n+α+1,2n+α+β+2,x)
     end
 end
+
 logjacobimoment(α::Real,β::Real,z) = logjacobimoment(α,β,0,z)
 
 
@@ -43,7 +44,7 @@ function logabsjacobimoment(α::Real,β::Real,n::Int,z)
     if n == 0
         2normalization(0,α,β)*(logabs(z-1)-real(dualpart(_₂F₁(dual(zero(α)+eps(α+β),one(β)),α+1,α+β+2,x))))
         # For testing purposes only, should be equivalent to above within radius of convergence
-        #2normalization(0,α,β)*(log(z-1)-(α+1)/(α+β+2)*x.*_₃F₂(α+2,α+β+3,x))
+        # 2normalization(0,α,β)*(logabs(z-1)-real((α+1)/(α+β+2)*x*_₃F₂(α+2,α+β+3,x)))
     else
         -2normalization(n,α,β)/n*real((-x)^n*_₂F₁(n,n+α+1,2n+α+β+2,x))
     end

test/logkerneltest.jl

@@ -36,7 +36,7 @@ import SingularIntegralEquations: testsies, testsieeval, stieltjesmoment, Direct
         @test (H*f)(3.5+0.5im) ≈ logkernel(f,3.5+0.5im)
     end
 
-    @testset "Circle" begin 
+    @testset "Circle" begin
         d=Circle(0.2,3.0)
         S=Fourier(d)
         ζ=Fun(d)
@@ -52,7 +52,7 @@ import SingularIntegralEquations: testsies, testsieeval, stieltjesmoment, Direct
         z=0.2+3im;@test (H*f)(z) ≈ logkernel(f,z)
     end
 
-    @testset "Legendre and Jacobi on intervel" begin
+    @testset "Legendre and Jacobi on interval" begin
         x=Fun()
         f=exp(x)