/*


  convolution using fftw

*/




//  void do_convolution( double* image, double* convolution, double* dest, int n );