Merge pull request #1 from hrobotron/mp-patch

hrobotron · web-flow · commit 5f833b1bafc3 · 2016-12-19T08:46:08.000+01:00
Enable complex and other float() types
diff --git a/src/Polynomials.jl b/src/Polynomials.jl
@@ -579,10 +579,10 @@ function polyfit(x, y, n::Int=length(x)-1, sym::Symbol=:x)
     1 <= n <= length(x) - 1 || throw(DomainError)
 
     #
-    # here unsure, whether zeros(Float64), or similar(x,...)
+    # here unsure, whether similar(float(x[1]),...), or similar(x,...)
     # however similar may yield unwanted surprise in case of e.g. x::Int
     #
-    A=zeros(Float64, length(x), n+1)
+    A=similar(float(x[1]), length(x), n+1)
     #
     # TODO: add support for poly coef bitmap
     # (i.e. polynomial with some terms fixed to zero)

Original file line number	Diff line number	Diff line change
`@@ -579,10 +579,10 @@ function polyfit(x, y, n::Int=length(x)-1, sym::Symbol=:x)`
`579`	`579`	`1 <= n <= length(x) - 1 \|\| throw(DomainError)`
`580`	`580`
`581`	`581`	`#`
`582`		`- # here unsure, whether zeros(Float64), or similar(x,...)`
	`582`	`+ # here unsure, whether similar(float(x[1]),...), or similar(x,...)`
`583`	`583`	`# however similar may yield unwanted surprise in case of e.g. x::Int`
`584`	`584`	`#`
`585`		`- A=zeros(Float64, length(x), n+1)`
	`585`	`+ A=similar(float(x[1]), length(x), n+1)`
`586`	`586`	`#`
`587`	`587`	`# TODO: add support for poly coef bitmap`
`588`	`588`	`# (i.e. polynomial with some terms fixed to zero)`