performance fix

Author: aquatiko

src/utils.jl

@@ -27,7 +27,7 @@ end
 Pads a given array and creates a border with edge elements.
 """
 function pad_edges(array_in::Matrix{T}) where {T}
-    image = zeros(size(array_in)[1] + 2, size(array_in)[2] + 2)
+    image = Matrix{T}(undef, size(array_in)[1] + 2, size(array_in)[2] + 2)
     image[2:end-1,2:end-1] = array_in
     image[2:end-1,1] = array_in[:,1]
     image[2:end-1,end] = array_in[:,end]

test/core.jl

@@ -25,10 +25,10 @@ end
     hdu1 = astropy.io.fits.open(joinpath("data", "gc_2mass_k.fits"))[1]
     hdu2 = astropy.io.fits.open(joinpath("data", "gc_msx_e.fits"))[1]
 
-    @test isapprox(reproject(imgin, imgout, order = 0)[1]', rp.reproject_interp(hdu2, hdu1.header, order = 0)[1], nans = true, atol = 1e-4)
-    @test isapprox(reproject(imgout, imgin, order = 0)[1]', rp.reproject_interp(hdu1, hdu2.header, order = 0)[1], nans = true, atol = 1e-3)
-    @test isapprox(reproject(imgin, imgout, order = 1)[1]', rp.reproject_interp(hdu2, hdu1.header, order = 1)[1], nans = true, atol = 1e-4)
-    @test isapprox(reproject(imgin, imgout, order = 2)[1]', rp.reproject_interp(hdu2, hdu1.header, order = 2)[1], nans = true, atol = 1e-4)
+    @test isapprox(reproject(imgin, imgout, order = 0)[1]', rp.reproject_interp(hdu2, hdu1.header, order = 0)[1], nans = true, rtol = 1e-4)
+    @test isapprox(reproject(imgout, imgin, order = 0)[1]', rp.reproject_interp(hdu1, hdu2.header, order = 0)[1], nans = true, rtol = 1e-3)
+    @test isapprox(reproject(imgin, imgout, order = 1)[1]', rp.reproject_interp(hdu2, hdu1.header, order = 1)[1], nans = true, rtol = 1e-4)
+    @test isapprox(reproject(imgin, imgout, order = 2)[1]', rp.reproject_interp(hdu2, hdu1.header, order = 2)[1], nans = true, rtol = 1e-4)
     @test isapprox(reproject(imgin[1], imgout[1], shape_out = (1000,1000))[1]',
             rp.reproject_interp(hdu2, astropy.wcs.WCS(hdu1.header), shape_out = (1000,1000))[1], nans = true, atol = 1e-4)