Additional fixes for Laplacian::solve

SPT and PDD implementations have their own solve(Field3D) functions,
which need to set the cell location of the return value.