Treat spaces as broadcast scalars (#417)

Author: jishnub

src/Space.jl

@@ -40,6 +40,8 @@ rangetype(::Space{D,R}) where {D,R} = R
 rangetype(::Type{Space{D,R}}) where {D,R} = R
 rangetype(::Type{FT}) where {FT<:Space} = rangetype(supertype(FT))
 
+Base.broadcastable(S::Space) = Ref(S)
+
 domaindimension(sp::Space) = dimension(domain(sp))
 """
     dimension(s::Space)

test/SpacesTest.jl

@@ -189,6 +189,15 @@ using LinearAlgebra
         @testset "call spaces" begin
             @test PointSpace(1:4)(3, 4) == 1
         end
+
+        @testset "Broadcasting" begin
+            f1 = x -> x
+            f2 = x -> x^2
+            sp = PointSpace(1:3)
+            fs = Fun.((f1, f2), sp)
+            @test fs[1] == Fun(f1, sp)
+            @test fs[2] == Fun(f2, sp)
+        end
     end
 
     @testset "DiracSpace" begin