Semi-Global Weighted Least Squares in Image Filtering
an introduction to the conjugate gradient method without the agonizing pain .