interpretor/libsndfile-1.0.25/tests/floating_point_test.tpl

   1 [+ AutoGen5 template c +]
   2 /*
   3 ** Copyright (C) 1999-2011 Erik de Castro Lopo <erikd@mega-nerd.com>
   4 **
   5 ** This program is free software; you can redistribute it and/or modify
   6 ** it under the terms of the GNU General Public License as published by
   7 ** the Free Software Foundation; either version 2 of the License, or
   8 ** (at your option) any later version.
   9 **
  10 ** This program is distributed in the hope that it will be useful,
  11 ** but WITHOUT ANY WARRANTY; without even the implied warranty of
  12 ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  13 ** GNU General Public License for more details.
  14 **
  15 ** You should have received a copy of the GNU General Public License
  16 ** along with this program; if not, write to the Free Software
  17 ** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
  18 */
  19
  20 #include "sfconfig.h"
  21
  22 #include <stdio.h>
  23 #include <stdlib.h>
  24 #include <string.h>
  25 #include <math.h>
  26
  27 #if HAVE_UNISTD_H
  28 #include <unistd.h>
  29 #endif
  30
  31 #include <sndfile.h>
  32
  33 #include "dft_cmp.h"
  34 #include "utils.h"
  35
  36 #define SAMPLE_RATE                     16000
  37
  38 static void     float_scaled_test       (const char *filename, int allow_exit, int replace_float, int filetype, double target_snr) ;
  39 static void     double_scaled_test      (const char *filename, int allow_exit, int replace_float, int filetype, double target_snr) ;
  40
  41 [+ FOR float_type +][+ FOR int_type +][+ FOR endian_type
  42 +]static void [+ (get "float_name") +]_[+ (get "int_name") +]_[+ (get "end_name") +]_test (const char * filename) ;
  43 [+ ENDFOR endian_type +][+ ENDFOR int_type +][+ ENDFOR float_type
  44 +]
  45
  46 static  double  double_data [DFT_DATA_LENGTH] ;
  47 static  double  double_test [DFT_DATA_LENGTH] ;
  48
  49 static float    float_data [DFT_DATA_LENGTH] ;
  50 static float    float_test [DFT_DATA_LENGTH] ;
  51
  52 static double   double_data [DFT_DATA_LENGTH] ;
  53 static short    short_data [DFT_DATA_LENGTH] ;
  54 static int              int_data [DFT_DATA_LENGTH] ;
  55
  56 int
  57 main (int argc, char *argv [])
  58 {       int allow_exit = 1 ;
  59
  60         if (argc == 2 && ! strstr (argv [1], "no-exit"))
  61                 allow_exit = 0 ;
  62
  63 #if ((HAVE_LRINTF == 0) && (HAVE_LRINT_REPLACEMENT == 0))
  64         puts ("*** Cannot run this test on this platform because it lacks lrintf().") ;
  65         exit (0) ;
  66 #endif
  67
  68         /* Float tests. */
  69         float_scaled_test       ("float.raw", allow_exit, SF_FALSE, SF_ENDIAN_LITTLE | SF_FORMAT_RAW | SF_FORMAT_FLOAT, -163.0) ;
  70
  71         /* Test both signed and unsigned 8 bit files. */
  72         float_scaled_test       ("pcm_s8.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_PCM_S8, -39.0) ;
  73         float_scaled_test       ("pcm_u8.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_PCM_U8, -39.0) ;
  74
  75         float_scaled_test       ("pcm_16.raw", allow_exit, SF_FALSE, SF_ENDIAN_BIG | SF_FORMAT_RAW | SF_FORMAT_PCM_16, -87.0) ;
  76         float_scaled_test       ("pcm_24.raw", allow_exit, SF_FALSE, SF_ENDIAN_LITTLE | SF_FORMAT_RAW | SF_FORMAT_PCM_24, -138.0) ;
  77         float_scaled_test       ("pcm_32.raw", allow_exit, SF_FALSE, SF_ENDIAN_BIG | SF_FORMAT_RAW | SF_FORMAT_PCM_32, -163.0) ;
  78
  79         float_scaled_test       ("ulaw.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_ULAW, -50.0) ;
  80         float_scaled_test       ("alaw.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_ALAW, -49.0) ;
  81
  82         float_scaled_test       ("ima_adpcm.wav", allow_exit, SF_FALSE, SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, -47.0) ;
  83         float_scaled_test       ("ms_adpcm.wav" , allow_exit, SF_FALSE, SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, -40.0) ;
  84         float_scaled_test       ("gsm610.raw"   , allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_GSM610, -33.0) ;
  85
  86         float_scaled_test       ("g721_32.au", allow_exit, SF_FALSE, SF_FORMAT_AU | SF_FORMAT_G721_32, -34.0) ;
  87         float_scaled_test       ("g723_24.au", allow_exit, SF_FALSE, SF_FORMAT_AU | SF_FORMAT_G723_24, -34.0) ;
  88         float_scaled_test       ("g723_40.au", allow_exit, SF_FALSE, SF_FORMAT_AU | SF_FORMAT_G723_40, -40.0) ;
  89
  90         /*      PAF files do not use the same encoding method for 24 bit PCM data as other file
  91         **      formats so we need to explicitly test it here.
  92         */
  93         float_scaled_test       ("le_paf_24.paf", allow_exit, SF_FALSE, SF_ENDIAN_LITTLE | SF_FORMAT_PAF | SF_FORMAT_PCM_24, -149.0) ;
  94         float_scaled_test       ("be_paf_24.paf", allow_exit, SF_FALSE, SF_ENDIAN_BIG | SF_FORMAT_PAF | SF_FORMAT_PCM_24, -149.0) ;
  95
  96         float_scaled_test       ("dwvw_12.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_DWVW_12, -64.0) ;
  97         float_scaled_test       ("dwvw_16.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_DWVW_16, -92.0) ;
  98         float_scaled_test       ("dwvw_24.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_DWVW_24, -151.0) ;
  99
 100         float_scaled_test       ("adpcm.vox", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_VOX_ADPCM, -40.0) ;
 101
 102         float_scaled_test       ("dpcm_16.xi", allow_exit, SF_FALSE, SF_FORMAT_XI | SF_FORMAT_DPCM_16, -90.0) ;
 103         float_scaled_test       ("dpcm_8.xi" , allow_exit, SF_FALSE, SF_FORMAT_XI | SF_FORMAT_DPCM_8 , -41.0) ;
 104
 105         float_scaled_test       ("pcm_s8.sds", allow_exit, SF_FALSE, SF_FORMAT_SDS | SF_FORMAT_PCM_S8, -90.0) ;
 106         float_scaled_test       ("pcm_16.sds", allow_exit, SF_FALSE, SF_FORMAT_SDS | SF_FORMAT_PCM_16, -140.0) ;
 107         float_scaled_test       ("pcm_24.sds", allow_exit, SF_FALSE, SF_FORMAT_SDS | SF_FORMAT_PCM_24, -170.0) ;
 108
 109 #if HAVE_EXTERNAL_LIBS
 110         float_scaled_test       ("flac_8.flac", allow_exit, SF_FALSE, SF_FORMAT_FLAC | SF_FORMAT_PCM_S8, -39.0) ;
 111         float_scaled_test       ("flac_16.flac", allow_exit, SF_FALSE, SF_FORMAT_FLAC | SF_FORMAT_PCM_16, -87.0) ;
 112         float_scaled_test       ("flac_24.flac", allow_exit, SF_FALSE, SF_FORMAT_FLAC | SF_FORMAT_PCM_24, -138.0) ;
 113
 114         float_scaled_test       ("vorbis.oga", allow_exit, SF_FALSE, SF_FORMAT_OGG | SF_FORMAT_VORBIS, -31.0) ;
 115 #endif
 116
 117         float_scaled_test       ("replace_float.raw", allow_exit, SF_TRUE, SF_ENDIAN_LITTLE | SF_FORMAT_RAW | SF_FORMAT_FLOAT, -163.0) ;
 118
 119         /*==============================================================================
 120         ** Double tests.
 121         */
 122
 123         double_scaled_test      ("double.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_DOUBLE, -300.0) ;
 124
 125         /* Test both signed (AIFF) and unsigned (WAV) 8 bit files. */
 126         double_scaled_test      ("pcm_s8.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_PCM_S8, -39.0) ;
 127         double_scaled_test      ("pcm_u8.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_PCM_U8, -39.0) ;
 128
 129         double_scaled_test      ("pcm_16.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_PCM_16, -87.0) ;
 130         double_scaled_test      ("pcm_24.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_PCM_24, -135.0) ;
 131         double_scaled_test      ("pcm_32.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_PCM_32, -184.0) ;
 132
 133         double_scaled_test      ("ulaw.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_ULAW, -50.0) ;
 134         double_scaled_test      ("alaw.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_ALAW, -49.0) ;
 135
 136         double_scaled_test      ("ima_adpcm.wav", allow_exit, SF_FALSE, SF_FORMAT_WAV | SF_FORMAT_IMA_ADPCM, -47.0) ;
 137         double_scaled_test      ("ms_adpcm.wav" , allow_exit, SF_FALSE, SF_FORMAT_WAV | SF_FORMAT_MS_ADPCM, -40.0) ;
 138         double_scaled_test      ("gsm610.raw"   , allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_GSM610, -33.0) ;
 139
 140         double_scaled_test      ("g721_32.au", allow_exit, SF_FALSE, SF_FORMAT_AU | SF_FORMAT_G721_32, -34.0) ;
 141         double_scaled_test      ("g723_24.au", allow_exit, SF_FALSE, SF_FORMAT_AU | SF_FORMAT_G723_24, -34.0) ;
 142         double_scaled_test      ("g723_40.au", allow_exit, SF_FALSE, SF_FORMAT_AU | SF_FORMAT_G723_40, -40.0) ;
 143
 144         /*      24 bit PCM PAF files tested here. */
 145         double_scaled_test      ("be_paf_24.paf", allow_exit, SF_FALSE, SF_ENDIAN_BIG | SF_FORMAT_PAF | SF_FORMAT_PCM_24, -151.0) ;
 146         double_scaled_test      ("le_paf_24.paf", allow_exit, SF_FALSE, SF_ENDIAN_LITTLE | SF_FORMAT_PAF | SF_FORMAT_PCM_24, -151.0) ;
 147
 148         double_scaled_test      ("dwvw_12.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_DWVW_12, -64.0) ;
 149         double_scaled_test      ("dwvw_16.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_DWVW_16, -92.0) ;
 150         double_scaled_test      ("dwvw_24.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_DWVW_24, -151.0) ;
 151
 152         double_scaled_test      ("adpcm.vox" , allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_VOX_ADPCM, -40.0) ;
 153
 154         double_scaled_test      ("dpcm_16.xi", allow_exit, SF_FALSE, SF_FORMAT_XI | SF_FORMAT_DPCM_16, -90.0) ;
 155         double_scaled_test      ("dpcm_8.xi" , allow_exit, SF_FALSE, SF_FORMAT_XI | SF_FORMAT_DPCM_8 , -42.0) ;
 156
 157         double_scaled_test      ("pcm_s8.sds", allow_exit, SF_FALSE, SF_FORMAT_SDS | SF_FORMAT_PCM_S8, -90.0) ;
 158         double_scaled_test      ("pcm_16.sds", allow_exit, SF_FALSE, SF_FORMAT_SDS | SF_FORMAT_PCM_16, -140.0) ;
 159         double_scaled_test      ("pcm_24.sds", allow_exit, SF_FALSE, SF_FORMAT_SDS | SF_FORMAT_PCM_24, -180.0) ;
 160
 161 #if HAVE_EXTERNAL_LIBS
 162         double_scaled_test      ("flac_8.flac", allow_exit, SF_FALSE, SF_FORMAT_FLAC | SF_FORMAT_PCM_S8, -39.0) ;
 163         double_scaled_test      ("flac_16.flac", allow_exit, SF_FALSE, SF_FORMAT_FLAC | SF_FORMAT_PCM_16, -87.0) ;
 164         double_scaled_test      ("flac_24.flac", allow_exit, SF_FALSE, SF_FORMAT_FLAC | SF_FORMAT_PCM_24, -138.0) ;
 165
 166         double_scaled_test      ("vorbis.oga", allow_exit, SF_FALSE, SF_FORMAT_OGG | SF_FORMAT_VORBIS, -29.0) ;
 167 #endif
 168
 169         double_scaled_test      ("replace_double.raw", allow_exit, SF_TRUE, SF_FORMAT_RAW | SF_FORMAT_DOUBLE, -300.0) ;
 170
 171         putchar ('\n') ;
 172         /* Float int tests. */
 173 [+ FOR float_type +][+ FOR int_type +][+ FOR endian_type
 174 +]      [+ (get "float_name") +]_[+ (get "int_name") +]_[+ (get "end_name") +]_test ("[+ (get "float_name") +]_[+ (get "int_name") +]_[+ (get "end_name") +].au") ;
 175 [+ ENDFOR endian_type +][+ ENDFOR int_type +][+ ENDFOR float_type
 176 +]
 177
 178         return 0 ;
 179 } /* main */
 180
 181 /*============================================================================================
 182  *      Here are the test functions.
 183  */
 184
 185 static void
 186 float_scaled_test (const char *filename, int allow_exit, int replace_float, int filetype, double target_snr)
 187 {       SNDFILE         *file ;
 188         SF_INFO         sfinfo ;
 189         double          snr ;
 190
 191         print_test_name ("float_scaled_test", filename) ;
 192
 193         gen_windowed_sine_float (float_data, DFT_DATA_LENGTH, 1.0) ;
 194
 195         sfinfo.samplerate       = SAMPLE_RATE ;
 196         sfinfo.frames           = DFT_DATA_LENGTH ;
 197         sfinfo.channels         = 1 ;
 198         sfinfo.format           = filetype ;
 199
 200         file = test_open_file_or_die (filename, SFM_WRITE, &sfinfo, SF_TRUE, __LINE__) ;
 201         sf_command (file, SFC_TEST_IEEE_FLOAT_REPLACE, NULL, replace_float) ;
 202
 203         test_write_float_or_die (file, 0, float_data, DFT_DATA_LENGTH, __LINE__) ;
 204
 205         sf_close (file) ;
 206
 207         memset (float_test, 0, sizeof (float_test)) ;
 208
 209         file = test_open_file_or_die (filename, SFM_READ, &sfinfo, SF_TRUE, __LINE__) ;
 210         sf_command (file, SFC_TEST_IEEE_FLOAT_REPLACE, NULL, replace_float) ;
 211
 212         exit_if_true (sfinfo.format != filetype, "\n\nLine %d: Returned format incorrect (0x%08X => 0x%08X).\n", __LINE__, filetype, sfinfo.format) ;
 213         exit_if_true (sfinfo.frames < DFT_DATA_LENGTH, "\n\nLine %d: Incorrect number of frames in file (too short). (%ld should be %d)\n", __LINE__, SF_COUNT_TO_LONG (sfinfo.frames), DFT_DATA_LENGTH) ;
 214         exit_if_true (sfinfo.channels != 1, "\n\nLine %d: Incorrect number of channels in file.\n", __LINE__) ;
 215
 216         check_log_buffer_or_die (file, __LINE__) ;
 217
 218         test_read_float_or_die (file, 0, float_test, DFT_DATA_LENGTH, __LINE__) ;
 219
 220         sf_close (file) ;
 221
 222         snr = dft_cmp_float (__LINE__, float_data, float_test, DFT_DATA_LENGTH, target_snr, allow_exit) ;
 223
 224         exit_if_true (snr > target_snr, "% 6.1fdB SNR\n\n    Error : should be better than % 6.1fdB\n\n", snr, target_snr) ;
 225
 226         printf ("% 6.1fdB SNR ... ok\n", snr) ;
 227
 228         unlink (filename) ;
 229
 230         return ;
 231 } /* float_scaled_test */
 232
 233 static void
 234 double_scaled_test (const char *filename, int allow_exit, int replace_float, int filetype, double target_snr)
 235 {       SNDFILE         *file ;
 236         SF_INFO         sfinfo ;
 237         double          snr ;
 238
 239         print_test_name ("double_scaled_test", filename) ;
 240
 241         gen_windowed_sine_double (double_data, DFT_DATA_LENGTH, 0.95) ;
 242
 243         sfinfo.samplerate       = SAMPLE_RATE ;
 244         sfinfo.frames           = DFT_DATA_LENGTH ;
 245         sfinfo.channels         = 1 ;
 246         sfinfo.format           = filetype ;
 247
 248         file = test_open_file_or_die (filename, SFM_WRITE, &sfinfo, SF_TRUE, __LINE__) ;
 249         sf_command (file, SFC_TEST_IEEE_FLOAT_REPLACE, NULL, replace_float) ;
 250
 251         test_write_double_or_die (file, 0, double_data, DFT_DATA_LENGTH, __LINE__) ;
 252
 253         sf_close (file) ;
 254
 255         memset (double_test, 0, sizeof (double_test)) ;
 256
 257         file = test_open_file_or_die (filename, SFM_READ, &sfinfo, SF_TRUE, __LINE__) ;
 258         sf_command (file, SFC_TEST_IEEE_FLOAT_REPLACE, NULL, replace_float) ;
 259
 260         exit_if_true (sfinfo.format != filetype, "\n\nLine %d: Returned format incorrect (0x%08X => 0x%08X).\n", __LINE__, filetype, sfinfo.format) ;
 261         exit_if_true (sfinfo.frames < DFT_DATA_LENGTH, "\n\nLine %d: Incorrect number of frames in file (too short). (%ld should be %d)\n", __LINE__, SF_COUNT_TO_LONG (sfinfo.frames), DFT_DATA_LENGTH) ;
 262         exit_if_true (sfinfo.channels != 1, "\n\nLine %d: Incorrect number of channels in file.\n", __LINE__) ;
 263
 264         check_log_buffer_or_die (file, __LINE__) ;
 265
 266         test_read_double_or_die (file, 0, double_test, DFT_DATA_LENGTH, __LINE__) ;
 267
 268         sf_close (file) ;
 269
 270         snr = dft_cmp_double (__LINE__, double_data, double_test, DFT_DATA_LENGTH, target_snr, allow_exit) ;
 271
 272         exit_if_true (snr > target_snr, "% 6.1fdB SNR\n\n    Error : should be better than % 6.1fdB\n\n", snr, target_snr) ;
 273
 274         printf ("% 6.1fdB SNR ... ok\n", snr) ;
 275
 276         unlink (filename) ;
 277
 278         return ;
 279 } /* double_scaled_test */
 280
 281 /*==============================================================================
 282 */
 283
 284 [+ FOR float_type +][+ FOR int_type +][+ FOR endian_type
 285 +]
 286 static void
 287 [+ (get "float_name") +]_[+ (get "int_name") +]_[+ (get "end_name") +]_test (const char * filename)
 288 {       SNDFILE         *file ;
 289         SF_INFO         sfinfo ;
 290         unsigned        k, max ;
 291
 292         print_test_name ("[+ (get "float_name") +]_[+ (get "int_name") +]_[+ (get "end_name") +]_test", filename) ;
 293
 294         gen_windowed_sine_[+ (get "float_name") +] ([+ (get "float_name") +]_data, ARRAY_LEN ([+ (get "float_name") +]_data), 0.98) ;
 295
 296         sfinfo.samplerate       = SAMPLE_RATE ;
 297         sfinfo.frames           = ARRAY_LEN ([+ (get "int_name") +]_data) ;
 298         sfinfo.channels         = 1 ;
 299         sfinfo.format           = [+ (get "end_type") +] | SF_FORMAT_AU | [+ (get "minor_type") +] ;
 300
 301         file = test_open_file_or_die (filename, SFM_WRITE, &sfinfo, SF_TRUE, __LINE__) ;
 302         test_write_[+ (get "float_name") +]_or_die (file, 0, [+ (get "float_name") +]_data, ARRAY_LEN ([+ (get "float_name") +]_data), __LINE__) ;
 303         sf_close (file) ;
 304
 305         file = test_open_file_or_die (filename, SFM_READ, &sfinfo, SF_TRUE, __LINE__) ;
 306
 307         if (sfinfo.frames != ARRAY_LEN ([+ (get "float_name") +]_data))
 308         {       printf ("\n\nLine %d: Incorrect number of frames in file (too short). (%ld should be %d)\n", __LINE__, SF_COUNT_TO_LONG (sfinfo.frames), DFT_DATA_LENGTH) ;
 309                 exit (1) ;
 310                 } ;
 311
 312         if (sfinfo.channels != 1)
 313         {       printf ("\n\nLine %d: Incorrect number of channels in file.\n", __LINE__) ;
 314                 exit (1) ;
 315                 } ;
 316
 317         sf_command (file, SFC_SET_SCALE_FLOAT_INT_READ, NULL, SF_TRUE) ;
 318
 319         test_read_[+ (get "int_name") +]_or_die (file, 0, [+ (get "int_name") +]_data, ARRAY_LEN ([+ (get "int_name") +]_data), __LINE__) ;
 320         sf_close (file) ;
 321
 322         max = 0 ;
 323         for (k = 0 ; k < ARRAY_LEN ([+ (get "int_name") +]_data) ; k++)
 324                 if ((unsigned) abs ([+ (get "int_name") +]_data [k]) > max)
 325                         max = abs ([+ (get "int_name") +]_data [k]) ;
 326
 327         if (1.0 * abs (max - [+ (get "int_max") +]) / [+ (get "int_max") +] > 0.01)
 328         {       printf ("\n\nLine %d: Bad maximum (%d should be %d).\n\n", __LINE__, max, [+ (get "int_max") +]) ;
 329                 exit (1) ;
 330                 } ;
 331
 332         unlink (filename) ;
 333         puts ("ok") ;
 334 } /* [+ (get "float_name") +]_[+ (get "int_name") +]_[+ (get "end_name") +]_test */
 335 [+ ENDFOR endian_type +][+ ENDFOR int_type +][+ ENDFOR float_type +]
 336