Block matrix multiplication with block memory storage works now. There was a bug in the block index computation when passing to do_block