Procedure from library linalg.lib (see linalg_lib).
linalg.lib
evnf(e[,m]); ideal e, intvec m
ncols(e)==size(m)
order eigenvalues e with multiplicities m