Added functions to create and smear sources.
Fixed bug in sparse matrix multiplication in smearing operator generation