Made the modifications to the block compressed sparse format so that it can deal with both BCSC and BCSR matrices. It passes some simple unit tests but I need to try it with nonsymmetric or rectangular matrices too.