ADDED: Created DIA sparse matrix-vector product

Since most Hamiltonians and Louvillians are banded matricies, taking
the spmv in DIA format may be faster.  Still need to check.