Has anyone written code to solve sparse linear systems (e.g. sparse cholesky ?) In general, how should this problem be tackeld in ublas, apart from running conjugate gradient?<br>Thanks for any insight,<br><br>John<br><br>