possible to increase speed of power function applied to a vector?

Question

ER2018 2018-7-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/408369-possible-to-increase-speed-of-power-function-applied-to-a-vector

评论： ER2018 2018-7-3

Is it possible to increase speed of power function applied to a vector?

 N=16e3;
 X=-2+4*rand(N,1);
 Y=10.^X;

I tried

Y=power(10,X);

But of course it gives me the same speed. Any ideas?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Yair Altman 2018-7-3

Cross-link to Stack Overflow

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jan 2018-7-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/408369-possible-to-increase-speed-of-power-function-applied-to-a-vector#answer_327124

No. power is expensive. There is no magic super power operator, which can do the same work in less time.

In some situations you can replace the power operation by accumulated multiplications:

a  = rand;
X1 = a .^ (1:10);
X2 = cumprod(repmat(a, 1, 10));

But remember that this accumulates the rounding errors also, so for large powers it gets less accurate.

Another common case, where the power operation can be avoided:

a = 1.5 * 10^8    % Expensive power and multiplication
b = 1.5e8         % Cheap constant

Especially if this appear in loops or in functions to be integrated, the constant can save much time.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

ER2018 2018-7-3

Thank you for the tips. Unfortunately there's no cheap constant I can use here, I really need the calculated values.

How would the solution with cumprod look in the example I provided?

请先登录，再进行评论。

Answer 2

Yair Altman 2018-7-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/408369-possible-to-increase-speed-of-power-function-applied-to-a-vector#answer_327251

In addition to what Jan said, perhaps your in your specific case it is possible to work at the exponent level, i.e. not Y=10^X but rather log10Y=X, and then continue processing using log10Y rather than with Y. Depending on what you do with the data later in your code, it might possibly be faster since it avoids the initial power operation.