X-Git-Url: https://git.stderr.nl/gitweb?a=blobdiff_plain;f=FFT_support.cpp;h=d3b5d09d28cc0917cd68106f5db53c7b7a53d916;hb=81823b7486c210898290292ed87142a89c496946;hp=23d07f08ca102e114e7428807701f364fcdcc566;hpb=a62dbd706365d44bdcac9f8ac4c7c6dd05484642;p=matthijs%2Fprojects%2Fmontium-fft.git

diff --git a/FFT_support.cpp b/FFT_support.cpp
index 23d07f0..d3b5d09 100644
--- a/FFT_support.cpp
+++ b/FFT_support.cpp
@@ -68,7 +68,9 @@ void pre_run()
 	
 	for (i=0;i<PARAM_N_t;i++)
 	{
-		int value = to_fixed(sin((float)i*2*M_PI/PARAM_N_t));
+		/* We take the sine from 0 to 2*2*Pi, ie two periods. We divide 
+		 * the value by PARAM_N_t to prevent overflow. */
+		int value = to_fixed(sin((float)i*2*2*M_PI/PARAM_N_t)/PARAM_N_t);
 
 		if (i<PARAM_N_t/2)
 		{