Why does a'*b yield a (slightly) different result than dot(a,b)?
2 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
I was working an building a matrix entry-wise and checking two different methods against each other and eventually found that dot(a,b) is giving me a different answer than a'*b, where a and b are both real 9948x1 column vectors. I can see that the implementation of dot uses conj(a)*b, but my vectors are real. In one example, dot(a,b) = -1.0351e-11 and a'*b = -8.6402e-12.
I think it's likely a floating point precision issue, but I'm curious about the implementation of the "*" operator and what might be leading to it. Can anyone shed some light on this?
3 Commenti
Risposta accettata
Walter Roberson
il 28 Nov 2017
For sufficiently large arrays, a'*b is going to invoke the high performance libraries, which are going to calculate the multiplications in multiple threads, doing partial summations, before eventually doing a final summation. Because the summations could be done in a different order than sum(conj(a).*b), the round-off can be different.
0 Commenti
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Matrices and Arrays in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!