|
| From: | John W. Eaton |
| Subject: | [Octave-bug-tracker] [bug #35021] Floating point exception on matrix multiplication |
| Date: | Fri, 09 Dec 2011 16:00:36 +0000 |
| User-agent: | Mozilla/5.0 (X11; Linux x86_64; rv:8.0) Gecko/20100101 Firefox/8.0 Iceweasel/8.0 |
Follow-up Comment #10, bug #35021 (project octave):
Octave recognizes
X'*X
as a special case and ultimately calls DSYRK to perform the calculation. For
(X+0)'*X
Octave does not recognize that X+0 is the same as X, so it calls DGEMM to do
Y'*X
in which Y is the result of actually computing X+0.
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?35021>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
| [Prev in Thread] | Current Thread | [Next in Thread] |