Function std::conv_fft
int conv_fft([]const double x, []const double h, []double out)FFT-accelerated linear convolution — same output as conv_direct (up to floating-point rounding), computed via zero-padded forward FFTs, pointwise complex multiply, and one inverse FFT. out must have space for (len_x + len_h - 1) samples.