Comment #2 on issue 1421 by jochen.voss: QRdecomposition fails for some
matrices
http://code.google.com/p/sympy/issues/detail?id=1421
Actully the three lines from the traceback, now that I look at them, show
what's
going on:
966 R[j,j] = tmp.norm()
967 Q[:,j] = tmp / R[j,j]
--> 968 assert Q[:,j].norm() == 1
The value "tmp" is divided by its norm and then we check whether the norm
of the
result is 1 (as a comparison between floats, I guess). This seems bound to
fail.
Instrumenting the function shows that the norm is indeed
1-1.1102230246e-16. Maybe
the assert should simply be removed? Or does it protect against something?
--
You received this message because you are listed in the owner
or CC fields of this issue, or because you starred this issue.
You may adjust your issue notification preferences at:
http://code.google.com/hosting/settings
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en
-~----------~----~----~----~------~----~------~--~---
