cannam@85: /* cannam@85: ** Copyright (C) 1999-2011 Erik de Castro Lopo cannam@85: ** cannam@85: ** This program is free software; you can redistribute it and/or modify cannam@85: ** it under the terms of the GNU General Public License as published by cannam@85: ** the Free Software Foundation; either version 2 of the License, or cannam@85: ** (at your option) any later version. cannam@85: ** cannam@85: ** This program is distributed in the hope that it will be useful, cannam@85: ** but WITHOUT ANY WARRANTY; without even the implied warranty of cannam@85: ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the cannam@85: ** GNU General Public License for more details. cannam@85: ** cannam@85: ** You should have received a copy of the GNU General Public License cannam@85: ** along with this program; if not, write to the Free Software cannam@85: ** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. cannam@85: */ cannam@85: cannam@85: #include "sfconfig.h" cannam@85: cannam@85: #include cannam@85: #include cannam@85: #include cannam@85: #include cannam@85: cannam@85: #if HAVE_UNISTD_H cannam@85: #include cannam@85: #endif cannam@85: cannam@85: #include cannam@85: cannam@85: #include "utils.h" cannam@85: cannam@85: #define BUFFER_SIZE (1<<14) /* Should be (1<<14) */ cannam@85: #define SAMPLE_RATE 11025 cannam@85: cannam@85: #ifndef M_PI cannam@85: #define M_PI 3.14159265358979323846264338 cannam@85: #endif cannam@85: cannam@85: #define LCT_MAX(x,y) ((x) > (y) ? (x) : (y)) cannam@85: cannam@85: static void lcomp_test_short (const char *filename, int filetype, int chan, double margin) ; cannam@85: static void lcomp_test_int (const char *filename, int filetype, int chan, double margin) ; cannam@85: static void lcomp_test_float (const char *filename, int filetype, int chan, double margin) ; cannam@85: static void lcomp_test_double (const char *filename, int filetype, int chan, double margin) ; cannam@85: cannam@85: static void sdlcomp_test_short (const char *filename, int filetype, int chan, double margin) ; cannam@85: static void sdlcomp_test_int (const char *filename, int filetype, int chan, double margin) ; cannam@85: static void sdlcomp_test_float (const char *filename, int filetype, int chan, double margin) ; cannam@85: static void sdlcomp_test_double (const char *filename, int filetype, int chan, double margin) ; cannam@85: cannam@85: static void read_raw_test (const char *filename, int filetype, int chan) ; cannam@85: cannam@85: static int error_function (double data, double orig, double margin) ; cannam@85: static int decay_response (int k) ; cannam@85: cannam@85: static void gen_signal_double (double *data, double scale, int channels, int datalen) ; cannam@85: cannam@85: static void smoothed_diff_short (short *data, unsigned int datalen) ; cannam@85: static void smoothed_diff_int (int *data, unsigned int datalen) ; cannam@85: static void smoothed_diff_float (float *data, unsigned int datalen) ; cannam@85: static void smoothed_diff_double (double *data, unsigned int datalen) ; cannam@85: cannam@85: static void check_comment (SNDFILE * file, int format, int lineno) ; cannam@85: cannam@85: static int is_lossy (int filetype) ; cannam@85: cannam@85: /* cannam@85: ** Force the start of these buffers to be double aligned. Sparc-solaris will cannam@85: ** choke if they are not. cannam@85: */ cannam@85: typedef union cannam@85: { double d [BUFFER_SIZE + 1] ; cannam@85: float f [BUFFER_SIZE + 1] ; cannam@85: int i [BUFFER_SIZE + 1] ; cannam@85: short s [BUFFER_SIZE + 1] ; cannam@85: char c [BUFFER_SIZE + 1] ; cannam@85: } BUFFER ; cannam@85: cannam@85: static BUFFER data_buffer ; cannam@85: static BUFFER orig_buffer ; cannam@85: static BUFFER smooth_buffer ; cannam@85: cannam@85: static const char *long_comment = cannam@85: "This is really quite a long comment. It is designed to be long enough " cannam@85: "to screw up the encoders and decoders if the file container format does " cannam@85: "not handle things correctly. If everything is working correctly, the " cannam@85: "decoder will only decode the actual audio data, and not this string at " cannam@85: "the end of the file." ; cannam@85: cannam@85: int cannam@85: main (int argc, char *argv []) cannam@85: { int do_all = 0 ; cannam@85: int test_count = 0 ; cannam@85: cannam@85: if (argc != 2) cannam@85: { printf ("Usage : %s \n", argv [0]) ; cannam@85: printf (" Where is one of the following:\n") ; cannam@85: printf (" wav_ima - test IMA ADPCM WAV file functions\n") ; cannam@85: printf (" wav_msadpcm - test MS ADPCM WAV file functions\n") ; cannam@85: printf (" wav_gsm610 - test GSM 6.10 WAV file functions\n") ; cannam@85: printf (" wav_ulaw - test u-law WAV file functions\n") ; cannam@85: printf (" wav_alaw - test A-law WAV file functions\n") ; cannam@85: printf (" wve - test Psion WVE file functions\n") ; cannam@85: printf (" all - perform all tests\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: do_all = ! strcmp (argv [1], "all") ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "wav_pcm") == 0) cannam@85: { /* This is just a sanity test for PCM encoding. */ cannam@85: lcomp_test_short ("pcm.wav", SF_FORMAT_WAV | SF_FORMAT_PCM_16, 2, 1e-50) ; cannam@85: lcomp_test_int ("pcm.wav", SF_FORMAT_WAV | SF_FORMAT_PCM_32, 2, 1e-50) ; cannam@85: lcomp_test_short ("pcm.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_PCM_16, 2, 1e-50) ; cannam@85: lcomp_test_int ("pcm.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_PCM_32, 2, 1e-50) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("pcm.wav", SF_FORMAT_WAV | SF_FORMAT_FLOAT, 2, 1e-50) ; cannam@85: lcomp_test_double ("pcm.wav", SF_FORMAT_WAV | SF_FORMAT_DOUBLE, 2, 1e-50) ; cannam@85: /* Lite remove end */ cannam@85: cannam@85: read_raw_test ("pcm.wav", SF_FORMAT_WAV | SF_FORMAT_PCM_U8, 2) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: /* For all the rest, if the file format supports more than 1 channel, use stereo. */ cannam@85: /* Lite remove start */ cannam@85: if (do_all || strcmp (argv [1], "wav_ima") == 0) cannam@85: { lcomp_test_short ("ima.wav", SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: lcomp_test_int ("ima.wav", SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, 2, 0.65) ; cannam@85: lcomp_test_float ("ima.wav", SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: lcomp_test_double ("ima.wav", SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: cannam@85: lcomp_test_short ("ima.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: lcomp_test_int ("ima.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: lcomp_test_float ("ima.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: lcomp_test_double ("ima.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: cannam@85: sdlcomp_test_short ("ima.wav", SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: sdlcomp_test_int ("ima.wav", SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: sdlcomp_test_float ("ima.wav", SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: sdlcomp_test_double ("ima.wav", SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "wav_msadpcm") == 0) cannam@85: { lcomp_test_short ("msadpcm.wav", SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: lcomp_test_int ("msadpcm.wav", SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: lcomp_test_float ("msadpcm.wav", SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: lcomp_test_double ("msadpcm.wav", SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: cannam@85: lcomp_test_short ("msadpcm.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: lcomp_test_int ("msadpcm.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: lcomp_test_float ("msadpcm.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: lcomp_test_double ("msadpcm.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: cannam@85: sdlcomp_test_short ("msadpcm.wav", SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: sdlcomp_test_int ("msadpcm.wav", SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: sdlcomp_test_float ("msadpcm.wav", SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: sdlcomp_test_double ("msadpcm.wav", SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "wav_g721") == 0) cannam@85: { printf ("**** Fix this later : error bound should be 0.06 ****\n") ; cannam@85: lcomp_test_short ("g721.wav", SF_FORMAT_WAV | SF_FORMAT_G721_32, 1, 0.7) ; cannam@85: lcomp_test_int ("g721.wav", SF_FORMAT_WAV | SF_FORMAT_G721_32, 1, 0.7) ; cannam@85: cannam@85: lcomp_test_short ("g721.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_G721_32, 1, 0.7) ; cannam@85: lcomp_test_int ("g721.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_G721_32, 1, 0.7) ; cannam@85: cannam@85: test_count++ ; cannam@85: } ; cannam@85: /* Lite remove end */ cannam@85: cannam@85: if (do_all || strcmp (argv [1], "wav_ulaw") == 0) cannam@85: { lcomp_test_short ("ulaw.wav", SF_FORMAT_WAV | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_int ("ulaw.wav", SF_FORMAT_WAV | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: cannam@85: lcomp_test_short ("ulaw.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_int ("ulaw.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("ulaw.wav", SF_FORMAT_WAV | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_double ("ulaw.wav", SF_FORMAT_WAV | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: /* Lite remove end */ cannam@85: cannam@85: read_raw_test ("ulaw.wav", SF_FORMAT_WAV | SF_FORMAT_ULAW, 2) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "wav_alaw") == 0) cannam@85: { lcomp_test_short ("alaw.wav", SF_FORMAT_WAV | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_int ("alaw.wav", SF_FORMAT_WAV | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("alaw.wav", SF_FORMAT_WAV | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_double ("alaw.wav", SF_FORMAT_WAV | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: /* Lite remove end */ cannam@85: cannam@85: read_raw_test ("alaw.wav", SF_FORMAT_WAV | SF_FORMAT_ALAW, 2) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "wav_gsm610") == 0) cannam@85: { /* Don't do lcomp_test_XXX as the errors are too big. */ cannam@85: sdlcomp_test_short ("gsm610.wav", SF_FORMAT_WAV | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: sdlcomp_test_int ("gsm610.wav", SF_FORMAT_WAV | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: cannam@85: sdlcomp_test_short ("gsm610.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: sdlcomp_test_int ("gsm610.rifx", SF_ENDIAN_BIG | SF_FORMAT_WAV | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: cannam@85: /* Lite remove start */ cannam@85: sdlcomp_test_float ("gsm610.wav", SF_FORMAT_WAV | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: sdlcomp_test_double ("gsm610.wav", SF_FORMAT_WAV | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: /* Lite remove end */ cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "aiff_ulaw") == 0) cannam@85: { lcomp_test_short ("ulaw.aiff", SF_FORMAT_AIFF | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_int ("ulaw.aiff", SF_FORMAT_AIFF | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("ulaw.aiff", SF_FORMAT_AIFF | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_double ("ulaw.aiff", SF_FORMAT_AIFF | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: /* Lite remove end */ cannam@85: cannam@85: read_raw_test ("ulaw.aiff", SF_FORMAT_AIFF | SF_FORMAT_ULAW, 2) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "aiff_alaw") == 0) cannam@85: { lcomp_test_short ("alaw.aiff", SF_FORMAT_AIFF | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_int ("alaw.aiff", SF_FORMAT_AIFF | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("alaw.aiff", SF_FORMAT_AIFF | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_double ("alaw.aiff", SF_FORMAT_AIFF | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: /* Lite remove end */ cannam@85: cannam@85: read_raw_test ("alaw.aiff", SF_FORMAT_AIFF | SF_FORMAT_ALAW, 2) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "aiff_gsm610") == 0) cannam@85: { /* Don't do lcomp_test_XXX as the errors are too big. */ cannam@85: sdlcomp_test_short ("gsm610.aiff", SF_FORMAT_AIFF | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: sdlcomp_test_int ("gsm610.aiff", SF_FORMAT_AIFF | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: /* Lite remove start */ cannam@85: sdlcomp_test_float ("gsm610.aiff", SF_FORMAT_AIFF | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: sdlcomp_test_double ("gsm610.aiff", SF_FORMAT_AIFF | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: /* Lite remove end */ cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (strcmp (argv [1], "aiff_ima") == 0) cannam@85: { lcomp_test_short ("ima.aiff", SF_FORMAT_AIFF | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: lcomp_test_int ("ima.aiff", SF_FORMAT_AIFF | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("ima.aiff", SF_FORMAT_AIFF | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: lcomp_test_double ("ima.aiff", SF_FORMAT_AIFF | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: /* Lite remove end */ cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "au_ulaw") == 0) cannam@85: { lcomp_test_short ("ulaw.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_int ("ulaw.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("ulaw.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_double ("ulaw.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: /* Lite remove end */ cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "au_alaw") == 0) cannam@85: { lcomp_test_short ("alaw.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_int ("alaw.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("alaw.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_double ("alaw.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: /* Lite remove end */ cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: /* Lite remove start */ cannam@85: if (do_all || strcmp (argv [1], "au_g721") == 0) cannam@85: { printf ("**** Fix this later : error bound should be 0.06 ****\n") ; cannam@85: lcomp_test_short ("g721.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_G721_32, 1, 0.7) ; cannam@85: lcomp_test_int ("g721.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_G721_32, 1, 0.7) ; cannam@85: lcomp_test_float ("g721.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_G721_32, 1, 0.7) ; cannam@85: lcomp_test_double ("g721.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_G721_32, 1, 0.7) ; cannam@85: cannam@85: /*- sdlcomp_test_short ("g721.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_G721_32, 1, 0.07) ; cannam@85: sdlcomp_test_int ("g721.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_G721_32, 1, 0.07) ; cannam@85: sdlcomp_test_float ("g721.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_G721_32, 1, 0.07) ; cannam@85: sdlcomp_test_double ("g721.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_G721_32, 1, 0.12) ; cannam@85: -*/ cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "au_g723") == 0) cannam@85: { printf ("**** Fix this later : error bound should be 0.16 ****\n") ; cannam@85: lcomp_test_short ("g723_24.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_G723_24, 1, 0.7) ; cannam@85: lcomp_test_int ("g723_24.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_G723_24, 1, 0.7) ; cannam@85: lcomp_test_float ("g723_24.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_G723_24, 1, 0.7) ; cannam@85: lcomp_test_double ("g723_24.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_G723_24, 1, 0.7) ; cannam@85: cannam@85: lcomp_test_short ("g723_40.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_G723_40, 1, 0.85) ; cannam@85: lcomp_test_int ("g723_40.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_G723_40, 1, 0.84) ; cannam@85: lcomp_test_float ("g723_40.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_G723_40, 1, 0.86) ; cannam@85: lcomp_test_double ("g723_40.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_G723_40, 1, 0.86) ; cannam@85: cannam@85: /*- sdlcomp_test_short ("g723.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_G723_24, 1, 0.15) ; cannam@85: sdlcomp_test_int ("g723.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_G723_24, 1, 0.15) ; cannam@85: sdlcomp_test_float ("g723.au", SF_ENDIAN_BIG | SF_FORMAT_AU | SF_FORMAT_G723_24, 1, 0.15) ; cannam@85: sdlcomp_test_double ("g723.au", SF_ENDIAN_LITTLE | SF_FORMAT_AU | SF_FORMAT_G723_24, 1, 0.15) ; cannam@85: -*/ cannam@85: test_count++ ; cannam@85: } ; cannam@85: /* Lite remove end */ cannam@85: cannam@85: if (do_all || strcmp (argv [1], "caf_ulaw") == 0) cannam@85: { lcomp_test_short ("ulaw.caf", SF_FORMAT_CAF | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_int ("ulaw.caf", SF_FORMAT_CAF | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("ulaw.caf", SF_FORMAT_CAF | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_double ("ulaw.caf", SF_FORMAT_CAF | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: /* Lite remove end */ cannam@85: cannam@85: read_raw_test ("ulaw.caf", SF_FORMAT_CAF | SF_FORMAT_ULAW, 2) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "caf_alaw") == 0) cannam@85: { lcomp_test_short ("alaw.caf", SF_FORMAT_CAF | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_int ("alaw.caf", SF_FORMAT_CAF | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("alaw.caf", SF_FORMAT_CAF | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_double ("alaw.caf", SF_FORMAT_CAF | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: /* Lite remove end */ cannam@85: cannam@85: read_raw_test ("alaw.caf", SF_FORMAT_CAF | SF_FORMAT_ALAW, 2) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: cannam@85: if (do_all || strcmp (argv [1], "raw_ulaw") == 0) cannam@85: { lcomp_test_short ("ulaw.raw", SF_ENDIAN_LITTLE | SF_FORMAT_RAW | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_int ("ulaw.raw", SF_ENDIAN_BIG | SF_FORMAT_RAW | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("ulaw.raw", SF_ENDIAN_LITTLE | SF_FORMAT_RAW | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_double ("ulaw.raw", SF_ENDIAN_BIG | SF_FORMAT_RAW | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: /* Lite remove end */ cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "raw_alaw") == 0) cannam@85: { lcomp_test_short ("alaw.raw", SF_ENDIAN_LITTLE | SF_FORMAT_RAW | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_int ("alaw.raw", SF_ENDIAN_BIG | SF_FORMAT_RAW | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("alaw.raw", SF_ENDIAN_LITTLE | SF_FORMAT_RAW | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_double ("alaw.raw", SF_ENDIAN_BIG | SF_FORMAT_RAW | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: /* Lite remove end */ cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "raw_gsm610") == 0) cannam@85: { /* Don't do lcomp_test_XXX as the errors are too big. */ cannam@85: sdlcomp_test_short ("raw.gsm", SF_FORMAT_RAW | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: sdlcomp_test_int ("raw.gsm", SF_FORMAT_RAW | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: sdlcomp_test_float ("raw.gsm", SF_FORMAT_RAW | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: sdlcomp_test_double ("raw.gsm", SF_FORMAT_RAW | SF_FORMAT_GSM610, 1, 0.24) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "ogg_vorbis") == 0) cannam@85: { if (HAVE_EXTERNAL_LIBS) cannam@85: { /* Don't do lcomp_test_XXX as the errors are too big. */ cannam@85: sdlcomp_test_short ("vorbis.oga", SF_FORMAT_OGG | SF_FORMAT_VORBIS, 1, 0.30) ; cannam@85: sdlcomp_test_int ("vorbis.oga", SF_FORMAT_OGG | SF_FORMAT_VORBIS, 1, 0.30) ; cannam@85: sdlcomp_test_float ("vorbis.oga", SF_FORMAT_OGG | SF_FORMAT_VORBIS, 1, 0.30) ; cannam@85: sdlcomp_test_double ("vorbis.oga", SF_FORMAT_OGG | SF_FORMAT_VORBIS, 1, 0.30) ; cannam@85: } cannam@85: else cannam@85: puts (" No Ogg/Vorbis tests because Ogg/Vorbis support was not compiled in.") ; cannam@85: cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: /* Lite remove start */ cannam@85: if (do_all || strcmp (argv [1], "ircam_ulaw") == 0) cannam@85: { lcomp_test_short ("ulaw.ircam", SF_ENDIAN_LITTLE | SF_FORMAT_IRCAM | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_int ("ulaw.ircam", SF_ENDIAN_BIG | SF_FORMAT_IRCAM | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_float ("ulaw.ircam", SF_ENDIAN_LITTLE | SF_FORMAT_IRCAM | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_double ("ulaw.ircam", SF_ENDIAN_BIG | SF_FORMAT_IRCAM | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "ircam_alaw") == 0) cannam@85: { lcomp_test_short ("alaw.ircam", SF_ENDIAN_LITTLE | SF_FORMAT_IRCAM | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_int ("alaw.ircam", SF_ENDIAN_BIG | SF_FORMAT_IRCAM | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_float ("alaw.ircam", SF_ENDIAN_LITTLE | SF_FORMAT_IRCAM | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_double ("alaw.ircam", SF_ENDIAN_BIG | SF_FORMAT_IRCAM | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "nist_ulaw") == 0) cannam@85: { lcomp_test_short ("ulaw.nist", SF_ENDIAN_LITTLE | SF_FORMAT_NIST | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_int ("ulaw.nist", SF_ENDIAN_BIG | SF_FORMAT_NIST | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_float ("ulaw.nist", SF_ENDIAN_LITTLE | SF_FORMAT_NIST | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_double ("ulaw.nist", SF_ENDIAN_BIG | SF_FORMAT_NIST | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "nist_alaw") == 0) cannam@85: { lcomp_test_short ("alaw.nist", SF_ENDIAN_LITTLE | SF_FORMAT_NIST | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_int ("alaw.nist", SF_ENDIAN_BIG | SF_FORMAT_NIST | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_float ("alaw.nist", SF_ENDIAN_LITTLE | SF_FORMAT_NIST | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_double ("alaw.nist", SF_ENDIAN_BIG | SF_FORMAT_NIST | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "voc_ulaw") == 0) cannam@85: { lcomp_test_short ("ulaw.voc", SF_FORMAT_VOC | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_int ("ulaw.voc", SF_FORMAT_VOC | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_float ("ulaw.voc", SF_FORMAT_VOC | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_double ("ulaw.voc", SF_FORMAT_VOC | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "voc_alaw") == 0) cannam@85: { lcomp_test_short ("alaw.voc", SF_FORMAT_VOC | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_int ("alaw.voc", SF_FORMAT_VOC | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_float ("alaw.voc", SF_FORMAT_VOC | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_double ("alaw.voc", SF_FORMAT_VOC | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: /* Lite remove end */ cannam@85: cannam@85: if (do_all || strcmp (argv [1], "w64_ulaw") == 0) cannam@85: { lcomp_test_short ("ulaw.w64", SF_FORMAT_W64 | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_int ("ulaw.w64", SF_FORMAT_W64 | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("ulaw.w64", SF_FORMAT_W64 | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: lcomp_test_double ("ulaw.w64", SF_FORMAT_W64 | SF_FORMAT_ULAW, 2, 0.04) ; cannam@85: /* Lite remove end */ cannam@85: cannam@85: read_raw_test ("ulaw.w64", SF_FORMAT_W64 | SF_FORMAT_ULAW, 2) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "w64_alaw") == 0) cannam@85: { lcomp_test_short ("alaw.w64", SF_FORMAT_W64 | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_int ("alaw.w64", SF_FORMAT_W64 | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("alaw.w64", SF_FORMAT_W64 | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: lcomp_test_double ("alaw.w64", SF_FORMAT_W64 | SF_FORMAT_ALAW, 2, 0.04) ; cannam@85: /* Lite remove end */ cannam@85: cannam@85: read_raw_test ("alaw.w64", SF_FORMAT_W64 | SF_FORMAT_ALAW, 2) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: /* Lite remove start */ cannam@85: if (do_all || strcmp (argv [1], "w64_ima") == 0) cannam@85: { lcomp_test_short ("ima.w64", SF_FORMAT_W64 | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: lcomp_test_int ("ima.w64", SF_FORMAT_W64 | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: lcomp_test_float ("ima.w64", SF_FORMAT_W64 | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: lcomp_test_double ("ima.w64", SF_FORMAT_W64 | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: cannam@85: sdlcomp_test_short ("ima.w64", SF_FORMAT_W64 | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: sdlcomp_test_int ("ima.w64", SF_FORMAT_W64 | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: sdlcomp_test_float ("ima.w64", SF_FORMAT_W64 | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: sdlcomp_test_double ("ima.w64", SF_FORMAT_W64 | SF_FORMAT_IMA_ADPCM, 2, 0.18) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "w64_msadpcm") == 0) cannam@85: { lcomp_test_short ("msadpcm.w64", SF_FORMAT_W64 | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: lcomp_test_int ("msadpcm.w64", SF_FORMAT_W64 | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: lcomp_test_float ("msadpcm.w64", SF_FORMAT_W64 | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: lcomp_test_double ("msadpcm.w64", SF_FORMAT_W64 | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: cannam@85: sdlcomp_test_short ("msadpcm.w64", SF_FORMAT_W64 | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: sdlcomp_test_int ("msadpcm.w64", SF_FORMAT_W64 | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: sdlcomp_test_float ("msadpcm.w64", SF_FORMAT_W64 | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: sdlcomp_test_double ("msadpcm.w64", SF_FORMAT_W64 | SF_FORMAT_MS_ADPCM, 2, 0.36) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "wve") == 0) cannam@85: { lcomp_test_short ("psion.wve", SF_FORMAT_WVE | SF_FORMAT_ALAW, 1, 0.04) ; cannam@85: lcomp_test_int ("psion.wve", SF_FORMAT_WVE | SF_FORMAT_ALAW, 1, 0.04) ; cannam@85: /* Lite remove start */ cannam@85: lcomp_test_float ("psion.wve", SF_FORMAT_WVE | SF_FORMAT_ALAW, 1, 0.04) ; cannam@85: lcomp_test_double ("psion.wve", SF_FORMAT_WVE | SF_FORMAT_ALAW, 1, 0.04) ; cannam@85: /* Lite remove end */ cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: /* Lite remove end */ cannam@85: cannam@85: if (do_all || strcmp (argv [1], "w64_gsm610") == 0) cannam@85: { /* Don't do lcomp_test_XXX as the errors are too big. */ cannam@85: sdlcomp_test_short ("gsm610.w64", SF_FORMAT_W64 | SF_FORMAT_GSM610, 1, 0.2) ; cannam@85: sdlcomp_test_int ("gsm610.w64", SF_FORMAT_W64 | SF_FORMAT_GSM610, 1, 0.2) ; cannam@85: /* Lite remove start */ cannam@85: sdlcomp_test_float ("gsm610.w64", SF_FORMAT_W64 | SF_FORMAT_GSM610, 1, 0.2) ; cannam@85: sdlcomp_test_double ("gsm610.w64", SF_FORMAT_W64 | SF_FORMAT_GSM610, 1, 0.2) ; cannam@85: /* Lite remove end */ cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: /* Lite remove start */ cannam@85: if (do_all || strcmp (argv [1], "vox_adpcm") == 0) cannam@85: { lcomp_test_short ("adpcm.vox", SF_FORMAT_RAW | SF_FORMAT_VOX_ADPCM, 1, 0.17) ; cannam@85: lcomp_test_int ("adpcm.vox", SF_FORMAT_RAW | SF_FORMAT_VOX_ADPCM, 1, 0.17) ; cannam@85: lcomp_test_float ("adpcm.vox", SF_FORMAT_RAW | SF_FORMAT_VOX_ADPCM, 1, 0.17) ; cannam@85: lcomp_test_double ("adpcm.vox", SF_FORMAT_RAW | SF_FORMAT_VOX_ADPCM, 1, 0.17) ; cannam@85: cannam@85: sdlcomp_test_short ("adpcm.vox", SF_FORMAT_RAW | SF_FORMAT_VOX_ADPCM, 1, 0.072) ; cannam@85: sdlcomp_test_int ("adpcm.vox", SF_FORMAT_RAW | SF_FORMAT_VOX_ADPCM, 1, 0.072) ; cannam@85: sdlcomp_test_float ("adpcm.vox", SF_FORMAT_RAW | SF_FORMAT_VOX_ADPCM, 1, 0.072) ; cannam@85: sdlcomp_test_double ("adpcm.vox", SF_FORMAT_RAW | SF_FORMAT_VOX_ADPCM, 1, 0.072) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: cannam@85: if (do_all || strcmp (argv [1], "xi_dpcm") == 0) cannam@85: { lcomp_test_short ("8bit.xi", SF_FORMAT_XI | SF_FORMAT_DPCM_8, 1, 0.25) ; cannam@85: lcomp_test_int ("8bit.xi", SF_FORMAT_XI | SF_FORMAT_DPCM_8, 1, 0.25) ; cannam@85: cannam@85: lcomp_test_short ("16bit.xi", SF_FORMAT_XI | SF_FORMAT_DPCM_16, 1, 0.002) ; cannam@85: lcomp_test_int ("16bit.xi", SF_FORMAT_XI | SF_FORMAT_DPCM_16, 1, 0.002) ; cannam@85: lcomp_test_float ("16bit.xi", SF_FORMAT_XI | SF_FORMAT_DPCM_16, 1, 0.002) ; cannam@85: lcomp_test_double ("16bit.xi", SF_FORMAT_XI | SF_FORMAT_DPCM_16, 1, 0.002) ; cannam@85: test_count++ ; cannam@85: } ; cannam@85: /* Lite remove end */ cannam@85: cannam@85: if (test_count == 0) cannam@85: { printf ("************************************\n") ; cannam@85: printf ("* No '%s' test defined.\n", argv [1]) ; cannam@85: printf ("************************************\n") ; cannam@85: return 1 ; cannam@85: } ; cannam@85: cannam@85: return 0 ; cannam@85: } /* main */ cannam@85: cannam@85: /*============================================================================================ cannam@85: ** Here are the test functions. cannam@85: */ cannam@85: cannam@85: static void cannam@85: lcomp_test_short (const char *filename, int filetype, int channels, double margin) cannam@85: { SNDFILE *file ; cannam@85: SF_INFO sfinfo ; cannam@85: int k, m, seekpos, half_max_abs ; cannam@85: long datalen ; cannam@85: short *orig, *data ; cannam@85: cannam@85: print_test_name ("lcomp_test_short", filename) ; cannam@85: cannam@85: datalen = BUFFER_SIZE / channels ; cannam@85: cannam@85: data = data_buffer.s ; cannam@85: orig = orig_buffer.s ; cannam@85: cannam@85: gen_signal_double (orig_buffer.d, 32000.0, channels, datalen) ; cannam@85: for (k = 0 ; k < channels * datalen ; k++) cannam@85: orig [k] = (short) (orig_buffer.d [k]) ; cannam@85: cannam@85: sfinfo.samplerate = SAMPLE_RATE ; cannam@85: sfinfo.frames = 123456789 ; /* Ridiculous value. */ cannam@85: sfinfo.channels = channels ; cannam@85: sfinfo.format = filetype ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_WRITE, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: test_writef_short_or_die (file, 0, orig, datalen, __LINE__) ; cannam@85: sf_set_string (file, SF_STR_COMMENT, long_comment) ; cannam@85: sf_close (file) ; cannam@85: cannam@85: memset (data, 0, datalen * sizeof (short)) ; cannam@85: cannam@85: if ((filetype & SF_FORMAT_TYPEMASK) != SF_FORMAT_RAW) cannam@85: memset (&sfinfo, 0, sizeof (sfinfo)) ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_READ, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: cannam@85: if ((sfinfo.format & (SF_FORMAT_TYPEMASK | SF_FORMAT_SUBMASK)) != (filetype & (SF_FORMAT_TYPEMASK | SF_FORMAT_SUBMASK))) cannam@85: { printf ("\n\nLine %d: Returned format incorrect (0x%08X => 0x%08X).\n", __LINE__, filetype, sfinfo.format) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames < datalen / channels) cannam@85: { printf ("Too few frames in file. (%ld should be a little more than %ld)\n", SF_COUNT_TO_LONG (sfinfo.frames), datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames > (datalen + datalen / 20)) cannam@85: { printf ("Too many frames in file. (%ld should be a little more than %ld)\n", SF_COUNT_TO_LONG (sfinfo.frames), datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.channels != channels) cannam@85: { printf ("Incorrect number of channels in file.\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: check_log_buffer_or_die (file, __LINE__) ; cannam@85: cannam@85: check_comment (file, filetype, __LINE__) ; cannam@85: cannam@85: test_readf_short_or_die (file, 0, data, datalen, __LINE__) ; cannam@85: cannam@85: half_max_abs = 0 ; cannam@85: for (k = 0 ; k < datalen ; k++) cannam@85: { if (error_function (data [k], orig [k], margin)) cannam@85: { printf ("\n\nLine %d: Incorrect sample A (#%d : %d should be %d).\n", __LINE__, k, data [k], orig [k]) ; cannam@85: oct_save_short (orig, data, datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: half_max_abs = LCT_MAX (half_max_abs, abs (data [k] / 2)) ; cannam@85: } ; cannam@85: cannam@85: if (half_max_abs < 1.0) cannam@85: { printf ("\n\nLine %d: Signal is all zeros.\n", __LINE__) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_readf_short (file, data, datalen)) != sfinfo.frames - datalen) cannam@85: { printf ("\n\nLine %d: Incorrect read length (%ld should be %d).\n", __LINE__, cannam@85: SF_COUNT_TO_LONG (channels * sfinfo.frames - datalen), k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* This check is only for block based encoders which must append silence cannam@85: ** to the end of a file so as to fill out a block. cannam@85: */ cannam@85: for (k = 0 ; k < sfinfo.frames - datalen ; k++) cannam@85: if (abs (data [channels * k]) > decay_response (channels * k)) cannam@85: { printf ("\n\nLine %d : Incorrect sample B (#%d : abs (%d) should be < %d).\n", __LINE__, channels * k, data [channels * k], decay_response (channels * k)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (! sfinfo.seekable) cannam@85: { sf_close (file) ; cannam@85: unlink (filename) ; cannam@85: printf ("ok\n") ; cannam@85: return ; cannam@85: } ; cannam@85: cannam@85: /* Now test sf_seek function. */ cannam@85: cannam@85: if ((k = sf_seek (file, 0, SEEK_SET)) != 0) cannam@85: { printf ("\n\nLine %d: Seek to start of file failed (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: for (m = 0 ; m < 3 ; m++) cannam@85: { test_readf_short_or_die (file, m, data, 11, __LINE__) ; cannam@85: cannam@85: for (k = 0 ; k < channels * 11 ; k++) cannam@85: if (error_function (1.0 * data [k], 1.0 * orig [k + channels * m * 11], margin)) cannam@85: { printf ("\n\nLine %d: Incorrect sample (m = %d) (#%d : %d => %d).\n", __LINE__, m, k + channels * m * 11, orig [k + channels * m * 11], data [k]) ; cannam@85: for (m = 0 ; m < channels ; m++) cannam@85: printf ("%d ", data [m]) ; cannam@85: printf ("\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: } ; cannam@85: cannam@85: seekpos = BUFFER_SIZE / 10 ; cannam@85: cannam@85: /* Check seek from start of file. */ cannam@85: if ((k = sf_seek (file, seekpos, SEEK_SET)) != seekpos) cannam@85: { printf ("Seek to start of file + %d failed (%d).\n", seekpos, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: test_readf_short_or_die (file, 0, data, 1, __LINE__) ; cannam@85: cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [seekpos * channels], margin)) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_SET) followed by sf_readf_short failed (%d, %d).\n", __LINE__, orig [1], data [0]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_seek (file, 0, SEEK_CUR)) != seekpos + 1) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_CUR) with 0 offset failed (%d should be %d)\n", __LINE__, k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) + BUFFER_SIZE / 5 ; cannam@85: k = sf_seek (file, BUFFER_SIZE / 5, SEEK_CUR) ; cannam@85: test_readf_short_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\n\nLine %d: sf_seek (forwards, SEEK_CUR) followed by sf_readf_short failed (%d, %d) (%d, %d).\n", __LINE__, data [0], orig [seekpos * channels], k, seekpos + 1) ; cannam@85: oct_save_short (orig, data, datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) - 20 ; cannam@85: /* Check seek backward from current position. */ cannam@85: k = sf_seek (file, -20, SEEK_CUR) ; cannam@85: test_readf_short_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (backwards, SEEK_CUR) followed by sf_readf_short failed (%d, %d) (%d, %d).\n", __LINE__, data [0], orig [seekpos * channels], k, seekpos) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check that read past end of file returns number of items. */ cannam@85: sf_seek (file, sfinfo.frames, SEEK_SET) ; cannam@85: cannam@85: if ((k = sf_readf_short (file, data, datalen)) != 0) cannam@85: { printf ("\n\nLine %d: Return value from sf_readf_short past end of file incorrect (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check seek backward from end. */ cannam@85: if ((k = sf_seek (file, 5 - sfinfo.frames, SEEK_END)) != 5) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_END) returned %d instead of %d.\n", __LINE__, k, 5) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: test_readf_short_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [5 * channels], margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_END) followed by sf_readf_short failed (%d should be %d).\n", __LINE__, data [0], orig [5 * channels]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: sf_close (file) ; cannam@85: cannam@85: unlink (filename) ; cannam@85: printf ("ok\n") ; cannam@85: } /* lcomp_test_short */ cannam@85: cannam@85: /*-------------------------------------------------------------------------------------------- cannam@85: */ cannam@85: cannam@85: static void cannam@85: lcomp_test_int (const char *filename, int filetype, int channels, double margin) cannam@85: { SNDFILE *file ; cannam@85: SF_INFO sfinfo ; cannam@85: int k, m, half_max_abs ; cannam@85: long datalen, seekpos ; cannam@85: double scale, max_val ; cannam@85: int *orig, *data ; cannam@85: cannam@85: print_test_name ("lcomp_test_int", filename) ; cannam@85: cannam@85: datalen = BUFFER_SIZE / channels ; cannam@85: cannam@85: if (is_lossy (filetype)) cannam@85: { scale = 1.0 * 0x10000 ; cannam@85: max_val = 32000.0 * scale ; cannam@85: } cannam@85: else cannam@85: { scale = 1.0 ; cannam@85: max_val = 0x7fffffff * scale ; cannam@85: } ; cannam@85: cannam@85: data = data_buffer.i ; cannam@85: orig = orig_buffer.i ; cannam@85: cannam@85: gen_signal_double (orig_buffer.d, max_val, channels, datalen) ; cannam@85: cannam@85: for (k = 0 ; k < channels * datalen ; k++) cannam@85: orig [k] = lrint (orig_buffer.d [k]) ; cannam@85: cannam@85: sfinfo.samplerate = SAMPLE_RATE ; cannam@85: sfinfo.frames = 123456789 ; /* Ridiculous value. */ cannam@85: sfinfo.channels = channels ; cannam@85: sfinfo.format = filetype ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_WRITE, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: test_writef_int_or_die (file, 0, orig, datalen, __LINE__) ; cannam@85: sf_set_string (file, SF_STR_COMMENT, long_comment) ; cannam@85: sf_close (file) ; cannam@85: cannam@85: memset (data, 0, datalen * sizeof (int)) ; cannam@85: cannam@85: if ((filetype & SF_FORMAT_TYPEMASK) != SF_FORMAT_RAW) cannam@85: memset (&sfinfo, 0, sizeof (sfinfo)) ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_READ, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: cannam@85: if ((sfinfo.format & (SF_FORMAT_TYPEMASK | SF_FORMAT_SUBMASK)) != (filetype & (SF_FORMAT_TYPEMASK | SF_FORMAT_SUBMASK))) cannam@85: { printf ("\n\nLine %d: Returned format incorrect (0x%08X => 0x%08X).\n", __LINE__, filetype, sfinfo.format) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames < datalen / channels) cannam@85: { printf ("Too few.frames in file. (%ld should be a little more than %ld)\n", datalen, SF_COUNT_TO_LONG (sfinfo.frames)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames > (datalen + datalen / 20)) cannam@85: { printf ("Too many.frames in file. (%ld should be a little more than %ld)\n", datalen, SF_COUNT_TO_LONG (sfinfo.frames)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.channels != channels) cannam@85: { printf ("Incorrect number of channels in file.\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: check_log_buffer_or_die (file, __LINE__) ; cannam@85: cannam@85: check_comment (file, filetype, __LINE__) ; cannam@85: cannam@85: test_readf_int_or_die (file, 0, data, datalen, __LINE__) ; cannam@85: cannam@85: half_max_abs = 0 ; cannam@85: for (k = 0 ; k < datalen ; k++) cannam@85: { if (error_function (data [k] / scale, orig [k] / scale, margin)) cannam@85: { printf ("\nLine %d: Incorrect sample (#%d : %f should be %f).\n", __LINE__, k, data [k] / scale, orig [k] / scale) ; cannam@85: oct_save_int (orig, data, datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: half_max_abs = LCT_MAX (half_max_abs, abs (data [k] / 2)) ; cannam@85: } ; cannam@85: cannam@85: if (half_max_abs < 1.0) cannam@85: { printf ("\n\nLine %d: Signal is all zeros (%d, 0x%X).\n", __LINE__, half_max_abs, half_max_abs) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_readf_int (file, data, datalen)) != sfinfo.frames - datalen) cannam@85: { printf ("\n\nLine %d: Incorrect read length (%ld should be %d).\n", __LINE__, cannam@85: SF_COUNT_TO_LONG (channels * sfinfo.frames - datalen), k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* This check is only for block based encoders which must append silence cannam@85: ** to the end of a file so as to fill out a block. cannam@85: */ cannam@85: if ((sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_MS_ADPCM) cannam@85: for (k = 0 ; k < sfinfo.frames - datalen ; k++) cannam@85: if (abs (data [channels * k] / scale) > decay_response (channels * k)) cannam@85: { printf ("\n\nLine %d : Incorrect sample B (#%d : abs (%d) should be < %d).\n", __LINE__, channels * k, data [channels * k], decay_response (channels * k)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (! sfinfo.seekable) cannam@85: { sf_close (file) ; cannam@85: unlink (filename) ; cannam@85: printf ("ok\n") ; cannam@85: return ; cannam@85: } ; cannam@85: cannam@85: /* Now test sf_seek function. */ cannam@85: cannam@85: if ((k = sf_seek (file, 0, SEEK_SET)) != 0) cannam@85: { printf ("\n\nLine %d: Seek to start of file failed (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: for (m = 0 ; m < 3 ; m++) cannam@85: { test_readf_int_or_die (file, m, data, 11, __LINE__) ; cannam@85: cannam@85: for (k = 0 ; k < channels * 11 ; k++) cannam@85: if (error_function (data [k] / scale, orig [k + channels * m * 11] / scale, margin)) cannam@85: { printf ("\nLine %d: Incorrect sample (m = %d) (#%d : %d => %d).\n", __LINE__, m, k + channels * m * 11, orig [k + channels * m * 11], data [k]) ; cannam@85: for (m = 0 ; m < channels ; m++) cannam@85: printf ("%d ", data [m]) ; cannam@85: printf ("\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: } ; cannam@85: cannam@85: seekpos = BUFFER_SIZE / 10 ; cannam@85: cannam@85: /* Check seek from start of file. */ cannam@85: if ((k = sf_seek (file, seekpos, SEEK_SET)) != seekpos) cannam@85: { printf ("Seek to start of file + %ld failed (%d).\n", seekpos, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: test_readf_int_or_die (file, 0, data, 1, __LINE__) ; cannam@85: cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [seekpos * channels], margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_SET) followed by sf_readf_int failed (%d, %d).\n", __LINE__, orig [1], data [0]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_seek (file, 0, SEEK_CUR)) != seekpos + 1) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_CUR) with 0 offset failed (%d should be %ld)\n", __LINE__, k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) + BUFFER_SIZE / 5 ; cannam@85: k = sf_seek (file, BUFFER_SIZE / 5, SEEK_CUR) ; cannam@85: test_readf_int_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (forwards, SEEK_CUR) followed by sf_readf_int failed (%d, %d) (%d, %ld).\n", __LINE__, data [0], orig [seekpos * channels], k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) - 20 ; cannam@85: /* Check seek backward from current position. */ cannam@85: k = sf_seek (file, -20, SEEK_CUR) ; cannam@85: test_readf_int_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (backwards, SEEK_CUR) followed by sf_readf_int failed (%d, %d) (%d, %ld).\n", __LINE__, data [0], orig [seekpos * channels], k, seekpos) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check that read past end of file returns number of items. */ cannam@85: sf_seek (file, sfinfo.frames, SEEK_SET) ; cannam@85: cannam@85: if ((k = sf_readf_int (file, data, datalen)) != 0) cannam@85: { printf ("\n\nLine %d: Return value from sf_readf_int past end of file incorrect (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check seek backward from end. */ cannam@85: if ((k = sf_seek (file, 5 - sfinfo.frames, SEEK_END)) != 5) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_END) returned %d instead of %d.\n", __LINE__, k, 5) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: test_readf_int_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (data [0] / scale, orig [5 * channels] / scale, margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_END) followed by sf_readf_short failed (%d should be %d).\n", __LINE__, data [0], orig [5]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: sf_close (file) ; cannam@85: cannam@85: unlink (filename) ; cannam@85: printf ("ok\n") ; cannam@85: } /* lcomp_test_int */ cannam@85: cannam@85: /*-------------------------------------------------------------------------------------------- cannam@85: */ cannam@85: cannam@85: static void cannam@85: lcomp_test_float (const char *filename, int filetype, int channels, double margin) cannam@85: { SNDFILE *file ; cannam@85: SF_INFO sfinfo ; cannam@85: int k, m, seekpos ; cannam@85: long datalen ; cannam@85: float *orig, *data ; cannam@85: double half_max_abs ; cannam@85: cannam@85: print_test_name ("lcomp_test_float", filename) ; cannam@85: cannam@85: datalen = BUFFER_SIZE / channels ; cannam@85: cannam@85: data = data_buffer.f ; cannam@85: orig = orig_buffer.f ; cannam@85: cannam@85: gen_signal_double (orig_buffer.d, 32000.0, channels, datalen) ; cannam@85: for (k = 0 ; k < channels * datalen ; k++) cannam@85: orig [k] = orig_buffer.d [k] ; cannam@85: cannam@85: sfinfo.samplerate = SAMPLE_RATE ; cannam@85: sfinfo.frames = 123456789 ; /* Ridiculous value. */ cannam@85: sfinfo.channels = channels ; cannam@85: sfinfo.format = filetype ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_WRITE, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: sf_command (file, SFC_SET_NORM_FLOAT, NULL, SF_FALSE) ; cannam@85: test_writef_float_or_die (file, 0, orig, datalen, __LINE__) ; cannam@85: sf_set_string (file, SF_STR_COMMENT, long_comment) ; cannam@85: sf_close (file) ; cannam@85: cannam@85: memset (data, 0, datalen * sizeof (float)) ; cannam@85: cannam@85: if ((filetype & SF_FORMAT_TYPEMASK) != SF_FORMAT_RAW) cannam@85: memset (&sfinfo, 0, sizeof (sfinfo)) ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_READ, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: cannam@85: if ((sfinfo.format & (SF_FORMAT_TYPEMASK | SF_FORMAT_SUBMASK)) != (filetype & (SF_FORMAT_TYPEMASK | SF_FORMAT_SUBMASK))) cannam@85: { printf ("\n\nLine %d: Returned format incorrect (0x%08X => 0x%08X).\n", __LINE__, filetype, sfinfo.format) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames < datalen / channels) cannam@85: { printf ("Too few.frames in file. (%ld should be a little more than %ld)\n", datalen, SF_COUNT_TO_LONG (sfinfo.frames)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames > (datalen + datalen / 20)) cannam@85: { printf ("Too many.frames in file. (%ld should be a little more than %ld)\n", datalen, SF_COUNT_TO_LONG (sfinfo.frames)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.channels != channels) cannam@85: { printf ("Incorrect number of channels in file.\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: check_comment (file, filetype, __LINE__) ; cannam@85: cannam@85: sf_command (file, SFC_SET_NORM_FLOAT, NULL, SF_FALSE) ; cannam@85: cannam@85: check_log_buffer_or_die (file, __LINE__) ; cannam@85: cannam@85: check_comment (file, filetype, __LINE__) ; cannam@85: cannam@85: sf_command (file, SFC_SET_NORM_FLOAT, NULL, SF_FALSE) ; cannam@85: cannam@85: test_readf_float_or_die (file, 0, data, datalen, __LINE__) ; cannam@85: cannam@85: half_max_abs = 0.0 ; cannam@85: for (k = 0 ; k < datalen ; k++) cannam@85: { if (error_function (data [k], orig [k], margin)) cannam@85: { printf ("\nLine %d: Incorrect sample A (#%d : %f should be %f).\n", __LINE__, k, data [k], orig [k]) ; cannam@85: oct_save_float (orig, data, datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: half_max_abs = LCT_MAX (half_max_abs, fabs (0.5 * data [k])) ; cannam@85: } ; cannam@85: cannam@85: if (half_max_abs < 1.0) cannam@85: { printf ("\n\nLine %d: Signal is all zeros.\n", __LINE__) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_readf_float (file, data, datalen)) != sfinfo.frames - datalen) cannam@85: { printf ("\n\nLine %d: Incorrect read length (%ld should be %d).\n", __LINE__, cannam@85: SF_COUNT_TO_LONG (channels * sfinfo.frames - datalen), k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* This check is only for block based encoders which must append silence cannam@85: ** to the end of a file so as to fill out a block. cannam@85: */ cannam@85: if ((sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_MS_ADPCM) cannam@85: for (k = 0 ; k < sfinfo.frames - datalen ; k++) cannam@85: if (abs (data [channels * k]) > decay_response (channels * k)) cannam@85: { printf ("\n\nLine %d : Incorrect sample B (#%d : abs (%f) should be < %d).\n", __LINE__, channels * k, data [channels * k], decay_response (channels * k)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (! sfinfo.seekable) cannam@85: { sf_close (file) ; cannam@85: unlink (filename) ; cannam@85: printf ("ok\n") ; cannam@85: return ; cannam@85: } ; cannam@85: cannam@85: /* Now test sf_seek function. */ cannam@85: cannam@85: if ((k = sf_seek (file, 0, SEEK_SET)) != 0) cannam@85: { printf ("\n\nLine %d: Seek to start of file failed (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: for (m = 0 ; m < 3 ; m++) cannam@85: { test_readf_float_or_die (file, 0, data, 11, __LINE__) ; cannam@85: cannam@85: for (k = 0 ; k < channels * 11 ; k++) cannam@85: if (error_function (data [k], orig [k + channels * m * 11], margin)) cannam@85: { printf ("\nLine %d: Incorrect sample (m = %d) (#%d : %f => %f).\n", __LINE__, m, k + channels * m * 11, orig [k + channels * m * 11], data [k]) ; cannam@85: for (m = 0 ; m < channels ; m++) cannam@85: printf ("%f ", data [m]) ; cannam@85: printf ("\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: } ; cannam@85: cannam@85: seekpos = BUFFER_SIZE / 10 ; cannam@85: cannam@85: /* Check seek from start of file. */ cannam@85: if ((k = sf_seek (file, seekpos, SEEK_SET)) != seekpos) cannam@85: { printf ("Seek to start of file + %d failed (%d).\n", seekpos, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: test_readf_float_or_die (file, 0, data, 1, __LINE__) ; cannam@85: cannam@85: if (error_function (data [0], orig [seekpos * channels], margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_SET) followed by sf_readf_float failed (%f, %f).\n", __LINE__, orig [1], data [0]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_seek (file, 0, SEEK_CUR)) != seekpos + 1) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_CUR) with 0 offset failed (%d should be %d)\n", __LINE__, k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) + BUFFER_SIZE / 5 ; cannam@85: k = sf_seek (file, BUFFER_SIZE / 5, SEEK_CUR) ; cannam@85: test_readf_float_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (data [0], orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (forwards, SEEK_CUR) followed by sf_readf_float failed (%f, %f) (%d, %d).\n", __LINE__, data [0], orig [seekpos * channels], k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) - 20 ; cannam@85: /* Check seek backward from current position. */ cannam@85: k = sf_seek (file, -20, SEEK_CUR) ; cannam@85: test_readf_float_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (data [0], orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (backwards, SEEK_CUR) followed by sf_readf_float failed (%f, %f) (%d, %d).\n", __LINE__, data [0], orig [seekpos * channels], k, seekpos) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check that read past end of file returns number of items. */ cannam@85: sf_seek (file, sfinfo.frames, SEEK_SET) ; cannam@85: cannam@85: if ((k = sf_readf_float (file, data, datalen)) != 0) cannam@85: { printf ("\n\nLine %d: Return value from sf_readf_float past end of file incorrect (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check seek backward from end. */ cannam@85: if ((k = sf_seek (file, 5 - sfinfo.frames, SEEK_END)) != 5) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_END) returned %d instead of %d.\n", __LINE__, k, 5) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: test_readf_float_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (data [0], orig [5 * channels], margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_END) followed by sf_readf_short failed (%f should be %f).\n", __LINE__, data [0], orig [5 * channels]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: sf_close (file) ; cannam@85: cannam@85: unlink (filename) ; cannam@85: printf ("ok\n") ; cannam@85: } /* lcomp_test_float */ cannam@85: cannam@85: /*-------------------------------------------------------------------------------------------- cannam@85: */ cannam@85: cannam@85: static void cannam@85: lcomp_test_double (const char *filename, int filetype, int channels, double margin) cannam@85: { SNDFILE *file ; cannam@85: SF_INFO sfinfo ; cannam@85: int k, m, seekpos ; cannam@85: long datalen ; cannam@85: double *orig, *data ; cannam@85: double half_max_abs ; cannam@85: cannam@85: print_test_name ("lcomp_test_double", filename) ; cannam@85: cannam@85: datalen = BUFFER_SIZE / channels ; cannam@85: cannam@85: data = data_buffer.d ; cannam@85: orig = orig_buffer.d ; cannam@85: cannam@85: gen_signal_double (orig_buffer.d, 32000.0, channels, datalen) ; cannam@85: for (k = 0 ; k < channels * datalen ; k++) cannam@85: orig [k] = orig_buffer.d [k] ; cannam@85: cannam@85: sfinfo.samplerate = SAMPLE_RATE ; cannam@85: sfinfo.frames = 123456789 ; /* Ridiculous value. */ cannam@85: sfinfo.channels = channels ; cannam@85: sfinfo.format = filetype ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_WRITE, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: sf_command (file, SFC_SET_NORM_DOUBLE, NULL, SF_FALSE) ; cannam@85: test_writef_double_or_die (file, 0, orig, datalen, __LINE__) ; cannam@85: sf_set_string (file, SF_STR_COMMENT, long_comment) ; cannam@85: sf_close (file) ; cannam@85: cannam@85: memset (data, 0, datalen * sizeof (double)) ; cannam@85: cannam@85: if ((filetype & SF_FORMAT_TYPEMASK) != SF_FORMAT_RAW) cannam@85: memset (&sfinfo, 0, sizeof (sfinfo)) ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_READ, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: cannam@85: if ((sfinfo.format & (SF_FORMAT_TYPEMASK | SF_FORMAT_SUBMASK)) != (filetype & (SF_FORMAT_TYPEMASK | SF_FORMAT_SUBMASK))) cannam@85: { printf ("\n\nLine %d: Returned format incorrect (0x%08X => 0x%08X).\n", __LINE__, filetype, sfinfo.format) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames < datalen / channels) cannam@85: { printf ("Too few.frames in file. (%ld should be a little more than %ld)\n", datalen, SF_COUNT_TO_LONG (sfinfo.frames)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames > (datalen + datalen / 20)) cannam@85: { printf ("Too many.frames in file. (%ld should be a little more than %ld)\n", datalen, SF_COUNT_TO_LONG (sfinfo.frames)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.channels != channels) cannam@85: { printf ("Incorrect number of channels in file.\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: check_comment (file, filetype, __LINE__) ; cannam@85: cannam@85: sf_command (file, SFC_SET_NORM_DOUBLE, NULL, SF_FALSE) ; cannam@85: cannam@85: check_log_buffer_or_die (file, __LINE__) ; cannam@85: cannam@85: check_comment (file, filetype, __LINE__) ; cannam@85: cannam@85: sf_command (file, SFC_SET_NORM_DOUBLE, NULL, SF_FALSE) ; cannam@85: cannam@85: test_readf_double_or_die (file, 0, data, datalen, __LINE__) ; cannam@85: cannam@85: half_max_abs = 0.0 ; cannam@85: for (k = 0 ; k < datalen ; k++) cannam@85: { if (error_function (data [k], orig [k], margin)) cannam@85: { printf ("\nLine %d: Incorrect sample A (#%d : %f should be %f).\n", __LINE__, k, data [k], orig [k]) ; cannam@85: oct_save_double (orig, data, datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: half_max_abs = LCT_MAX (half_max_abs, abs (0.5 * data [k])) ; cannam@85: } ; cannam@85: cannam@85: if (half_max_abs < 1.0) cannam@85: { printf ("\n\nLine %d: Signal is all zeros.\n", __LINE__) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_readf_double (file, data, datalen)) != sfinfo.frames - datalen) cannam@85: { printf ("\n\nLine %d: Incorrect read length (%ld should be %d).\n", __LINE__, cannam@85: SF_COUNT_TO_LONG (channels * sfinfo.frames - datalen), k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* This check is only for block based encoders which must append silence cannam@85: ** to the end of a file so as to fill out a block. cannam@85: */ cannam@85: if ((sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_MS_ADPCM) cannam@85: for (k = 0 ; k < sfinfo.frames - datalen ; k++) cannam@85: if (abs (data [channels * k]) > decay_response (channels * k)) cannam@85: { printf ("\n\nLine %d : Incorrect sample B (#%d : abs (%f) should be < %d).\n", __LINE__, channels * k, data [channels * k], decay_response (channels * k)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (! sfinfo.seekable) cannam@85: { sf_close (file) ; cannam@85: unlink (filename) ; cannam@85: printf ("ok\n") ; cannam@85: return ; cannam@85: } ; cannam@85: cannam@85: /* Now test sf_seek function. */ cannam@85: cannam@85: if ((k = sf_seek (file, 0, SEEK_SET)) != 0) cannam@85: { printf ("\n\nLine %d: Seek to start of file failed (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: for (m = 0 ; m < 3 ; m++) cannam@85: { test_readf_double_or_die (file, m, data, 11, __LINE__) ; cannam@85: cannam@85: for (k = 0 ; k < channels * 11 ; k++) cannam@85: if (error_function (data [k], orig [k + channels * m * 11], margin)) cannam@85: { printf ("\nLine %d: Incorrect sample (m = %d) (#%d : %f => %f).\n", __LINE__, m, k + channels * m * 11, orig [k + channels * m * 11], data [k]) ; cannam@85: for (m = 0 ; m < channels ; m++) cannam@85: printf ("%f ", data [m]) ; cannam@85: printf ("\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: } ; cannam@85: cannam@85: seekpos = BUFFER_SIZE / 10 ; cannam@85: cannam@85: /* Check seek from start of file. */ cannam@85: if ((k = sf_seek (file, seekpos, SEEK_SET)) != seekpos) cannam@85: { printf ("Seek to start of file + %d failed (%d).\n", seekpos, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: test_readf_double_or_die (file, 0, data, 1, __LINE__) ; cannam@85: cannam@85: if (error_function (data [0], orig [seekpos * channels], margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_SET) followed by sf_readf_double failed (%f, %f).\n", __LINE__, orig [1], data [0]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_seek (file, 0, SEEK_CUR)) != seekpos + 1) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_CUR) with 0 offset failed (%d should be %d)\n", __LINE__, k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) + BUFFER_SIZE / 5 ; cannam@85: k = sf_seek (file, BUFFER_SIZE / 5, SEEK_CUR) ; cannam@85: test_readf_double_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (data [0], orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (forwards, SEEK_CUR) followed by sf_readf_double failed (%f, %f) (%d, %d).\n", __LINE__, data [0], orig [seekpos * channels], k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) - 20 ; cannam@85: /* Check seek backward from current position. */ cannam@85: k = sf_seek (file, -20, SEEK_CUR) ; cannam@85: test_readf_double_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (data [0], orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (backwards, SEEK_CUR) followed by sf_readf_double failed (%f, %f) (%d, %d).\n", __LINE__, data [0], orig [seekpos * channels], k, seekpos) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check that read past end of file returns number of items. */ cannam@85: sf_seek (file, sfinfo.frames, SEEK_SET) ; cannam@85: cannam@85: if ((k = sf_readf_double (file, data, datalen)) != 0) cannam@85: { printf ("\n\nLine %d: Return value from sf_readf_double past end of file incorrect (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check seek backward from end. */ cannam@85: if ((k = sf_seek (file, 5 - sfinfo.frames, SEEK_END)) != 5) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_END) returned %d instead of %d.\n", __LINE__, k, 5) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: test_readf_double_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (data [0], orig [5 * channels], margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_END) followed by sf_readf_short failed (%f should be %f).\n", __LINE__, data [0], orig [5 * channels]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: sf_close (file) ; cannam@85: cannam@85: unlink (filename) ; cannam@85: printf ("ok\n") ; cannam@85: } /* lcomp_test_double */ cannam@85: cannam@85: /*======================================================================================== cannam@85: ** Smoothed differential loss compression tests. cannam@85: */ cannam@85: cannam@85: static void cannam@85: sdlcomp_test_short (const char *filename, int filetype, int channels, double margin) cannam@85: { SNDFILE *file ; cannam@85: SF_INFO sfinfo ; cannam@85: int k, m, seekpos, half_max_abs ; cannam@85: long datalen ; cannam@85: short *orig, *data, *smooth ; cannam@85: cannam@85: channels = 1 ; cannam@85: print_test_name ("sdlcomp_test_short", filename) ; cannam@85: cannam@85: datalen = BUFFER_SIZE ; cannam@85: cannam@85: orig = orig_buffer.s ; cannam@85: data = data_buffer.s ; cannam@85: smooth = smooth_buffer.s ; cannam@85: cannam@85: gen_signal_double (orig_buffer.d, 32000.0, channels, datalen) ; cannam@85: for (k = 0 ; k < datalen ; k++) cannam@85: orig [k] = lrint (orig_buffer.d [k]) ; cannam@85: cannam@85: sfinfo.samplerate = SAMPLE_RATE ; cannam@85: sfinfo.frames = 123456789 ; /* Ridiculous value. */ cannam@85: sfinfo.channels = channels ; cannam@85: sfinfo.format = filetype ; cannam@85: cannam@85: /* The Vorbis encoder has a bug on PowerPC and X86-64 with sample rates cannam@85: ** <= 22050. Increasing the sample rate to 32000 avoids triggering it. cannam@85: ** See https://trac.xiph.org/ticket/1229 cannam@85: */ cannam@85: if ((file = sf_open (filename, SFM_WRITE, &sfinfo)) == NULL) cannam@85: { const char * errstr ; cannam@85: cannam@85: errstr = sf_strerror (NULL) ; cannam@85: if (strstr (errstr, "Sample rate chosen is known to trigger a Vorbis") == NULL) cannam@85: { printf ("Line %d: sf_open_fd (SFM_WRITE) failed : %s\n", __LINE__, errstr) ; cannam@85: dump_log_buffer (NULL) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: printf ("\n Sample rate -> 32kHz ") ; cannam@85: sfinfo.samplerate = 32000 ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_WRITE, &sfinfo, SF_TRUE, __LINE__) ; cannam@85: } ; cannam@85: cannam@85: test_write_short_or_die (file, 0, orig, datalen, __LINE__) ; cannam@85: sf_set_string (file, SF_STR_COMMENT, long_comment) ; cannam@85: sf_close (file) ; cannam@85: cannam@85: memset (data, 0, datalen * sizeof (short)) ; cannam@85: cannam@85: if ((filetype & SF_FORMAT_TYPEMASK) != SF_FORMAT_RAW) cannam@85: memset (&sfinfo, 0, sizeof (sfinfo)) ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_READ, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: cannam@85: if (sfinfo.format != filetype) cannam@85: { printf ("\n\nLine %d: Returned format incorrect (0x%08X => 0x%08X).\n", __LINE__, filetype, sfinfo.format) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames < datalen / channels) cannam@85: { printf ("Too few.frames in file. (%ld should be a little more than %ld)\n", datalen, SF_COUNT_TO_LONG (sfinfo.frames)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames > (datalen + 400)) cannam@85: { printf ("Too many.frames in file. (%ld should be a little more than %ld)\n", SF_COUNT_TO_LONG (sfinfo.frames), datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.channels != channels) cannam@85: { printf ("Incorrect number of channels in file.\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: check_comment (file, filetype, __LINE__) ; cannam@85: cannam@85: sf_command (file, SFC_SET_NORM_FLOAT, NULL, SF_FALSE) ; cannam@85: cannam@85: check_log_buffer_or_die (file, __LINE__) ; cannam@85: cannam@85: test_readf_short_or_die (file, 0, data, datalen, __LINE__) ; cannam@85: cannam@85: memcpy (smooth, orig, datalen * sizeof (short)) ; cannam@85: smoothed_diff_short (data, datalen) ; cannam@85: smoothed_diff_short (smooth, datalen) ; cannam@85: cannam@85: half_max_abs = 0.0 ; cannam@85: for (k = 0 ; k < datalen ; k++) cannam@85: { if (error_function (1.0 * data [k], 1.0 * smooth [k], margin)) cannam@85: { printf ("\nLine %d: Incorrect sample (#%d : %d should be %d).\n", __LINE__, k, data [k], smooth [k]) ; cannam@85: oct_save_short (orig, smooth, datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: half_max_abs = LCT_MAX (half_max_abs, abs (0.5 * data [k])) ; cannam@85: } ; cannam@85: cannam@85: if (half_max_abs < 1) cannam@85: { printf ("\n\nLine %d: Signal is all zeros.\n", __LINE__) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_read_short (file, data, datalen)) != sfinfo.frames - datalen) cannam@85: { printf ("\n\nLine %d: Incorrect read length (%d should be %ld).\n", __LINE__, k, SF_COUNT_TO_LONG (sfinfo.frames - datalen)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_MS_ADPCM && cannam@85: (sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_GSM610) cannam@85: for (k = 0 ; k < sfinfo.frames - datalen ; k++) cannam@85: if (abs (data [k]) > decay_response (k)) cannam@85: { printf ("\n\nLine %d: Incorrect sample (#%ld : abs (%d) should be < %d).\n", __LINE__, datalen + k, data [k], decay_response (k)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Now test sf_seek function. */ cannam@85: if (sfinfo.seekable) cannam@85: { if ((k = sf_seek (file, 0, SEEK_SET)) != 0) cannam@85: { printf ("\n\nLine %d: Seek to start of file failed (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: for (m = 0 ; m < 3 ; m++) cannam@85: { test_readf_short_or_die (file, m, data, datalen / 7, __LINE__) ; cannam@85: cannam@85: smoothed_diff_short (data, datalen / 7) ; cannam@85: memcpy (smooth, orig + m * datalen / 7, datalen / 7 * sizeof (short)) ; cannam@85: smoothed_diff_short (smooth, datalen / 7) ; cannam@85: cannam@85: for (k = 0 ; k < datalen / 7 ; k++) cannam@85: if (error_function (1.0 * data [k], 1.0 * smooth [k], margin)) cannam@85: { printf ("\nLine %d: Incorrect sample C (#%d (%ld) : %d => %d).\n", __LINE__, k, k + m * (datalen / 7), smooth [k], data [k]) ; cannam@85: for (m = 0 ; m < 10 ; m++) cannam@85: printf ("%d ", data [k]) ; cannam@85: printf ("\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: } ; /* for (m = 0 ; m < 3 ; m++) */ cannam@85: cannam@85: seekpos = BUFFER_SIZE / 10 ; cannam@85: cannam@85: /* Check seek from start of file. */ cannam@85: if ((k = sf_seek (file, seekpos, SEEK_SET)) != seekpos) cannam@85: { printf ("Seek to start of file + %d failed (%d).\n", seekpos, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: test_readf_short_or_die (file, 0, data, 1, __LINE__) ; cannam@85: cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [seekpos * channels], margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_SET) followed by sf_read_short failed (%d, %d).\n", __LINE__, orig [1], data [0]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_seek (file, 0, SEEK_CUR)) != seekpos + 1) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_CUR) with 0 offset failed (%d should be %d)\n", __LINE__, k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) + BUFFER_SIZE / 5 ; cannam@85: k = sf_seek (file, BUFFER_SIZE / 5, SEEK_CUR) ; cannam@85: test_readf_short_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (forwards, SEEK_CUR) followed by sf_read_short failed (%d, %d) (%d, %d).\n", __LINE__, data [0], orig [seekpos * channels], k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) - 20 ; cannam@85: /* Check seek backward from current position. */ cannam@85: k = sf_seek (file, -20, SEEK_CUR) ; cannam@85: test_readf_short_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (backwards, SEEK_CUR) followed by sf_read_short failed (%d, %d) (%d, %d).\n", __LINE__, data [0], orig [seekpos * channels], k, seekpos) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check that read past end of file returns number of items. */ cannam@85: sf_seek (file, sfinfo.frames, SEEK_SET) ; cannam@85: cannam@85: if ((k = sf_read_short (file, data, datalen)) != 0) cannam@85: { printf ("\n\nLine %d: Return value from sf_read_short past end of file incorrect (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check seek backward from end. */ cannam@85: cannam@85: if ((k = sf_seek (file, 5 - sfinfo.frames, SEEK_END)) != 5) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_END) returned %d instead of %d.\n", __LINE__, k, 5) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: test_read_short_or_die (file, 0, data, channels, __LINE__) ; cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [5 * channels], margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_END) followed by sf_read_short failed (%d should be %d).\n", __LINE__, data [0], orig [5 * channels]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: } /* if (sfinfo.seekable) */ cannam@85: cannam@85: sf_close (file) ; cannam@85: cannam@85: unlink (filename) ; cannam@85: printf ("ok\n") ; cannam@85: } /* sdlcomp_test_short */ cannam@85: cannam@85: static void cannam@85: sdlcomp_test_int (const char *filename, int filetype, int channels, double margin) cannam@85: { SNDFILE *file ; cannam@85: SF_INFO sfinfo ; cannam@85: int k, m, seekpos, half_max_abs ; cannam@85: long datalen ; cannam@85: int *orig, *data, *smooth ; cannam@85: double scale ; cannam@85: cannam@85: channels = 1 ; cannam@85: cannam@85: print_test_name ("sdlcomp_test_int", filename) ; cannam@85: cannam@85: datalen = BUFFER_SIZE ; cannam@85: scale = 1.0 * 0x10000 ; cannam@85: cannam@85: orig = orig_buffer.i ; cannam@85: data = data_buffer.i ; cannam@85: smooth = smooth_buffer.i ; cannam@85: cannam@85: gen_signal_double (orig_buffer.d, 32000.0 * scale, channels, datalen) ; cannam@85: for (k = 0 ; k < datalen ; k++) cannam@85: orig [k] = lrint (orig_buffer.d [k]) ; cannam@85: cannam@85: sfinfo.samplerate = SAMPLE_RATE ; cannam@85: sfinfo.frames = 123456789 ; /* Ridiculous value. */ cannam@85: sfinfo.channels = channels ; cannam@85: sfinfo.format = filetype ; cannam@85: cannam@85: /* The Vorbis encoder has a bug on PowerPC and X86-64 with sample rates cannam@85: ** <= 22050. Increasing the sample rate to 32000 avoids triggering it. cannam@85: ** See https://trac.xiph.org/ticket/1229 cannam@85: */ cannam@85: if ((file = sf_open (filename, SFM_WRITE, &sfinfo)) == NULL) cannam@85: { const char * errstr ; cannam@85: cannam@85: errstr = sf_strerror (NULL) ; cannam@85: if (strstr (errstr, "Sample rate chosen is known to trigger a Vorbis") == NULL) cannam@85: { printf ("Line %d: sf_open_fd (SFM_WRITE) failed : %s\n", __LINE__, errstr) ; cannam@85: dump_log_buffer (NULL) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: printf ("\n Sample rate -> 32kHz ") ; cannam@85: sfinfo.samplerate = 32000 ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_WRITE, &sfinfo, SF_TRUE, __LINE__) ; cannam@85: } ; cannam@85: cannam@85: test_writef_int_or_die (file, 0, orig, datalen, __LINE__) ; cannam@85: sf_set_string (file, SF_STR_COMMENT, long_comment) ; cannam@85: sf_close (file) ; cannam@85: cannam@85: memset (data, 0, datalen * sizeof (int)) ; cannam@85: cannam@85: if ((filetype & SF_FORMAT_TYPEMASK) != SF_FORMAT_RAW) cannam@85: memset (&sfinfo, 0, sizeof (sfinfo)) ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_READ, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: cannam@85: if (sfinfo.format != filetype) cannam@85: { printf ("Returned format incorrect (0x%08X => 0x%08X).\n", filetype, sfinfo.format) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames < datalen / channels) cannam@85: { printf ("Too few.frames in file. (%ld should be a little more than %ld)\n", datalen, SF_COUNT_TO_LONG (sfinfo.frames)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames > (datalen + 400)) cannam@85: { printf ("Too many.frames in file. (%ld should be a little more than %ld)\n", SF_COUNT_TO_LONG (sfinfo.frames), datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.channels != channels) cannam@85: { printf ("Incorrect number of channels in file.\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: check_log_buffer_or_die (file, __LINE__) ; cannam@85: cannam@85: test_readf_int_or_die (file, 0, data, datalen, __LINE__) ; cannam@85: cannam@85: memcpy (smooth, orig, datalen * sizeof (int)) ; cannam@85: smoothed_diff_int (data, datalen) ; cannam@85: smoothed_diff_int (smooth, datalen) ; cannam@85: cannam@85: half_max_abs = abs (data [0] >> 16) ; cannam@85: for (k = 1 ; k < datalen ; k++) cannam@85: { if (error_function (data [k] / scale, smooth [k] / scale, margin)) cannam@85: { printf ("\nLine %d: Incorrect sample (#%d : %d should be %d).\n", __LINE__, k, data [k], smooth [k]) ; cannam@85: oct_save_int (orig, smooth, datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: half_max_abs = LCT_MAX (half_max_abs, abs (data [k] / 2)) ; cannam@85: } ; cannam@85: cannam@85: if (half_max_abs < 1) cannam@85: { printf ("\n\nLine %d: Signal is all zeros.\n", __LINE__) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_readf_int (file, data, datalen)) != sfinfo.frames - datalen) cannam@85: { printf ("\n\nLine %d: Incorrect read length (%d should be %ld).\n", __LINE__, k, SF_COUNT_TO_LONG (sfinfo.frames - datalen)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_IMA_ADPCM && cannam@85: (sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_MS_ADPCM && cannam@85: (sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_GSM610 && cannam@85: (sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_G721_32 && cannam@85: (sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_G723_24) cannam@85: for (k = 0 ; k < sfinfo.frames - datalen ; k++) cannam@85: if (abs (data [k]) > decay_response (k)) cannam@85: { printf ("\n\nLine %d: Incorrect sample (#%ld : abs (%d) should be < %d).\n", __LINE__, datalen + k, data [k], decay_response (k)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Now test sf_seek function. */ cannam@85: if (sfinfo.seekable) cannam@85: { if ((k = sf_seek (file, 0, SEEK_SET)) != 0) cannam@85: { printf ("\n\nLine %d: Seek to start of file failed (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: for (m = 0 ; m < 3 ; m++) cannam@85: { test_readf_int_or_die (file, m, data, datalen / 7, __LINE__) ; cannam@85: cannam@85: smoothed_diff_int (data, datalen / 7) ; cannam@85: memcpy (smooth, orig + m * datalen / 7, datalen / 7 * sizeof (int)) ; cannam@85: smoothed_diff_int (smooth, datalen / 7) ; cannam@85: cannam@85: for (k = 0 ; k < datalen / 7 ; k++) cannam@85: if (error_function (data [k] / scale, smooth [k] / scale, margin)) cannam@85: { printf ("\nLine %d: Incorrect sample (#%d (%ld) : %d => %d).\n", __LINE__, k, k + m * (datalen / 7), smooth [k], data [k]) ; cannam@85: for (m = 0 ; m < 10 ; m++) cannam@85: printf ("%d ", data [k]) ; cannam@85: printf ("\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: } ; /* for (m = 0 ; m < 3 ; m++) */ cannam@85: cannam@85: seekpos = BUFFER_SIZE / 10 ; cannam@85: cannam@85: /* Check seek from start of file. */ cannam@85: if ((k = sf_seek (file, seekpos, SEEK_SET)) != seekpos) cannam@85: { printf ("Seek to start of file + %d failed (%d).\n", seekpos, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: test_readf_int_or_die (file, 0, data, 1, __LINE__) ; cannam@85: cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [seekpos * channels], margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_SET) followed by sf_readf_int failed (%d, %d).\n", __LINE__, orig [1], data [0]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_seek (file, 0, SEEK_CUR)) != seekpos + 1) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_CUR) with 0 offset failed (%d should be %d)\n", __LINE__, k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) + BUFFER_SIZE / 5 ; cannam@85: k = sf_seek (file, BUFFER_SIZE / 5, SEEK_CUR) ; cannam@85: test_readf_int_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (forwards, SEEK_CUR) followed by sf_readf_int failed (%d, %d) (%d, %d).\n", __LINE__, data [0], orig [seekpos * channels], k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) - 20 ; cannam@85: /* Check seek backward from current position. */ cannam@85: k = sf_seek (file, -20, SEEK_CUR) ; cannam@85: test_readf_int_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (1.0 * data [0], 1.0 * orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (backwards, SEEK_CUR) followed by sf_readf_int failed (%d, %d) (%d, %d).\n", __LINE__, data [0], orig [seekpos * channels], k, seekpos) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check that read past end of file returns number of items. */ cannam@85: sf_seek (file, sfinfo.frames, SEEK_SET) ; cannam@85: cannam@85: if ((k = sf_readf_int (file, data, datalen)) != 0) cannam@85: { printf ("\n\nLine %d: Return value from sf_readf_int past end of file incorrect (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check seek backward from end. */ cannam@85: cannam@85: if ((k = sf_seek (file, 5 - sfinfo.frames, SEEK_END)) != 5) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_END) returned %d instead of %d.\n", __LINE__, k, 5) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: test_readf_int_or_die (file, 0, data, 1, __LINE__) ; cannam@85: if (error_function (data [0] / scale, orig [5] / scale, margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_END) followed by sf_readf_int failed (%d should be %d).\n", __LINE__, data [0], orig [5]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: } /* if (sfinfo.seekable) */ cannam@85: cannam@85: sf_close (file) ; cannam@85: cannam@85: unlink (filename) ; cannam@85: printf ("ok\n") ; cannam@85: } /* sdlcomp_test_int */ cannam@85: cannam@85: static void cannam@85: sdlcomp_test_float (const char *filename, int filetype, int channels, double margin) cannam@85: { SNDFILE *file ; cannam@85: SF_INFO sfinfo ; cannam@85: int k, m, seekpos ; cannam@85: long datalen ; cannam@85: float *orig, *data, *smooth ; cannam@85: double half_max_abs ; cannam@85: cannam@85: channels = 1 ; cannam@85: cannam@85: print_test_name ("sdlcomp_test_float", filename) ; cannam@85: cannam@85: if ((filetype & SF_FORMAT_SUBMASK) == SF_FORMAT_VORBIS) cannam@85: { puts ("Not working for this format.") ; cannam@85: return ; cannam@85: } ; cannam@85: cannam@85: printf ("** fix this ** ") ; cannam@85: cannam@85: datalen = BUFFER_SIZE ; cannam@85: cannam@85: orig = orig_buffer.f ; cannam@85: data = data_buffer.f ; cannam@85: smooth = smooth_buffer.f ; cannam@85: cannam@85: gen_signal_double (orig_buffer.d, 32000.0, channels, datalen) ; cannam@85: for (k = 0 ; k < datalen ; k++) cannam@85: orig [k] = lrint (orig_buffer.d [k]) ; cannam@85: cannam@85: sfinfo.samplerate = SAMPLE_RATE ; cannam@85: sfinfo.frames = 123456789 ; /* Ridiculous value. */ cannam@85: sfinfo.channels = channels ; cannam@85: sfinfo.format = filetype ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_WRITE, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: sf_command (file, SFC_SET_NORM_FLOAT, NULL, SF_FALSE) ; cannam@85: test_write_float_or_die (file, 0, orig, datalen, __LINE__) ; cannam@85: sf_set_string (file, SF_STR_COMMENT, long_comment) ; cannam@85: sf_close (file) ; cannam@85: cannam@85: memset (data, 0, datalen * sizeof (float)) ; cannam@85: cannam@85: if ((filetype & SF_FORMAT_TYPEMASK) != SF_FORMAT_RAW) cannam@85: memset (&sfinfo, 0, sizeof (sfinfo)) ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_READ, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: cannam@85: if ((sfinfo.format & (SF_FORMAT_TYPEMASK | SF_FORMAT_SUBMASK)) != (filetype & (SF_FORMAT_TYPEMASK | SF_FORMAT_SUBMASK))) cannam@85: { printf ("\n\nLine %d: Returned format incorrect (0x%08X => 0x%08X).\n", __LINE__, filetype, sfinfo.format) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames < datalen / channels) cannam@85: { printf ("Too few.frames in file. (%ld should be a little more than %ld)\n", datalen, SF_COUNT_TO_LONG (sfinfo.frames)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames > (datalen + 400)) cannam@85: { printf ("Too many.frames in file. (%ld should be a little more than %ld)\n", SF_COUNT_TO_LONG (sfinfo.frames), datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.channels != channels) cannam@85: { printf ("Incorrect number of channels in file.\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: check_comment (file, filetype, __LINE__) ; cannam@85: cannam@85: sf_command (file, SFC_SET_NORM_FLOAT, NULL, SF_FALSE) ; cannam@85: cannam@85: check_log_buffer_or_die (file, __LINE__) ; cannam@85: cannam@85: test_read_float_or_die (file, 0, data, datalen, __LINE__) ; cannam@85: cannam@85: memcpy (smooth, orig, datalen * sizeof (float)) ; cannam@85: smoothed_diff_float (data, datalen) ; cannam@85: smoothed_diff_float (smooth, datalen) ; cannam@85: cannam@85: half_max_abs = fabs (data [0]) ; cannam@85: for (k = 1 ; k < datalen ; k++) cannam@85: { if (error_function (data [k], smooth [k], margin)) cannam@85: { printf ("\nLine %d: Incorrect sample (#%d : %d should be %d).\n", __LINE__, k, (int) data [k], (int) smooth [k]) ; cannam@85: oct_save_float (orig, smooth, datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: half_max_abs = LCT_MAX (half_max_abs, abs (0.5 * data [k])) ; cannam@85: } ; cannam@85: cannam@85: if (half_max_abs <= 0.0) cannam@85: { printf ("\n\nLine %d: Signal is all zeros.\n", __LINE__) ; cannam@85: printf ("half_max_abs : % 10.6f\n", half_max_abs) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_read_float (file, data, datalen)) != sfinfo.frames - datalen) cannam@85: { printf ("\n\nLine %d: Incorrect read length (%d should be %ld).\n", __LINE__, k, SF_COUNT_TO_LONG (sfinfo.frames - datalen)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_MS_ADPCM && cannam@85: (sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_GSM610) cannam@85: for (k = 0 ; k < sfinfo.frames - datalen ; k++) cannam@85: if (abs (data [k]) > decay_response (k)) cannam@85: { printf ("\n\nLine %d: Incorrect sample (#%ld : abs (%d) should be < %d).\n", __LINE__, datalen + k, (int) data [k], (int) decay_response (k)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Now test sf_seek function. */ cannam@85: if (sfinfo.seekable) cannam@85: { if ((k = sf_seek (file, 0, SEEK_SET)) != 0) cannam@85: { printf ("\n\nLine %d: Seek to start of file failed (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: for (m = 0 ; m < 3 ; m++) cannam@85: { test_read_float_or_die (file, 0, data, datalen / 7, __LINE__) ; cannam@85: cannam@85: smoothed_diff_float (data, datalen / 7) ; cannam@85: memcpy (smooth, orig + m * datalen / 7, datalen / 7 * sizeof (float)) ; cannam@85: smoothed_diff_float (smooth, datalen / 7) ; cannam@85: cannam@85: for (k = 0 ; k < datalen / 7 ; k++) cannam@85: if (error_function (data [k], smooth [k], margin)) cannam@85: { printf ("\nLine %d: Incorrect sample C (#%d (%ld) : %d => %d).\n", __LINE__, k, k + m * (datalen / 7), (int) smooth [k], (int) data [k]) ; cannam@85: for (m = 0 ; m < 10 ; m++) cannam@85: printf ("%d ", (int) data [k]) ; cannam@85: printf ("\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: } ; /* for (m = 0 ; m < 3 ; m++) */ cannam@85: cannam@85: seekpos = BUFFER_SIZE / 10 ; cannam@85: cannam@85: /* Check seek from start of file. */ cannam@85: if ((k = sf_seek (file, seekpos, SEEK_SET)) != seekpos) cannam@85: { printf ("Seek to start of file + %d failed (%d).\n", seekpos, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: test_read_float_or_die (file, 0, data, channels, __LINE__) ; cannam@85: cannam@85: if (error_function (data [0], orig [seekpos * channels], margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_SET) followed by sf_read_float failed (%d, %d).\n", __LINE__, (int) orig [1], (int) data [0]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_seek (file, 0, SEEK_CUR)) != seekpos + 1) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_CUR) with 0 offset failed (%d should be %d)\n", __LINE__, k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) + BUFFER_SIZE / 5 ; cannam@85: k = sf_seek (file, BUFFER_SIZE / 5, SEEK_CUR) ; cannam@85: test_read_float_or_die (file, 0, data, channels, __LINE__) ; cannam@85: if (error_function (data [0], orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (forwards, SEEK_CUR) followed by sf_read_float failed (%d, %d) (%d, %d).\n", __LINE__, (int) data [0], (int) orig [seekpos * channels], k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) - 20 ; cannam@85: /* Check seek backward from current position. */ cannam@85: k = sf_seek (file, -20, SEEK_CUR) ; cannam@85: test_read_float_or_die (file, 0, data, channels, __LINE__) ; cannam@85: if (error_function (data [0], orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (backwards, SEEK_CUR) followed by sf_read_float failed (%d, %d) (%d, %d).\n", __LINE__, (int) data [0], (int) orig [seekpos * channels], k, seekpos) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check that read past end of file returns number of items. */ cannam@85: sf_seek (file, sfinfo.frames, SEEK_SET) ; cannam@85: cannam@85: if ((k = sf_read_float (file, data, datalen)) != 0) cannam@85: { printf ("\n\nLine %d: Return value from sf_read_float past end of file incorrect (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check seek backward from end. */ cannam@85: cannam@85: if ((k = sf_seek (file, 5 - sfinfo.frames, SEEK_END)) != 5) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_END) returned %d instead of %d.\n", __LINE__, k, 5) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: test_read_float_or_die (file, 0, data, channels, __LINE__) ; cannam@85: if (error_function (data [0], orig [5 * channels], margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_END) followed by sf_read_float failed (%f should be %f).\n", __LINE__, data [0], orig [5 * channels]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: } /* if (sfinfo.seekable) */ cannam@85: cannam@85: sf_close (file) ; cannam@85: cannam@85: unlink (filename) ; cannam@85: printf ("ok\n") ; cannam@85: } /* sdlcomp_test_float */ cannam@85: cannam@85: static void cannam@85: sdlcomp_test_double (const char *filename, int filetype, int channels, double margin) cannam@85: { SNDFILE *file ; cannam@85: SF_INFO sfinfo ; cannam@85: int k, m, seekpos ; cannam@85: long datalen ; cannam@85: double *orig, *data, *smooth, half_max_abs ; cannam@85: cannam@85: channels = 1 ; cannam@85: print_test_name ("sdlcomp_test_double", filename) ; cannam@85: cannam@85: if ((filetype & SF_FORMAT_SUBMASK) == SF_FORMAT_VORBIS) cannam@85: { puts ("Not working for this format.") ; cannam@85: return ; cannam@85: } ; cannam@85: cannam@85: datalen = BUFFER_SIZE ; cannam@85: cannam@85: orig = orig_buffer.d ; cannam@85: data = data_buffer.d ; cannam@85: smooth = smooth_buffer.d ; cannam@85: cannam@85: gen_signal_double (orig_buffer.d, 32000.0, channels, datalen) ; cannam@85: cannam@85: sfinfo.samplerate = SAMPLE_RATE ; cannam@85: sfinfo.frames = 123456789 ; /* Ridiculous value. */ cannam@85: sfinfo.channels = channels ; cannam@85: sfinfo.format = filetype ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_WRITE, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: sf_command (file, SFC_SET_NORM_DOUBLE, NULL, SF_FALSE) ; cannam@85: test_write_double_or_die (file, 0, orig, datalen, __LINE__) ; cannam@85: sf_set_string (file, SF_STR_COMMENT, long_comment) ; cannam@85: sf_close (file) ; cannam@85: cannam@85: memset (data, 0, datalen * sizeof (double)) ; cannam@85: cannam@85: if ((filetype & SF_FORMAT_TYPEMASK) != SF_FORMAT_RAW) cannam@85: memset (&sfinfo, 0, sizeof (sfinfo)) ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_READ, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: cannam@85: if (sfinfo.format != filetype) cannam@85: { printf ("Returned format incorrect (0x%08X => 0x%08X).\n", filetype, sfinfo.format) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames < datalen / channels) cannam@85: { printf ("Too few.frames in file. (%ld should be a little more than %ld)\n", datalen, SF_COUNT_TO_LONG (sfinfo.frames)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames > (datalen + 400)) cannam@85: { printf ("Too many.frames in file. (%ld should be a little more than %ld)\n", SF_COUNT_TO_LONG (sfinfo.frames), datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.channels != channels) cannam@85: { printf ("Incorrect number of channels in file.\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: check_comment (file, filetype, __LINE__) ; cannam@85: cannam@85: check_comment (file, filetype, __LINE__) ; cannam@85: cannam@85: sf_command (file, SFC_SET_NORM_DOUBLE, NULL, SF_FALSE) ; cannam@85: cannam@85: check_log_buffer_or_die (file, __LINE__) ; cannam@85: cannam@85: test_read_double_or_die (file, 0, data, datalen, __LINE__) ; cannam@85: cannam@85: memcpy (smooth, orig, datalen * sizeof (double)) ; cannam@85: smoothed_diff_double (data, datalen) ; cannam@85: smoothed_diff_double (smooth, datalen) ; cannam@85: cannam@85: half_max_abs = 0.0 ; cannam@85: for (k = 0 ; k < datalen ; k++) cannam@85: { if (error_function (data [k], smooth [k], margin)) cannam@85: { printf ("\n\nLine %d: Incorrect sample (#%d : %d should be %d).\n", __LINE__, k, (int) data [k], (int) smooth [k]) ; cannam@85: oct_save_double (orig, smooth, datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: half_max_abs = LCT_MAX (half_max_abs, 0.5 * fabs (data [k])) ; cannam@85: } ; cannam@85: cannam@85: if (half_max_abs < 1.0) cannam@85: { printf ("\n\nLine %d: Signal is all zeros.\n", __LINE__) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_read_double (file, data, datalen)) != sfinfo.frames - datalen) cannam@85: { printf ("\n\nLine %d: Incorrect read length (%d should be %ld).\n", __LINE__, k, SF_COUNT_TO_LONG (sfinfo.frames - datalen)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_MS_ADPCM && cannam@85: (sfinfo.format & SF_FORMAT_SUBMASK) != SF_FORMAT_GSM610) cannam@85: for (k = 0 ; k < sfinfo.frames - datalen ; k++) cannam@85: if (abs (data [k]) > decay_response (k)) cannam@85: { printf ("\n\nLine %d: Incorrect sample (#%ld : abs (%d) should be < %d).\n", __LINE__, datalen + k, (int) data [k], (int) decay_response (k)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Now test sf_seek function. */ cannam@85: if (sfinfo.seekable) cannam@85: { if ((k = sf_seek (file, 0, SEEK_SET)) != 0) cannam@85: { printf ("\n\nLine %d: Seek to start of file failed (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: for (m = 0 ; m < 3 ; m++) cannam@85: { test_read_double_or_die (file, m, data, datalen / 7, __LINE__) ; cannam@85: cannam@85: smoothed_diff_double (data, datalen / 7) ; cannam@85: memcpy (smooth, orig + m * datalen / 7, datalen / 7 * sizeof (double)) ; cannam@85: smoothed_diff_double (smooth, datalen / 7) ; cannam@85: cannam@85: for (k = 0 ; k < datalen / 7 ; k++) cannam@85: if (error_function (data [k], smooth [k], margin)) cannam@85: { printf ("\nLine %d: Incorrect sample C (#%d (%ld) : %d => %d).\n", __LINE__, k, k + m * (datalen / 7), (int) smooth [k], (int) data [k]) ; cannam@85: for (m = 0 ; m < 10 ; m++) cannam@85: printf ("%d ", (int) data [k]) ; cannam@85: printf ("\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: } ; /* for (m = 0 ; m < 3 ; m++) */ cannam@85: cannam@85: seekpos = BUFFER_SIZE / 10 ; cannam@85: cannam@85: /* Check seek from start of file. */ cannam@85: if ((k = sf_seek (file, seekpos, SEEK_SET)) != seekpos) cannam@85: { printf ("Seek to start of file + %d failed (%d).\n", seekpos, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: test_read_double_or_die (file, 0, data, channels, __LINE__) ; cannam@85: cannam@85: if (error_function (data [0], orig [seekpos * channels], margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_SET) followed by sf_read_double failed (%d, %d).\n", __LINE__, (int) orig [1], (int) data [0]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if ((k = sf_seek (file, 0, SEEK_CUR)) != seekpos + 1) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_CUR) with 0 offset failed (%d should be %d)\n", __LINE__, k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) + BUFFER_SIZE / 5 ; cannam@85: k = sf_seek (file, BUFFER_SIZE / 5, SEEK_CUR) ; cannam@85: test_read_double_or_die (file, 0, data, channels, __LINE__) ; cannam@85: if (error_function (data [0], orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (forwards, SEEK_CUR) followed by sf_read_double failed (%d, %d) (%d, %d).\n", __LINE__, (int) data [0], (int) orig [seekpos * channels], k, seekpos + 1) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: seekpos = sf_seek (file, 0, SEEK_CUR) - 20 ; cannam@85: /* Check seek backward from current position. */ cannam@85: k = sf_seek (file, -20, SEEK_CUR) ; cannam@85: test_read_double_or_die (file, 0, data, channels, __LINE__) ; cannam@85: if (error_function (data [0], orig [seekpos * channels], margin) || k != seekpos) cannam@85: { printf ("\nLine %d: sf_seek (backwards, SEEK_CUR) followed by sf_read_double failed (%d, %d) (%d, %d).\n", __LINE__, (int) data [0], (int) orig [seekpos * channels], k, seekpos) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check that read past end of file returns number of items. */ cannam@85: sf_seek (file, sfinfo.frames, SEEK_SET) ; cannam@85: cannam@85: if ((k = sf_read_double (file, data, datalen)) != 0) cannam@85: { printf ("\n\nLine %d: Return value from sf_read_double past end of file incorrect (%d).\n", __LINE__, k) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: /* Check seek backward from end. */ cannam@85: cannam@85: if ((k = sf_seek (file, 5 - sfinfo.frames, SEEK_END)) != 5) cannam@85: { printf ("\n\nLine %d: sf_seek (SEEK_END) returned %d instead of %d.\n", __LINE__, k, 5) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: test_read_double_or_die (file, 0, data, channels, __LINE__) ; cannam@85: if (error_function (data [0], orig [5 * channels], margin)) cannam@85: { printf ("\nLine %d: sf_seek (SEEK_END) followed by sf_read_double failed (%f should be %f).\n", __LINE__, data [0], orig [5 * channels]) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: } /* if (sfinfo.seekable) */ cannam@85: cannam@85: sf_close (file) ; cannam@85: cannam@85: unlink (filename) ; cannam@85: printf ("ok\n") ; cannam@85: } /* sdlcomp_test_double */ cannam@85: cannam@85: static void cannam@85: read_raw_test (const char *filename, int filetype, int channels) cannam@85: { SNDFILE *file ; cannam@85: SF_INFO sfinfo ; cannam@85: sf_count_t count ; cannam@85: long datalen ; cannam@85: short *orig, *data ; cannam@85: int k ; cannam@85: cannam@85: print_test_name ("read_raw_test", filename) ; cannam@85: cannam@85: datalen = ARRAY_LEN (orig_buffer.s) / 2 ; cannam@85: cannam@85: orig = orig_buffer.s ; cannam@85: data = data_buffer.s ; cannam@85: cannam@85: gen_signal_double (orig_buffer.d, 32000.0, channels, datalen) ; cannam@85: for (k = 0 ; k < datalen ; k++) cannam@85: orig [k] = lrint (orig_buffer.d [k]) ; cannam@85: cannam@85: sfinfo.samplerate = SAMPLE_RATE ; cannam@85: sfinfo.frames = 123456789 ; /* Ridiculous value. */ cannam@85: sfinfo.channels = channels ; cannam@85: sfinfo.format = filetype ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_WRITE, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: test_write_short_or_die (file, 0, orig, datalen, __LINE__) ; cannam@85: sf_set_string (file, SF_STR_COMMENT, long_comment) ; cannam@85: sf_close (file) ; cannam@85: cannam@85: memset (data, 0, datalen * sizeof (double)) ; cannam@85: cannam@85: if ((filetype & SF_FORMAT_TYPEMASK) != SF_FORMAT_RAW) cannam@85: memset (&sfinfo, 0, sizeof (sfinfo)) ; cannam@85: cannam@85: file = test_open_file_or_die (filename, SFM_READ, &sfinfo, SF_FALSE, __LINE__) ; cannam@85: cannam@85: if (sfinfo.format != filetype) cannam@85: { printf ("Returned format incorrect (0x%08X => 0x%08X).\n", filetype, sfinfo.format) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames < datalen / channels) cannam@85: { printf ("Too few.frames in file. (%ld should be a little more than %ld)\n", datalen, SF_COUNT_TO_LONG (sfinfo.frames)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.frames > (datalen + 400)) cannam@85: { printf ("Too many.frames in file. (%ld should be a little more than %ld)\n", SF_COUNT_TO_LONG (sfinfo.frames), datalen) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (sfinfo.channels != channels) cannam@85: { printf ("Incorrect number of channels in file.\n") ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: check_comment (file, filetype, __LINE__) ; cannam@85: cannam@85: count = sf_read_raw (file, orig_buffer.c, datalen + 5 * channels) ; cannam@85: if (count != sfinfo.channels * sfinfo.frames) cannam@85: { printf ("\nLine %d : sf_read_raw returned %ld should be %ld\n", __LINE__, SF_COUNT_TO_LONG (count), sfinfo.channels * SF_COUNT_TO_LONG (sfinfo.frames)) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: sf_close (file) ; cannam@85: cannam@85: unlink (filename) ; cannam@85: printf ("ok\n") ; cannam@85: } /* read_raw_test */ cannam@85: cannam@85: /*======================================================================================== cannam@85: ** Auxiliary functions cannam@85: */ cannam@85: cannam@85: #define SIGNAL_MAXVAL 30000.0 cannam@85: #define DECAY_COUNT 1000 cannam@85: cannam@85: static int cannam@85: decay_response (int k) cannam@85: { if (k < 1) cannam@85: return (int) (1.2 * SIGNAL_MAXVAL) ; cannam@85: if (k > DECAY_COUNT) cannam@85: return 0 ; cannam@85: return (int) (1.2 * SIGNAL_MAXVAL * (DECAY_COUNT - k) / (1.0 * DECAY_COUNT)) ; cannam@85: } /* decay_response */ cannam@85: cannam@85: static void cannam@85: gen_signal_double (double *data, double scale, int channels, int datalen) cannam@85: { int k, ramplen ; cannam@85: double amp = 0.0 ; cannam@85: cannam@85: ramplen = DECAY_COUNT ; cannam@85: cannam@85: if (channels == 1) cannam@85: { for (k = 0 ; k < datalen ; k++) cannam@85: { if (k <= ramplen) cannam@85: amp = scale * k / ((double) ramplen) ; cannam@85: else if (k > datalen - ramplen) cannam@85: amp = scale * (datalen - k) / ((double) ramplen) ; cannam@85: cannam@85: /*-printf ("%3d : %g\n", k, amp) ;-*/ cannam@85: cannam@85: data [k] = amp * (0.4 * sin (33.3 * 2.0 * M_PI * ((double) (k+1)) / ((double) SAMPLE_RATE)) cannam@85: + 0.3 * cos (201.1 * 2.0 * M_PI * ((double) (k+1)) / ((double) SAMPLE_RATE))) ; cannam@85: } ; cannam@85: } cannam@85: else cannam@85: { for (k = 0 ; k < datalen ; k ++) cannam@85: { if (k <= ramplen) cannam@85: amp = scale * k / ((double) ramplen) ; cannam@85: else if (k > datalen - ramplen) cannam@85: amp = scale * (datalen - k) / ((double) ramplen) ; cannam@85: cannam@85: data [2 * k] = amp * (0.4 * sin (33.3 * 2.0 * M_PI * ((double) (k+1)) / ((double) SAMPLE_RATE)) cannam@85: + 0.3 * cos (201.1 * 2.0 * M_PI * ((double) (k+1)) / ((double) SAMPLE_RATE))) ; cannam@85: data [2 * k + 1] = amp * (0.4 * sin (55.5 * 2.0 * M_PI * ((double) (k+1)) / ((double) SAMPLE_RATE)) cannam@85: + 0.3 * cos (201.1 * 2.0 * M_PI * ((double) (k+1)) / ((double) SAMPLE_RATE))) ; cannam@85: } ; cannam@85: } ; cannam@85: cannam@85: return ; cannam@85: } /* gen_signal_double */ cannam@85: cannam@85: static int cannam@85: error_function (double data, double orig, double margin) cannam@85: { double error ; cannam@85: cannam@85: if (fabs (orig) <= 500.0) cannam@85: error = fabs (fabs (data) - fabs (orig)) / 2000.0 ; cannam@85: else if (fabs (orig) <= 1000.0) cannam@85: error = fabs (data - orig) / 3000.0 ; cannam@85: else cannam@85: error = fabs (data - orig) / fabs (orig) ; cannam@85: cannam@85: if (error > margin) cannam@85: { printf ("\n\nerror_function (data = %f, orig = %f, margin = %f) -> %f\n", data, orig, margin, error) ; cannam@85: return 1 ; cannam@85: } ; cannam@85: return 0 ; cannam@85: } /* error_function */ cannam@85: cannam@85: static void cannam@85: smoothed_diff_short (short *data, unsigned int datalen) cannam@85: { unsigned int k ; cannam@85: double memory = 0.0 ; cannam@85: cannam@85: /* Calculate the smoothed sample-to-sample difference. */ cannam@85: for (k = 0 ; k < datalen - 1 ; k++) cannam@85: { memory = 0.7 * memory + (1 - 0.7) * (double) (data [k+1] - data [k]) ; cannam@85: data [k] = (short) memory ; cannam@85: } ; cannam@85: data [datalen-1] = data [datalen-2] ; cannam@85: cannam@85: } /* smoothed_diff_short */ cannam@85: cannam@85: static void cannam@85: smoothed_diff_int (int *data, unsigned int datalen) cannam@85: { unsigned int k ; cannam@85: double memory = 0.0 ; cannam@85: cannam@85: /* Calculate the smoothed sample-to-sample difference. */ cannam@85: for (k = 0 ; k < datalen - 1 ; k++) cannam@85: { memory = 0.7 * memory + (1 - 0.7) * (double) (data [k+1] - data [k]) ; cannam@85: data [k] = (int) memory ; cannam@85: } ; cannam@85: data [datalen-1] = data [datalen-2] ; cannam@85: cannam@85: } /* smoothed_diff_int */ cannam@85: cannam@85: static void cannam@85: smoothed_diff_float (float *data, unsigned int datalen) cannam@85: { unsigned int k ; cannam@85: float memory = 0.0 ; cannam@85: cannam@85: /* Calculate the smoothed sample-to-sample difference. */ cannam@85: for (k = 0 ; k < datalen - 1 ; k++) cannam@85: { memory = 0.7 * memory + (1 - 0.7) * (data [k+1] - data [k]) ; cannam@85: data [k] = memory ; cannam@85: } ; cannam@85: data [datalen-1] = data [datalen-2] ; cannam@85: cannam@85: } /* smoothed_diff_float */ cannam@85: cannam@85: static void cannam@85: smoothed_diff_double (double *data, unsigned int datalen) cannam@85: { unsigned int k ; cannam@85: double memory = 0.0 ; cannam@85: cannam@85: /* Calculate the smoothed sample-to-sample difference. */ cannam@85: for (k = 0 ; k < datalen - 1 ; k++) cannam@85: { memory = 0.7 * memory + (1 - 0.7) * (data [k+1] - data [k]) ; cannam@85: data [k] = memory ; cannam@85: } ; cannam@85: data [datalen-1] = data [datalen-2] ; cannam@85: cannam@85: } /* smoothed_diff_double */ cannam@85: cannam@85: static void cannam@85: check_comment (SNDFILE * file, int format, int lineno) cannam@85: { const char *comment ; cannam@85: cannam@85: switch (format & SF_FORMAT_TYPEMASK) cannam@85: { case SF_FORMAT_AIFF : cannam@85: case SF_FORMAT_WAV : cannam@85: case SF_FORMAT_WAVEX : cannam@85: break ; cannam@85: default : cannam@85: return ; cannam@85: } ; cannam@85: cannam@85: comment = sf_get_string (file, SF_STR_COMMENT) ; cannam@85: if (comment == NULL) cannam@85: { printf ("\n\nLine %d : File does not contain a comment string.\n\n", lineno) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: if (strcmp (comment, long_comment) != 0) cannam@85: { printf ("\n\nLine %d : File comment does not match comment written.\n\n", lineno) ; cannam@85: exit (1) ; cannam@85: } ; cannam@85: cannam@85: return ; cannam@85: } /* check_comment */ cannam@85: cannam@85: static int cannam@85: is_lossy (int filetype) cannam@85: { cannam@85: switch (SF_FORMAT_SUBMASK & filetype) cannam@85: { case SF_FORMAT_PCM_U8 : cannam@85: case SF_FORMAT_PCM_S8 : cannam@85: case SF_FORMAT_PCM_16 : cannam@85: case SF_FORMAT_PCM_24 : cannam@85: case SF_FORMAT_PCM_32 : cannam@85: case SF_FORMAT_FLOAT : cannam@85: case SF_FORMAT_DOUBLE : cannam@85: return 0 ; cannam@85: cannam@85: default : cannam@85: break ; cannam@85: } ; cannam@85: cannam@85: return 1 ; cannam@85: } /* is_lossy */ cannam@85: