Conjugate Gradient iterative solver

std::vector<double> cgsolve(const CooMatrix<double> &A, const std::vector<double> &b)