[SciPy-user] Solution y <- Ax, where A is sparse

[Dinesh B Vadhia]

Hello!  Co-incidentally, I too am trying (or rather struggling) to write Python code to solve the problem:

y <- Ax

where A is an m (column) by n (row) sparse matrix (with m>>n) and the matrix entries are either 0's or 1's (with no descernible pattern ie. assume random entries), x is a floating point vector and so is the result vector y.  

I was looking at the SciPy sparse matrix functions (see http://www.rexx.com/~dkuhlman/scipy_course_01.html#sparse) and will look at PySparse.  

I was hoping to find in the Python world something as simple to use as the SparseLib++ for C++ library (see http://math.nist.gov/sparselib++/) which we have used successfully for many years.

In the meantime, any help and pointers would be truly appreciated.

Dinesh


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20071009/4cfac119/attachment.html>