Minor debugging modifications
Included a standalone solver instead of Eigen

Classical Jacobi preconditioned conjugate gradient
Minimum number of iterations specified in order to amortize fixed costs
Relative improvement as termination condition
Two types of matrices (CSR/Ellpack)