home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c017 / 7.ddi / GPPLIB.ZIP / VEC.CC < prev    next >
Encoding:
C/C++ Source or Header  |  1988-08-17  |  7.4 KB  |  373 lines
  1. // This may look like C code, but it is really -*- C++ -*-
  2. /* 
  3. Copyright (C) 1988 Free Software Foundation
  4.     written by Doug Lea (dl@rocky.oswego.edu)
  5.  
  6. This file is part of GNU CC.
  7.  
  8. GNU CC is distributed in the hope that it will be useful,
  9. but WITHOUT ANY WARRANTY.  No author or distributor
  10. accepts responsibility to anyone for the consequences of using it
  11. or for whether it serves any particular purpose or works at all,
  12. unless he says so in writing.  Refer to the GNU CC General Public
  13. License for full details.
  14.  
  15. Everyone is granted permission to copy, modify and redistribute
  16. GNU CC, but only under the conditions described in the
  17. GNU CC General Public License.   A copy of this license is
  18. supposed to have been given to you along with GNU CC so you
  19. can know your rights and responsibilities.  It should be in a
  20. file named COPYING.  Among other things, the copyright notice
  21. and this notice must be preserved on all copies.  
  22. */
  23.  
  24. #include <stream.h>
  25. #include "<T>Vec.h"
  26.  
  27. // error handling
  28.  
  29.  
  30. void default_<T>Vec_error_handler(char* msg)
  31. {
  32.   cerr << "Fatal <T>Vec error. " << msg << "\n";
  33.   exit(1);
  34. }
  35.  
  36. one_arg_error_handler_t <T>Vec_error_handler = default_<T>Vec_error_handler;
  37.  
  38. one_arg_error_handler_t set_<T>Vec_error_handler(one_arg_error_handler_t f)
  39. {
  40.   one_arg_error_handler_t old = <T>Vec_error_handler;
  41.   <T>Vec_error_handler = f;
  42.   return old;
  43. }
  44.  
  45. void <T>Vec::error(const char* msg)
  46. {
  47.   (*<T>Vec_error_handler)(msg);
  48. }
  49.  
  50. void <T>Vec::range_error()
  51. {
  52.   (*<T>Vec_error_handler)("Index out of range.");
  53. }
  54.  
  55.  
  56. // can't just realloc since there may be need for constructors/destructors
  57. void <T>Vec::resize(int newl)
  58. {
  59.   <T>* news = new <T> [newl];
  60.   <T>* p = news;
  61.   int minl = len <? newl;
  62.   <T>* top = &(s[minl]);
  63.   <T>* t = s;
  64.   while (t < top) *p++ = *t++;
  65.   delete [len] s;
  66.   s = news;
  67.   len = newl;
  68. }
  69.  
  70. <T>Vec concat(<T>Vec & a, <T>Vec & b)
  71. {
  72.   int newl = a.len + b.len;
  73.   <T>* news = new <T> [newl];
  74.   <T>* p = news;
  75.   <T>* top = &(a.s[a.len]);
  76.   <T>* t = a.s;
  77.   while (t < top) *p++ = *t++;
  78.   top = &(b.s[b.len]);
  79.   t = b.s;
  80.   while (t < top) *p++ = *t++;
  81.   return <T>Vec(newl, news);
  82. }
  83.  
  84.  
  85. <T>Vec combine(<T>Combiner f, <T>Vec& a, <T>Vec& b)
  86. {
  87.   int newl = a.len <? b.len;
  88.   <T>* news = new <T> [newl];
  89.   <T>* p = news;
  90.   <T>* top = &(a.s[newl]);
  91.   <T>* t = a.s;
  92.   <T>* u = b.s;
  93.   while (t < top) *p++ = (*f)(*t++, *u++);
  94.   return <T>Vec(newl, news);
  95. }
  96.  
  97. <T> <T>Vec::reduce(<T>Combiner f, <T&> base)
  98. {
  99.   <T> r = base;
  100.   <T>* top = &(s[len]);
  101.   <T>* t = s;
  102.   while (t < top) r = (*f)(r, *t++);
  103.   return r;
  104. }
  105.  
  106. <T>Vec reverse(<T>Vec& a)
  107. {
  108.   <T>* news = new <T> [a.len];
  109.   if (a.len != 0)
  110.   {
  111.     <T>* lo = news;
  112.     <T>* hi = &(news[a.len - 1]);
  113.     while (lo < hi)
  114.     {
  115.       <T> tmp = *lo;
  116.       *lo++ = *hi;
  117.       *hi-- = tmp;
  118.     }
  119.   }
  120.   return <T>Vec(a.len, news);
  121. }
  122.  
  123. void <T>Vec::reverse()
  124. {
  125.   if (len != 0)
  126.   {
  127.     <T>* lo = s;
  128.     <T>* hi = &(s[len - 1]);
  129.     while (lo < hi)
  130.     {
  131.       <T> tmp = *lo;
  132.       *lo++ = *hi;
  133.       *hi-- = tmp;
  134.     }
  135.   }
  136. }
  137.  
  138. int <T>Vec::index(<T&> targ)
  139. {
  140.   for (int i = 0; i < len; ++i) if (targ == s[i]) return i;
  141.   return -1;
  142. }
  143.  
  144. <T>Vec map(<T>Mapper f, <T>Vec& a)
  145. {
  146.   <T>* news = new <T> [a.len];
  147.   <T>* p = news;
  148.   <T>* top = &(a.s[a.len]);
  149.   <T>* t = a.s;
  150.   while(t < top) *p++ = (*f)(*t++);
  151.   return <T>Vec(a.len, news);
  152. }
  153.  
  154. int operator == (<T>Vec& a, <T>Vec& b)
  155. {
  156.   if (a.len != b.len)
  157.     return 0;
  158.   <T>* top = &(a.s[a.len]);
  159.   <T>* t = a.s;
  160.   <T>* u = b.s;
  161.   while (t < top) if (*t++ != *u++) return 0;
  162.   return 1;
  163. }
  164.  
  165. void <T>Vec::fill(<T&> val, int from = 0, int n = -1)
  166. {
  167.   int to;
  168.   if (n < 0)
  169.     to = len - 1;
  170.   else
  171.     to = from + n - 1;
  172.   if ((unsigned)from > to)
  173.     range_error();
  174.   <T>* t = &(s[from]);
  175.   <T>* top = &(s[to]);
  176.   while (t <= top) *t++ = val;
  177. }
  178.  
  179. <T>Vec <T>Vec::at(int from = 0, int n = -1)
  180. {
  181.   int to;
  182.   if (n < 0)
  183.   {
  184.     n = len - from;
  185.     to = len - 1;
  186.   }
  187.   else
  188.     to = from + n - 1;
  189.   if ((unsigned)from > to)
  190.     range_error();
  191.   <T>* news = new <T> [n];
  192.   <T>* p = news;
  193.   <T>* t = &(s[from]);
  194.   <T>* top = &(s[to]);
  195.   while (t <= top) *p++ = *t++;
  196.   return <T>Vec(n, news);
  197. }
  198.  
  199. <T>Vec merge(<T>Vec & a, <T>Vec & b, <T>Comparator f)
  200. {
  201.   int newl = a.len + b.len;
  202.   <T>* news = new <T> [newl];
  203.   <T>* p = news;
  204.   <T>* topa = &(a.s[a.len]);
  205.   <T>* as = a.s;
  206.   <T>* topb = &(b.s[b.len]);
  207.   <T>* bs = b.s;
  208.  
  209.   for (;;)
  210.   {
  211.     if (as >= topa)
  212.     {
  213.       while (bs < topb) *p++ = *bs++;
  214.       break;
  215.     }
  216.     else if (bs >= topb)
  217.     {
  218.       while (as < topa) *p++ = *as++;
  219.       break;
  220.     }
  221.     else if ((*f)(*as, *bs) <= 0)
  222.       *p++ = *as++;
  223.     else
  224.       *p++ = *bs++;
  225.   }
  226.   return <T>Vec(newl, news);
  227. }
  228.  
  229. /*
  230.  * a lightly edited version of emacs-18.51 etc/qsort.c
  231.  * The THRESHold below is the insertion sort threshold, and has been adjusted
  232.  * for records of size 48 bytes.
  233.  * The MTHREShold is where we stop finding a better median.
  234.  */
  235.  
  236.  
  237. #define        THRESH        4        /* threshold for insertion */
  238.  
  239. #define        MTHRESH        6        /* threshold for median */
  240.  
  241. static  <T>Comparator qcmp;        /* the comparison routine */
  242.  
  243. static void qst (<T>* base, <T>* max)
  244. {
  245.   <T> *i, *j, *jj, *mid;
  246.   <T> *tmp;
  247.   int lo, hi;
  248.  
  249.   lo = max - base;
  250.   do
  251.   {
  252.     mid = i = base + (lo >> 1);
  253.     if (lo >= MTHRESH)
  254.     {
  255.       j = ((*qcmp) (*(jj = base), *i) > 0 ? jj : i);
  256.       if ((*qcmp) (*j, *(tmp = max - 1)) > 0)
  257.       {
  258.         j = (j == jj ? i : jj);   /* switch to first loser */
  259.         if ((*qcmp) (*j, *tmp) < 0)
  260.           j = tmp;
  261.       }
  262.       if (j != i)
  263.       {
  264.         <T> c = *i;
  265.         *i++ = *j;
  266.         *j++ = c;
  267.       }
  268.     }
  269.     for (i = base, j = max - 1; ;)
  270.     {
  271.       while (i < mid && (*qcmp) (*i, *mid) <= 0)
  272.         i += 1;
  273.       while (j > mid)
  274.       {
  275.         if ((*qcmp) (*mid, *j) <= 0)
  276.         {
  277.           j -= 1;
  278.           continue;
  279.         }
  280.         tmp = i + 1;
  281.         if (i == mid)
  282.         { 
  283.           mid = jj = j;
  284.         }
  285.         else
  286.         { 
  287.           jj = j;
  288.           j -= 1;
  289.         }
  290.         goto  swap;
  291.       }
  292.       if (i == mid)
  293.       {
  294.         break;
  295.       }
  296.       else
  297.       {               /* i <-> mid, new mid is i */
  298.         jj = mid;
  299.         tmp = mid = i;        /* value of i after swap */
  300.         j -= 1;
  301.       }
  302.     swap:
  303.       <T> c = *i;
  304.       *i++ = *jj;
  305.       *jj++ = c;
  306.       i = tmp;
  307.     }
  308.     i = (j = mid) + 1;
  309.     if ((lo = j - base) <= (hi = max - i))
  310.     {
  311.       if (lo >= THRESH)
  312.         qst (base, j);
  313.       base = i;
  314.       lo = hi;
  315.     }
  316.     else
  317.     {
  318.       if (hi >= THRESH)
  319.         qst (i, max);
  320.       max = j;
  321.     }
  322.   }
  323.   while (lo >= THRESH);
  324. }
  325.  
  326. void <T>Vec::sort (<T>Comparator compar)
  327. {
  328.   <T>* base = s;
  329.   <T> *i, *j, *lo, *hi, *min;
  330.  
  331.   if (len <= 1)  return;
  332.   qcmp = compar;
  333.  
  334.   <T>* max = base + len;
  335.  
  336.   if (len >= THRESH)
  337.   {
  338.     qst (base, max);
  339.     hi = base + THRESH;
  340.   }
  341.   else
  342.   {
  343.     hi = max;
  344.   }
  345.   for (j = lo = base; (lo += 1) < hi; )
  346.   {
  347.     if ((*qcmp) (*j, *lo) > 0)
  348.       j = lo;
  349.   }
  350.   if (j != base)
  351.   { 
  352.     <T> c = *j;
  353.     *j++ = *base;
  354.     *i  = c;
  355.   }
  356.   for (min = base; (hi = min += 1) < max;)
  357.   {
  358.     while ( (*qcmp) (*(hi -= 1), *min) > 0);
  359.     if ((hi += 1) != min)
  360.     {
  361.       for (lo = min + 1; --lo >= min;)
  362.       {
  363.         <T> c = *lo;
  364.         for (i = j = lo; (j -= 1) >= hi; i = j)
  365.           *i = *j;
  366.         *i = c;
  367.       }
  368.     }
  369.   }
  370. }
  371.  
  372.  
  373.