func setDiagonal(M *matrix.FloatMatrix, srow, scol, erow, ecol int, val float64) { for i := srow; i < erow; i++ { if i < ecol { M.SetAt(i, i, val) } } }