Matrix multiplication gives different result than manual dot products with each column

Why should they? The algorithm can perform diffrent sum orders, different accumulation scheme, different thread numbers. One should not expect they to be exactly equal.

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

Bruno Luong 2023-3-2

编辑：Bruno Luong 2023-3-2

When one works with numerical calculation on computer, one should be awared about the imperfect of finite precision aritmetics (round off error, non associativity, non commutivity, ...).

The "*" MATLAB operation is performed by Blas and Lapack library and can be differently implemented by platform and version. They are not documented by TMW because they reserved the freedom to change it (and they did serveral times in the pass).

The proof is that your code returns identical result on Linux machine R2022b (see @KSSV answer), but not on Windows (for example my laptop or your computer). It can also change with HW such as the number of cores of the CPU.

Bottomline is that use should not numerical calculation of the same mathematical expression returns the exact same answer if they call different function calls to achieve the same expression.

请先登录，再进行评论。

Answer 2

KSSV 2023-3-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1921680-matrix-multiplication-gives-different-result-than-manual-dot-products-with-each-column#answer_1183510

As you are comparing floatting point number, you should not use isequal. You should proceed like shown below:

N = 1000;
B = 10;
C = 2;
x = randn(N,B);
b = randn(B,C);
matMult = x * b;
dotProd = [x*b(:,1), x*b(:,2)];
tol = 10^-5 ; % can be changed 
all(abs(matMult-dotProd)<tol,'all')
ans = logical
   1