#include "openeye.h"
#include <iomanip>
#include <iostream>
#include <string>
#include "oesystem.h"

using namespace std;
using namespace OESystem;

string bitstring(const OEBitVector bv)
{
  int i;
  string s("");
  for (i=0;i<bv.GetSize();++i)
  {
    if (bv.IsBitOn(i)) s+='1';
    else               s+='0';
  }
  return s;
}

int main()
{
  const char *hex1="408713C1E87C000BD18187000F70000000FFB807";
  const char *hex2="78EF93C1EE7FFF48018000005FF840BD9C0F9807";
  string bitstr="00000000111111110000111100000001";
  string hexstr;
  OEBitVector bv1,bv2;
  unsigned int i,count;
  float t;

  bv1.FromHexString(hex1);
  cout << "bv1 size=" << bv1.GetSize();
  cout << "...bitcount=" << bv1.CountBits();
  count = 0;
  for (i=0; i<bv1.GetSize(); ++i)
    if (bv1.IsBitOn(i))
      ++count;
  cout << "...recount=" << count << endl;

  bv2.FromHexString(hex2);
  cout << "bv2 size=" << bv2.GetSize();
  cout << "...bitcount=" << bv2.CountBits();
  count = 0;
  for (i=0; i<bv2.GetSize(); ++i)
    if (bv2.IsBitOn(i))
      ++count;
  cout << "...recount=" << count << endl;

  std::string str;
  bv1.ToHexString(str);
  cout << "FP1 is " << str << endl;
  cout << "orig=  " << hex1 << endl;
  bv2.ToHexString(str);
  cout << "FP2 is " << str << endl;
  cout << "orig=  " << hex2 << endl;

  t = OETanimoto(bv1,bv2);

  cout << "OETanimoto(bv1,bv2)=" << setprecision(3) << t << endl;

  bv1.SetSize(bitstr.size());

  for (i=0;i<bv1.GetSize();++i) {
    if (bitstr[i]=='1') bv1.SetBitOn(i);
  }
  bv1.ToHexString(hexstr);
  cout << bitstr << "->" << hexstr << "->" << endl;

  bv2.FromHexString(hexstr);
  cout << bitstring(bv2) << endl;

  hex1="FFFF0000";
  bv1.FromHexString(hex1);
  cout << "bv1 = " << hex1 << " = " << bitstring(bv1) << endl;
  cout << "bv1[0]=" << bv1[0]?1:0;
  cout << " bv1[1]=" << bv1[1]?1:0;
  cout << " bv1[15]=" << bv1[15]?1:0;
  cout << " bv1[16]=" << bv1[16]?1:0;
  cout << " bv1[31]=" << bv1[31]?1:0;
  cout << " bv1[32]=" << (bv1[32]?1:0) << endl;

  OEBitVector bv3;
  bv3 = ~bv1;
  bv3.ToHexString(hexstr);
  cout << "~bv1 = " << hexstr << " = " << bitstring(bv3) << endl;

  hex2="00FF0011";
  bv2.FromHexString(hex2);
  cout << "bv2 = " << hex2 << " = " << bitstring(bv2) << endl;

  bv3 = bv1 & bv2;
  bv3.ToHexString(hexstr);
  cout << "bv1 & bv2 = " << hexstr << " = " << bitstring(bv3) << endl;

  bv3 = bv1 | bv2;
  bv3.ToHexString(hexstr);
  cout << "bv1 | bv2 = " << hexstr << " = " << bitstring(bv3) << endl;

  bv3 = OEBitVectorXor(bv1,bv2);
  bv3.ToHexString(hexstr);
  cout << "OEBitVectorXor(bv1,bv2) = " << hexstr << " = " << bitstring(bv3) << endl;

  bv3 = OEBitVectorSub(bv1,bv2);
  bv3.ToHexString(hexstr);
  cout << "OEBitVectorSub(bv1,bv2) = " << hexstr << " = " << bitstring(bv3) << endl;

  // one way to tell if bv1 is subset of bv2
  bv1.FromHexString("00F00000");
  cout << "bv1 = " << bitstring(bv1) << endl;
  bool ok = OEBitVectorSub(bv1,bv2).IsEmpty();
  cout << "OEBitVectorSub(bv1,bv2).IsEmpty()=" << ok;
  cout << " (bv1 subset of bv2)" << endl;

  return 0;
}