use simple one-state-at-a-time-in-specified-order occupation scheme for complex scaling.  Add inputfile keywords for penalizing energies when ordering, and to rotate the spectrum when running the arpack eigensolver so that states can be ordered in other ways than smallest real, largest imaginary and so on.  Make half-hartree exchange available as a separate functional (only with complex scaling so far)
correct size of arrays for arnoldivectors in parallel.  Fiddle more with penalization of positive eigenvalues.  Various minor cleanups
merge git branches and resolve lots of conflicts