On Tue, Apr 20, 2010 at 9:34 AM, Jason Grout <jason-s...@creativetrax.com> wrote: > > Do you mind posting this example? Was it truly a bug (i.e., unintentional > wrong behavior), or a result of double precision computation rounding things > off? In other words, were the returned evals almost real? What was the > condition number of the matrix? >
The condition number wasn't large (all the numbers were of the same order). It certainly was a bug which I reported to Mathworks. Last time I tried it (which was some time ago), it still wasn't fixed. Plus, for real SPD matrices, they should be using one of the routines that guarantees real eigenvalues. The code is on one of my backup CDs so it would be some work to dig it out. Cheers, Tim. --- Tim Lahey PhD Candidate, Systems Design Engineering University of Waterloo http://www.linkedin.com/in/timlahey -- To post to this group, send an email to sage-devel@googlegroups.com To unsubscribe from this group, send an email to sage-devel+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URL: http://www.sagemath.org
