[SciPy-user] Solution y <- Ax, where A is sparse
Dinesh B Vadhia
Tue Oct 9 15:48:07 CDT 2007
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.
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the SciPy-user