Add filtering for phase and magnitude in phase retrieval function

[david-hoffman]

This will need phase unwrapping at each iteration, because magnitude and phase need to be smoothed independently, I think.

Needs some testing to see if raw complex pupil can be smoothed directly ...