Category

Development Platform
Home > SourceCode/Document > Mathimatics-Numerical algorithms > View Code
_int_set.c
Package: leda.tar.gz　[view]
Upload User: gzelex
Upload Date: 2007-01-07
Package Size: 707k
Code Size: 2k
Category:

Mathimatics-Numerical algorithms

Development Platform:

MultiPlatform

_int_set.c：Code Content
  1. /*******************************************************************************
  2. +
  3. +  LEDA-R  3.2.3
  4. +
  5. +  _int_set.c
  6. +
  7. +  Copyright (c) 1995  by  Max-Planck-Institut fuer Informatik
  8. +  Im Stadtwald, 66123 Saarbruecken, Germany     
  9. +  All rights reserved.
  10. + 
  11. *******************************************************************************/
  13. #include <LEDA/int_set.h>
  15. typedef unsigned long  word;
  17. int_set::int_set(int n)
  18. { size = n; 
  19.   low = 0;
  20.   register int i = 1+size/SIZE_OF_ULONG;
  21.   if ((V=new word[i]) == 0) error_handler(1,"int_set: out of memory"); 
  22.   while (i--) V[i]=0;
  23.  } 
  25. int_set::int_set(int a, int b)
  26. { size = b-a+1; 
  27.   low = a;
  28.   register int i = 1+size/SIZE_OF_ULONG;
  29.   if ((V=new word[i]) == 0) error_handler(1,"int_set: out of memory"); 
  30.   while (i--) V[i]=0;
  31.  } 
  33. int_set::int_set(const int_set& b)
  34. { size = b.size;
  35.   low  = b.low;
  36.   register int n = 1+size/SIZE_OF_ULONG;
  37.   V = new word[n];
  38.   while (n--) V[n] = b.V[n];
  39. }
  41. int_set& int_set::operator=(const int_set& b)
  42. { if (this == &b) return *this;
  43.   delete V;
  44.   size = b.size;
  45.   low  = b.low;
  46.   register int n = 1+size/SIZE_OF_ULONG;
  47.   V = new word[n];
  48.   while (n--) V[n] = b.V[n];
  49.   return *this;
  50. }
  52. void int_set::clear()
  53. { register int i = 1+size/SIZE_OF_ULONG;
  54.   while (i--) V[i]=0;
  55.  }
  56.   
  58. int_set& int_set::join(const int_set& b) 
  59. { word* stop = V+size/SIZE_OF_ULONG +1;
  60.   word* p;
  61.   word* q;
  62.   for(p = V, q = b.V; p<stop; p++, q++) *p |= *q;
  63.   return *this;
  64.  }
  66. int_set& int_set::intersect(const int_set& b) 
  67. { word* stop = V+size/SIZE_OF_ULONG +1;
  68.   word* p;
  69.   word* q;
  70.   for(p = V, q = b.V; p<stop; p++, q++) *p &= *q;
  71.   return *this;
  72.  }
  74. int_set& int_set::complement() 
  75. { word* stop = V+size/SIZE_OF_ULONG +1;
  76.   for(word* p = V; p<stop; p++) *p = ~(*p);
  77.   return *this;
  78.  }
  80. int_set  int_set::operator|(const int_set& b) 
  81. { int_set res(*this); 
  82.   return res.join(b); 
  83.  }
  85. int_set  int_set::operator&(const int_set& b) 
  86. { int_set res(*this); 
  87.   return res.intersect(b); 
  88.  }
  90. int_set  int_set::operator~()   
  91. { int_set res(*this); 
  92.   return res.complement(); 
  93.  }