revert despecialization of type in zero/one (#162)

Author: jishnub

src/constructors.jl

@@ -144,11 +144,11 @@ Fun(f::Type, d::Space) = error("Not implemented")
 
 # special case constructors
 zeros(S::Space) = zeros(Float64, S)
-zeros(T::Type, S::Space) = Fun(S,zeros(T,1))
+zeros(::Type{T}, S::Space) where {T<:Number} = Fun(S,zeros(T,1))
 
 # catch all
 ones(S::Space) = ones(Float64, S)
-ones(T::Type, S::Space) = Fun(x->one(T),S)
+ones(::Type{T}, S::Space) where {T<:Number} = Fun(x->one(T),S)
 
 function Fun(::typeof(identity), d::Domain)
     cd=canonicaldomain(d)

test/SpacesTest.jl

@@ -80,6 +80,12 @@ using ApproxFunOrthogonalPolynomials
             end
         end
 
+        f = @inferred (T -> ones(T, PointSpace(1:3)))(Float64)
+        @test f == Fun(PointSpace(1:3), [1.0, 1.0, 1.0])
+        f = @inferred (T -> zeros(T, PointSpace(1:3)))(Float64)
+        @test space(f) == PointSpace(1:3)
+        @test all(iszero(coefficients(f)))
+
         M = Multiplication(Fun(PointSpace(1:3), [1:3;]), PointSpace(1:3))
         @test (@inferred size(M)) == (3,3)
         M2 = M + M