8

我写了一个代码，它在数字上使用了最高 n 阶的勒让德多项式。例如：

....
case 8 
p = (6435*x.^8-12012*x.^6+6930*x.^4-1260*x.^2+35)/128; return
case 9 
...

如果向量x很长，这可能会变慢。我看到说x.^4和之间存在性能差异，x.*x.*x.*x并认为我可以使用它来改进我的代码。我已经使用timeit并发现：

x=linspace(0,10,1e6);
f1= @() power(x,4)
f2= @() x.4;
f3= @() x.^2.^2
f4= @() x.*x.*x.*x

f4比其他人快2 倍。但是，当我去时，和x.^6之间几乎没有区别 （而所有其他选项都较慢）。(x.*x.*x).^2x.*x.*x.*x.*x.*x

有没有告诉我们什么是最有效的方式来获取向量的力量？你能解释一下为什么会有这么大的性能差异吗？

4

3 回答 3

8

这不完全是您问题的答案，但它可能会解决您的问题：

x2 = x.*x; % or x.^2 or power(x,2), whichever is most efficient
p = ((((6435*x2-12012)*x2+6930)*x2-1260)*x2+35)/128

这样你只做一次幂，而且只用指数 2。这个技巧可以应用于所有勒让德多项式（在奇数多项式中，一个x2被替换为x）。

于 2013-09-24T23:48:26.890 回答
1

以下是一些想法：

power(x,4)并且x.^4是等效的（只需阅读文档）。

x.*x.*x.*x可能已优化为类似x.^2.^2

x.^2.^2可能被评估为：取每个元素的平方（快速），然后再次取其平方（再次快速）。

x.^4大概直接评价为：取每个元素的四次方（慢）。

看到 2 次快速操作比 1 次慢速操作花费的时间更少，这并不奇怪。太糟糕了，在幂 4 的情况下没有执行优化，但也许它并不总是有效或需要付出代价（输入检查，内存？）。

关于时间安排：实际上有比因素 2 更多的差异！

当您现在在函数中调用它们时，在每种情况下都会增加函数开销，从而使相对差异更小：

y=x;tic,power(x,4);toc
y=x;tic,x.^4;toc
y=x;tic,x.^2.^2;toc
y=x;tic,x.*x.*x.*x;toc

会给：

Elapsed time is 0.034826 seconds.
Elapsed time is 0.029186 seconds.
Elapsed time is 0.003891 seconds.
Elapsed time is 0.003840 seconds.

因此，这几乎是 10 倍的差异。但是，请注意，以秒为单位的时间差异仍然很小，因此对于大多数实际应用程序，我只会使用简单的语法。

于 2013-09-25T14:21:21.587 回答
1

似乎 Mathworks 在其幂函数中具有特殊的大小写正方形（不幸的是，它都是我们看不到的内置封闭源代码）。在我对 R2013b 的测试中，它看起来好像.^,power和realpow使用相同的算法。对于正方形，我相信他们已经将其特例化为x.*x.

1.0x (4.4ms):   @()x.^2
1.0x (4.4ms):   @()power(x,2)
1.0x (4.5ms):   @()x.*x
1.0x (4.5ms):   @()realpow(x,2)
6.1x (27.1ms):  @()exp(2*log(x))

对于立方体，情况就不同了。它们不再是特殊情况。再一次，.^,power和realpow所有都是相似的，但这次要慢得多：

1.0x (4.5ms):   @()x.*x.*x
1.0x (4.6ms):   @()x.*x.^2
5.9x (26.9ms):  @()exp(3*log(x))
13.8x (62.3ms): @()power(x,3)
14.0x (63.2ms): @()x.^3
14.1x (63.7ms): @()realpow(x,3)

让我们跳到 16 次方，看看这些算法是如何扩展的：

1.0x (8.1ms):   @()x.*x.*x.*x.*x.*x.*x.*x.*x.*x.*x.*x.*x.*x.*x.*x
2.2x (17.4ms):  @()x.^2.^2.^2.^2
3.5x (27.9ms):  @()exp(16*log(x))
7.9x (63.8ms):  @()power(x,16)
7.9x (63.9ms):  @()realpow(x,16)
8.3x (66.9ms):  @()x.^16

所以: .^，power并且realpow所有关于指数的运行时间都是恒定的，除非它是特殊情况（-1 也似乎是特殊情况）。使用这个exp(n*log(x))技巧也是关于指数的恒定时间，而且速度更快。唯一的结果我不太明白为什么重复平方比乘法慢。

正如预期的那样，将 的大小x增加 100 倍同样会增加所有算法的时间。

那么，这个故事的寓意是什么？使用标量整数指数时，请始终自己进行乘法运算。和朋友们有很多聪明才智power（指数可以是浮点数、向量等）。唯一的例外是 Mathworks 为您完成优化的地方。在 2013b 中，它似乎是x^2和x^(-1). 希望随着时间的推移他们会增加更多。但是，一般来说，求幂很困难，而乘法很容易。在对性能敏感的代码中，我认为你不会因为总是输入而出错x.*x.*x.*x。（当然，在您的情况下，请遵循 Luis 的建议，并利用每个学期的中间结果！）

function powerTest(x)

f{1} = @() x.*x.*x.*x.*x.*x.*x.*x.*x.*x.*x.*x.*x.*x.*x.*x;
f{2} = @() x.^2.^2.^2.^2;
f{3} = @() exp(16.*log(x));
f{4} = @() x.^16;
f{5} = @() power(x,16);
f{6} = @() realpow(x,16);

for i = 1:length(f)
    t(i) = timeit(f{i});
end

[t,idxs] = sort(t);
fcns = f(idxs);

for i = 1:length(fcns)
    fprintf('%.1fx (%.1fms):\t%s\n',t(i)/t(1),t(i)*1e3,func2str(fcns{i}));
end
于 2013-09-25T14:52:06.563 回答