#7 zeros were added to north and south diags
#8 residual norm computing were improved