Added a feature to form explicitly the matrix.
Fix a bug in BlockLinearOperator for block matrices
bigger than 2x2.
Add a test.