* https://github.com/kangwonlee/nmisp/commit/aebaad178bd306e5d29f723a5e9632cacb1789e1 * Sequential method may need many iterations * Can parallelization help? CUDA? Numba.CUDA ?