1*a58d3d2aSXin Li /* Copyright (c) 2011 Xiph.Org Foundation
2*a58d3d2aSXin Li Written by Gregory Maxwell */
3*a58d3d2aSXin Li /*
4*a58d3d2aSXin Li Redistribution and use in source and binary forms, with or without
5*a58d3d2aSXin Li modification, are permitted provided that the following conditions
6*a58d3d2aSXin Li are met:
7*a58d3d2aSXin Li
8*a58d3d2aSXin Li - Redistributions of source code must retain the above copyright
9*a58d3d2aSXin Li notice, this list of conditions and the following disclaimer.
10*a58d3d2aSXin Li
11*a58d3d2aSXin Li - Redistributions in binary form must reproduce the above copyright
12*a58d3d2aSXin Li notice, this list of conditions and the following disclaimer in the
13*a58d3d2aSXin Li documentation and/or other materials provided with the distribution.
14*a58d3d2aSXin Li
15*a58d3d2aSXin Li THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
16*a58d3d2aSXin Li ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
17*a58d3d2aSXin Li LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
18*a58d3d2aSXin Li A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER
19*a58d3d2aSXin Li OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
20*a58d3d2aSXin Li EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
21*a58d3d2aSXin Li PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
22*a58d3d2aSXin Li PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
23*a58d3d2aSXin Li LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
24*a58d3d2aSXin Li NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
25*a58d3d2aSXin Li SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
26*a58d3d2aSXin Li */
27*a58d3d2aSXin Li
deb2_impl(unsigned char * _t,unsigned char ** _p,int _k,int _x,int _y)28*a58d3d2aSXin Li static OPUS_INLINE void deb2_impl(unsigned char *_t,unsigned char **_p,int _k,int _x,int _y)
29*a58d3d2aSXin Li {
30*a58d3d2aSXin Li int i;
31*a58d3d2aSXin Li if(_x>2){
32*a58d3d2aSXin Li if(_y<3)for(i=0;i<_y;i++)*(--*_p)=_t[i+1];
33*a58d3d2aSXin Li }else{
34*a58d3d2aSXin Li _t[_x]=_t[_x-_y];
35*a58d3d2aSXin Li deb2_impl(_t,_p,_k,_x+1,_y);
36*a58d3d2aSXin Li for(i=_t[_x-_y]+1;i<_k;i++){
37*a58d3d2aSXin Li _t[_x]=i;
38*a58d3d2aSXin Li deb2_impl(_t,_p,_k,_x+1,_x);
39*a58d3d2aSXin Li }
40*a58d3d2aSXin Li }
41*a58d3d2aSXin Li }
42*a58d3d2aSXin Li
43*a58d3d2aSXin Li /*Generates a De Bruijn sequence (k,2) with length k^2*/
debruijn2(int _k,unsigned char * _res)44*a58d3d2aSXin Li static OPUS_INLINE void debruijn2(int _k, unsigned char *_res)
45*a58d3d2aSXin Li {
46*a58d3d2aSXin Li unsigned char *p;
47*a58d3d2aSXin Li unsigned char *t;
48*a58d3d2aSXin Li t=malloc(sizeof(unsigned char)*_k*2);
49*a58d3d2aSXin Li memset(t,0,sizeof(unsigned char)*_k*2);
50*a58d3d2aSXin Li p=&_res[_k*_k];
51*a58d3d2aSXin Li deb2_impl(t,&p,_k,1,1);
52*a58d3d2aSXin Li free(t);
53*a58d3d2aSXin Li }
54*a58d3d2aSXin Li
55*a58d3d2aSXin Li /*MWC RNG of George Marsaglia*/
56*a58d3d2aSXin Li static opus_uint32 Rz, Rw;
fast_rand(void)57*a58d3d2aSXin Li static OPUS_INLINE opus_uint32 fast_rand(void)
58*a58d3d2aSXin Li {
59*a58d3d2aSXin Li Rz=36969*(Rz&65535)+(Rz>>16);
60*a58d3d2aSXin Li Rw=18000*(Rw&65535)+(Rw>>16);
61*a58d3d2aSXin Li return (Rz<<16)+Rw;
62*a58d3d2aSXin Li }
63*a58d3d2aSXin Li static opus_uint32 iseed;
64*a58d3d2aSXin Li
65*a58d3d2aSXin Li #ifdef __GNUC__
66*a58d3d2aSXin Li __attribute__((noreturn))
67*a58d3d2aSXin Li #elif defined(_MSC_VER)
68*a58d3d2aSXin Li __declspec(noreturn)
69*a58d3d2aSXin Li #endif
_test_failed(const char * file,int line)70*a58d3d2aSXin Li static OPUS_INLINE void _test_failed(const char *file, int line)
71*a58d3d2aSXin Li {
72*a58d3d2aSXin Li fprintf(stderr,"\n ***************************************************\n");
73*a58d3d2aSXin Li fprintf(stderr," *** A fatal error was detected. ***\n");
74*a58d3d2aSXin Li fprintf(stderr," ***************************************************\n");
75*a58d3d2aSXin Li fprintf(stderr,"Please report this failure and include\n");
76*a58d3d2aSXin Li fprintf(stderr,"'make check SEED=%u fails %s at line %d for %s'\n",iseed,file,line,opus_get_version_string());
77*a58d3d2aSXin Li fprintf(stderr,"and any relevant details about your system.\n\n");
78*a58d3d2aSXin Li #if defined(_MSC_VER)
79*a58d3d2aSXin Li _set_abort_behavior( 0, _WRITE_ABORT_MSG);
80*a58d3d2aSXin Li #endif
81*a58d3d2aSXin Li abort();
82*a58d3d2aSXin Li }
83*a58d3d2aSXin Li #define test_failed() _test_failed(__FILE__, __LINE__);
84*a58d3d2aSXin Li #define opus_test_assert(cond) {if (!(cond)) {test_failed();}}
85*a58d3d2aSXin Li #define expect_true(cond, msg) {if (!(cond)) {fprintf(stderr, "FAIL - %s\n", msg); test_failed();}}
86*a58d3d2aSXin Li void regression_test(void);
87