Miscellaneous progress

* fixed a bug in the Sparse(SparseMap&) constructor where the last line
wasn't taken over properly
* many other changes preparing the solution of the Helmholtz equation,
we're not there yet but better commit WIP code than being sorry for data
loss later!