window: fabs > fabsf and cos > cosf where it makes sense
authorRosen Penev <rosenp@gmail.com>
Thu, 29 Aug 2019 21:40:57 +0000 (14:40 -0700)
committerErik de Castro Lopo <erikd@mega-nerd.com>
Sun, 15 Sep 2019 01:34:08 +0000 (11:34 +1000)
The f functions are faster. It seems the arguments are all floats.

src/libFLAC/window.c

index 7698328..94495f3 100644 (file)
@@ -68,7 +68,7 @@ void FLAC__window_bartlett_hann(FLAC__real *window, const FLAC__int32 L)
        FLAC__int32 n;
 
        for (n = 0; n < L; n++)
-               window[n] = (FLAC__real)(0.62f - 0.48f * fabs((float)n/(float)N-0.5f) - 0.38f * cos(2.0f * M_PI * ((float)n/(float)N)));
+               window[n] = (FLAC__real)(0.62f - 0.48f * fabsf((float)n/(float)N-0.5f) - 0.38f * cosf(2.0f * M_PI * ((float)n/(float)N)));
 }
 
 void FLAC__window_blackman(FLAC__real *window, const FLAC__int32 L)
@@ -77,7 +77,7 @@ void FLAC__window_blackman(FLAC__real *window, const FLAC__int32 L)
        FLAC__int32 n;
 
        for (n = 0; n < L; n++)
-               window[n] = (FLAC__real)(0.42f - 0.5f * cos(2.0f * M_PI * n / N) + 0.08f * cos(4.0f * M_PI * n / N));
+               window[n] = (FLAC__real)(0.42f - 0.5f * cosf(2.0f * M_PI * n / N) + 0.08f * cosf(4.0f * M_PI * n / N));
 }
 
 /* 4-term -92dB side-lobe */
@@ -87,7 +87,7 @@ void FLAC__window_blackman_harris_4term_92db_sidelobe(FLAC__real *window, const
        FLAC__int32 n;
 
        for (n = 0; n <= N; n++)
-               window[n] = (FLAC__real)(0.35875f - 0.48829f * cos(2.0f * M_PI * n / N) + 0.14128f * cos(4.0f * M_PI * n / N) - 0.01168f * cos(6.0f * M_PI * n / N));
+               window[n] = (FLAC__real)(0.35875f - 0.48829f * cosf(2.0f * M_PI * n / N) + 0.14128f * cosf(4.0f * M_PI * n / N) - 0.01168f * cosf(6.0f * M_PI * n / N));
 }
 
 void FLAC__window_connes(FLAC__real *window, const FLAC__int32 L)
@@ -109,7 +109,7 @@ void FLAC__window_flattop(FLAC__real *window, const FLAC__int32 L)
        FLAC__int32 n;
 
        for (n = 0; n < L; n++)
-               window[n] = (FLAC__real)(0.21557895f - 0.41663158f * cos(2.0f * M_PI * n / N) + 0.277263158f * cos(4.0f * M_PI * n / N) - 0.083578947f * cos(6.0f * M_PI * n / N) + 0.006947368f * cos(8.0f * M_PI * n / N));
+               window[n] = (FLAC__real)(0.21557895f - 0.41663158f * cosf(2.0f * M_PI * n / N) + 0.277263158f * cosf(4.0f * M_PI * n / N) - 0.083578947f * cosf(6.0f * M_PI * n / N) + 0.006947368f * cosf(8.0f * M_PI * n / N));
 }
 
 void FLAC__window_gauss(FLAC__real *window, const FLAC__int32 L, const FLAC__real stddev)
@@ -130,7 +130,7 @@ void FLAC__window_hamming(FLAC__real *window, const FLAC__int32 L)
        FLAC__int32 n;
 
        for (n = 0; n < L; n++)
-               window[n] = (FLAC__real)(0.54f - 0.46f * cos(2.0f * M_PI * n / N));
+               window[n] = (FLAC__real)(0.54f - 0.46f * cosf(2.0f * M_PI * n / N));
 }
 
 void FLAC__window_hann(FLAC__real *window, const FLAC__int32 L)
@@ -139,7 +139,7 @@ void FLAC__window_hann(FLAC__real *window, const FLAC__int32 L)
        FLAC__int32 n;
 
        for (n = 0; n < L; n++)
-               window[n] = (FLAC__real)(0.5f - 0.5f * cos(2.0f * M_PI * n / N));
+               window[n] = (FLAC__real)(0.5f - 0.5f * cosf(2.0f * M_PI * n / N));
 }
 
 void FLAC__window_kaiser_bessel(FLAC__real *window, const FLAC__int32 L)
@@ -148,7 +148,7 @@ void FLAC__window_kaiser_bessel(FLAC__real *window, const FLAC__int32 L)
        FLAC__int32 n;
 
        for (n = 0; n < L; n++)
-               window[n] = (FLAC__real)(0.402f - 0.498f * cos(2.0f * M_PI * n / N) + 0.098f * cos(4.0f * M_PI * n / N) - 0.001f * cos(6.0f * M_PI * n / N));
+               window[n] = (FLAC__real)(0.402f - 0.498f * cosf(2.0f * M_PI * n / N) + 0.098f * cosf(4.0f * M_PI * n / N) - 0.001f * cosf(6.0f * M_PI * n / N));
 }
 
 void FLAC__window_nuttall(FLAC__real *window, const FLAC__int32 L)
@@ -157,7 +157,7 @@ void FLAC__window_nuttall(FLAC__real *window, const FLAC__int32 L)
        FLAC__int32 n;
 
        for (n = 0; n < L; n++)
-               window[n] = (FLAC__real)(0.3635819f - 0.4891775f*cos(2.0f*M_PI*n/N) + 0.1365995f*cos(4.0f*M_PI*n/N) - 0.0106411f*cos(6.0f*M_PI*n/N));
+               window[n] = (FLAC__real)(0.3635819f - 0.4891775f*cosf(2.0f*M_PI*n/N) + 0.1365995f*cosf(4.0f*M_PI*n/N) - 0.0106411f*cosf(6.0f*M_PI*n/N));
 }
 
 void FLAC__window_rectangle(FLAC__real *window, const FLAC__int32 L)
@@ -200,8 +200,8 @@ void FLAC__window_tukey(FLAC__real *window, const FLAC__int32 L, const FLAC__rea
                /* ...replace ends with hann */
                if (Np > 0) {
                        for (n = 0; n <= Np; n++) {
-                               window[n] = (FLAC__real)(0.5f - 0.5f * cos(M_PI * n / Np));
-                               window[L-Np-1+n] = (FLAC__real)(0.5f - 0.5f * cos(M_PI * (n+Np) / Np));
+                               window[n] = (FLAC__real)(0.5f - 0.5f * cosf(M_PI * n / Np));
+                               window[L-Np-1+n] = (FLAC__real)(0.5f - 0.5f * cosf(M_PI * (n+Np) / Np));
                        }
                }
        }
@@ -225,11 +225,11 @@ void FLAC__window_partial_tukey(FLAC__real *window, const FLAC__int32 L, const F
                for (n = 0; n < start_n && n < L; n++)
                        window[n] = 0.0f;
                for (i = 1; n < (start_n+Np) && n < L; n++, i++)
-                       window[n] = (FLAC__real)(0.5f - 0.5f * cos(M_PI * i / Np));
+                       window[n] = (FLAC__real)(0.5f - 0.5f * cosf(M_PI * i / Np));
                for (; n < (end_n-Np) && n < L; n++)
                        window[n] = 1.0f;
                for (i = Np; n < end_n && n < L; n++, i--)
-                       window[n] = (FLAC__real)(0.5f - 0.5f * cos(M_PI * i / Np));
+                       window[n] = (FLAC__real)(0.5f - 0.5f * cosf(M_PI * i / Np));
                for (; n < L; n++)
                        window[n] = 0.0f;
        }
@@ -251,19 +251,19 @@ void FLAC__window_punchout_tukey(FLAC__real *window, const FLAC__int32 L, const
                Ne = (FLAC__int32)(p / 2.0f * (L - end_n));
 
                for (n = 0, i = 1; n < Ns && n < L; n++, i++)
-                       window[n] = (FLAC__real)(0.5f - 0.5f * cos(M_PI * i / Ns));
+                       window[n] = (FLAC__real)(0.5f - 0.5f * cosf(M_PI * i / Ns));
                for (; n < start_n-Ns && n < L; n++)
                        window[n] = 1.0f;
                for (i = Ns; n < start_n && n < L; n++, i--)
-                       window[n] = (FLAC__real)(0.5f - 0.5f * cos(M_PI * i / Ns));
+                       window[n] = (FLAC__real)(0.5f - 0.5f * cosf(M_PI * i / Ns));
                for (; n < end_n && n < L; n++)
                        window[n] = 0.0f;
                for (i = 1; n < end_n+Ne && n < L; n++, i++)
-                       window[n] = (FLAC__real)(0.5f - 0.5f * cos(M_PI * i / Ne));
+                       window[n] = (FLAC__real)(0.5f - 0.5f * cosf(M_PI * i / Ne));
                for (; n < L - (Ne) && n < L; n++)
                        window[n] = 1.0f;
                for (i = Ne; n < L; n++, i--)
-                       window[n] = (FLAC__real)(0.5f - 0.5f * cos(M_PI * i / Ne));
+                       window[n] = (FLAC__real)(0.5f - 0.5f * cosf(M_PI * i / Ne));
        }
 }