yaffs_qsort.c 4.6 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163
  1. /*
  2. * Copyright (c) 1992, 1993
  3. * The Regents of the University of California. All rights reserved.
  4. *
  5. * Redistribution and use in source and binary forms, with or without
  6. * modification, are permitted provided that the following conditions
  7. * are met:
  8. * 1. Redistributions of source code must retain the above copyright
  9. * notice, this list of conditions and the following disclaimer.
  10. * 2. Redistributions in binary form must reproduce the above copyright
  11. * notice, this list of conditions and the following disclaimer in the
  12. * documentation and/or other materials provided with the distribution.
  13. * 3. Neither the name of the University nor the names of its contributors
  14. * may be used to endorse or promote products derived from this software
  15. * without specific prior written permission.
  16. *
  17. * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
  18. * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  19. * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  20. * ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
  21. * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  22. * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  23. * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  24. * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  25. * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  26. * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  27. * SUCH DAMAGE.
  28. */
  29. /* XXX U-BOOT XXX */
  30. #include <common.h>
  31. #include "yportenv.h"
  32. //#include <linux/string.h>
  33. /*
  34. * Qsort routine from Bentley & McIlroy's "Engineering a Sort Function".
  35. */
  36. #define swapcode(TYPE, parmi, parmj, n) { \
  37. long i = (n) / sizeof (TYPE); \
  38. register TYPE *pi = (TYPE *) (parmi); \
  39. register TYPE *pj = (TYPE *) (parmj); \
  40. do { \
  41. register TYPE t = *pi; \
  42. *pi++ = *pj; \
  43. *pj++ = t; \
  44. } while (--i > 0); \
  45. }
  46. #define SWAPINIT(a, es) swaptype = ((char *)a - (char *)0) % sizeof(long) || \
  47. es % sizeof(long) ? 2 : es == sizeof(long)? 0 : 1;
  48. static __inline void
  49. swapfunc(char *a, char *b, int n, int swaptype)
  50. {
  51. if (swaptype <= 1)
  52. swapcode(long, a, b, n)
  53. else
  54. swapcode(char, a, b, n)
  55. }
  56. #define swap(a, b) \
  57. if (swaptype == 0) { \
  58. long t = *(long *)(a); \
  59. *(long *)(a) = *(long *)(b); \
  60. *(long *)(b) = t; \
  61. } else \
  62. swapfunc(a, b, es, swaptype)
  63. #define vecswap(a, b, n) if ((n) > 0) swapfunc(a, b, n, swaptype)
  64. static __inline char *
  65. med3(char *a, char *b, char *c, int (*cmp)(const void *, const void *))
  66. {
  67. return cmp(a, b) < 0 ?
  68. (cmp(b, c) < 0 ? b : (cmp(a, c) < 0 ? c : a ))
  69. :(cmp(b, c) > 0 ? b : (cmp(a, c) < 0 ? a : c ));
  70. }
  71. #ifndef min
  72. #define min(a,b) (((a) < (b)) ? (a) : (b))
  73. #endif
  74. void
  75. yaffs_qsort(void *aa, size_t n, size_t es,
  76. int (*cmp)(const void *, const void *))
  77. {
  78. char *pa, *pb, *pc, *pd, *pl, *pm, *pn;
  79. int d, r, swaptype, swap_cnt;
  80. register char *a = aa;
  81. loop: SWAPINIT(a, es);
  82. swap_cnt = 0;
  83. if (n < 7) {
  84. for (pm = (char *)a + es; pm < (char *) a + n * es; pm += es)
  85. for (pl = pm; pl > (char *) a && cmp(pl - es, pl) > 0;
  86. pl -= es)
  87. swap(pl, pl - es);
  88. return;
  89. }
  90. pm = (char *)a + (n / 2) * es;
  91. if (n > 7) {
  92. pl = (char *)a;
  93. pn = (char *)a + (n - 1) * es;
  94. if (n > 40) {
  95. d = (n / 8) * es;
  96. pl = med3(pl, pl + d, pl + 2 * d, cmp);
  97. pm = med3(pm - d, pm, pm + d, cmp);
  98. pn = med3(pn - 2 * d, pn - d, pn, cmp);
  99. }
  100. pm = med3(pl, pm, pn, cmp);
  101. }
  102. swap(a, pm);
  103. pa = pb = (char *)a + es;
  104. pc = pd = (char *)a + (n - 1) * es;
  105. for (;;) {
  106. while (pb <= pc && (r = cmp(pb, a)) <= 0) {
  107. if (r == 0) {
  108. swap_cnt = 1;
  109. swap(pa, pb);
  110. pa += es;
  111. }
  112. pb += es;
  113. }
  114. while (pb <= pc && (r = cmp(pc, a)) >= 0) {
  115. if (r == 0) {
  116. swap_cnt = 1;
  117. swap(pc, pd);
  118. pd -= es;
  119. }
  120. pc -= es;
  121. }
  122. if (pb > pc)
  123. break;
  124. swap(pb, pc);
  125. swap_cnt = 1;
  126. pb += es;
  127. pc -= es;
  128. }
  129. if (swap_cnt == 0) { /* Switch to insertion sort */
  130. for (pm = (char *) a + es; pm < (char *) a + n * es; pm += es)
  131. for (pl = pm; pl > (char *) a && cmp(pl - es, pl) > 0;
  132. pl -= es)
  133. swap(pl, pl - es);
  134. return;
  135. }
  136. pn = (char *)a + n * es;
  137. r = min(pa - (char *)a, pb - pa);
  138. vecswap(a, pb - r, r);
  139. r = min((long)(pd - pc), (long)(pn - pd - es));
  140. vecswap(pb, pn - r, r);
  141. if ((r = pb - pa) > es)
  142. yaffs_qsort(a, r / es, es, cmp);
  143. if ((r = pd - pc) > es) {
  144. /* Iterate rather than recurse to save stack space */
  145. a = pn - r;
  146. n = r / es;
  147. goto loop;
  148. }
  149. /* yaffs_qsort(pn - r, r / es, es, cmp);*/
  150. }