It would be nice if SPARC supported some way to start a calculation from the results of the previous, for example when using an ASE dynamics algorithm (MD, relaxation, NEB, whatever). The simplest approach would be to save and reload wavefunctions, but another (arguably better) way would be to have the sparc executable work in a sort of client-server mode, where it waited for new input positions/cell passed from the ASE calculator somehow. I've done things like that with file based I/O for VASP, so if that seems useful, I'd be happy to help if I can.