compare for h = ndplex(5, edges(3)) - h.diag2((0,0,0),(1,1)) works - h.diag2((0,0),(1,1,1)) does not work
compare for h = ndplex(5, edges(3))