Merge branch 'libsndfile'
[Faustine.git] / interpretor / lib / src / libsndfile-1.0.25 / tests / fix_this.c
1 /*
2 ** Copyright (C) 1999-2011 Erik de Castro Lopo <erikd@mega-nerd.com>
3 **
4 ** This program is free software; you can redistribute it and/or modify
5 ** it under the terms of the GNU General Public License as published by
6 ** the Free Software Foundation; either version 2 of the License, or
7 ** (at your option) any later version.
8 **
9 ** This program is distributed in the hope that it will be useful,
10 ** but WITHOUT ANY WARRANTY; without even the implied warranty of
11 ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 ** GNU General Public License for more details.
13 **
14 ** You should have received a copy of the GNU General Public License
15 ** along with this program; if not, write to the Free Software
16 ** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
17 */
18
19 #include <stdio.h>
20 #include <stdlib.h>
21 #include <string.h>
22 #include <unistd.h>
23 #include <math.h>
24
25 #include <sndfile.h>
26
27 #include "utils.h"
28
29 #define BUFFER_SIZE (1<<14) /* Should be (1<<14) */
30 #define SAMPLE_RATE (11025)
31
32 #ifndef M_PI
33 #define M_PI 3.14159265358979323846264338
34 #endif
35
36 static void lcomp_test_int (const char *str, const char *filename, int filetype, double margin) ;
37
38 static int error_function (double data, double orig, double margin) ;
39 static int decay_response (int k) ;
40
41 static void gen_signal_double (double *data, double scale, int datalen) ;
42
43 /* Force the start of these buffers to be double aligned. Sparc-solaris will
44 ** choke if they are not.
45 */
46
47 typedef union
48 { double d [BUFFER_SIZE + 1] ;
49 int i [BUFFER_SIZE + 1] ;
50 } BUFFER ;
51
52 static BUFFER data_buffer ;
53 static BUFFER orig_buffer ;
54
55 int
56 main (void)
57 { const char *filename = "test.au" ;
58
59 lcomp_test_int ("au_g721", filename, SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_G721_32, 0.06) ;
60
61 return 0 ;
62 } /* main */
63
64 /*============================================================================================
65 ** Here are the test functions.
66 */
67
68 static void
69 lcomp_test_int (const char *str, const char *filename, int filetype, double margin)
70 { SNDFILE *file ;
71 SF_INFO sfinfo ;
72 int k, m, *orig, *data, sum_abs ;
73 long datalen, seekpos ;
74 double scale ;
75
76 printf ("\nThis is program is not part of the libsndfile test suite.\n\n") ;
77
78 printf (" lcomp_test_int : %s ... ", str) ;
79 fflush (stdout) ;
80
81 datalen = BUFFER_SIZE ;
82
83 scale = 1.0 * 0x10000 ;
84
85 data = data_buffer.i ;
86 orig = orig_buffer.i ;
87
88 gen_signal_double (orig_buffer.d, 32000.0 * scale, datalen) ;
89 for (k = 0 ; k < datalen ; k++)
90 orig [k] = orig_buffer.d [k] ;
91
92
93 sfinfo.samplerate = SAMPLE_RATE ;
94 sfinfo.frames = 123456789 ; /* Ridiculous value. */
95 sfinfo.channels = 1 ;
96 sfinfo.format = filetype ;
97
98 if (! (file = sf_open (filename, SFM_WRITE, &sfinfo)))
99 { printf ("sf_open_write failed with error : ") ;
100 puts (sf_strerror (NULL)) ;
101 exit (1) ;
102 } ;
103
104 if ((k = sf_writef_int (file, orig, datalen)) != datalen)
105 { printf ("sf_writef_int failed with short write (%ld => %d).\n", datalen, k) ;
106 exit (1) ;
107 } ;
108 sf_close (file) ;
109
110 memset (data, 0, datalen * sizeof (int)) ;
111
112 if ((filetype & SF_FORMAT_TYPEMASK) != SF_FORMAT_RAW)
113 memset (&sfinfo, 0, sizeof (sfinfo)) ;
114
115 if (! (file = sf_open (filename, SFM_READ, &sfinfo)))
116 { printf ("sf_open_read failed with error : ") ;
117 puts (sf_strerror (NULL)) ;
118 exit (1) ;
119 } ;
120
121 if ((sfinfo.format & (SF_FORMAT_TYPEMASK | SF_FORMAT_SUBMASK)) != (filetype & (SF_FORMAT_TYPEMASK | SF_FORMAT_SUBMASK)))
122 { printf ("Line %d: Returned format incorrect (0x%08X => 0x%08X).\n", __LINE__, filetype, sfinfo.format) ;
123 exit (1) ;
124 } ;
125
126 if (sfinfo.frames < datalen)
127 { printf ("Too few.frames in file. (%ld should be a little more than %ld)\n", datalen, SF_COUNT_TO_LONG (sfinfo.frames)) ;
128 exit (1) ;
129 } ;
130
131 if (sfinfo.frames > (datalen + datalen / 2))
132 { printf ("Too many.frames in file. (%ld should be a little more than %ld)\n", datalen, SF_COUNT_TO_LONG (sfinfo.frames)) ;
133 exit (1) ;
134 } ;
135
136 if (sfinfo.channels != 1)
137 { printf ("Incorrect number of channels in file.\n") ;
138 exit (1) ;
139 } ;
140
141 check_log_buffer_or_die (file, __LINE__) ;
142
143 if ((k = sf_readf_int (file, data, datalen)) != datalen)
144 { printf ("Line %d: short read (%d should be %ld).\n", __LINE__, k, datalen) ;
145 exit (1) ;
146 } ;
147
148 sum_abs = 0 ;
149 for (k = 0 ; k < datalen ; k++)
150 { if (error_function (data [k] / scale, orig [k] / scale, margin))
151 { printf ("Line %d: Incorrect sample (#%d : %f should be %f).\n", __LINE__, k, data [k] / scale, orig [k] / scale) ;
152 oct_save_int (orig, data, datalen) ;
153 exit (1) ;
154 } ;
155 sum_abs = abs (sum_abs + abs (data [k])) ;
156 } ;
157
158 if (sum_abs < 1.0)
159 { printf ("Line %d: Signal is all zeros.\n", __LINE__) ;
160 exit (1) ;
161 } ;
162
163 if ((k = sf_readf_int (file, data, datalen)) != sfinfo.frames - datalen)
164 { printf ("Line %d: Incorrect read length (%ld should be %d).\n", __LINE__, SF_COUNT_TO_LONG (sfinfo.frames - datalen), k) ;
165 exit (1) ;
166 } ;
167
168 /* This check is only for block based encoders which must append silence
169 ** to the end of a file so as to fill out a block.
170 */
171 if ((sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_MS_ADPCM)
172 for (k = 0 ; k < sfinfo.frames - datalen ; k++)
173 if (abs (data [k] / scale) > decay_response (k))
174 { printf ("Line %d : Incorrect sample B (#%d : abs (%d) should be < %d).\n", __LINE__, k, data [k], decay_response (k)) ;
175 exit (1) ;
176 } ;
177
178 if (! sfinfo.seekable)
179 { printf ("ok\n") ;
180 return ;
181 } ;
182
183 /* Now test sf_seek function. */
184
185 if ((k = sf_seek (file, 0, SEEK_SET)) != 0)
186 { printf ("Line %d: Seek to start of file failed (%d).\n", __LINE__, k) ;
187 exit (1) ;
188 } ;
189
190 for (m = 0 ; m < 3 ; m++)
191 { int n ;
192
193 if ((k = sf_readf_int (file, data, 11)) != 11)
194 { printf ("Line %d: Incorrect read length (11 => %d).\n", __LINE__, k) ;
195 exit (1) ;
196 } ;
197
198 for (k = 0 ; k < 11 ; k++)
199 if (error_function (data [k] / scale, orig [k + m * 11] / scale, margin))
200 { printf ("Line %d: Incorrect sample (m = %d) (#%d : %d => %d).\n", __LINE__, m, k + m * 11, orig [k + m * 11], data [k]) ;
201 for (n = 0 ; n < 1 ; n++)
202 printf ("%d ", data [n]) ;
203 printf ("\n") ;
204 exit (1) ;
205 } ;
206 } ;
207
208 seekpos = BUFFER_SIZE / 10 ;
209
210 /* Check seek from start of file. */
211 if ((k = sf_seek (file, seekpos, SEEK_SET)) != seekpos)
212 { printf ("Seek to start of file + %ld failed (%d).\n", seekpos, k) ;
213 exit (1) ;
214 } ;
215
216 if ((k = sf_readf_int (file, data, 1)) != 1)
217 { printf ("Line %d: sf_readf_int (file, data, 1) returned %d.\n", __LINE__, k) ;
218 exit (1) ;
219 } ;
220
221 if (error_function ((double) data [0], (double) orig [seekpos], margin))
222 { printf ("Line %d: sf_seek (SEEK_SET) followed by sf_readf_int failed (%d, %d).\n", __LINE__, orig [1], data [0]) ;
223 exit (1) ;
224 } ;
225
226 if ((k = sf_seek (file, 0, SEEK_CUR)) != seekpos + 1)
227 { printf ("Line %d: sf_seek (SEEK_CUR) with 0 offset failed (%d should be %ld)\n", __LINE__, k, seekpos + 1) ;
228 exit (1) ;
229 } ;
230
231 seekpos = sf_seek (file, 0, SEEK_CUR) + BUFFER_SIZE / 5 ;
232 k = sf_seek (file, BUFFER_SIZE / 5, SEEK_CUR) ;
233 sf_readf_int (file, data, 1) ;
234 if (error_function ((double) data [0], (double) orig [seekpos], margin) || k != seekpos)
235 { printf ("Line %d: sf_seek (forwards, SEEK_CUR) followed by sf_readf_int failed (%d, %d) (%d, %ld).\n", __LINE__, data [0], orig [seekpos], k, seekpos + 1) ;
236 exit (1) ;
237 } ;
238
239 seekpos = sf_seek (file, 0, SEEK_CUR) - 20 ;
240 /* Check seek backward from current position. */
241 k = sf_seek (file, -20, SEEK_CUR) ;
242 sf_readf_int (file, data, 1) ;
243 if (error_function ((double) data [0], (double) orig [seekpos], margin) || k != seekpos)
244 { printf ("sf_seek (backwards, SEEK_CUR) followed by sf_readf_int failed (%d, %d) (%d, %ld).\n", data [0], orig [seekpos], k, seekpos) ;
245 exit (1) ;
246 } ;
247
248 /* Check that read past end of file returns number of items. */
249 sf_seek (file, (int) sfinfo.frames, SEEK_SET) ;
250
251 if ((k = sf_readf_int (file, data, datalen)) != 0)
252 { printf ("Line %d: Return value from sf_readf_int past end of file incorrect (%d).\n", __LINE__, k) ;
253 exit (1) ;
254 } ;
255
256 /* Check seek backward from end. */
257 if ((k = sf_seek (file, 5 - (int) sfinfo.frames, SEEK_END)) != 5)
258 { printf ("sf_seek (SEEK_END) returned %d instead of %d.\n", k, 5) ;
259 exit (1) ;
260 } ;
261
262 sf_readf_int (file, data, 1) ;
263 if (error_function (data [0] / scale, orig [5] / scale, margin))
264 { printf ("Line %d: sf_seek (SEEK_END) followed by sf_readf_short failed (%d should be %d).\n", __LINE__, data [0], orig [5]) ;
265 exit (1) ;
266 } ;
267
268 sf_close (file) ;
269
270 printf ("ok\n") ;
271 } /* lcomp_test_int */
272
273 /*========================================================================================
274 ** Auxiliary functions
275 */
276
277 #define SIGNAL_MAXVAL 30000.0
278 #define DECAY_COUNT 800
279
280 static int
281 decay_response (int k)
282 { if (k < 1)
283 return (int) (1.2 * SIGNAL_MAXVAL) ;
284 if (k > DECAY_COUNT)
285 return 0 ;
286 return (int) (1.2 * SIGNAL_MAXVAL * (DECAY_COUNT - k) / (1.0 * DECAY_COUNT)) ;
287 } /* decay_response */
288
289 static void
290 gen_signal_double (double *data, double scale, int datalen)
291 { int k, ramplen ;
292 double amp = 0.0 ;
293
294 ramplen = datalen / 18 ;
295
296 for (k = 0 ; k < datalen ; k++)
297 { if (k <= ramplen)
298 amp = scale * k / ((double) ramplen) ;
299 else if (k > datalen - ramplen)
300 amp = scale * (datalen - k) / ((double) ramplen) ;
301
302 data [k] = amp * (0.4 * sin (33.3 * 2.0 * M_PI * ((double) (k + 1)) / ((double) SAMPLE_RATE))
303 + 0.3 * cos (201.1 * 2.0 * M_PI * ((double) (k + 1)) / ((double) SAMPLE_RATE))) ;
304 } ;
305
306 return ;
307 } /* gen_signal_double */
308
309 static int
310 error_function (double data, double orig, double margin)
311 { double error ;
312
313 if (fabs (orig) <= 500.0)
314 error = fabs (fabs (data) - fabs (orig)) / 2000.0 ;
315 else if (fabs (orig) <= 1000.0)
316 error = fabs (data - orig) / 3000.0 ;
317 else
318 error = fabs (data - orig) / fabs (orig) ;
319
320 if (error > margin)
321 { printf ("\n\n*******************\nError : %f\n", error) ;
322 return 1 ;
323 } ;
324 return 0 ;
325 } /* error_function */
326