Free cookie consent management tool by TermsFeed Policy Generator

source: branches/HeuristicLab.Problems.GaussianProcessTuning/ILNumerics.2.14.4735.573/Data/ILBTree.cs @ 12547

Last change on this file since 12547 was 9102, checked in by gkronber, 12 years ago

#1967: ILNumerics source for experimentation

File size: 18.4 KB
RevLine 
[9102]1///
2///    This file is part of ILNumerics Community Edition.
3///
4///    ILNumerics Community Edition - high performance computing for applications.
5///    Copyright (C) 2006 - 2012 Haymo Kutschbach, http://ilnumerics.net
6///
7///    ILNumerics Community Edition is free software: you can redistribute it and/or modify
8///    it under the terms of the GNU General Public License version 3 as published by
9///    the Free Software Foundation.
10///
11///    ILNumerics Community Edition 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
14///    GNU General Public License for more details.
15///
16///    You should have received a copy of the GNU General Public License
17///    along with ILNumerics Community Edition. See the file License.txt in the root
18///    of your distribution package. If not, see <http://www.gnu.org/licenses/>.
19///
20///    In addition this software uses the following components and/or licenses:
21///
22///    =================================================================================
23///    The Open Toolkit Library License
24///   
25///    Copyright (c) 2006 - 2009 the Open Toolkit library.
26///   
27///    Permission is hereby granted, free of charge, to any person obtaining a copy
28///    of this software and associated documentation files (the "Software"), to deal
29///    in the Software without restriction, including without limitation the rights to
30///    use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of
31///    the Software, and to permit persons to whom the Software is furnished to do
32///    so, subject to the following conditions:
33///
34///    The above copyright notice and this permission notice shall be included in all
35///    copies or substantial portions of the Software.
36///
37///    =================================================================================
38///   
39
40using System;
41using System.Collections.Generic;
42using System.Linq;
43using System.Text;
44using keytype = System.Int32;
45
46namespace ILNumerics.Data {
47    /// <summary>
48    /// incomplete BTree implementation (it lags of a remove method)
49    /// </summary>
50    /// <typeparam name="T">value element type</typeparam>
51    public class ILBTree<T> {
52
53        //private static readonly int ORDER = 4;   // max number of childs
54        //private static readonly int MAX_KEYS = 3;  // uneven ! (order - 1)
55        //private static readonly int MIN_CHILDS = 2; 
56        //private static readonly int MIN_KEYS = 1;  // uneven ! (order / 2 - 1)       
57        private static int m_order = 24;
58        public static int ORDER  {
59            get { return m_order; }
60            set {
61                if (value / 2 * 2 != value) {
62                    value++;
63                }
64                m_order = value;
65                MAX_KEYS = value -1;
66                MIN_KEYS = value / 2 - 1;
67            }
68        }   // max number of childs
69
70        public static int MAX_KEYS = 15;  // uneven ! (order - 1)
71        public static int MIN_KEYS = 7;  // uneven ! (order / 2 - 1)
72        public static int DEFAULT_LIST_CAPACITY = 8;
73
74        private int m_nodeCount = 0;
75        private int m_levelCount = 1;
76       
77        private nodetype m_root;
78
79        //private class itemtype {
80        //    public keytype Key;
81        //    public List<T> Values; 
82        //    //private itemtype () {
83        //    //    Key = -1;
84        //    //    Values = new LinkedList<T>();
85        //    //}
86        //    public itemtype (keytype key, T value) {
87        //        Key = key;
88        //        Values = new List<T>(8);
89        //        Values.Add(value); 
90        //    }
91        //    public override string ToString() {
92        //        return String.Format("K:{0} Count:{1}",Key,Values.Count);
93        //    }
94        //}
95        private class nodetype {
96            //public itemtype[] Items;
97            public int[] Keys;
98            public List<T>[] Values;
99            public nodetype[] Childs;
100            /// <summary>
101            /// number of keys in the node
102            /// </summary>
103            public int Count;
104
105            public nodetype() {
106                Count = 0;
107                //Items = new itemtype[ MAX_KEYS ];
108                Keys = new int[MAX_KEYS];
109                Values = new List<T>[MAX_KEYS];
110                Childs = new nodetype[ MAX_KEYS + 1 ];
111            }
112            public override string ToString() {
113                string itemDesc = "---";
114                if (Count > 0) {
115                    itemDesc = String.Format("{0} - {1}",Keys[0],Keys[Count-1]);
116                }
117                return String.Format("{0} Items:[{1}]",Count,itemDesc);
118            }
119
120        }
121
122        #region constructor
123        public ILBTree() {
124            m_root = new nodetype();
125            m_nodeCount = 1;
126        }
127        #endregion
128
129        #region public interface
130
131        public int NodeCount {
132            get { return m_nodeCount; }
133        }
134        public int LevelCount {
135            get { return m_levelCount; }
136        }
137
138
139        public void Add(keytype key, T value) {
140            if (m_root.Count >= MAX_KEYS) {
141                nodetype newright = new nodetype();
142                keytype upkey;
143                List<T> upval;
144                SplitNode(m_root,newright, out upkey, out upval);
145                nodetype newRoot = new nodetype();
146                newRoot.Childs[0] = m_root;
147                newRoot.Childs[1] = newright;
148                newRoot.Keys[0] = upkey;
149                newRoot.Values[0] = upval;
150                newRoot.Count = 1;
151                m_root = newRoot;
152                m_nodeCount += 2;
153                m_levelCount++;
154            }
155            AddValue(key, value, m_root);
156        }
157        #endregion
158       
159        #region private helpers
160        /// <summary>
161        ///  is the key contained in curNode ?
162        /// </summary>
163        /// <param name="curNode"></param>
164        /// <param name="key"></param>
165        /// <param name="location"></param>
166        /// <returns></returns>
167        private unsafe static bool SearchNode(nodetype curNode, keytype key, ref int location) {
168            int loc = 0, max = curNode.Count; int[] curNodeKeys = curNode.Keys;
169            for (; loc < max; loc++) {
170                int curKey = curNodeKeys[loc];
171                if (curKey == key) {
172                    location = loc;
173                    return true;
174                } else if (curKey > key) {
175                    location = loc;
176                    return false;
177                }
178            }
179            location = loc;
180            return false;
181
182            //int loc = 0, max = curNode.Count;
183            //foreach (keytype curKey in curNode.Keys) {
184            //    if (curKey == key) {
185            //        location = loc;
186            //        return true;
187            //    } else if (curKey > key) {
188            //        location = loc;
189            //        return false;
190            //    }
191            //    if (++loc >= max) break;
192            //}
193            //location = loc;
194            //return false;
195
196            //int loc = curNode.Count - 1;
197            //fixed (int* keys = curNode.Keys) {
198            //    int* curKey = keys + loc;
199            //    do {
200            //        if (*curKey == key) {
201            //            location = loc;
202            //            return true;
203            //        } else if (*curKey < key) {
204            //            location = ++loc;
205            //            return false;
206            //        }
207            //        curKey--;
208            //    } while (loc --> 0);
209            //}
210            //location = 0;
211            //return false;
212        }
213        private void AddValue(keytype key, T value, nodetype curNode) {
214            System.Diagnostics.Debug.Assert(curNode.Count < MAX_KEYS);
215            while (true) {
216                // search node (inlined from func)
217                short max = (short)curNode.Count;
218                short location = 0;
219                int[] curNodeKeys = curNode.Keys;
220                for (; location < max; location++) {
221                    int curKey = curNodeKeys[location];
222                    if (curKey == key) {
223                        curNode.Values[location].Add(value);
224                        return;
225                    }
226                    if (curKey > key) {
227                        break;
228                    }
229                }
230                nodetype childNode = curNode.Childs[location];
231                if (childNode == null) {
232                    // insert right here
233                    var newVals = new List<T>(8);
234                    newVals.Add(value);
235                    InsertItemIntoNode(curNode,key,newVals,location,null); 
236                    break;
237                } else {
238                    // check if the child is filled up allready
239                    if (childNode.Count >= MAX_KEYS) {
240                        nodetype newRightChild = new nodetype();
241                        int upkey;
242                        List<T> upval;
243                        SplitNode(childNode, newRightChild, out upkey, out upval);
244                        m_nodeCount++;
245                        InsertItemIntoNode(curNode,upkey, upval, location, newRightChild);  // cur node should always have space left!
246                        // in which subtree to go?
247                        if (key == upkey) {
248                            curNode.Values[location].Add(value);
249                            return;
250                        }
251                        if (key > upkey)
252                            location++;
253                    }
254                    //AddValue(key, value, curNode.Childs[location]);
255                    curNode = curNode.Childs[location];
256                }
257            }
258        }
259
260        //private void AddValue(keytype key, T value, nodetype curNode) {
261        //    System.Diagnostics.Debug.Assert(curNode.Count < MAX_KEYS);
262        //    int location;
263        //    if (SearchNode(curNode, key, out location)) {
264        //        curNode.Values[location].Add(value);
265        //    } else {
266        //        nodetype childNode = curNode.Childs[location];
267        //        if (childNode == null) {
268        //            // insert right here
269        //            var newVals = new List<T>(8);
270        //            newVals.Add(value);
271        //            InsertItemIntoNode(curNode,key,newVals,location,null); 
272        //        } else {
273        //            // check if the child is filled up allready
274        //            if (childNode.Count >= MAX_KEYS) {
275        //                nodetype newRightChild = new nodetype();
276        //                int upkey;
277        //                List<T> upval;
278        //                SplitNode(childNode, newRightChild, out upkey, out upval);
279        //                m_nodeCount++;
280        //                InsertItemIntoNode(curNode,upkey, upval, location, newRightChild);  // cur node should always have space left!
281        //                // in which subtree to go?
282        //                if (key == upkey) {
283        //                    curNode.Values[location].Add(value);
284        //                    return;
285        //                }
286        //                if (key > upkey)
287        //                    location++;
288        //            }
289        //            AddValue(key, value, curNode.Childs[location]);
290        //        }
291        //    }
292        //}
293
294        #region break
295        #endregion
296
297        //private void AddValue(keytype key, T value, nodetype curNode, ref itemtype upcomingItem, ref nodetype rightchild) {
298        //    if (curNode == null) {
299        //        // we came from a leaf! let the parent handle the insertion
300        //        upcomingItem = new itemtype(key,value);
301        //        rightchild = null;
302        //        return;
303        //    }
304        //    System.Diagnostics.Debug.Assert(curNode.Count < MAX_KEYS);
305        //    int location;
306        //    if (SearchNode(curNode, key, out location)) {
307        //        curNode.Items[location].Values.AddFirst(value);
308        //        upcomingItem = null;
309        //    } else {
310        //        nodetype childNode = curNode.Childs[location];
311        //        if (childNode.Count >= MAX_KEYS) {
312        //            nodetype newRightChild = new nodetype();
313        //            itemtype newUpItem;
314        //            SplitNode(childNode, newRightChild, out newUpItem);
315        //            InsertItemIntoNode(curNode,newUpItem,location,newRightChild);  // cur node should always have space left!
316        //            // in which subtree to go?
317        //            if (key > newUpItem.Key)
318        //                location++;
319        //            if (key == newUpItem.Key) {
320        //                curNode.Items[location].Values.AddFirst(value);
321        //                upcomingItem = null;
322        //                return;
323        //            }
324        //        }
325        //        // location = index of child/ insertion point
326        //        AddValue(key, value, curNode.Childs[location], ref upcomingItem, ref rightchild);
327        //        // handle upcoming node
328        //        if (upcomingItem != null) {
329        //            if (curNode.Count < MAX_KEYS) {
330        //                // insert right here
331        //                InsertItemIntoNode(curNode, upcomingItem, location, rightchild);
332        //                upcomingItem = null;
333        //            } else {
334        //                // split the node and move median up
335        //                SplitNode(curNode, upcomingItem, location, ref upcomingItem, ref rightchild);
336        //            }
337        //        }
338        //    }
339        //}
340
341        private void SplitNode(nodetype curNode, nodetype newRight, out keytype upKey, out List<T> upval) {
342            System.Diagnostics.Debug.Assert(Math.Round(curNode.Count / 2.0) != curNode.Count / 2.0);
343            System.Diagnostics.Debug.Assert(curNode.Count == (MIN_KEYS * 2 + 1));
344            System.Diagnostics.Debug.Assert(curNode.Count == MAX_KEYS);
345            System.Diagnostics.Debug.Assert(newRight != null && newRight.Childs != null && newRight.Keys != null);
346            System.Diagnostics.Debug.Assert(newRight.Childs.Length == ORDER && newRight.Keys.Length == MAX_KEYS && newRight.Values.Length == MAX_KEYS);
347            for (int i = MIN_KEYS + 1, ni = 0; ni < MIN_KEYS; ni++, i++) {
348                newRight.Keys[ni] = curNode.Keys[i];
349                newRight.Values[ni] = curNode.Values[i];
350                newRight.Childs[ni+1] = curNode.Childs[i+1];
351                // clear old data
352                curNode.Keys[i] = -1;
353                curNode.Values[i] = null;
354                curNode.Childs[i+1] = null;
355            }
356            curNode.Count = MIN_KEYS;
357            newRight.Count = MIN_KEYS;
358            newRight.Childs[0] = curNode.Childs[MIN_KEYS + 1];
359            upKey = curNode.Keys[MIN_KEYS];
360            upval = curNode.Values[MIN_KEYS];
361            // clear old data
362            curNode.Keys[MIN_KEYS] = -1;
363            curNode.Values[MIN_KEYS] = null;
364            curNode.Childs[MIN_KEYS + 1] = null;
365        }
366
367        //private void SplitNode(nodetype curNode, itemtype newItem, int location, ref itemtype upcomingItem, ref nodetype upcomingrightnode) {
368        //    var newnode = new nodetype();
369        //    if (location <= MIN_KEYS) {
370        //        // copy right node 
371        //        for (int ci = MIN_KEYS, ni = 0; ci < MAX_KEYS; ci++, ni++ ) {
372        //            newnode.Items[ni] = curNode.Items[ci];
373        //            newnode.Childs[ni+1] = curNode.Childs[ci];
374        //        }
375        //        // copy curNode
376        //        for (int ci = MIN_KEYS; ci --> location; ) {
377        //            curNode.Items[ci+1] = curNode.Items[ci];
378        //            curNode.Childs[ci+1] = curNode.Childs[ci];
379        //        }
380        //        curNode.Count = MIN_KEYS + 1;
381        //        newnode.Count = MAX_KEYS - MIN_KEYS;
382        //        // insert new item 
383        //        curNode.Items[location] = newItem;
384        //        curNode.Childs[location] = upcomingrightnode;
385        //        newnode.Childs[0] = curNode.Childs[MIN_KEYS];
386        //        // give the median back upwards
387        //        upcomingItem = curNode.Items[MIN_KEYS];
388        //        upcomingrightnode = newnode;
389        //    } else {
390        //        // copy right node -> up to insert location
391        //        int i = MIN_KEYS, ni = 0;
392        //        for (; i < location; i++, ni++) {
393        //            newnode.Items[ni] = curNode.Items[i];
394        //            newnode.Childs[ni+1] = curNode.Childs[i];
395        //        }
396        //        // insert new item
397        //        newnode.Items[ni] = newItem;
398        //        newnode.Childs[ni++] = upcomingrightnode;
399        //        // copy rest of node to right side
400        //        for (; i < curNode.Count; i++, ni++) {
401        //            newnode.Items[ni] = curNode.Items[i];
402        //            newnode.Childs[ni+1] = curNode.Childs[i];
403        //        }
404        //        // update curNode
405        //        curNode.Count = MIN_KEYS + 1;
406        //        newnode.Count = ni;
407        //        newnode.Childs[0] = curNode.Childs[MIN_KEYS + 1];
408        //        upcomingItem = curNode.Items[MIN_KEYS + 1];
409        //        upcomingrightnode = newnode;
410        //    }
411        //    m_nodeCount++;
412        //}
413
414        private void InsertItemIntoNode(nodetype curNode, keytype newKey,List<T> newvalue, int location, nodetype rigthchild) {
415            int a = curNode.Count;
416            for (; a --> location; ) {
417                curNode.Keys[a+1] = curNode.Keys[a];
418                curNode.Values[a+1] = curNode.Values[a];
419                curNode.Childs[a+2] = curNode.Childs[a+1];
420            }
421            // insert a new linked list at freed position
422            curNode.Keys[location] = newKey;
423            curNode.Values[location] = newvalue;
424            curNode.Childs[location+1] = rigthchild;
425            curNode.Count++;
426        }
427
428        #endregion
429
430    }
431}
Note: See TracBrowser for help on using the repository browser.