constructor fix

Author: jverzani

src/Polynomials.jl

@@ -58,22 +58,39 @@ q = Poly([1, 2, 3], :s)
 p + q                  # ERROR: Polynomials must have same variable.
 ```
 
-"""
-immutable Poly{T}
-  a::Vector{T}
-  var::Symbol
-  @compat function (::Type{Poly}){T<:Number}(a::Vector{T}, var::SymbolLike = :x)
-    # if a == [] we replace it with a = [0]
-    if length(a) == 0
-      return new{T}(zeros(T,1), @compat Symbol(var))
-    else
-      # determine the last nonzero element and truncate a accordingly
-      a_last = max(1,findlast(x->x!=zero(T), a))
-      new{T}(a[1:a_last], @compat Symbol(var))
+    """
+
+immutable Poly{T<:Number}
+    a::Vector{T}
+    var::Symbol
+    function Poly(a::Vector{T}, var::SymbolLike=:x)
+        # if a == [] we replace it with a = [0]
+        if length(a) == 0
+            return new{T}(zeros(T,1), @compat Symbol(var))
+        else
+            # determine the last nonzero element and truncate a accordingly
+            a_last = max(1,findlast(x->x!=zero(T), a))
+            new{T}(a[1:a_last], @compat Symbol(var))
+        end
     end
-  end
 end
 
+# immutable Poly{T}
+#   a::Vector{T}
+#   var::Symbol
+#   @compat function (::Type{Poly}){T<:Number}(a::Vector{T}, var::SymbolLike = :x)
+#     # if a == [] we replace it with a = [0]
+#     if length(a) == 0
+#       return new{T}(zeros(T,1), @compat Symbol(var))
+#     else
+#       # determine the last nonzero element and truncate a accordingly
+#       a_last = max(1,findlast(x->x!=zero(T), a))
+#       new{T}(a[1:a_last], @compat Symbol(var))
+#     end
+#   end
+# end
+
+Poly{T <: Number}(x::Vector{T}, var::SymbolLike=:x) = Poly{T}(x, var)
 Poly(n::Number, var::SymbolLike = :x) = Poly([n], var)
 
 

test/runtests.jl

@@ -13,6 +13,10 @@ p5 = Poly([1,4,6,4,1,0,0,0,0,0,0,0,0,0,0,0,0,0])
 pN = Poly([276,3,87,15,24,0])
 pR = Poly([3//4, -2//1, 1//1])
 X = Poly([0.0, 1.0])
+T = Int64
+Poly{T}([zero(T), one(T)])
+Poly{T}([zero(T), one(T)], :y)
+
 p1000 = Poly(randn(1000))
 
 @test length(pNULL) == 1
@@ -279,3 +283,4 @@ pint  = polyint(p)
 @test isnan(p(1))                 # p(1) evaluates to NaN
 @test isequal(pder, Poly([NaN]))  # derivative will give Poly([NaN])
 @test isequal(pint, Poly([NaN]))  # integral will give Poly([NaN])
+