Reuse memory in RHS

Author: dlfivefifty

src/logkernelsquare.jl

@@ -156,8 +156,8 @@ function logkernelsquare!(A::AbstractMatrix{T}, z, F_1, F_2) where T
     logkernelsquare_populatefirstcolumn!(A, z, F_1, F_2)
     logkernelsquare_populatefirstrow!(A, z, F_1, F_2)
 
-    # F = F_1 # reuse the memory
-    F = F_2 .- F_1
+    F = F_1 # reuse the memory
+    F .= F_2 .- F_1
 
     # 2nd row/column
 

test/runtests.jl

@@ -41,7 +41,7 @@ end
 
     @testset "BigFloat" begin
         setprecision(2000) do
-            for z in (big(2.0), big(2.0im), big(2.0+2im), big(2.0-2im), big(-2.0 - 2im))
+            for z in big.(Z̃)
                 P = Legendre{BigFloat}()
                 n = 100
                 M = Diagonal((P'P)[1:n,1:n])