mat: triangular matrices don't respect the .mat.Diag field
Created by: kortschak
The .mat.Diag
field of the triangular types is used by nearly none of the operations in mat (one internal method checks it to panic and the new diagonal extraction method on Diagonal
allows unit diagonals to be obtained (FWIW). However indexing a TriDense
for example does not check t.mat.Diag
. This should probably be done (although there is currently no way for a user to set the state to unit, that may happen).