Implement custom broadcasting #35
Description
In order to make element-wise operations faster between BlockDiagonal matrices with the same block structure, a custom broadcasting interface should be implemented (see https://docs.julialang.org/en/v1/manual/interfaces/#man-interfaces-broadcasting-1). I would like to help, but I don't think I have the necessary knowledge of how this interface works.