special cases

Author: dlfivefifty

src/MultivariateSingularIntegrals.jl

@@ -2,7 +2,7 @@ module MultivariateSingularIntegrals
 using LinearAlgebra, ClassicalOrthogonalPolynomials, SingularIntegrals, FillArrays
 import Base: size
 
-export logkernelsquare, stieltjessquare
+export logkernelsquare, stieltjessquare, logkernelsquare!, stieltjessquare!
 
 
 function imlogkernel_vec(n, z)

src/logkernelsquare.jl

@@ -150,9 +150,11 @@ end
 function logkernelsquare!(A::AbstractMatrix{T}, z, F_1, F_2) where T
     m,n = size(A)
     @assert m  == n
+    n == 0 && return A
+    A[1,1] = L₀₀(z)
+    n == 1 && return A
     rec_rhs_1!(F_1, z)
     rec_rhs_2!(F_2, z)
-    A[1,1] = L₀₀(z)
     logkernelsquare_populatefirstcolumn!(A, z, F_1, F_2)
     logkernelsquare_populatefirstrow!(A, z, F_1, F_2)
 
@@ -169,7 +171,7 @@ function logkernelsquare!(A::AbstractMatrix{T}, z, F_1, F_2) where T
         A[2,j+1] = F[1,j+1] + im*(A[1,j+2] - A[1,j])/(2j+1)
     end
 
-    for ℓ = 1:((n-1)÷2-1)
+    @inbounds for ℓ = 1:((n-1)÷2-1)
         for k = ℓ+1:n-(ℓ+2)
             A[k+1,ℓ+2] = (2ℓ+1)*im*(F[k+1,ℓ+1] + (A[k,ℓ+1] - A[k+2,ℓ+1])/(2k+1)) + A[k+1,ℓ]
         end