TODO

   1 Fix BUGS -- see BUGS file
   2
   3 Convert pol to C++ object.
   4
   5 Add divergent filtering to frequency/direct filtering
   6
   7 FFT filtering requires FFTW library. Add a public domain FFT routine.
   8