00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040 #include "defs.h"
00041 #include "mempool.h"
00042 #include "cxx_memory.h"
00043 #include "errors.h"
00044
00045 #include "sparse_bv.h"
00046
00047 void
00048 SPARSE_BV::Set_Bit (UINT32 index)
00049 {
00050 if (buffer == NULL) {
00051 buffer = CXX_NEW (SPARSE_BV_BUFFER, Malloc_Mem_Pool);
00052 data = CXX_NEW (SPARSE_BV_DATA, Malloc_Mem_Pool);
00053 data->push_back (0LL);
00054 }
00055
00056 if (index >= buffer->size () * mUINT64_BIT) {
00057 INT n = ((index + mUINT64_BIT) / mUINT64_BIT) - buffer->size ();
00058 buffer->insert (buffer->end (), n, 0);
00059 }
00060
00061 UINT32 idx = index / mUINT64_BIT;
00062
00063 Is_True (idx < buffer->size (), ("Bit vector index out of range"));
00064
00065 if ((*buffer)[idx] == 0 ) {
00066 (*buffer)[idx] = data->size();
00067 data->push_back (1LL << (index % mUINT64_BIT));
00068 } else
00069 (*data)[(*buffer)[idx]] |= 1LL << (index % mUINT64_BIT);
00070
00071 #ifdef Is_True_On
00072
00073 #endif
00074 }
00075
00076
00077 void
00078 SPARSE_BV::Union_2_Core (const SPARSE_BV& vec, BUFFER_ITER vec_iter,
00079 BUFFER_ITER iter)
00080 {
00081 while (vec_iter != vec.buffer->end ()) {
00082 Is_True (iter != buffer->end (), ("bit vector buffer too short"));
00083
00084 if (*vec_iter != 0) {
00085 if (*iter == 0) {
00086
00087 *iter = data->size ();
00088 data->push_back ((*vec.data)[*vec_iter]);
00089 } else {
00090 (*data)[*iter] |= (*vec.data)[*vec_iter];
00091 }
00092 }
00093
00094 ++vec_iter;
00095 ++iter;
00096 }
00097
00098 #ifdef Is_True_On
00099
00100 #endif
00101 }
00102
00103
00104 void
00105 SPARSE_BV::Union_2 (const SPARSE_BV& vec)
00106 {
00107 if (vec.buffer == NULL)
00108 return;
00109
00110 if (buffer == NULL) {
00111 buffer = CXX_NEW (SPARSE_BV_BUFFER, Malloc_Mem_Pool);
00112 buffer->reserve (vec.buffer->size ());
00113 data = CXX_NEW (SPARSE_BV_DATA, Malloc_Mem_Pool);
00114 data->reserve (vec.data->size ());
00115 data->push_back (0LL);
00116 }
00117
00118
00119 if (vec.buffer->size () > buffer->size ()) {
00120 buffer->insert (buffer->end (),
00121 vec.buffer->size() - buffer->size (), 0);
00122 }
00123
00124 Union_2_Core (vec, vec.buffer->begin (), buffer->begin ());
00125
00126 #ifdef Is_True_On
00127
00128 #endif
00129
00130 }
00131
00132
00133 BOOL
00134 SPARSE_BV::Union_2_Diff_Core (const SPARSE_BV& vec)
00135 {
00136 BOOL modified = FALSE;
00137
00138
00139 if (vec.buffer->size () > buffer->size ()) {
00140 buffer->insert (buffer->end (),
00141 vec.buffer->size() - buffer->size (), 0);
00142 modified = TRUE;
00143 }
00144
00145 BUFFER_ITER vec_iter (vec.buffer->begin ());
00146 BUFFER_ITER iter (buffer->begin ());
00147
00148 if (!modified) {
00149
00150
00151 while (!modified && vec_iter != vec.buffer->end ()) {
00152 Is_True (iter != buffer->end (),
00153 ("Inconsisten bit vector comparison"));
00154 if (*vec_iter != 0) {
00155 if (*iter == 0) {
00156 modified = TRUE;
00157 break;
00158 }
00159 mUINT64& data_word = (*data)[*iter];
00160 mUINT64 xword = data_word | (*vec.data)[*vec_iter];
00161 if (xword != data_word) {
00162 modified = TRUE;
00163 data_word = xword;
00164 }
00165 }
00166
00167 ++vec_iter;
00168 ++iter;
00169 }
00170
00171 }
00172
00173 if (modified)
00174 Union_2_Core (vec, vec_iter, iter);
00175
00176 return modified;
00177
00178 }
00179
00180
00181 #ifdef Is_True_On
00182
00183 void
00184 SPARSE_BV::Verify () const
00185 {
00186 if (buffer == NULL)
00187 return;
00188
00189 Is_True (buffer->size () != 0, ("SPARSE_BV: empty map"));
00190 Is_True (buffer->back () != 0, ("SPARSE_BV: empty trailing slot"));
00191
00192 for (BUFFER_ITER first (buffer->begin ()); first != buffer->end ();
00193 ++first) {
00194
00195 if (*first != 0) {
00196 Is_True (*first < data->size (), ("SPARSE_BV: invalid index"));
00197 Is_True ((*data)[*first] != 0LL, ("SPARSE_BV: empty buffer"));
00198 }
00199 }
00200
00201 }
00202 #endif // Is_True_On