Use spdiags instead of diags

Diags is a new scipy function?  Not found in the version on travis for testing.