copy returns Fill

Author: dlfivefifty

src/FillArrays.jl

@@ -3,7 +3,8 @@ module FillArrays
 using LinearAlgebra, SparseArrays
 import Base: size, getindex, setindex!, IndexStyle, checkbounds, convert,
     +, -, *, /, \, diff, sum, cumsum, maximum, minimum, sort, sort!,
-    any, all, axes, isone, iterate, unique, allunique, permutedims, inv
+    any, all, axes, isone, iterate, unique, allunique, permutedims, inv,
+    copy
 
 import LinearAlgebra: rank, svdvals!, tril, triu, tril!, triu!, diag
 
@@ -74,6 +75,8 @@ convert(::Type{AbstractFill}, F::AbstractFill) = F
 convert(::Type{AbstractFill{T}}, F::AbstractFill) where T = convert(AbstractArray{T}, F)
 convert(::Type{AbstractFill{T,N}}, F::AbstractFill) where {T,N} = convert(AbstractArray{T,N}, F)
 
+copy(F::Fill) = Fill(copy(F.value), F.axes)
+
 """ Throws an error if `arr` does not contain one and only one unique value. """
 function unique_value(arr::AbstractArray)
     if isempty(arr) error("Cannot convert empty array to Fill") end
@@ -174,6 +177,8 @@ for (Typ, funcs, func) in ((:Zeros, :zeros, :zero), (:Ones, :ones, :one))
         convert(::Type{AbstractArray{T}}, F::$Typ) where T = AbstractArray{T}(F)
         convert(::Type{AbstractArray{T,N}}, F::$Typ) where {T,N} = AbstractArray{T,N}(F)
 
+        copy(F::$Typ) = F
+
         getindex(F::$Typ{T,0}) where T = getindex_value(F)
         function getindex(F::$Typ{T}, kj::Vararg{AbstractVector{II},N}) where {T,II<:Integer,N}
             checkbounds(F, kj...)

test/runtests.jl

@@ -129,6 +129,17 @@ import FillArrays: AbstractFill, RectDiagonal
         @test Zeros{Bool}(5) + x ≡ x
         @test -x ≡ Fill(-1,5)
     end
+
+    @testset "copy should return Fill" begin
+        x = Fill(1.0,10)
+        @test copy(x) ≡ x
+        x = Zeros(10)
+        @test copy(x) ≡ x
+        x = Fill([1.,2.],10)
+        @test copy(x) == x
+        @test copy(x) !== x
+        @test copy(x) isa Fill
+    end
 end
 
 @testset "RectDiagonal" begin