add a dune.istl.BCRSMatrix method for arbitrary local block size
A BCRSMatrix can be constructed using
Matrix = dune.istl.BCRSMatrix( (brow,bcol) )
matrix = Matrix( (10,10),2,Matrix.random )
The solvers are still missing for more general BCRSMatrix types but this can be added later.
Edited by Andreas Dedner