We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 1995031 commit d22d8f1Copy full SHA for d22d8f1
base/math.jl
@@ -118,7 +118,8 @@ function _evalpoly(x, p)
118
return s
119
end
120
121
-function evalpoly(z::Complex, p::Tuple)
+# Goertzel-like algorithm from Knuth, TAOCP vol. 2, section 4.6.4:
122
+function evalpoly(z::Complex, p::Tuple{Any, Any, Vararg})
123
if @generated
124
N = length(p.parameters)
125
a = :(p[end])
@@ -143,10 +144,6 @@ function evalpoly(z::Complex, p::Tuple)
143
144
_evalpoly(z, p)
145
146
-evalpoly(z::Complex, p::Tuple{<:Any}) = p[1]
147
-evalpoly(z::Complex, ::Tuple{}) = zero(one(z)) # dimensionless zero, i.e. 0 * z^0
148
-
149
-evalpoly(z::Complex, p::AbstractVector) = _evalpoly(z, p)
150
151
function _evalpoly(z::Complex, p)
152
Base.require_one_based_indexing(p)
0 commit comments