Test with DoubleFloat

Author: dlfivefifty

Project.toml

@@ -11,6 +11,7 @@ SingularIntegrals = "d7440221-8b5e-42fc-909c-0567823f424a"
 
 [compat]
 ClassicalOrthogonalPolynomials = "0.15.1"
+DoubleFloats = "1"
 FillArrays = "1.13.0"
 LinearAlgebra = "1"
 QuadGK = "2.11.2"
@@ -19,8 +20,9 @@ Test = "1"
 julia = "1.10"
 
 [extras]
+DoubleFloats = "497a8b3b-efae-58df-a0af-a86822472b78"
 QuadGK = "1fd47b50-473d-5c70-9696-f719f8f3bcdc"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [targets]
-test = ["Test", "QuadGK"]
+test = ["Test", "DoubleFloats", "QuadGK"]

test/runtests.jl

@@ -1,4 +1,4 @@
-using MultivariateSingularIntegrals, QuadGK, ClassicalOrthogonalPolynomials, LinearAlgebra, Test
+using MultivariateSingularIntegrals, QuadGK, ClassicalOrthogonalPolynomials, LinearAlgebra, DoubleFloats, Test
 
 
 Z̃ = (2.0, 2.0im, 2.0+2im, 2.0-2im, -2.0 - 2im)
@@ -47,6 +47,8 @@ end
                 M = Diagonal((P'P)[1:n,1:n])
                 @test P[big(0.), 1:n]' *  inv(M)  * logkernelsquare(z, n) * inv(M) * P[big(0.), 1:n] ≈ log(z) atol=1E-40
             end
+
+            @test logkernelsquare(Double64(0) + 0im, 10) ≈ logkernelsquare(big(0.0) + 0im, 10) atol=1E-30
         end
     end
 end