The time series is padded with zeroes to avoid problems with
circular versus linear convolution. This does not affect usage,
as the matrix returned has the added columns removed. (JML Sep 29, 2021).
The output contains the (complex) values of the wavelet transform of the
input signal. The format of the output can be
2D array (signal_size x nb_scales)
3D array (signal_size x noctave x nvoice)
Since Morlet's wavelet is not strictly speaking a wavelet (it is not of
vanishing integral), artifacts may occur for certain signals.