ref: 7aa0ed40c1b0ae4edccb6fb6bafafdc9fcdef073
parent: c4b146c4ba38de86f6970a90463ea35c320a78b6
author: Gregory Maxwell <[email protected]>
date: Thu Feb 19 16:28:33 EST 2009
Make ectest use random random numbers. Add a new test "tandem-test" which tests the encoder and decoder in many modes.
--- a/tests/Makefile.am
+++ b/tests/Makefile.am
@@ -1,9 +1,9 @@
INCLUDES = -I$(top_srcdir)/libcelt
METASOURCES = AUTO
-TESTS = type-test ectest cwrs32-test real-fft-test dft-test laplace-test mdct-test mathops-test
+TESTS = type-test ectest cwrs32-test real-fft-test dft-test laplace-test mdct-test mathops-test tandem-test
-noinst_PROGRAMS = type-test ectest cwrs32-test real-fft-test dft-test laplace-test mdct-test mathops-test
+noinst_PROGRAMS = type-test ectest cwrs32-test real-fft-test dft-test laplace-test mdct-test mathops-test tandem-test
type_test_SOURCES = type-test.c
ectest_SOURCES = ectest.c
@@ -14,3 +14,5 @@
mdct_test_SOURCES = mdct-test.c
#rotation_test_SOURCES = rotation-test.c
mathops_test_SOURCES = mathops-test.c
+tandem_test_SOURCES = tandem-test.c
+tandem_test_LDADD = $(top_builddir)/libcelt/libcelt.la
--- a/tests/ectest.c
+++ b/tests/ectest.c
@@ -5,6 +5,9 @@
#include <stdlib.h>
#include <stdio.h>
#include <math.h>
+#include <sys/types.h>
+#include <unistd.h>
+#include <time.h>
#include "entcode.h"
#include "entenc.h"
#include "entdec.h"
@@ -32,8 +35,18 @@
int sz;
int i;
int ret;
+ unsigned int seed;
ret=0;
entropy=0;
+
+ if (_argc > 2) {+ fprintf(stderr, "Usage: %s [<seed>]\n", _argv[0]);
+ return 1;
+ }
+ if (_argc > 1)
+ seed = atoi(_argv[1]);
+ else
+ seed = (time(NULL) ^ (getpid()%(1<<16) << 16));
/*Testing encoding of raw bit values.*/
ec_byte_writeinit(&buf);
ec_enc_init(&enc,&buf);
@@ -91,8 +104,8 @@
ret=-1;
}
ec_byte_writeclear(&buf);
- fprintf(stderr,"Testing random streams...\n");
- srand(0);
+ srand(seed);
+ fprintf(stderr,"Testing random streams... Random seed: %u (%.4X)\n", seed, rand() % 65536);
for(i=0;i<409600;i++){unsigned *data;
int j;
@@ -118,8 +131,8 @@
ec_enc_done(&enc);
if ((tell_bits+7)/8 < ec_byte_bytes(&buf))
{- fprintf (stderr, "tell() lied, there's %li bytes instead of %d\n",
- ec_byte_bytes(&buf), (tell_bits+7)/8);
+ fprintf (stderr, "tell() lied, there's %li bytes instead of %d (Random seed: %u)\n",
+ ec_byte_bytes(&buf), (tell_bits+7)/8,seed);
ret=-1;
}
tell_bits -= 8*ec_byte_bytes(&buf);
@@ -129,8 +142,8 @@
sym=ec_dec_uint(&dec,ft);
if(sym!=data[j]){fprintf(stderr,
- "Decoded %i instead of %i with ft of %i at position %i of %i.\n",
- sym,data[j],ft,j,sz);
+ "Decoded %i instead of %i with ft of %i at position %i of %i (Random seed: %u).\n",
+ sym,data[j],ft,j,sz,seed);
ret=-1;
}
}
--- /dev/null
+++ b/tests/tandem-test.c
@@ -1,0 +1,191 @@
+/* (C) 2009 Gregory Maxwell
+
+ This test runs pink noise through the encoder and decoder many times
+ while feeding the output back into the input. It checks that after
+ a number of cycles the energy has not increased or decreased by too
+ large an amount.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions
+ are met:
+
+ - Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ - Redistributions in binary form must reproduce the above copyright
+ notice, this list of conditions and the following disclaimer in the
+ documentation and/or other materials provided with the distribution.
+
+ - Neither the name of the Xiph.org Foundation nor the names of its
+ contributors may be used to endorse or promote products derived from
+ this software without specific prior written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR
+ CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
+ EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
+ PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
+ PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
+ LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
+ NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
+ SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+
+#ifdef HAVE_CONFIG_H
+#include "config.h"
+#endif
+
+#include "celt.h"
+#include "arch.h"
+#include <stdio.h>
+#include <stdlib.h>
+#include <sys/types.h>
+#include <time.h>
+#include <unistd.h>
+#include <math.h>
+#include <string.h>
+
+
+int async_tandem(int rate, int frame_size, int channels, int bitrate_min,
+ int bitrate_max)
+{+ unsigned char data[250];
+ CELTMode *mode = NULL;
+ CELTEncoder *enc;
+ short carry[2];
+ short pcm[512 * 2];
+ CELTDecoder *dec;
+ int bmin, bmax;
+ float ms;
+ int ret, i, j, bytes_per_frame;
+ int increment = 1;
+
+ bmin = floor((bitrate_min / (rate / (float) frame_size)) / 8.0);
+ bmax = ceil((bitrate_max / (rate / (float) frame_size)) / 8.0);
+ if (bmin < 12)
+ bmin = 12;
+ if (bmax > 250)
+ bmax = 250;
+ if (bmin >= bmax)
+ bmax = bmin + 8;
+
+ /*increment += (bmax - bmin) / 64; */
+
+ printf ("Testing asynchronous tandeming (%dHz, %dch, %d samples, %d - %d bytes).\n",+ rate, channels, frame_size, bmin, bmax);
+
+ mode = celt_mode_create(rate, channels, frame_size, NULL);
+ if (mode == NULL) {+ fprintf(stderr, "Error: failed to create a mode\n");
+ exit(1);
+ }
+
+ dec = celt_decoder_create(mode);
+ enc = celt_encoder_create(mode);
+
+ for (j = 0; j < frame_size * channels; j++)
+ pcm[j] = 0;
+
+ for (bytes_per_frame = bmin; bytes_per_frame <= bmax;
+ bytes_per_frame += increment) {+
+ /*Prime the encoder and decoder */
+ for (i = 0; i < (1024 + (frame_size >> 1)) / frame_size + 2; i++) {+
+ for (j = 0; j < channels; j++)
+ pcm[j] = pcm[frame_size * channels - (channels - j + 1)];
+ for (j = channels; j < frame_size * channels - 1; j++)
+ pcm[j] = ((rand() % 4096) - 2048) + .9 * pcm[j - channels];
+
+ ret = celt_encode(enc, pcm, NULL, data, bytes_per_frame);
+ if (ret != bytes_per_frame) {+ fprintf(stderr, "Error: during init celt_encode returned %d\n", ret);
+ exit(1);
+ }
+
+ ret = celt_decode(dec, data, ret, pcm);
+ if (ret != CELT_OK) {+ fprintf(stderr, "Error: during init celt_decode returned %d\n", ret);
+ }
+ }
+
+ for (j = 0; j < channels; j++)
+ pcm[j] = pcm[frame_size * channels - (channels - j)];
+ for (j = channels; j < frame_size * channels - 1; j++)
+ pcm[j] = ((rand() % 4096) - 2048) + .9 * pcm[j - channels];
+
+ for (i = 0; i < 8; i++) {+ for (j = 0; j < channels; j++)
+ carry[j] = pcm[frame_size * channels - (channels - j)];
+ memmove(pcm + channels, pcm, sizeof(short) * frame_size * channels);
+ for (j = 0; j < channels; j++)
+ pcm[j] = carry[j];
+
+ ret = celt_encode(enc, pcm, NULL, data, bytes_per_frame);
+ if (ret != bytes_per_frame) {+ fprintf(stderr, "Error: at %d bytes_per_frame celt_encode returned %d\n",
+ bytes_per_frame, ret);
+ exit(1);
+ }
+
+ ret = celt_decode(dec, data, ret, pcm);
+ if (ret != CELT_OK) {+ fprintf(stderr, "Error: at %d bytes_per_frame celt_decode returned %d\n",
+ bytes_per_frame, ret);
+ exit(1);
+ }
+ }
+ ms = 0;
+ for (j = 0; j < frame_size * channels; j++)
+ ms += pcm[j] * pcm[j];
+ ms = sqrt(ms / (frame_size * channels));
+ if (ms > 7000 || ms < 1000) {+ fprintf(stderr, "Error: Signal out of expected range. %d %d %d %d %f\n",
+ rate, channels, frame_size, bytes_per_frame, ms);
+ exit(1);
+ }
+ }
+
+ celt_encoder_destroy(enc);
+ celt_decoder_destroy(dec);
+ celt_mode_destroy(mode);
+
+ return 0;
+}
+
+int main(int argc, char *argv[])
+{+ unsigned int seed;
+ int ch, n;
+
+ if (argc > 2) {+ fprintf(stderr, "Usage: %s [<seed>]\n", argv[0]);
+ return 1;
+ }
+
+ if (argc > 1)
+ seed = atoi(argv[1]);
+ else
+ seed = (time(NULL) ^ ((getpid() % 65536) << 16));
+
+ srand(seed);
+ printf("CELT codec tests. Random seed: %u (%.4X)\n", seed, rand() % 65536);+
+ for (n = 6; n < 10; n++) {+ for (ch = 1; ch <= 2; ch++) {+ async_tandem(44100, 1 << n, ch, 47000 * ch, 77000 * ch);
+ async_tandem(48000, 1 << n, ch, 47000 * ch, 77000 * ch);
+ async_tandem(32000, 1 << n, ch, 31000 * ch, 65000 * ch);
+ }
+ }
+
+ for (ch = 1; ch <= 2; ch++)
+ async_tandem(32000, 320, ch, 31000 * ch, 65000 * ch);
+
+ for (ch = 1; ch <= 2; ch++)
+ async_tandem(48000, 480, ch, 31000 * ch, 77000 * ch);
+
+ return 0;
+}