Function sparse_identity 

pub fn sparse_identity(n: usize) -> (Vec<usize>, Vec<usize>, Vec<f64>)

Construct the CSR representation of an n x n identity matrix.

Returns (row_offsets, col_indices, values).