[Scipy-tickets] [SciPy] #1456: scipy.sparse.linalg.lobpcg fails for large sparse matrices

SciPy Trac scipy-tickets@scipy....
Fri Jun 10 14:44:43 CDT 2011


#1456: scipy.sparse.linalg.lobpcg fails for large sparse matrices
---------------------------------+------------------------------------------
 Reporter:  AndrewMacLean        |       Owner:  wnbell      
     Type:  defect               |      Status:  needs_review
 Priority:  normal               |   Milestone:              
Component:  scipy.sparse.linalg  |     Version:  0.7.2       
 Keywords:                       |  
---------------------------------+------------------------------------------
Changes (by AndrewMacLean):

 * cc: andrew.maclean@… (added)
  * status:  needs_info => needs_review


Comment:

 Thanks for these answers. It appears that slow convergence like this is
 normal for matrices of this size for this algorithm:
 [http://www.mathworks.com/matlabcentral/answers/9063-lobpcg-returning-
 incorrect-results-for-large-sparse-symmetric-matrices]

 So, there is nothing wrong with Scipy's implementation for this case.
 Sending a pre-conditioner to speed up convergence, as described in the
 response in the link above, should solve my problem.

-- 
Ticket URL: <http://projects.scipy.org/scipy/ticket/1456#comment:4>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list