X-Git-Url: http://git.nguyen.vg/gitweb/?a=blobdiff_plain;f=libcds%2Fsrc%2Fstatic_sequence%2Fstatic_sequence.cpp;h=05b77538d6a3511b7fff21ff7d75827ad91f2dba;hb=6a9117a97674437c9fe00b099c40ccb61d45d704;hp=0ba5c2a15940a281f74dac73624369b897f22d38;hpb=450ba3c9c74665094fb8f6821d6cc92d2bf23011;p=SXSI%2FXMLTree.git

diff --git a/libcds/src/static_sequence/static_sequence.cpp b/libcds/src/static_sequence/static_sequence.cpp
index 0ba5c2a..05b7753 100644
--- a/libcds/src/static_sequence/static_sequence.cpp
+++ b/libcds/src/static_sequence/static_sequence.cpp
@@ -25,12 +25,58 @@ static_sequence::static_sequence() {}
 static_sequence::~static_sequence() {}
 uint static_sequence::length() { return len; }
 
+uint static_sequence::count(uint s) {
+  return rank(s,len-1);
+}
+
 static_sequence * static_sequence::load(FILE * fp) {
   uint rd;
   if(fread(&rd,sizeof(uint),1,fp)!=1) return NULL;
   fseek(fp,-sizeof(uint),SEEK_CUR);
   switch(rd) {
     case WVTREE_HDR: return static_sequence_wvtree::load(fp);
+    case GMR_CHUNK_HDR: return static_sequence_gmr_chunk::load(fp);
+    case GMR_HDR: return static_sequence_gmr::load(fp);
+    case WVTREE_NOPTRS_HDR: return static_sequence_wvtree_noptrs::load(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;
+}
+