[Scipy-tickets] [SciPy] #1307: Sparse matrices are limited to 2^32 non-zero elements

SciPy Trac scipy-tickets@scipy....
Wed Oct 13 20:19:32 CDT 2010


#1307: Sparse matrices are limited to 2^32 non-zero elements
--------------------------------------+-------------------------------------
 Reporter:  peb                       |       Owner:  wnbell     
     Type:  enhancement               |      Status:  new        
 Priority:  normal                    |   Milestone:  Unscheduled
Component:  scipy.sparse              |     Version:  0.8.0      
 Keywords:  compressed sparse matrix  |  
--------------------------------------+-------------------------------------
Changes (by wnbell):

  * type:  defect => enhancement
  * milestone:  0.9.0 => Unscheduled


Comment:

 Extending sparsetools to 64-bit indices should be straightforward.  All
 you need to do is tell SWIG to instantiate the template functions with
 64-bit integer indices here [1].  A portable type (e.g. int64) should be
 used instead of long or long long since those vary in size from platform
 to platform.

 The larger issue is (1) modifying the Python code to accept both int32 and
 int64 (2) ensuring that external libraries (ARPACK, SuperLU) continue to
 receive 32-bit integers (if that is what they require) and (3) silencing
 the complains that arise when sparsetools compile times double :)

 [1]
 http://projects.scipy.org/scipy/browser/trunk/scipy/sparse/sparsetools/sparsetools.i#L167

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


More information about the Scipy-tickets mailing list