Simplify this: https://github.com/poldracklab/fmriprep/blob/9c57dded04c7c0f26768264afaaee1d4aa158b93/fmriprep/workflows/fieldmap/base.py#L202-L241 At this time, ``init_fmap_wf`` and ``init_phdiff_wf`` are very similar. The idea would be to implement: - [x] Magnitude registration; and - [x] regularization and displacements field calculation only once on ``init_fmap_wf`` from a fieldmap in Hz. Phase1/2 and phasediff would then focus on phase unwrapping and conversion to Hz.