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