Function ls_solver::gs_solve::get_gs_p

pub fn get_gs_p(a: &CsrMatrix<f64>, omega: f64) -> CscMatrix<f64>