blob: c709064bf9ebae0a2e98f99ec1f509480cb1abf3 [file] [log] [blame]
using System;
using System.Diagnostics;
/*
* Licensed to the Apache Software Foundation (ASF) under one or more
* contributor license agreements. See the NOTICE file distributed with
* this work for additional information regarding copyright ownership.
* The ASF licenses this file to You under the Apache License, Version 2.0
* (the "License"); you may not use this file except in compliance with
* the License. You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
namespace Lucene.Net.Util
{
using Lucene.Net.Support;
using DocIdSet = Lucene.Net.Search.DocIdSet;
using DocIdSetIterator = Lucene.Net.Search.DocIdSetIterator;
/// <summary>
/// An "open" BitSet implementation that allows direct access to the array of words
/// storing the bits.
/// <p/>
/// Unlike java.util.bitset, the fact that bits are packed into an array of longs
/// is part of the interface. this allows efficient implementation of other algorithms
/// by someone other than the author. It also allows one to efficiently implement
/// alternate serialization or interchange formats.
/// <p/>
/// <code>OpenBitSet</code> is faster than <code>java.util.BitSet</code> in most operations
/// and *much* faster at calculating cardinality of sets and results of set operations.
/// It can also handle sets of larger cardinality (up to 64 * 2**32-1)
/// <p/>
/// The goals of <code>OpenBitSet</code> are the fastest implementation possible, and
/// maximum code reuse. Extra safety and encapsulation
/// may always be built on top, but if that's built in, the cost can never be removed (and
/// hence people re-implement their own version in order to get better performance).
/// If you want a "safe", totally encapsulated (and slower and limited) BitSet
/// class, use <code>java.util.BitSet</code>.
/// <p/>
/// <h3>Performance Results</h3>
///
/// Test system: Pentium 4, Sun Java 1.5_06 -server -Xbatch -Xmx64M
/// <br/>BitSet size = 1,000,000
/// <br/>Results are java.util.BitSet time divided by OpenBitSet time.
/// <table border="1">
/// <tr>
/// <th></th> <th>cardinality</th> <th>intersect_count</th> <th>union</th> <th>nextSetBit</th> <th>get</th> <th>iterator</th>
/// </tr>
/// <tr>
/// <th>50% full</th> <td>3.36</td> <td>3.96</td> <td>1.44</td> <td>1.46</td> <td>1.99</td> <td>1.58</td>
/// </tr>
/// <tr>
/// <th>1% full</th> <td>3.31</td> <td>3.90</td> <td>&nbsp;</td> <td>1.04</td> <td>&nbsp;</td> <td>0.99</td>
/// </tr>
/// </table>
/// <br/>
/// Test system: AMD Opteron, 64 bit linux, Sun Java 1.5_06 -server -Xbatch -Xmx64M
/// <br/>BitSet size = 1,000,000
/// <br/>Results are java.util.BitSet time divided by OpenBitSet time.
/// <table border="1">
/// <tr>
/// <th></th> <th>cardinality</th> <th>intersect_count</th> <th>union</th> <th>nextSetBit</th> <th>get</th> <th>iterator</th>
/// </tr>
/// <tr>
/// <th>50% full</th> <td>2.50</td> <td>3.50</td> <td>1.00</td> <td>1.03</td> <td>1.12</td> <td>1.25</td>
/// </tr>
/// <tr>
/// <th>1% full</th> <td>2.51</td> <td>3.49</td> <td>&nbsp;</td> <td>1.00</td> <td>&nbsp;</td> <td>1.02</td>
/// </tr>
/// </table>
/// </summary>
public class OpenBitSet : DocIdSet, Bits, ICloneable
{
protected internal long[] bits;
protected internal int Wlen; // number of words (elements) used in the array
// Used only for assert:
private long NumBits;
/// <summary>
/// Constructs an OpenBitSet large enough to hold {@code numBits}. </summary>
public OpenBitSet(long numBits)
{
this.NumBits = numBits;
bits = new long[Bits2words(numBits)];
Wlen = bits.Length;
}
/// <summary>
/// Constructor: allocates enough space for 64 bits. </summary>
public OpenBitSet()
: this(64)
{
}
/// <summary>
/// Constructs an OpenBitSet from an existing long[].
/// <p>
/// The first 64 bits are in long[0], with bit index 0 at the least significant
/// bit, and bit index 63 at the most significant. Given a bit index, the word
/// containing it is long[index/64], and it is at bit number index%64 within
/// that word.
/// <p>
/// numWords are the number of elements in the array that contain set bits
/// (non-zero longs). numWords should be &lt= bits.length, and any existing
/// words in the array at position &gt= numWords should be zero.
///
/// </summary>
public OpenBitSet(long[] bits, int numWords)
{
if (numWords > bits.Length)
{
throw new System.ArgumentException("numWords cannot exceed bits.length");
}
this.bits = bits;
this.Wlen = numWords;
this.NumBits = Wlen * 64;
}
public override DocIdSetIterator GetIterator()
{
return new OpenBitSetIterator(bits, Wlen);
}
public override Bits GetBits()
{
return this;
}
/// <summary>
/// this DocIdSet implementation is cacheable. </summary>
public override bool Cacheable
{
get
{
return true;
}
}
/// <summary>
/// Returns the current capacity in bits (1 greater than the index of the last bit) </summary>
public virtual long Capacity()
{
return bits.Length << 6;
}
/// <summary>
/// Returns the current capacity of this set. Included for
/// compatibility. this is *not* equal to <seealso cref="#cardinality"/>
/// </summary>
public virtual long Size()
{
return Capacity();
}
public virtual int Length()
{
return bits.Length << 6;
}
/// <summary>
/// Returns true if there are no set bits </summary>
public virtual bool Empty
{
get
{
return Cardinality() == 0;
}
}
/// <summary>
/// Expert: returns the long[] storing the bits </summary>
public virtual long[] Bits
{
get
{
return bits;
}
}
/// <summary>
/// Expert: gets the number of longs in the array that are in use </summary>
public virtual int NumWords
{
get
{
return Wlen;
}
}
/// <summary>
/// Returns true or false for the specified bit index. </summary>
public virtual bool Get(int index)
{
int i = index >> 6; // div 64
// signed shift will keep a negative index and force an
// array-index-out-of-bounds-exception, removing the need for an explicit check.
if (i >= bits.Length)
{
return false;
}
int bit = index & 0x3f; // mod 64
long bitmask = 1L << bit;
return (bits[i] & bitmask) != 0;
}
/// <summary>
/// Returns true or false for the specified bit index.
/// The index should be less than the OpenBitSet size
/// </summary>
public virtual bool FastGet(int index)
{
Debug.Assert(index >= 0 && index < NumBits);
int i = index >> 6; // div 64
// signed shift will keep a negative index and force an
// array-index-out-of-bounds-exception, removing the need for an explicit check.
int bit = index & 0x3f; // mod 64
long bitmask = 1L << bit;
return (bits[i] & bitmask) != 0;
}
/// <summary>
/// Returns true or false for the specified bit index
/// </summary>
public virtual bool Get(long index)
{
int i = (int)(index >> 6); // div 64
if (i >= bits.Length)
{
return false;
}
int bit = (int)index & 0x3f; // mod 64
long bitmask = 1L << bit;
return (bits[i] & bitmask) != 0;
}
/// <summary>
/// Returns true or false for the specified bit index.
/// The index should be less than the OpenBitSet size.
/// </summary>
public virtual bool FastGet(long index)
{
Debug.Assert(index >= 0 && index < NumBits);
int i = (int)(index >> 6); // div 64
int bit = (int)index & 0x3f; // mod 64
long bitmask = 1L << bit;
return (bits[i] & bitmask) != 0;
}
/*
// alternate implementation of get()
public boolean get1(int index) {
int i = index >> 6; // div 64
int bit = index & 0x3f; // mod 64
return ((bits[i]>>>bit) & 0x01) != 0;
// this does a long shift and a bittest (on x86) vs
// a long shift, and a long AND, (the test for zero is prob a no-op)
// testing on a P4 indicates this is slower than (bits[i] & bitmask) != 0;
}
*/
/// <summary>
/// returns 1 if the bit is set, 0 if not.
/// The index should be less than the OpenBitSet size
/// </summary>
public virtual int GetBit(int index)
{
Debug.Assert(index >= 0 && index < NumBits);
int i = index >> 6; // div 64
int bit = index & 0x3f; // mod 64
return ((int)((long)((ulong)bits[i] >> bit))) & 0x01;
}
/*
public boolean get2(int index) {
int word = index >> 6; // div 64
int bit = index & 0x0000003f; // mod 64
return (bits[word] << bit) < 0; // hmmm, this would work if bit order were reversed
// we could right shift and check for parity bit, if it was available to us.
}
*/
/// <summary>
/// sets a bit, expanding the set size if necessary </summary>
public virtual void Set(long index)
{
int wordNum = ExpandingWordNum(index);
int bit = (int)index & 0x3f;
long bitmask = 1L << bit;
bits[wordNum] |= bitmask;
}
/// <summary>
/// Sets the bit at the specified index.
/// The index should be less than the OpenBitSet size.
/// </summary>
public virtual void FastSet(int index)
{
Debug.Assert(index >= 0 && index < NumBits);
int wordNum = index >> 6; // div 64
int bit = index & 0x3f; // mod 64
long bitmask = 1L << bit;
bits[wordNum] |= bitmask;
}
/// <summary>
/// Sets the bit at the specified index.
/// The index should be less than the OpenBitSet size.
/// </summary>
public virtual void FastSet(long index)
{
Debug.Assert(index >= 0 && index < NumBits);
int wordNum = (int)(index >> 6);
int bit = (int)index & 0x3f;
long bitmask = 1L << bit;
bits[wordNum] |= bitmask;
}
/// <summary>
/// Sets a range of bits, expanding the set size if necessary
/// </summary>
/// <param name="startIndex"> lower index </param>
/// <param name="endIndex"> one-past the last bit to set </param>
public virtual void Set(long startIndex, long endIndex)
{
if (endIndex <= startIndex)
{
return;
}
int startWord = (int)(startIndex >> 6);
// since endIndex is one past the end, this is index of the last
// word to be changed.
int endWord = ExpandingWordNum(endIndex - 1);
long startmask = -1L << (int)startIndex;
long endmask = (long)(0xffffffffffffffffUL >> (int)-endIndex); // 64-(endIndex&0x3f) is the same as -endIndex due to wrap
if (startWord == endWord)
{
bits[startWord] |= (startmask & endmask);
return;
}
bits[startWord] |= startmask;
Arrays.Fill(bits, startWord + 1, endWord, -1L);
bits[endWord] |= endmask;
}
protected internal virtual int ExpandingWordNum(long index)
{
int wordNum = (int)(index >> 6);
if (wordNum >= Wlen)
{
EnsureCapacity(index + 1);
}
return wordNum;
}
/// <summary>
/// clears a bit.
/// The index should be less than the OpenBitSet size.
/// </summary>
public virtual void FastClear(int index)
{
Debug.Assert(index >= 0 && index < NumBits);
int wordNum = index >> 6;
int bit = index & 0x03f;
long bitmask = 1L << bit;
bits[wordNum] &= ~bitmask;
// hmmm, it takes one more instruction to clear than it does to set... any
// way to work around this? If there were only 63 bits per word, we could
// use a right shift of 10111111...111 in binary to position the 0 in the
// correct place (using sign extension).
// Could also use Long.rotateRight() or rotateLeft() *if* they were converted
// by the JVM into a native instruction.
// bits[word] &= Long.rotateLeft(0xfffffffe,bit);
}
/// <summary>
/// clears a bit.
/// The index should be less than the OpenBitSet size.
/// </summary>
public virtual void FastClear(long index)
{
Debug.Assert(index >= 0 && index < NumBits);
int wordNum = (int)(index >> 6); // div 64
int bit = (int)index & 0x3f; // mod 64
long bitmask = 1L << bit;
bits[wordNum] &= ~bitmask;
}
/// <summary>
/// clears a bit, allowing access beyond the current set size without changing the size. </summary>
public virtual void Clear(long index)
{
int wordNum = (int)(index >> 6); // div 64
if (wordNum >= Wlen)
{
return;
}
int bit = (int)index & 0x3f; // mod 64
long bitmask = 1L << bit;
bits[wordNum] &= ~bitmask;
}
/// <summary>
/// Clears a range of bits. Clearing past the end does not change the size of the set.
/// </summary>
/// <param name="startIndex"> lower index </param>
/// <param name="endIndex"> one-past the last bit to clear </param>
public virtual void Clear(int startIndex, int endIndex)
{
if (endIndex <= startIndex)
{
return;
}
int startWord = (startIndex >> 6);
if (startWord >= Wlen)
{
return;
}
// since endIndex is one past the end, this is index of the last
// word to be changed.
int endWord = ((endIndex - 1) >> 6);
long startmask = (-1L) << startIndex; // -1 << (startIndex mod 64)
long endmask = (-1L) << endIndex; // -1 << (endIndex mod 64)
if ((endIndex & 0x3f) == 0)
{
endmask = 0;
}
startmask = ~startmask;
if (startWord == endWord)
{
bits[startWord] &= (startmask | endmask);
return;
}
bits[startWord] &= startmask;
int middle = Math.Min(Wlen, endWord);
Arrays.Fill(bits, startWord + 1, middle, 0L);
if (endWord < Wlen)
{
bits[endWord] &= endmask;
}
}
/// <summary>
/// Clears a range of bits. Clearing past the end does not change the size of the set.
/// </summary>
/// <param name="startIndex"> lower index </param>
/// <param name="endIndex"> one-past the last bit to clear </param>
public virtual void Clear(long startIndex, long endIndex)
{
if (endIndex <= startIndex)
{
return;
}
int startWord = (int)(startIndex >> 6);
if (startWord >= Wlen)
{
return;
}
// since endIndex is one past the end, this is index of the last
// word to be changed.
int endWord = (int)((endIndex - 1) >> 6);
//LUCENE TO-DO
long startmask = -1L << (int)startIndex;
long endmask = -(int)((uint)1L >> (int)-endIndex); // 64-(endIndex&0x3f) is the same as -endIndex due to wrap
// invert masks since we are clearing
startmask = ~startmask;
endmask = ~endmask;
if (startWord == endWord)
{
bits[startWord] &= (startmask | endmask);
return;
}
bits[startWord] &= startmask;
int middle = Math.Min(Wlen, endWord);
Arrays.Fill(bits, startWord + 1, middle, 0L);
if (endWord < Wlen)
{
bits[endWord] &= endmask;
}
}
/// <summary>
/// Sets a bit and returns the previous value.
/// The index should be less than the OpenBitSet size.
/// </summary>
public virtual bool GetAndSet(int index)
{
Debug.Assert(index >= 0 && index < NumBits);
int wordNum = index >> 6; // div 64
int bit = index & 0x3f; // mod 64
long bitmask = 1L << bit;
bool val = (bits[wordNum] & bitmask) != 0;
bits[wordNum] |= bitmask;
return val;
}
/// <summary>
/// Sets a bit and returns the previous value.
/// The index should be less than the OpenBitSet size.
/// </summary>
public virtual bool GetAndSet(long index)
{
Debug.Assert(index >= 0 && index < NumBits);
int wordNum = (int)(index >> 6); // div 64
int bit = (int)index & 0x3f; // mod 64
long bitmask = 1L << bit;
bool val = (bits[wordNum] & bitmask) != 0;
bits[wordNum] |= bitmask;
return val;
}
/// <summary>
/// flips a bit.
/// The index should be less than the OpenBitSet size.
/// </summary>
public virtual void FastFlip(int index)
{
Debug.Assert(index >= 0 && index < NumBits);
int wordNum = index >> 6; // div 64
int bit = index & 0x3f; // mod 64
long bitmask = 1L << bit;
bits[wordNum] ^= bitmask;
}
/// <summary>
/// flips a bit.
/// The index should be less than the OpenBitSet size.
/// </summary>
public virtual void FastFlip(long index)
{
Debug.Assert(index >= 0 && index < NumBits);
int wordNum = (int)(index >> 6); // div 64
int bit = (int)index & 0x3f; // mod 64
long bitmask = 1L << bit;
bits[wordNum] ^= bitmask;
}
/// <summary>
/// flips a bit, expanding the set size if necessary </summary>
public virtual void Flip(long index)
{
int wordNum = ExpandingWordNum(index);
int bit = (int)index & 0x3f; // mod 64
long bitmask = 1L << bit;
bits[wordNum] ^= bitmask;
}
/// <summary>
/// flips a bit and returns the resulting bit value.
/// The index should be less than the OpenBitSet size.
/// </summary>
public virtual bool FlipAndGet(int index)
{
Debug.Assert(index >= 0 && index < NumBits);
int wordNum = index >> 6; // div 64
int bit = index & 0x3f; // mod 64
long bitmask = 1L << bit;
bits[wordNum] ^= bitmask;
return (bits[wordNum] & bitmask) != 0;
}
/// <summary>
/// flips a bit and returns the resulting bit value.
/// The index should be less than the OpenBitSet size.
/// </summary>
public virtual bool FlipAndGet(long index)
{
Debug.Assert(index >= 0 && index < NumBits);
int wordNum = (int)(index >> 6); // div 64
int bit = (int)index & 0x3f; // mod 64
long bitmask = 1L << bit;
bits[wordNum] ^= bitmask;
return (bits[wordNum] & bitmask) != 0;
}
/// <summary>
/// Flips a range of bits, expanding the set size if necessary
/// </summary>
/// <param name="startIndex"> lower index </param>
/// <param name="endIndex"> one-past the last bit to flip </param>
public virtual void Flip(long startIndex, long endIndex)
{
if (endIndex <= startIndex)
{
return;
}
int startWord = (int)(startIndex >> 6);
// since endIndex is one past the end, this is index of the last
// word to be changed.
int endWord = ExpandingWordNum(endIndex - 1);
/// <summary>
///* Grrr, java shifting wraps around so -1L>>>64 == -1
/// for that reason, make sure not to use endmask if the bits to flip will
/// be zero in the last word (redefine endWord to be the last changed...)
/// long startmask = -1L << (startIndex & 0x3f); // example: 11111...111000
/// long endmask = -1L >>> (64-(endIndex & 0x3f)); // example: 00111...111111
/// **
/// </summary>
//LUCENE TO-DO
long startmask = -1L << (int)startIndex;
long endmask = (long)(0xffffffffffffffffUL >> (int)-endIndex); // 64-(endIndex&0x3f) is the same as -endIndex due to wrap
if (startWord == endWord)
{
bits[startWord] ^= (startmask & endmask);
return;
}
bits[startWord] ^= startmask;
for (int i = startWord + 1; i < endWord; i++)
{
bits[i] = ~bits[i];
}
bits[endWord] ^= endmask;
}
/*
public static int pop(long v0, long v1, long v2, long v3) {
// derived from pop_array by setting last four elems to 0.
// exchanges one pop() call for 10 elementary operations
// saving about 7 instructions... is there a better way?
long twosA=v0 & v1;
long ones=v0^v1;
long u2=ones^v2;
long twosB =(ones&v2)|(u2&v3);
ones=u2^v3;
long fours=(twosA&twosB);
long twos=twosA^twosB;
return (pop(fours)<<2)
+ (pop(twos)<<1)
+ pop(ones);
}
*/
/// <returns> the number of set bits </returns>
public virtual long Cardinality()
{
return BitUtil.Pop_array(bits, 0, Wlen);
}
/// <summary>
/// Returns the popcount or cardinality of the intersection of the two sets.
/// Neither set is modified.
/// </summary>
public static long IntersectionCount(OpenBitSet a, OpenBitSet b)
{
return BitUtil.Pop_intersect(a.bits, b.bits, 0, Math.Min(a.Wlen, b.Wlen));
}
/// <summary>
/// Returns the popcount or cardinality of the union of the two sets.
/// Neither set is modified.
/// </summary>
public static long UnionCount(OpenBitSet a, OpenBitSet b)
{
long tot = BitUtil.Pop_union(a.bits, b.bits, 0, Math.Min(a.Wlen, b.Wlen));
if (a.Wlen < b.Wlen)
{
tot += BitUtil.Pop_array(b.bits, a.Wlen, b.Wlen - a.Wlen);
}
else if (a.Wlen > b.Wlen)
{
tot += BitUtil.Pop_array(a.bits, b.Wlen, a.Wlen - b.Wlen);
}
return tot;
}
/// <summary>
/// Returns the popcount or cardinality of "a and not b"
/// or "intersection(a, not(b))".
/// Neither set is modified.
/// </summary>
public static long AndNotCount(OpenBitSet a, OpenBitSet b)
{
long tot = BitUtil.Pop_andnot(a.bits, b.bits, 0, Math.Min(a.Wlen, b.Wlen));
if (a.Wlen > b.Wlen)
{
tot += BitUtil.Pop_array(a.bits, b.Wlen, a.Wlen - b.Wlen);
}
return tot;
}
/// <summary>
/// Returns the popcount or cardinality of the exclusive-or of the two sets.
/// Neither set is modified.
/// </summary>
public static long XorCount(OpenBitSet a, OpenBitSet b)
{
long tot = BitUtil.Pop_xor(a.bits, b.bits, 0, Math.Min(a.Wlen, b.Wlen));
if (a.Wlen < b.Wlen)
{
tot += BitUtil.Pop_array(b.bits, a.Wlen, b.Wlen - a.Wlen);
}
else if (a.Wlen > b.Wlen)
{
tot += BitUtil.Pop_array(a.bits, b.Wlen, a.Wlen - b.Wlen);
}
return tot;
}
/// <summary>
/// Returns the index of the first set bit starting at the index specified.
/// -1 is returned if there are no more set bits.
/// </summary>
public virtual int NextSetBit(int index)
{
int i = index >> 6;
if (i >= Wlen)
{
return -1;
}
int subIndex = index & 0x3f; // index within the word
long word = bits[i] >> subIndex; // skip all the bits to the right of index
if (word != 0)
{
return (i << 6) + subIndex + Number.NumberOfTrailingZeros(word);
}
while (++i < Wlen)
{
word = bits[i];
if (word != 0)
{
return (i << 6) + Number.NumberOfTrailingZeros(word);
}
}
return -1;
}
/// <summary>
/// Returns the index of the first set bit starting at the index specified.
/// -1 is returned if there are no more set bits.
/// </summary>
public virtual long NextSetBit(long index)
{
int i = (int)((long)((ulong)index >> 6));
if (i >= Wlen)
{
return -1;
}
int subIndex = (int)index & 0x3f; // index within the word
long word = (long)((ulong)bits[i] >> subIndex); // skip all the bits to the right of index
if (word != 0)
{
return (((long)i) << 6) + (subIndex + Number.NumberOfTrailingZeros(word));
}
while (++i < Wlen)
{
word = bits[i];
if (word != 0)
{
return (((long)i) << 6) + Number.NumberOfTrailingZeros(word);
}
}
return -1;
}
/// <summary>
/// Returns the index of the first set bit starting downwards at
/// the index specified.
/// -1 is returned if there are no more set bits.
/// </summary>
public virtual int PrevSetBit(int index)
{
int i = index >> 6;
int subIndex;
long word;
if (i >= Wlen)
{
i = Wlen - 1;
if (i < 0)
{
return -1;
}
subIndex = 63; // last possible bit
word = bits[i];
}
else
{
if (i < 0)
{
return -1;
}
subIndex = index & 0x3f; // index within the word
word = (bits[i] << (63 - subIndex)); // skip all the bits to the left of index
}
if (word != 0)
{
return (i << 6) + subIndex - Number.NumberOfLeadingZeros(word); // See LUCENE-3197
}
while (--i >= 0)
{
word = bits[i];
if (word != 0)
{
return (i << 6) + 63 - Number.NumberOfLeadingZeros(word);
}
}
return -1;
}
/// <summary>
/// Returns the index of the first set bit starting downwards at
/// the index specified.
/// -1 is returned if there are no more set bits.
/// </summary>
public virtual long PrevSetBit(long index)
{
int i = (int)(index >> 6);
int subIndex;
long word;
if (i >= Wlen)
{
i = Wlen - 1;
if (i < 0)
{
return -1;
}
subIndex = 63; // last possible bit
word = bits[i];
}
else
{
if (i < 0)
{
return -1;
}
subIndex = (int)index & 0x3f; // index within the word
word = (bits[i] << (63 - subIndex)); // skip all the bits to the left of index
}
if (word != 0)
{
return (((long)i) << 6) + subIndex - Number.NumberOfLeadingZeros(word); // See LUCENE-3197
}
while (--i >= 0)
{
word = bits[i];
if (word != 0)
{
return (((long)i) << 6) + 63 - Number.NumberOfLeadingZeros(word);
}
}
return -1;
}
public object Clone()
{
try
{
//OpenBitSet obs = (OpenBitSet)base.Clone();
//obs.bits = (long[])obs.bits.Clone(); // hopefully an array clone is as fast(er) than arraycopy
OpenBitSet obs = new OpenBitSet((long[])bits.Clone(), Wlen);
return obs;
}
catch (Exception e)
{
throw new SystemException(e.Message, e);
}
}
/// <summary>
/// this = this AND other </summary>
public virtual void Intersect(OpenBitSet other)
{
int newLen = Math.Min(this.Wlen, other.Wlen);
long[] thisArr = this.bits;
long[] otherArr = other.bits;
// testing against zero can be more efficient
int pos = newLen;
while (--pos >= 0)
{
thisArr[pos] &= otherArr[pos];
}
if (this.Wlen > newLen)
{
// fill zeros from the new shorter length to the old length
Arrays.Fill(bits, newLen, this.Wlen, 0);
}
this.Wlen = newLen;
}
/// <summary>
/// this = this OR other </summary>
public virtual void Union(OpenBitSet other)
{
int newLen = Math.Max(Wlen, other.Wlen);
EnsureCapacityWords(newLen);
Debug.Assert((NumBits = Math.Max(other.NumBits, NumBits)) >= 0);
long[] thisArr = this.bits;
long[] otherArr = other.bits;
int pos = Math.Min(Wlen, other.Wlen);
while (--pos >= 0)
{
thisArr[pos] |= otherArr[pos];
}
if (this.Wlen < newLen)
{
Array.Copy(otherArr, this.Wlen, thisArr, this.Wlen, newLen - this.Wlen);
}
this.Wlen = newLen;
}
/// <summary>
/// Remove all elements set in other. this = this AND_NOT other </summary>
public virtual void Remove(OpenBitSet other)
{
int idx = Math.Min(Wlen, other.Wlen);
long[] thisArr = this.bits;
long[] otherArr = other.bits;
while (--idx >= 0)
{
thisArr[idx] &= ~otherArr[idx];
}
}
/// <summary>
/// this = this XOR other </summary>
public virtual void Xor(OpenBitSet other)
{
int newLen = Math.Max(Wlen, other.Wlen);
EnsureCapacityWords(newLen);
Debug.Assert((NumBits = Math.Max(other.NumBits, NumBits)) >= 0);
long[] thisArr = this.bits;
long[] otherArr = other.bits;
int pos = Math.Min(Wlen, other.Wlen);
while (--pos >= 0)
{
thisArr[pos] ^= otherArr[pos];
}
if (this.Wlen < newLen)
{
Array.Copy(otherArr, this.Wlen, thisArr, this.Wlen, newLen - this.Wlen);
}
this.Wlen = newLen;
}
// some BitSet compatability methods
//** see {@link intersect} */
public virtual void And(OpenBitSet other)
{
Intersect(other);
}
//** see {@link union} */
public virtual void Or(OpenBitSet other)
{
Union(other);
}
//** see {@link andNot} */
public virtual void AndNot(OpenBitSet other)
{
Remove(other);
}
/// <summary>
/// returns true if the sets have any elements in common </summary>
public virtual bool Intersects(OpenBitSet other)
{
int pos = Math.Min(this.Wlen, other.Wlen);
long[] thisArr = this.bits;
long[] otherArr = other.bits;
while (--pos >= 0)
{
if ((thisArr[pos] & otherArr[pos]) != 0)
{
return true;
}
}
return false;
}
/// <summary>
/// Expand the long[] with the size given as a number of words (64 bit longs). </summary>
public virtual void EnsureCapacityWords(int numWords)
{
bits = ArrayUtil.Grow(bits, numWords);
Wlen = numWords;
Debug.Assert((this.NumBits = Math.Max(this.NumBits, numWords << 6)) >= 0);
}
/// <summary>
/// Ensure that the long[] is big enough to hold numBits, expanding it if
/// necessary.
/// </summary>
public virtual void EnsureCapacity(long numBits)
{
EnsureCapacityWords(Bits2words(numBits));
// ensureCapacityWords sets numBits to a multiple of 64, but we want to set
// it to exactly what the app asked.
Debug.Assert((this.NumBits = Math.Max(this.NumBits, numBits)) >= 0);
}
/// <summary>
/// Lowers numWords, the number of words in use,
/// by checking for trailing zero words.
/// </summary>
public virtual void TrimTrailingZeros()
{
int idx = Wlen - 1;
while (idx >= 0 && bits[idx] == 0)
{
idx--;
}
Wlen = idx + 1;
}
/// <summary>
/// returns the number of 64 bit words it would take to hold numBits </summary>
public static int Bits2words(long numBits)
{
return (int)(((numBits - 1) >> 6) + 1);
}
/// <summary>
/// returns true if both sets have the same bits set </summary>
public override bool Equals(object o)
{
if (this == o)
{
return true;
}
if (!(o is OpenBitSet))
{
return false;
}
OpenBitSet a;
OpenBitSet b = (OpenBitSet)o;
// make a the larger set.
if (b.Wlen > this.Wlen)
{
a = b;
b = this;
}
else
{
a = this;
}
// check for any set bits out of the range of b
for (int i = a.Wlen - 1; i >= b.Wlen; i--)
{
if (a.bits[i] != 0)
{
return false;
}
}
for (int i = b.Wlen - 1; i >= 0; i--)
{
if (a.bits[i] != b.bits[i])
{
return false;
}
}
return true;
}
public override int GetHashCode()
{
// Start with a zero hash and use a mix that results in zero if the input is zero.
// this effectively truncates trailing zeros without an explicit check.
long h = 0;
for (int i = bits.Length; --i >= 0; )
{
h ^= bits[i];
h = (h << 1) | ((long)((ulong)h >> 63)); // rotate left
}
// fold leftmost bits into right and add a constant to prevent
// empty sets from returning 0, which is too common.
return (int)((h >> 32) ^ h) + unchecked((int)0x98761234);
}
}
}