Dear numpy users, I have a 3x3 matrix which I want to repeat 50 times along a diagonal, thus creating a 150x150 block diagonal matrix.
I know of a method usin scipy.linalg.block_diag, but I don't know if this is the best one: a = random.randn(3,3) b = a.reshape(1,3,3).repeat(50,axis=0) scipy.linalg.block_diag( *b ) Jose
_______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion
