patches

[SXSI/XMLTree.git] / libcds / src / static_sequence / static_sequence.cpp

diff --git a/libcds/src/static_sequence/static_sequence.cpp b/libcds/src/static_sequence/static_sequence.cpp
index e97a112..05b7753 100644 (file)
--- a/libcds/src/static_sequence/static_sequence.cpp
+++ b/libcds/src/static_sequence/static_sequence.cpp
@@ -41,3 +41,42 @@ static_sequence * static_sequence::load(FILE * fp) {
   }
   return NULL;
 }
+
+bool static_sequence::test(uint * seq, uint n) {
+       uint sigma = 0;
+       for(uint i=0;i<n;i++)
+               sigma = max(sigma,seq[i]);
+       uint * occ = new uint[sigma+1];
+       for(uint i=0;i<=sigma;i++)
+               occ[i] = 0;
+       for(uint i=0;i<n;i++) {
+               occ[seq[i]]++;
+               if(rank(seq[i],i)!=occ[seq[i]]) {
+                       cout << "rank failed!" << endl;
+                       cout << "rank("<<seq[i]<<","<<i<<")="<<rank(seq[i],i)<<endl;
+                       cout << "expected result: " << occ[seq[i]] << endl;
+                       delete [] occ;
+                       return false;
+               }
+               if(i>0 && rank(seq[i],i-1)!=occ[seq[i]]-1) {
+                       cout << "rank-1 failed!" << endl;
+                       delete [] occ;
+                       return false;
+               }
+               if(select(seq[i],occ[seq[i]])!=i) {
+                       cout << "select failed!" << endl;
+                       cout << "select(" << seq[i] << "," << occ[seq[i]] << ")="<<select(seq[i],occ[seq[i]]) << endl;
+                       cout << "i=" << i << "  rank(" << seq[i] << ",i)=" << rank(seq[i],i) << endl;
+                       delete [] occ;
+                       return false;
+               }
+               if(access(i)!=seq[i]) {
+                       cout << "access failed!" << endl;
+                       delete [] occ;
+                       return false;
+               }
+       }
+       delete [] occ;
+       return true;
+}
+