pub fn jacobi_scale(a: &CsrMatrix) -> CsrMatrixExpand description
Diagonal (Jacobi) scaling: return D^{-1} A where D = diag(A).
Rows with zero diagonal are left unchanged (scale factor = 1). This is the simplest preconditioner and is used to equilibrate the system before iterative solves.