BlockLinearOperator for sizes more than 2x2

Hello,

I found a bug in the BlockLinearOperator class for matrices bigger than 2x2.
Here is a fix.

I also added a feature to LinearOperator which allows us to form the actual matrix.

Sylvain