Just a guess, but I think the problem is just that usually when people
solve sparse equations the matrix is sparse, but the RHS vector is
dense.   I don't know if UMFPACK or SuperLU can actually handle a
sparse RHS vector, but if they can, then it sounds like a SciPy bug.

I guess it should work if you convert your 'f' to a dense array before
making the calls to the solver.

Regards,
--bb
_______________________________________________
SciPy-user mailing list
SciPy-user@[...].org
http://projects.scipy.org/mailman/listinfo/scipy-user