Fig. REF  presents a comparison of fps between the inherently sequential CPU and the massively parallel sIM for the factorization problem. The CPU calculation is done in MATLAB by an iterative solution of Eq. (REF )-(REF ) using standard Gibbs sampling [1]} with simulated annealing. Even though optimization techniques including graph coloring by multiple threads could improve our MATLAB implementation, these should not change our conclusions since fps between the CPU and sIM differ by 6-orders of magnitude.
