pub fn frobenius_norm(a: &CsrMatrix) -> f64
Compute the Frobenius norm of a sparse matrix: sqrt(sum(a_ij^2)).
sqrt(sum(a_ij^2))