User Portlet
| Discussions |
|---|
| In case it's helpful, the code below seems to trigger the same issue. L = 5; n = 4; id = ConstantArray[1, L] // SparseArray // DiagonalMatrix; D0 = SparseArray[{{i_, j_} /; Mod[j - i, L] == 1 :> 1, {i_, j_} /; Mod[j - i, L] == 0... |