this question was already asked here. However, that was 10 years ago; I'd like to bring it back, to see if there are new ways to store/retrieve/work with a 3D sparse matrix.
I know scipy.sparse still doesn't support 3D matrices, but maybe some other library/tool does it.