Fix BUGS -- see BUGS file

Convert pol to C++ object.

Add divergent filtering to frequency/direct filtering

FFT filtering requires FFTW library. Add a public domain FFT routine.