Functionality for working efficiently with block diagonal matrices. - View it on GitHub