[SciPy-dev] Sparse feature suggestion: zero-sized sparse matrices
Mon Feb 11 11:09:49 CST 2008
On Feb 11, 2008 10:50 AM, Neilen Marais <firstname.lastname@example.org> wrote:
> This kind of ties in to the sparse block matrix assembly I brought up
> last week. I build up the sparse blocks from different subsystems, but in
> some cases some of the subsystems fall away. It would be very convenient
> if I could represent the "null" systems simply as zero-sized sparse
> matrices that will then make no contribution to the global matrix.
> It seems normal numpy arrays support the idea of zero-dimension 2-D
> arrays. Would anyone object to this feature in sparse matrices? If not, I
> could try to cook up a patch.
You're welcome to try, but I think it would be a nightmare to implement.
I would only object if the necessary changes introduced a large number
of statements of the form
if shape == 0:
do strange thing
do normal thing
because that sort of thing is difficult to maintain.
On the other hand, there may be a comprehensive solution that resolves
such matters in a simple and consistent way.
Nathan Bell email@example.com
More information about the Scipy-dev