improve MAD SNR estimation

[n-west]

MAD (Median Absolute Deviation) provides a good alternative to SNR estimation. Right now it works on CPU for non-masked SNR. This is largely due to the median implementation being a relatively simple C++-container based method. Implement median for the general case and get masked medians + cuda version running