Is the following post-processing of the trace enough to fix this problem?
For each value of matrix \mathcal{A} in the trace,
- Calculate
indices = np.argsort(np.sum(A, axis=0)) - Reorder the matrix by
A = A[indices, :][:, indices]
Is the following post-processing of the trace enough to fix this problem?
For each value of matrix \mathcal{A} in the trace,
indices = np.argsort(np.sum(A, axis=0))
A = A[indices, :][:, indices]