TODO

   1 Convert pol to C++ object.
   2
   3 Add divergent filtering to frequency/direct filtering
   4
   5 FFT filtering requires FFTW library. Add a basic FFT routine.
   6