bugfix: invlaid memory access in dm_update

Src/m_ts_dm_update.F90

@@ -510,11 +510,13 @@ contains
        dE => val(spEDM)
 
        ! the actual size of the shift
-       jo = nnzs(spDM)
+         jo = nnzs(spDM)
 
        ! As we have shifted the fermi-level up to 0, we need to shift the
        ! energy-density matrix back
-       call daxpy(jo,Ef,dD(1,1),1,dE(1,1),1)
+       if (jo > 0) then
+         call daxpy(jo,Ef,dD(1,1),1,dE(1,1),1)
+       end if
 
     end if