CG: Pre-allocation of A*x and using A_mul_B! for sparse matrices for speedup
CG: Trying to use Blas axpy! and scal! routines for speed improvement