[SciPy-user] Sparse indexing workarounds
Tue Apr 10 04:11:51 CDT 2007
You can use pysparse instead of sparse. The update_add_mask function will do
On 4/10/07, Robert Cimrman <email@example.com> wrote:
> Bill Baxter wrote:
> > Does anyone have an easy (and efficient way) to update a submatrix of
> > a big sparse matrix?
> > K = scipy.sparse.lil_matrix(bigN,bigN)
> > ...
> > conn = [1,4,11,12]
> > K[ix_(conn,conn)] += elemK
> > where elemK is a 4x4 dense matrix.
> > This kind of thing is very commonly needed in FEM codes for assembling
> > the global stiffness matrix. But sparse doesn't seem to support
> > either += or indexing with the open grid of indices returned by ix_.
> You may have a look at
> it uses CSR matrix though...
> SciPy-user mailing list
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the SciPy-user