fix broadccasting

Author: mileslucas

src/core.jl

@@ -29,11 +29,12 @@ abstract type PSFModel{T} <: AbstractMatrix{T} end
 # getindex just calls model with reversed indices
 Base.getindex(model::PSFModel, idx::Vararg{<:Integer,2}) = model(reverse(idx))
 # always inbounds
-Base.checkbounds(::Type{Bool}, ::PSFModel, idx::Vararg{Int,2}) = true
+Base.checkbounds(::Type{Bool}, ::PSFModel, idx) = true
+Base.checkbounds(::Type{Bool}, ::PSFModel, idx...) = true
 
 # broadcasting hack to slurp other axes (doesn't work for numbers)
-Broadcast.combine_axes(psf::PSFModel, other::AbstractMatrix) = axes(other)#map(intersect, axes(other), axes(psf))
-Broadcast.combine_axes(other::AbstractMatrix, psf::PSFModel) = map(intersect, axes(other), axes(psf))
+Broadcast.combine_axes(psf::PSFModel, other::AbstractMatrix) = axes(other)
+Broadcast.combine_axes(other::AbstractMatrix, psf::PSFModel) = axes(other)
 
 ## get the position depending on the keyword inputs
 function _position(nt::NamedTuple{(:r, :theta)})

test/runtests.jl

@@ -21,8 +21,8 @@ function test_model_interface(K)
 
     # broadcasting to an array that's smaller
     @test m .* ones(10, 10) ≈ ones(10, 10) .* m ≈ m[1:10, 1:10]
-    # broadcasting to a larger array cuts out a stamp
-    @test m .* ones(30, 30) ≈ ones(30, 30) .* m ≈ m[1:15, 1:15]
+    # broadcasting to a larger array
+    @test m .* ones(30, 30) ≈ ones(30, 30) .* m ≈ m[1:30, 1:30]
 
     # test new position
     m = K(x=12, y=13, fwhm=10)