[SciPy-User] isn't it a bug in scipy.sparse? + some questions

Dmitrey tmp50@ukr....
Wed Nov 11 12:28:24 CST 2009


So, anyone doesn't know answers to the questions about scipy.sparse module mentioned below?  
As for the bug mentioned, I have installed latest numpy & scipy svn snapshots (1.4.0.dev7726 and 0.8.0.dev6096), the bug still exist.  
D.  
  
--- Исходное сообщение ---  
От кого: "Dmitrey" <tmp50@ukr.net>  
Кому: scipy-user@scipy.org  
Дата: 8 ноября, 13:22:34  
Тема: [SciPy-User] isn't it a bug in scipy.sparse? + some questions  
  
  Hi scipy.sparse developers and all other scipy users,  
I'm trying to take benefits for solving SLEs in FuncDesigner via involving scipy.sparse.  
Some examples are here  
http://openopt.org/FuncDesignerDoc#Solving_systems_of_linear_equations  
and example for sparse SLEs is here  
http://trac.openopt.org/openopt/browser/PythonPackages/FuncDesigner/FuncDesigner/examples/sparseSLE.py  
It already works faster than using dense matrices, but I want to speedup it even more, so I have some questions and seems like bug report (scipy.__version__ 0.7.0):  
  
from scipy import sparse  
from numpy import *  
a=sparse.lil_matrix((3,1))  
a[0:3,:] = ones(3)  
print a.todense()  
#prints  
[[ 1.]  
 [ 0.]  
 [ 0.]]  
while I expect all-ones  
  
Questions:  
1) Seems like a[some_ind,:]=something works very, very slow for lil. I have implemented a workaround, but can I use a[some_ind,:] for another format than lil? (seems like all other ones doesn't support it).  
2) What is current situation with matmat and matvec functions? They say "deprecated" but no alternative is mentioned.  
3) What is current situation with scipy.sparse.linalg.spsolve? It says  
/usr/lib/python2.6/dist-packages/scipy/sparse/linalg/dsolve/linsolve.py:78: DeprecationWarning: scipy.sparse.linalg.dsolve.umfpack will be removed, install scikits.umfpack instead  
  ' install scikits.umfpack instead', DeprecationWarning )  
But I don't want my code to be dependent on a scikits module. Are there another default/autoselect solver for sparse SLEs?  
If no, which one would you recommend me to use as default for sparse SLEs - bicg, gmres, something else?  
  
Thank you in advance, D.  
  
_______________________________________________  
SciPy-User mailing list  
SciPy-User@scipy.org  
http://mail.scipy.org/mailman/listinfo/scipy-user  
  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/scipy-user/attachments/20091111/4a0aac53/attachment.html 


More information about the SciPy-User mailing list