octave-maintainers
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

sparse is_symmetric functions


From: John W. Eaton
Subject: sparse is_symmetric functions
Date: Wed, 6 Dec 2006 13:11:34 -0500

On  4-Dec-2006, David Bateman wrote:

| While starting to write svds I found that I need to check if a matrix
| passed to eigs is symmetric, rather than symmetric with a positive
| diagonal as the cheap test for positive definite matrices does. At the
| moment eigs uses arpack's dseupd only for real positive definite
| matrices, where in fact it could use and should use it for the case of
| [sparse(size(A,1),size(A,1)),A;A',sparse(size(A,2),size(A,2))] as
| created by svds. I therefore need the is_symmetric function of the
| Matrix and SparseMatrix classes.
| 
| However, the current sparse is_symmetric is very inefficient. The
| attached patch makes this method of the Sparse classes much better.

OK, please check in these changes.

Thanks,

jwe


reply via email to

[Prev in Thread] Current Thread [Next in Thread]