带有 varargs 的 Julia @evalpoly 宏

我正在尝试使用 Julia 的@evalpoly宏。当我手动提供系数时它可以工作，但我一直无法弄清楚如何通过数组提供这些系数

julia> VERSION
v"0.3.5"

julia> @evalpoly 0.5 1 2 3 4
3.25

julia> c = [1, 2, 3, 4]
4-element Array{Int64,1}: 
 1
 2
 3
 4

julia> @evalpoly 0.5 c
ERROR: BoundsError()

julia> @evalpoly 0.5 c...
ERROR: BoundsError()

julia> @evalpoly(0.5, c...)
ERROR: BoundsError()

有人能指出我在这方面的正确方向吗？

看到这个问题的精彩答案后添加

直到我研究其中一些答案之前，我才发现其中的一个微妙之处。这z论证@evalpoly可以是变量，但系数应为文字

julia> z = 0.5
0.5

julia> @evalpoly z 1 2 3 4
3.25

julia> @evalpoly z c[1] c[2] c[3] c[4]
ERROR: c not defined

查看最后一个命令的扩展的输出，可以看到，确实是 z 被分配给扩展中的变量，但系数被逐字插入到代码中。

julia> macroexpand(:@evalpoly z c[1] c[2] c[3] c[4])
:(if Base.Math.isa(z,Base.Math.Complex)
        #291#t = z
        #292#x = Base.Math.real(#291#t)
        #293#y = Base.Math.imag(#291#t)
        #294#r = Base.Math.+(#292#x,#292#x)
        #295#s = Base.Math.+(Base.Math.*(#292#x,#292#x),Base.Math.*(#293#y,#293#y))
        #296#a2 = c[4]
        #297#a1 = Base.Math.+(c[3],Base.Math.*(#294#r,#296#a2))
        #298#a0 = Base.Math.+(Base.Math.-(c[2],Base.Math.*(#295#s,#296#a2)),Base.Math.*(#294#r,#297#a1))
        Base.Math.+(Base.Math.*(#298#a0,#291#t),Base.Math.-(c[1],Base.Math.*(#295#s,#297#a1)))
    else 
        #299#t = z
        Base.Math.+(Base.Math.c[1],Base.Math.*(#299#t,Base.Math.+(Base.Math.c[2],Base.Math.*(#299#t,Base.Math.+(Base.Math.c[3],Base.Math.*(#299#t,Base.Math.c[4]))))))
    end)

我不相信你想做的事是可能的，因为@evalpoly是一个宏 - 这意味着它在编译时生成代码。它生成的是霍纳方法的非常有效的实现（在实数情况下），但要做到这一点，它需要知道多项式的次数。长度为c在编译时未知，因此它不起作用（也不能），而当您直接提供系数时，它就拥有所需的一切。

不过，错误消息不太好，所以如果可以的话，您可以在 Julia Github 页面上提交问题吗？

更新：为了回答问题的更新，是的，第一个参数可以是变量。你可以这样想：

function dostuff()
  z = 0.0
  # Do some stuff to z
  # Time to evaluate a polynomial!
  y = @evalpoly z 1 2 3 4
  return y
end

正变成

function dostuff()
  z = 0.0
  # Do some stuff to z
  # Time to evaluate a polynomial!
  y = z + 2z^2 + 3z^3 + 4z^4
  return y
end

除了，不是那样，因为它使用霍纳规则，但无论如何。问题是，它无法生成该表达式编译时间不知道系数的数量。但它不需要知道什么z根本就是。