Spiral readout trajectories have become popular in several applications of 2D and 3D MR imaging applications because of their ability to acquire more data in less time than their traditional Cartesian counterparts. 3D sampling is especially desirable for high-resolution imaging, which results in the combined burden of the (i) complexity inherent to non-Cartesian reconstruction, (ii) large data matrices for high resolution, and (iii) the need to process the entire 3D dataset in one go (instead of slice-by-slice as in 2D). The consequences are high computational cost and memory requirements during the reconstruction process, especially in combination with parallel imaging and off-resonance correction.