bump version number to 1.1.3
[flac.git] / src / libFLAC / window.c
1 /* libFLAC - Free Lossless Audio Codec library
2  * Copyright (C) 2006  Josh Coalson
3  *
4  * Redistribution and use in source and binary forms, with or without
5  * modification, are permitted provided that the following conditions
6  * are met:
7  *
8  * - Redistributions of source code must retain the above copyright
9  * notice, this list of conditions and the following disclaimer.
10  *
11  * - Redistributions in binary form must reproduce the above copyright
12  * notice, this list of conditions and the following disclaimer in the
13  * documentation and/or other materials provided with the distribution.
14  *
15  * - Neither the name of the Xiph.org Foundation nor the names of its
16  * contributors may be used to endorse or promote products derived from
17  * this software without specific prior written permission.
18  *
19  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
20  * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
21  * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
22  * A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR
23  * CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
24  * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
25  * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
26  * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
27  * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
28  * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
29  * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
30  */
31
32 #if HAVE_CONFIG_H
33 #  include <config.h>
34 #endif
35
36 #include <math.h>
37 #include "FLAC/assert.h"
38 #include "FLAC/format.h"
39 #include "private/window.h"
40
41 #ifndef FLAC__INTEGER_ONLY_LIBRARY
42
43 #ifndef M_PI
44 /* math.h in VC++ doesn't seem to have this (how Microsoft is that?) */
45 #define M_PI 3.14159265358979323846
46 #endif
47
48
49 void FLAC__window_bartlett(FLAC__real *window, const FLAC__int32 L)
50 {
51         const FLAC__int32 N = L - 1;
52         FLAC__int32 n;
53
54         if (L & 1) {
55                 for (n = 0; n <= N/2; n++)
56                         window[n] = 2.0f * n / (float)N;
57                 for (; n <= N; n++)
58                         window[n] = 2.0f - 2.0f * n / (float)N;
59         }
60         else {
61                 for (n = 0; n <= L/2-1; n++)
62                         window[n] = 2.0f * n / (float)N;
63                 for (; n <= N; n++)
64                         window[n] = 2.0f - 2.0f * (N-n) / (float)N;
65         }
66 }
67
68 void FLAC__window_bartlett_hann(FLAC__real *window, const FLAC__int32 L)
69 {
70         const FLAC__int32 N = L - 1;
71         FLAC__int32 n;
72
73         for (n = 0; n < L; n++)
74                 window[n] = 0.62f - 0.48f * fabs((float)n/(float)N+0.5f) + 0.38f * cos(2.0f * M_PI * ((float)n/(float)N+0.5f));
75 }
76
77 void FLAC__window_blackman(FLAC__real *window, const FLAC__int32 L)
78 {
79         const FLAC__int32 N = L - 1;
80         FLAC__int32 n;
81
82         for (n = 0; n < L; n++)
83                 window[n] = 0.42f - 0.5f * cos(2.0f * M_PI * n / N) + 0.08f * cos(4.0f * M_PI * n / N);
84 }
85
86 /* 4-term -92dB side-lobe */
87 void FLAC__window_blackman_harris_4term_92db_sidelobe(FLAC__real *window, const FLAC__int32 L)
88 {
89         const FLAC__int32 N = L - 1;
90         FLAC__int32 n;
91
92         for (n = 0; n <= N; n++)
93                 window[n] = 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);
94 }
95
96 void FLAC__window_connes(FLAC__real *window, const FLAC__int32 L)
97 {
98         const FLAC__int32 N = L - 1;
99         const double N2 = (double)N / 2.;
100         FLAC__int32 n;
101
102         for (n = 0; n <= N; n++) {
103                 double k = ((double)n - N2) / N2;
104                 k = 1.0f - k * k;
105                 window[n] = k * k;
106         }
107 }
108
109 void FLAC__window_flattop(FLAC__real *window, const FLAC__int32 L)
110 {
111         const FLAC__int32 N = L - 1;
112         FLAC__int32 n;
113
114         for (n = 0; n < L; n++)
115                 window[n] = 1.0f - 1.93f * cos(2.0f * M_PI * n / N) + 1.29f * cos(4.0f * M_PI * n / N) - 0.388f * cos(6.0f * M_PI * n / N) + 0.0322f * cos(8.0f * M_PI * n / N);
116 }
117
118 void FLAC__window_gauss(FLAC__real *window, const FLAC__int32 L, const FLAC__real stddev)
119 {
120         const FLAC__int32 N = L - 1;
121         const double N2 = (double)N / 2.;
122         FLAC__int32 n;
123
124         for (n = 0; n <= N; n++) {
125                 const double k = ((double)n - N2) / (stddev * N2);
126                 window[n] = exp(-0.5f * k * k);
127         }
128 }
129
130 void FLAC__window_hamming(FLAC__real *window, const FLAC__int32 L)
131 {
132         const FLAC__int32 N = L - 1;
133         FLAC__int32 n;
134
135         for (n = 0; n < L; n++)
136                 window[n] = 0.54f - 0.46f * cos(2.0f * M_PI * n / N);
137 }
138
139 void FLAC__window_hann(FLAC__real *window, const FLAC__int32 L)
140 {
141         const FLAC__int32 N = L - 1;
142         FLAC__int32 n;
143
144         for (n = 0; n < L; n++)
145                 window[n] = 0.5f - 0.5f * cos(2.0f * M_PI * n / N);
146 }
147
148 void FLAC__window_kaiser_bessel(FLAC__real *window, const FLAC__int32 L)
149 {
150         const FLAC__int32 N = L - 1;
151         FLAC__int32 n;
152
153         for (n = 0; n < L; n++)
154                 window[n] = 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);
155 }
156
157 void FLAC__window_nuttall(FLAC__real *window, const FLAC__int32 L)
158 {
159         const FLAC__int32 N = L - 1;
160         FLAC__int32 n;
161
162         for (n = 0; n < L; n++)
163                 window[n] = 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);
164 }
165
166 void FLAC__window_rectangle(FLAC__real *window, const FLAC__int32 L)
167 {
168         FLAC__int32 n;
169
170         for (n = 0; n < L; n++)
171                 window[n] = 1.0f;
172 }
173
174 void FLAC__window_triangle(FLAC__real *window, const FLAC__int32 L)
175 {
176         FLAC__int32 n;
177
178         if (L & 1) {
179                 for (n = 1; n <= L+1/2; n++)
180                         window[n-1] = 2.0f * n / ((float)L + 1.0f);
181                 for (; n <= L; n++)
182                         window[n-1] = - (float)(2 * (L - n + 1)) / ((float)L + 1.0f);
183         }
184         else {
185                 for (n = 1; n <= L/2; n++)
186                         window[n-1] = 2.0f * n / (float)L;
187                 for (; n <= L; n++)
188                         window[n-1] = ((float)(2 * (L - n)) + 1.0f) / (float)L;
189         }
190 }
191
192 void FLAC__window_tukey(FLAC__real *window, const FLAC__int32 L, const FLAC__real p)
193 {
194         if (p <= 0.0)
195                 FLAC__window_rectangle(window, L);
196         else if (p >= 1.0)
197                 FLAC__window_hann(window, L);
198         else {
199                 const FLAC__int32 Np = (FLAC__int32)(p / 2.0f * L) - 1;
200                 FLAC__int32 n;
201                 /* start with rectangle... */
202                 FLAC__window_rectangle(window, L);
203                 /* ...replace ends with hann */
204                 if (Np > 0) {
205                         for (n = 0; n <= Np; n++) {
206                                 window[n] = 0.5f - 0.5f * cos(M_PI * n / Np);
207                                 window[L-Np-1+n] = 0.5f - 0.5f * cos(M_PI * (n+Np) / Np);
208                         }
209                 }
210         }
211 }
212
213 void FLAC__window_welch(FLAC__real *window, const FLAC__int32 L)
214 {
215         const FLAC__int32 N = L - 1;
216         const double N2 = (double)N / 2.;
217         FLAC__int32 n;
218
219         for (n = 0; n <= N; n++) {
220                 const double k = ((double)n - N2) / N2;
221                 window[n] = 1.0f - k * k;
222         }
223 }
224
225 #endif /* !defined FLAC__INTEGER_ONLY_LIBRARY */