SIAM Journal o Scientific Computing, 37 (2015), B962-B979 http://dx.doi.org/10.1137/15M1022562