-
Notifications
You must be signed in to change notification settings - Fork 14
/
bitboard.cpp
176 lines (141 loc) · 3.69 KB
/
bitboard.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
// BitBoard.cpp: implementation of the BitBoard class.
//
//////////////////////////////////////////////////////////////////////
#include "bitboard.h"
#include <iostream>
using namespace std;
//ostream & operator<<(ostream& o, const BITBOARD& bb_data){
// BITBOARD bb=bb_data;
// if(bb){
// int nBit=EMPTY_ELEM;
// while(bb)
// {
// nBit=BitBoard::lsb64_de_Bruijn(bb);
// if(nBit==EMPTY_ELEM) break;
// o<<nBit<<" ";
// bb^=Tables::mask[nBit];
// }
// }
//
// o<<" ["<<BitBoard::popc64(bb_data)<<"]"<<endl;
//
//return o;
//}
int BitBoard::lsb64_lup (const BITBOARD bb){
U16 bb16;
if(bb){
bb16=(U16)bb;
if(bb16) return (Tables::lsba[0][bb16]);
bb16=(U16)(bb>>16);
if(bb16) return (Tables::lsba[1][bb16]);
bb16=(U16)(bb>>32);
if(bb16) return (Tables::lsba[2][bb16]);
bb16=(U16)(bb>>48);
if(bb16) return (Tables::lsba[3][bb16]);
}
return EMPTY_ELEM; //should not occur
}
int BitBoard::lsb64_pc(const BITBOARD bb_dato) {
/////////////////////
// isolani-1
// Comments: useful if you have anchitecturewith a fast popcount
if(bb_dato){
return BitBoard::popc64( (bb_dato & -bb_dato) - 1 );
}
return EMPTY_ELEM;
}
int BitBoard::lsb64_mod(const BITBOARD bb_dato){
//////////////////
// Modulus implementation of perfect hashing
// Use T_64[67] to convert the class of BITBOARDS 000...1...000in the position of the bit 1
// COMMENTS: Modulus % x reduces efficiency
if(bb_dato){
return(Tables::T_64[(bb_dato&(~bb_dato+1))%67]);
}
return EMPTY_ELEM;
}
int BitBoard::lsb64_lup_1(const BITBOARD bb_dato){
//////////////////////////
//Look up table implementation
//
// COMMENTS
// 1-Best implementation for 32 bits on x86 on average. Is is worse than LSB_32x only for sparse BB with
// 1bits in the final BB
union x {
U16 c[4];
BITBOARD b;
}val; /*register*/
val.b = bb_dato; //union load
//Control
if(bb_dato){
if(val.c[0])
return (Tables::lsb[val.c[0]]);
else if (val.c[1])
return (Tables::lsb[val.c[1]] + 16);
else if (val.c[2])
return (Tables::lsb[val.c[2]] + 32);
else return (Tables::lsb[val.c[3]] + 48); //nonzero by (1)
}
return EMPTY_ELEM; //Should not reach
}
//////////////////////////////
//
// PopCount
//
////////////////////////////
int BitBoard::popc64_lup_1(const BITBOARD bb_dato){
//uses 16 bit look up table
return (Tables::pc[(U16)bb_dato]+Tables::pc[(U16)(bb_dato>>16)]+Tables::pc[(U16)(bb_dato>>32)]+Tables::pc[(U16)(bb_dato>>48)]);
}
/////////////////////
//
// MSB_64 (...)
//
// No bit shuffling implememation seems close to LOOK UP for MSB(64 bits)
//
/////////////////////
int BitBoard::msb64_lup(const BITBOARD bb){
//////////////////////
// lookup table variant
union x {
U16 c[4];
BITBOARD b;
} val;
val.b = bb;
//if(bb==0) return -1; //for sparse data
if(val.b){
if(val.c[3]) return (Tables::msba[3][val.c[3]]);
if(val.c[2]) return (Tables::msba[2][val.c[2]]);
if(val.c[1]) return (Tables::msba[1][val.c[1]]);
if(val.c[0]) return (Tables::msba[0][val.c[0]]);
}
/* if (val.c[3]) //access table msb[65536](valores:0-15) in blocks of 16 bits
return (msb[val.c[3]] + 48);
else if (val.c[2])
return (msb[val.c[2]] + 32);
else if (val.c[1])
return (msb[val.c[1]] + 16);
else return (msb[val.c[0]]); //should be nonzero by (1)*/
return EMPTY_ELEM; //should not reach here
}
///////////////////////
//
// VIEW
//
///////////////////////
void BitBoard::print(const BITBOARD bb_data, std::ostream& o){
///////////////////////
// Shows bits and popcount
BITBOARD bb=bb_data;
if(bb){
int nBit=EMPTY_ELEM;
while(bb)
{
nBit=BitBoard::lsb64_de_Bruijn(bb);
if(nBit==EMPTY_ELEM) break;
o<<nBit<<" ";
bb^=Tables::mask[nBit];
}
}
o<<"["<<BitBoard::popc64(bb_data)<<"]"<<endl;
}