길이 n인 real wavelet1과
길이 m인 real wavelet2를
각각 complex로 만들어준 후 길이 m+n-1이 되도록
Zero padding을 해준다.

각각의 wavelet들을 FFT한 후에

k=m+n-1
output(1:k)=cwave1(1:k)*cwave2(1:k) !곱하기

그리고 나서 inverse Fourier transform 하면
Time domain에서의 convolution 결과와 같게 나온다.

01:   L=m+n-1
02:   cw1(1:n)=cmplx(wave1(1:n),0.)
03:   cw1(n+1:L)=cmplx(0.,0.)
04:   
05:   cw2(1:m)=cmplx(wave2(1:m),0.)
06:   cw2(m+1:L)=cmplx(0.,0.)
07:   
08:   call fft(cw1,L)
09:   call fft(cw2,L)
10:   
11:   conv(1:L)=cw1(1:L)*cw2(1:L)
12:   
13:   call ift(conv,L)