开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > Audio > 查看源码
floating_point_test.tpl
资源名称：libsndfile-1.0.21.tar.gz　[点击查看]
上传用户：shw771010
上传日期：2022-01-05
资源大小：991k
文件大小：15k
源码类别：

Audio

开发平台：

Unix_Linux

floating_point_test.tpl：源码内容
  1. [+ AutoGen5 template c +]
  2. /*
  3. ** Copyright (C) 1999-2009 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. */
  20. #include "sfconfig.h"
  22. #include <stdio.h>
  23. #include <stdlib.h>
  24. #include <string.h>
  25. #include <math.h>
  27. #if HAVE_UNISTD_H
  28. #include <unistd.h>
  29. #endif
  31. #include <sndfile.h>
  33. #include "dft_cmp.h"
  34. #include "utils.h"
  36. #define SAMPLE_RATE 16000
  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) ;
  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. +]
  46. static double double_data [DFT_DATA_LENGTH] ;
  47. static double double_test [DFT_DATA_LENGTH] ;
  49. static float float_data [DFT_DATA_LENGTH] ;
  50. static float float_test [DFT_DATA_LENGTH] ;
  52. static double double_data [DFT_DATA_LENGTH] ;
  53. static short short_data [DFT_DATA_LENGTH] ;
  54. static int int_data [DFT_DATA_LENGTH] ;
  56. int
  57. main (int argc, char *argv [])
  58. { int allow_exit = 1 ;
  60. if (argc == 2 && ! strstr (argv [1], "no-exit"))
  61. allow_exit = 0 ;
  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
  68. /* Float tests. */
  69. float_scaled_test ("float.raw", allow_exit, SF_FALSE, SF_ENDIAN_LITTLE | SF_FORMAT_RAW | SF_FORMAT_FLOAT, -163.0) ;
  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) ;
  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) ;
  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) ;
  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) ;
  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) ;
  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) ;
  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) ;
  100. float_scaled_test ("adpcm.vox", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_VOX_ADPCM, -40.0) ;
  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) ;
  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) ;
  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) ;
  114. float_scaled_test ("vorbis.oga", allow_exit, SF_FALSE, SF_FORMAT_OGG | SF_FORMAT_VORBIS, -31.0) ;
  115. #endif
  117. float_scaled_test ("replace_float.raw", allow_exit, SF_TRUE, SF_ENDIAN_LITTLE | SF_FORMAT_RAW | SF_FORMAT_FLOAT, -163.0) ;
  119. /*==============================================================================
  120. ** Double tests.
  121. */
  123. double_scaled_test ("double.raw", allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_DOUBLE, -300.0) ;
  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) ;
  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) ;
  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) ;
  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) ;
  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) ;
  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) ;
  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) ;
  152. double_scaled_test ("adpcm.vox" , allow_exit, SF_FALSE, SF_FORMAT_RAW | SF_FORMAT_VOX_ADPCM, -40.0) ;
  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) ;
  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) ;
  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) ;
  166. double_scaled_test ("vorbis.oga", allow_exit, SF_FALSE, SF_FORMAT_OGG | SF_FORMAT_VORBIS, -29.0) ;
  167. #endif
  169. double_scaled_test ("replace_double.raw", allow_exit, SF_TRUE, SF_FORMAT_RAW | SF_FORMAT_DOUBLE, -300.0) ;
  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. +]
  178. return 0 ;
  179. } /* main */
  181. /*============================================================================================
  182.  * Here are the test functions.
  183.  */
  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 ;
  191. print_test_name ("float_scaled_test", filename) ;
  193. gen_windowed_sine_float (float_data, DFT_DATA_LENGTH, 1.0) ;
  195. sfinfo.samplerate = SAMPLE_RATE ;
  196. sfinfo.frames = DFT_DATA_LENGTH ;
  197. sfinfo.channels = 1 ;
  198. sfinfo.format = filetype ;
  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) ;
  203. test_write_float_or_die (file, 0, float_data, DFT_DATA_LENGTH, __LINE__) ;
  205. sf_close (file) ;
  207. memset (float_test, 0, sizeof (float_test)) ;
  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) ;
  212. exit_if_true (sfinfo.format != filetype, "nnLine %d: Returned format incorrect (0x%08X => 0x%08X).n", __LINE__, filetype, sfinfo.format) ;
  213. exit_if_true (sfinfo.frames < DFT_DATA_LENGTH, "nnLine %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, "nnLine %d: Incorrect number of channels in file.n", __LINE__) ;
  216. check_log_buffer_or_die (file, __LINE__) ;
  218. test_read_float_or_die (file, 0, float_test, DFT_DATA_LENGTH, __LINE__) ;
  220. sf_close (file) ;
  222. snr = dft_cmp_float (__LINE__, float_data, float_test, DFT_DATA_LENGTH, target_snr, allow_exit) ;
  224. exit_if_true (snr > target_snr, "% 6.1fdB SNRnn    Error : should be better than % 6.1fdBnn", snr, target_snr) ;
  226. printf ("% 6.1fdB SNR ... okn", snr) ;
  228. unlink (filename) ;
  230. return ;
  231. } /* float_scaled_test */
  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 ;
  239. print_test_name ("double_scaled_test", filename) ;
  241. gen_windowed_sine_double (double_data, DFT_DATA_LENGTH, 0.95) ;
  243. sfinfo.samplerate = SAMPLE_RATE ;
  244. sfinfo.frames = DFT_DATA_LENGTH ;
  245. sfinfo.channels = 1 ;
  246. sfinfo.format = filetype ;
  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) ;
  251. test_write_double_or_die (file, 0, double_data, DFT_DATA_LENGTH, __LINE__) ;
  253. sf_close (file) ;
  255. memset (double_test, 0, sizeof (double_test)) ;
  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) ;
  260. exit_if_true (sfinfo.format != filetype, "nnLine %d: Returned format incorrect (0x%08X => 0x%08X).n", __LINE__, filetype, sfinfo.format) ;
  261. exit_if_true (sfinfo.frames < DFT_DATA_LENGTH, "nnLine %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, "nnLine %d: Incorrect number of channels in file.n", __LINE__) ;
  264. check_log_buffer_or_die (file, __LINE__) ;
  266. test_read_double_or_die (file, 0, double_test, DFT_DATA_LENGTH, __LINE__) ;
  268. sf_close (file) ;
  270. snr = dft_cmp_double (__LINE__, double_data, double_test, DFT_DATA_LENGTH, target_snr, allow_exit) ;
  272. exit_if_true (snr > target_snr, "% 6.1fdB SNRnn    Error : should be better than % 6.1fdBnn", snr, target_snr) ;
  274. printf ("% 6.1fdB SNR ... okn", snr) ;
  276. unlink (filename) ;
  278. return ;
  279. } /* double_scaled_test */
  281. /*==============================================================================
  282. */
  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 ;
  292. print_test_name ("[+ (get "float_name") +]_[+ (get "int_name") +]_[+ (get "end_name") +]_test", filename) ;
  294. gen_windowed_sine_[+ (get "float_name") +] ([+ (get "float_name") +]_data, ARRAY_LEN ([+ (get "float_name") +]_data), 0.98) ;
  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") +] ;
  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) ;
  305. file = test_open_file_or_die (filename, SFM_READ, &sfinfo, SF_TRUE, __LINE__) ;
  307. if (sfinfo.frames != ARRAY_LEN ([+ (get "float_name") +]_data))
  308. { printf ("nnLine %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. } ;
  312. if (sfinfo.channels != 1)
  313. { printf ("nnLine %d: Incorrect number of channels in file.n", __LINE__) ;
  314. exit (1) ;
  315. } ;
  317. sf_command (file, SFC_SET_SCALE_FLOAT_INT_READ, NULL, SF_TRUE) ;
  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) ;
  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]) ;
  327. if (1.0 * abs (max - [+ (get "int_max") +]) / [+ (get "int_max") +] > 0.01)
  328. { printf ("nnLine %d: Bad maximum (%d should be %d).nn", __LINE__, max, [+ (get "int_max") +]) ;
  329. exit (1) ;
  330. } ;
  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 +]