Just wondering whether there have emerged any packages for this in the meantime. I need a rank 3 sparse matrix to implement a retarded potential integral equation solver. The structure will have non zero entries for all values of the first two indices and a varying but fixed length set for the third index:
A[i,j,k] != 0 for k0(i,j) <= k <= k1(i,j)
