Function laplacian_unwrap_fft_mirrored 

fn laplacian_unwrap_fft_mirrored(phase: &[f64], grid: &Grid) -> Vec<f64>

Schofield & Zhu sin/cos unwrapping on the volume mirrored to 2n along each axis (lapunwrap.m). Returns the whole volume, unmasked.