Sparse Matrix Multiplications

[Peter Otten]

Peter Otten wrote:

> # use at your own risk
> import numpy
> 
> N = 10**4 # I get a MemoryError for a bigger exponent
> b = numpy.array(range(N))
> a = numpy.zeros((N, N)) + b
> a *= a.transpose()
> a[0,0] = (b*b).sum()
> print a

Sorry, this is nonsense.