Implement "inhomogeneous arrays" for `ldos` and `wavefunc`
The point is to have ldos
or wavefunc
return objects that behave like 2d arrays, where the site index and orbital index are separate. However, this must work in the general case where the number of orbitals per site varies.
Depends on #7 (closed) and #6 (closed).