Statistics
| Branch: | Revision:

ffmpeg / libavutil / tree.c @ d7fb5a18

History | View | Annotate | Download (6.3 KB)

1 9eb88fde Michael Niedermayer
/*
2
 * copyright (c) 2006 Michael Niedermayer <michaelni@gmx.at>
3
 *
4
 * This file is part of FFmpeg.
5
 *
6
 * FFmpeg is free software; you can redistribute it and/or
7
 * modify it under the terms of the GNU Lesser General Public
8
 * License as published by the Free Software Foundation; either
9
 * version 2.1 of the License, or (at your option) any later version.
10
 *
11
 * FFmpeg is distributed in the hope that it will be useful,
12
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14
 * Lesser General Public License for more details.
15
 *
16
 * You should have received a copy of the GNU Lesser General Public
17
 * License along with FFmpeg; if not, write to the Free Software
18
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19
 */
20
21
#include "common.h"
22
#include "log.h"
23
#include "tree.h"
24
25
typedef struct AVTreeNode{
26
    struct AVTreeNode *child[2];
27
    void *elem;
28
    int state;
29
}AVTreeNode;
30
31 6e8b982b Michael Niedermayer
const int av_tree_node_size = sizeof(AVTreeNode);
32
33 9eb88fde Michael Niedermayer
void *av_tree_find(const AVTreeNode *t, void *key, int (*cmp)(void *key, const void *b), void *next[2]){
34
    if(t){
35 2e1d2873 Michael Niedermayer
        unsigned int v= cmp(key, t->elem);
36 9eb88fde Michael Niedermayer
        if(v){
37 2e1d2873 Michael Niedermayer
            if(next) next[v>>31]= t->elem;
38
            return av_tree_find(t->child[(v>>31)^1], key, cmp, next);
39 9eb88fde Michael Niedermayer
        }else{
40 116d15cc Michael Niedermayer
            if(next){
41
                av_tree_find(t->child[0], key, cmp, next);
42
                av_tree_find(t->child[1], key, cmp, next);
43
            }
44 9eb88fde Michael Niedermayer
            return t->elem;
45
        }
46
    }
47
    return NULL;
48
}
49
50 6e8b982b Michael Niedermayer
void *av_tree_insert(AVTreeNode **tp, void *key, int (*cmp)(void *key, const void *b), AVTreeNode **next){
51 9eb88fde Michael Niedermayer
    AVTreeNode *t= *tp;
52
    if(t){
53
        unsigned int v= cmp(t->elem, key);
54 f05dda3b Michael Niedermayer
        void *ret;
55
        if(!v){
56
            if(*next)
57
                return t->elem;
58
            else if(t->child[0]||t->child[1]){
59
                int i= !t->child[0];
60
                void *next_elem[2];
61 3f161c7e Michael Niedermayer
                av_tree_find(t->child[i], key, cmp, next_elem);
62 f05dda3b Michael Niedermayer
                key= t->elem= next_elem[i];
63
                v= -i;
64
            }else{
65
                *next= t;
66
                *tp=NULL;
67
                return NULL;
68
            }
69
        }
70 a35bf971 Michael Niedermayer
        ret= av_tree_insert(&t->child[v>>31], key, cmp, next);
71
        if(!ret){
72
            int i= (v>>31) ^ !!*next;
73
            AVTreeNode **child= &t->child[i];
74
            t->state += 2*i - 1;
75 f05dda3b Michael Niedermayer
76 a35bf971 Michael Niedermayer
            if(!(t->state&1)){
77
                if(t->state){
78 51198a87 Michael Niedermayer
                    /* The following code is equivalent to
79
                    if((*child)->state*2 == -t->state)
80
                        rotate(child, i^1);
81
                    rotate(tp, i);
82

83
                    with rotate():
84
                    static void rotate(AVTreeNode **tp, int i){
85
                        AVTreeNode *t= *tp;
86

87
                        *tp= t->child[i];
88
                        t->child[i]= t->child[i]->child[i^1];
89
                        (*tp)->child[i^1]= t;
90
                        i= 4*t->state + 2*(*tp)->state + 12;
91
                          t  ->state=                     ((0x614586 >> i) & 3)-1;
92
                        (*tp)->state= ((*tp)->state>>1) + ((0x400EEA >> i) & 3)-1;
93
                    }
94
                    but such a rotate function is both bigger and slower
95
                    */
96 a35bf971 Michael Niedermayer
                    if((*child)->state*2 == -t->state){
97
                        *tp= (*child)->child[i^1];
98
                        (*child)->child[i^1]= (*tp)->child[i];
99
                        (*tp)->child[i]= *child;
100
                        *child= (*tp)->child[i^1];
101
                        (*tp)->child[i^1]= t;
102 9eb88fde Michael Niedermayer
103 a35bf971 Michael Niedermayer
                        (*tp)->child[0]->state= -((*tp)->state>0);
104
                        (*tp)->child[1]->state=   (*tp)->state<0 ;
105
                        (*tp)->state=0;
106
                    }else{
107
                        *tp= *child;
108
                        *child= (*child)->child[i^1];
109
                        (*tp)->child[i^1]= t;
110
                        if((*tp)->state) t->state  = 0;
111
                        else             t->state>>= 1;
112
                        (*tp)->state= -t->state;
113 9eb88fde Michael Niedermayer
                    }
114
                }
115
            }
116 a35bf971 Michael Niedermayer
            if(!(*tp)->state ^ !!*next)
117
                return key;
118
        }
119
        return ret;
120 9eb88fde Michael Niedermayer
    }else{
121 6e8b982b Michael Niedermayer
        *tp= *next; *next= NULL;
122 eed36075 Michael Niedermayer
        if(*tp){
123
            (*tp)->elem= key;
124
            return NULL;
125
        }else
126
            return key;
127 9eb88fde Michael Niedermayer
    }
128
}
129
130
void av_tree_destroy(AVTreeNode *t){
131
    av_tree_destroy(t->child[0]);
132
    av_tree_destroy(t->child[1]);
133
    av_free(t);
134
}
135
136
#if 0
137
void av_tree_enumerate(AVTreeNode *t, void *opaque, int (*f)(void *opaque, void *elem)){
138 33206876 Michael Niedermayer
    int v= f(opaque, t->elem);
139
    if(v>=0) av_tree_enumerate(t->child[0], opaque, f);
140
    if(v<=0) av_tree_enumerate(t->child[1], opaque, f);
141 9eb88fde Michael Niedermayer
}
142
#endif
143
144
#ifdef TEST
145 d5cb5fe8 Michael Niedermayer
#undef random
146 9eb88fde Michael Niedermayer
static int check(AVTreeNode *t){
147
    if(t){
148
        int left= check(t->child[0]);
149
        int right= check(t->child[1]);
150
151
        if(left>999 || right>999)
152
            return 1000;
153
        if(right - left != t->state)
154
            return 1000;
155
        if(t->state>1 || t->state<-1)
156
            return 1000;
157
        return FFMAX(left, right)+1;
158
    }
159
    return 0;
160
}
161
162
static void print(AVTreeNode *t, int depth){
163
    int i;
164
    for(i=0; i<depth*4; i++) av_log(NULL, AV_LOG_ERROR, " ");
165
    if(t){
166
        av_log(NULL, AV_LOG_ERROR, "Node %p %2d %4d\n", t, t->state, t->elem);
167
        print(t->child[0], depth+1);
168
        print(t->child[1], depth+1);
169
    }else
170
        av_log(NULL, AV_LOG_ERROR, "NULL\n");
171
}
172
173
int cmp(const void *a, const void *b){
174
    return a-b;
175
}
176
177 f8a80fd6 Diego Biurrun
int main(void){
178 a005768d Diego Biurrun
    int i,k;
179 f05dda3b Michael Niedermayer
    AVTreeNode *root= NULL, *node=NULL;
180 9eb88fde Michael Niedermayer
181
    for(i=0; i<10000; i++){
182 f05dda3b Michael Niedermayer
        int j= (random()%86294);
183 9eb88fde Michael Niedermayer
        if(check(root) > 999){
184
            av_log(NULL, AV_LOG_ERROR, "FATAL error %d\n", i);
185
        print(root, 0);
186
            return -1;
187
        }
188
        av_log(NULL, AV_LOG_ERROR, "inserting %4d\n", j);
189 f05dda3b Michael Niedermayer
        if(!node)
190
            node= av_mallocz(av_tree_node_size);
191
        av_tree_insert(&root, (void*)(j+1), cmp, &node);
192
193
        j= (random()%86294);
194 eed36075 Michael Niedermayer
        {
195 f05dda3b Michael Niedermayer
            AVTreeNode *node2=NULL;
196 5d1e3d22 Michael Niedermayer
            av_log(NULL, AV_LOG_ERROR, "removing %4d\n", j);
197 f05dda3b Michael Niedermayer
            av_tree_insert(&root, (void*)(j+1), cmp, &node2);
198
            k= av_tree_find(root, (void*)(j+1), cmp, NULL);
199
            if(k)
200
                av_log(NULL, AV_LOG_ERROR, "removial failure %d\n", i);
201
        }
202 9eb88fde Michael Niedermayer
    }
203
    return 0;
204
}
205
#endif