John Travers <[email protected]> writes: > On 28 Aug 2013, at 14:58, Jed Brown <[email protected]> wrote: >> csr=(A.indptr[rstart:rend] - A.indptr[rstart], >> A.indices[A.indptr[rstart]:A.indptr[rend]], >> A.data[A.indptr[rstart]:A.indptr[rend]]) >> > > Thanks, this works (except it should be A.indptr[rstart:rend+1] in the > first line I think).
Yes, good catch. > Well, not in the long run, but at the moment that code is debugged and > working and I only wanted to try the slepc solvers first. Great, let us know if you have any problems.
pgp0p67tIvZbQ.pgp
Description: PGP signature
