why square takes longer than multiplication

Question

Jerry 2011-12-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/23164-why-square-takes-longer-than-multiplication

My question is why in matlab square takes more time than multiplication. I read an article, say, square actually can be done a little bit quicker than multiplication. But the test code below doesn't support this. So how does Matlab implement square(and other integer power, which I'm also interested in)? The test code is following:

a=1; tic; for k=1:10000 b=a*a; end toc;

Then, without any changes except "b=a*a" is changed into "b=a^2":

a=1; tic; for k=1:10000 b=a^2; end toc;

Consequently, multiplication takes about 77 microseconds while square takes 500+.

Thanks in advance!! Jerry

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2011-12-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/23164-why-square-takes-longer-than-multiplication#answer_30421

The .^ operator is implemented in terms of pow(), which can involve logs. .^ does not appear to optimize .^2

If I recall correctly, there is a File Exchange contribution that decomposes integer powers to create the optimal multiplication series.

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

Jerry 2011-12-6

Thank you very much. Do you know if there is anyway to find how Matlab writes its built-in function?

Walter Roberson 2011-12-6

The way to find the code for a built-in MATLAB function is to get a code development job with MathWorks.

I did at one time write code (for my work) that decomposed integer powers into multiplications. Eventually, though, I removed that code again, as I was able to show that the result was lower precision than .^ was able to get.

请先登录，再进行评论。

Answer 2

Jan 2011-12-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/23164-why-square-takes-longer-than-multiplication#answer_30442

在 MATLAB Online 中打开

Squaring is cheaper than the multiplication if it is performed in C and if it is implemented using the multiplication (!):

double a = 3.14159265, b;
b = a * a;  // Standard 
a *= a;     // Slightly faster inplace squaring

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

Jerry 2011-12-6

Thank you! I hope to accept your answer as well, but unexpectedly I can accept just one.

请先登录，再进行评论。