This fixes the error issues found by fast_complementary_polynomial by using the algorithm described in Berntson and Sunderhauf (2024) (Algorithm 2 based off the code in Figure 1). This method takes the target precision as an input parameter and automatically determines the number of FFT modes needed to guarantee the given precision.
Added fft_qsp module to address Issue 1076
This fixes the error issues found by fast_complementary_polynomial by using the algorithm described in Berntson and Sunderhauf (2024) (Algorithm 2 based off the code in Figure 1). This method takes the target precision as an input parameter and automatically determines the number of FFT modes needed to guarantee the given precision.