Just as I was making up an example for the block diagonal matrix question, I remembered the old problem I had with <br>consistency of nympy functions. <br><br>If you want to generate a random number matrix, you can make the same call as with matlab: <br>
<br>rand(2,2) for 2x2 matrix,<br>randn(1,5) for 1x5 etc. <br><br>but if you want to generate ones or zeros matrices, you cannot say<br>ones(3,3), you have to write ones([3,3]) or zeros([3,3]) (note the extra brackets). <br>
<br>It is not a big deal, but it seems a bit inconsistent for me. <br>