New input option `restartdir`:  the program will save Hd coefficients to `$restartdir/hd.data.$iter` every iteration if this option is nonempty
- Fixed a problem where the point indices are not correctly printed in `surfgen.out` when some ab initio are absent