The problem is here. My latest attempt is here.

I've implemented a version of this algorithm. It works for the smaller tests, however in the worst case (100000x100000 matrix, nothing blurry) my implementation still takes about 13 seconds.

I've optimized it as far as I'm able to. Is there something I'm missing? Or another approach entirely? I'm stuck.