starting to work on Jacobian and psisq. However, scatters don't seem to work the way I think they do, as illustrated by the redundancies in VecCollect. Study them more before you continue