media/libopus/silk/sort.c

Thu, 22 Jan 2015 13:21:57 +0100

author
Michael Schloh von Bennewitz <michael@schloh.com>
date
Thu, 22 Jan 2015 13:21:57 +0100
branch
TOR_BUG_9701
changeset 15
b8a032363ba2
permissions
-rw-r--r--

Incorporate requested changes from Mozilla in review:
https://bugzilla.mozilla.org/show_bug.cgi?id=1123480#c6

michael@0 1 /***********************************************************************
michael@0 2 Copyright (c) 2006-2011, Skype Limited. All rights reserved.
michael@0 3 Redistribution and use in source and binary forms, with or without
michael@0 4 modification, are permitted provided that the following conditions
michael@0 5 are met:
michael@0 6 - Redistributions of source code must retain the above copyright notice,
michael@0 7 this list of conditions and the following disclaimer.
michael@0 8 - Redistributions in binary form must reproduce the above copyright
michael@0 9 notice, this list of conditions and the following disclaimer in the
michael@0 10 documentation and/or other materials provided with the distribution.
michael@0 11 - Neither the name of Internet Society, IETF or IETF Trust, nor the
michael@0 12 names of specific contributors, may be used to endorse or promote
michael@0 13 products derived from this software without specific prior written
michael@0 14 permission.
michael@0 15 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
michael@0 16 AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
michael@0 17 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
michael@0 18 ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
michael@0 19 LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
michael@0 20 CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
michael@0 21 SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
michael@0 22 INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
michael@0 23 CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
michael@0 24 ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
michael@0 25 POSSIBILITY OF SUCH DAMAGE.
michael@0 26 ***********************************************************************/
michael@0 27
michael@0 28 #ifdef HAVE_CONFIG_H
michael@0 29 #include "config.h"
michael@0 30 #endif
michael@0 31
michael@0 32 /* Insertion sort (fast for already almost sorted arrays): */
michael@0 33 /* Best case: O(n) for an already sorted array */
michael@0 34 /* Worst case: O(n^2) for an inversely sorted array */
michael@0 35 /* */
michael@0 36 /* Shell short: http://en.wikipedia.org/wiki/Shell_sort */
michael@0 37
michael@0 38 #include "SigProc_FIX.h"
michael@0 39
michael@0 40 void silk_insertion_sort_increasing(
michael@0 41 opus_int32 *a, /* I/O Unsorted / Sorted vector */
michael@0 42 opus_int *idx, /* O Index vector for the sorted elements */
michael@0 43 const opus_int L, /* I Vector length */
michael@0 44 const opus_int K /* I Number of correctly sorted positions */
michael@0 45 )
michael@0 46 {
michael@0 47 opus_int32 value;
michael@0 48 opus_int i, j;
michael@0 49
michael@0 50 /* Safety checks */
michael@0 51 silk_assert( K > 0 );
michael@0 52 silk_assert( L > 0 );
michael@0 53 silk_assert( L >= K );
michael@0 54
michael@0 55 /* Write start indices in index vector */
michael@0 56 for( i = 0; i < K; i++ ) {
michael@0 57 idx[ i ] = i;
michael@0 58 }
michael@0 59
michael@0 60 /* Sort vector elements by value, increasing order */
michael@0 61 for( i = 1; i < K; i++ ) {
michael@0 62 value = a[ i ];
michael@0 63 for( j = i - 1; ( j >= 0 ) && ( value < a[ j ] ); j-- ) {
michael@0 64 a[ j + 1 ] = a[ j ]; /* Shift value */
michael@0 65 idx[ j + 1 ] = idx[ j ]; /* Shift index */
michael@0 66 }
michael@0 67 a[ j + 1 ] = value; /* Write value */
michael@0 68 idx[ j + 1 ] = i; /* Write index */
michael@0 69 }
michael@0 70
michael@0 71 /* If less than L values are asked for, check the remaining values, */
michael@0 72 /* but only spend CPU to ensure that the K first values are correct */
michael@0 73 for( i = K; i < L; i++ ) {
michael@0 74 value = a[ i ];
michael@0 75 if( value < a[ K - 1 ] ) {
michael@0 76 for( j = K - 2; ( j >= 0 ) && ( value < a[ j ] ); j-- ) {
michael@0 77 a[ j + 1 ] = a[ j ]; /* Shift value */
michael@0 78 idx[ j + 1 ] = idx[ j ]; /* Shift index */
michael@0 79 }
michael@0 80 a[ j + 1 ] = value; /* Write value */
michael@0 81 idx[ j + 1 ] = i; /* Write index */
michael@0 82 }
michael@0 83 }
michael@0 84 }
michael@0 85
michael@0 86 #ifdef FIXED_POINT
michael@0 87 /* This function is only used by the fixed-point build */
michael@0 88 void silk_insertion_sort_decreasing_int16(
michael@0 89 opus_int16 *a, /* I/O Unsorted / Sorted vector */
michael@0 90 opus_int *idx, /* O Index vector for the sorted elements */
michael@0 91 const opus_int L, /* I Vector length */
michael@0 92 const opus_int K /* I Number of correctly sorted positions */
michael@0 93 )
michael@0 94 {
michael@0 95 opus_int i, j;
michael@0 96 opus_int value;
michael@0 97
michael@0 98 /* Safety checks */
michael@0 99 silk_assert( K > 0 );
michael@0 100 silk_assert( L > 0 );
michael@0 101 silk_assert( L >= K );
michael@0 102
michael@0 103 /* Write start indices in index vector */
michael@0 104 for( i = 0; i < K; i++ ) {
michael@0 105 idx[ i ] = i;
michael@0 106 }
michael@0 107
michael@0 108 /* Sort vector elements by value, decreasing order */
michael@0 109 for( i = 1; i < K; i++ ) {
michael@0 110 value = a[ i ];
michael@0 111 for( j = i - 1; ( j >= 0 ) && ( value > a[ j ] ); j-- ) {
michael@0 112 a[ j + 1 ] = a[ j ]; /* Shift value */
michael@0 113 idx[ j + 1 ] = idx[ j ]; /* Shift index */
michael@0 114 }
michael@0 115 a[ j + 1 ] = value; /* Write value */
michael@0 116 idx[ j + 1 ] = i; /* Write index */
michael@0 117 }
michael@0 118
michael@0 119 /* If less than L values are asked for, check the remaining values, */
michael@0 120 /* but only spend CPU to ensure that the K first values are correct */
michael@0 121 for( i = K; i < L; i++ ) {
michael@0 122 value = a[ i ];
michael@0 123 if( value > a[ K - 1 ] ) {
michael@0 124 for( j = K - 2; ( j >= 0 ) && ( value > a[ j ] ); j-- ) {
michael@0 125 a[ j + 1 ] = a[ j ]; /* Shift value */
michael@0 126 idx[ j + 1 ] = idx[ j ]; /* Shift index */
michael@0 127 }
michael@0 128 a[ j + 1 ] = value; /* Write value */
michael@0 129 idx[ j + 1 ] = i; /* Write index */
michael@0 130 }
michael@0 131 }
michael@0 132 }
michael@0 133 #endif
michael@0 134
michael@0 135 void silk_insertion_sort_increasing_all_values_int16(
michael@0 136 opus_int16 *a, /* I/O Unsorted / Sorted vector */
michael@0 137 const opus_int L /* I Vector length */
michael@0 138 )
michael@0 139 {
michael@0 140 opus_int value;
michael@0 141 opus_int i, j;
michael@0 142
michael@0 143 /* Safety checks */
michael@0 144 silk_assert( L > 0 );
michael@0 145
michael@0 146 /* Sort vector elements by value, increasing order */
michael@0 147 for( i = 1; i < L; i++ ) {
michael@0 148 value = a[ i ];
michael@0 149 for( j = i - 1; ( j >= 0 ) && ( value < a[ j ] ); j-- ) {
michael@0 150 a[ j + 1 ] = a[ j ]; /* Shift value */
michael@0 151 }
michael@0 152 a[ j + 1 ] = value; /* Write value */
michael@0 153 }
michael@0 154 }

mercurial