13cf59944a961c75f3b30afc2e7fd7c95eb4dbf2
[stockfish] / src / position.cpp
1 /*
2   Stockfish, a UCI chess playing engine derived from Glaurung 2.1
3   Copyright (C) 2004-2008 Tord Romstad (Glaurung author)
4   Copyright (C) 2008-2010 Marco Costalba, Joona Kiiski, Tord Romstad
5
6   Stockfish is free software: you can redistribute it and/or modify
7   it under the terms of the GNU General Public License as published by
8   the Free Software Foundation, either version 3 of the License, or
9   (at your option) any later version.
10
11   Stockfish is distributed in the hope that it will be useful,
12   but WITHOUT ANY WARRANTY; without even the implied warranty of
13   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14   GNU General Public License for more details.
15
16   You should have received a copy of the GNU General Public License
17   along with this program.  If not, see <http://www.gnu.org/licenses/>.
18 */
19
20 #include <cassert>
21 #include <cstring>
22 #include <fstream>
23 #include <iostream>
24 #include <sstream>
25
26 #include "bitcount.h"
27 #include "movegen.h"
28 #include "position.h"
29 #include "psqtab.h"
30 #include "rkiss.h"
31 #include "thread.h"
32 #include "tt.h"
33
34 using std::string;
35 using std::cout;
36 using std::endl;
37
38 Key Position::zobrist[2][8][64];
39 Key Position::zobEp[64];
40 Key Position::zobCastle[16];
41 Key Position::zobSideToMove;
42 Key Position::zobExclusion;
43
44 Score Position::pieceSquareTable[16][64];
45
46 // Material values arrays, indexed by Piece
47 const Value PieceValueMidgame[17] = {
48   VALUE_ZERO,
49   PawnValueMidgame, KnightValueMidgame, BishopValueMidgame,
50   RookValueMidgame, QueenValueMidgame,
51   VALUE_ZERO, VALUE_ZERO, VALUE_ZERO,
52   PawnValueMidgame, KnightValueMidgame, BishopValueMidgame,
53   RookValueMidgame, QueenValueMidgame
54 };
55
56 const Value PieceValueEndgame[17] = {
57   VALUE_ZERO,
58   PawnValueEndgame, KnightValueEndgame, BishopValueEndgame,
59   RookValueEndgame, QueenValueEndgame,
60   VALUE_ZERO, VALUE_ZERO, VALUE_ZERO,
61   PawnValueEndgame, KnightValueEndgame, BishopValueEndgame,
62   RookValueEndgame, QueenValueEndgame
63 };
64
65
66 namespace {
67
68   // Bonus for having the side to move (modified by Joona Kiiski)
69   const Score TempoValue = make_score(48, 22);
70
71   // To convert a Piece to and from a FEN char
72   const string PieceToChar(".PNBRQK  pnbrqk  ");
73 }
74
75
76 /// CheckInfo c'tor
77
78 CheckInfo::CheckInfo(const Position& pos) {
79
80   Color them = flip(pos.side_to_move());
81   Square ksq = pos.king_square(them);
82
83   pinned = pos.pinned_pieces();
84   dcCandidates = pos.discovered_check_candidates();
85
86   checkSq[PAWN]   = pos.attacks_from<PAWN>(ksq, them);
87   checkSq[KNIGHT] = pos.attacks_from<KNIGHT>(ksq);
88   checkSq[BISHOP] = pos.attacks_from<BISHOP>(ksq);
89   checkSq[ROOK]   = pos.attacks_from<ROOK>(ksq);
90   checkSq[QUEEN]  = checkSq[BISHOP] | checkSq[ROOK];
91   checkSq[KING]   = EmptyBoardBB;
92 }
93
94
95 /// Position c'tors. Here we always create a copy of the original position
96 /// or the FEN string, we want the new born Position object do not depend
97 /// on any external data so we detach state pointer from the source one.
98
99 Position::Position(const Position& pos, int th) {
100
101   memcpy(this, &pos, sizeof(Position));
102   threadID = th;
103   nodes = 0;
104
105   assert(pos_is_ok());
106 }
107
108 Position::Position(const string& fen, bool isChess960, int th) {
109
110   from_fen(fen, isChess960);
111   threadID = th;
112 }
113
114
115 /// Position::from_fen() initializes the position object with the given FEN
116 /// string. This function is not very robust - make sure that input FENs are
117 /// correct (this is assumed to be the responsibility of the GUI).
118
119 void Position::from_fen(const string& fenStr, bool isChess960) {
120 /*
121    A FEN string defines a particular position using only the ASCII character set.
122
123    A FEN string contains six fields. The separator between fields is a space. The fields are:
124
125    1) Piece placement (from white's perspective). Each rank is described, starting with rank 8 and ending
126       with rank 1; within each rank, the contents of each square are described from file A through file H.
127       Following the Standard Algebraic Notation (SAN), each piece is identified by a single letter taken
128       from the standard English names. White pieces are designated using upper-case letters ("PNBRQK")
129       while Black take lowercase ("pnbrqk"). Blank squares are noted using digits 1 through 8 (the number
130       of blank squares), and "/" separate ranks.
131
132    2) Active color. "w" means white moves next, "b" means black.
133
134    3) Castling availability. If neither side can castle, this is "-". Otherwise, this has one or more
135       letters: "K" (White can castle kingside), "Q" (White can castle queenside), "k" (Black can castle
136       kingside), and/or "q" (Black can castle queenside).
137
138    4) En passant target square in algebraic notation. If there's no en passant target square, this is "-".
139       If a pawn has just made a 2-square move, this is the position "behind" the pawn. This is recorded
140       regardless of whether there is a pawn in position to make an en passant capture.
141
142    5) Halfmove clock: This is the number of halfmoves since the last pawn advance or capture. This is used
143       to determine if a draw can be claimed under the fifty-move rule.
144
145    6) Fullmove number: The number of the full move. It starts at 1, and is incremented after Black's move.
146 */
147
148   char col, row, token;
149   size_t p;
150   Square sq = SQ_A8;
151   std::istringstream fen(fenStr);
152
153   clear();
154   fen >> std::noskipws;
155
156   // 1. Piece placement
157   while ((fen >> token) && !isspace(token))
158   {
159       if (token == '/')
160           sq -= Square(16); // Jump back of 2 rows
161
162       else if (isdigit(token))
163           sq += Square(token - '0'); // Skip the given number of files
164
165       else if ((p = PieceToChar.find(token)) != string::npos)
166       {
167           put_piece(Piece(p), sq);
168           sq++;
169       }
170   }
171
172   // 2. Active color
173   fen >> token;
174   sideToMove = (token == 'w' ? WHITE : BLACK);
175   fen >> token;
176
177   // 3. Castling availability. Compatible with 3 standards: Normal FEN standard,
178   // Shredder-FEN that uses the letters of the columns on which the rooks began
179   // the game instead of KQkq and also X-FEN standard that, in case of Chess960,
180   // if an inner rook is associated with the castling right, the castling tag is
181   // replaced by the file letter of the involved rook, as for the Shredder-FEN.
182   while ((fen >> token) && !isspace(token))
183   {
184       Square rsq;
185       Color c = islower(token) ? BLACK : WHITE;
186       Piece rook = make_piece(c, ROOK);
187
188       token = char(toupper(token));
189
190       if (token == 'K')
191           for (rsq = relative_square(c, SQ_H1); piece_on(rsq) != rook; rsq--) {}
192
193       else if (token == 'Q')
194           for (rsq = relative_square(c, SQ_A1); piece_on(rsq) != rook; rsq++) {}
195
196       else if (token >= 'A' && token <= 'H')
197           rsq = make_square(File(token - 'A'), relative_rank(c, RANK_1));
198
199       else
200           continue;
201
202       set_castle_right(king_square(c), rsq);
203   }
204
205   // 4. En passant square. Ignore if no pawn capture is possible
206   if (   ((fen >> col) && (col >= 'a' && col <= 'h'))
207       && ((fen >> row) && (row == '3' || row == '6')))
208   {
209       st->epSquare = make_square(File(col - 'a'), Rank(row - '1'));
210
211       if (!(attackers_to(st->epSquare) & pieces(PAWN, sideToMove)))
212           st->epSquare = SQ_NONE;
213   }
214
215   // 5-6. Halfmove clock and fullmove number
216   fen >> std::skipws >> st->rule50 >> startPosPly;
217
218   // Convert from fullmove starting from 1 to ply starting from 0,
219   // handle also common incorrect FEN with fullmove = 0.
220   startPosPly = Max(2 * (startPosPly - 1), 0) + int(sideToMove == BLACK);
221
222   st->key = compute_key();
223   st->pawnKey = compute_pawn_key();
224   st->materialKey = compute_material_key();
225   st->value = compute_value();
226   st->npMaterial[WHITE] = compute_non_pawn_material(WHITE);
227   st->npMaterial[BLACK] = compute_non_pawn_material(BLACK);
228   st->checkersBB = attackers_to(king_square(sideToMove)) & pieces(flip(sideToMove));
229   chess960 = isChess960;
230
231   assert(pos_is_ok());
232 }
233
234
235 /// Position::set_castle_right() is an helper function used to set castling
236 /// rights given the corresponding king and rook starting squares.
237
238 void Position::set_castle_right(Square ksq, Square rsq) {
239
240   int f = (rsq < ksq ? WHITE_OOO : WHITE_OO) << color_of(piece_on(ksq));
241
242   st->castleRights |= f;
243   castleRightsMask[ksq] ^= f;
244   castleRightsMask[rsq] ^= f;
245   castleRookSquare[f] = rsq;
246 }
247
248
249 /// Position::to_fen() returns a FEN representation of the position. In case
250 /// of Chess960 the Shredder-FEN notation is used. Mainly a debugging function.
251
252 const string Position::to_fen() const {
253
254   std::ostringstream fen;
255   Square sq;
256   int emptyCnt;
257
258   for (Rank rank = RANK_8; rank >= RANK_1; rank--)
259   {
260       emptyCnt = 0;
261
262       for (File file = FILE_A; file <= FILE_H; file++)
263       {
264           sq = make_square(file, rank);
265
266           if (!square_is_empty(sq))
267           {
268               if (emptyCnt)
269               {
270                   fen << emptyCnt;
271                   emptyCnt = 0;
272               }
273               fen << PieceToChar[piece_on(sq)];
274           }
275           else
276               emptyCnt++;
277       }
278
279       if (emptyCnt)
280           fen << emptyCnt;
281
282       if (rank > RANK_1)
283           fen << '/';
284   }
285
286   fen << (sideToMove == WHITE ? " w " : " b ");
287
288   if (st->castleRights != CASTLES_NONE)
289   {
290       if (can_castle(WHITE_OO))
291           fen << (chess960 ? char(toupper(file_to_char(file_of(castle_rook_square(WHITE_OO))))) : 'K');
292
293       if (can_castle(WHITE_OOO))
294           fen << (chess960 ? char(toupper(file_to_char(file_of(castle_rook_square(WHITE_OOO))))) : 'Q');
295
296       if (can_castle(BLACK_OO))
297           fen << (chess960 ? file_to_char(file_of(castle_rook_square(BLACK_OO))) : 'k');
298
299       if (can_castle(BLACK_OOO))
300           fen << (chess960 ? file_to_char(file_of(castle_rook_square(BLACK_OOO))) : 'q');
301   } else
302       fen << '-';
303
304   fen << (ep_square() == SQ_NONE ? " -" : " " + square_to_string(ep_square()))
305       << " " << st->rule50 << " " << 1 + (startPosPly - int(sideToMove == BLACK)) / 2;
306
307   return fen.str();
308 }
309
310
311 /// Position::print() prints an ASCII representation of the position to
312 /// the standard output. If a move is given then also the san is printed.
313
314 void Position::print(Move move) const {
315
316   const char* dottedLine = "\n+---+---+---+---+---+---+---+---+\n";
317
318   if (move)
319   {
320       Position p(*this, thread());
321       string dd = (sideToMove == BLACK ? ".." : "");
322       cout << "\nMove is: " << dd << move_to_san(p, move);
323   }
324
325   for (Rank rank = RANK_8; rank >= RANK_1; rank--)
326   {
327       cout << dottedLine << '|';
328       for (File file = FILE_A; file <= FILE_H; file++)
329       {
330           Square sq = make_square(file, rank);
331           Piece piece = piece_on(sq);
332
333           if (piece == PIECE_NONE && color_of(sq) == DARK)
334               piece = PIECE_NONE_DARK_SQ;
335
336           char c = (color_of(piece_on(sq)) == BLACK ? '=' : ' ');
337           cout << c << PieceToChar[piece] << c << '|';
338       }
339   }
340   cout << dottedLine << "Fen is: " << to_fen() << "\nKey is: " << st->key << endl;
341 }
342
343
344 /// Position:hidden_checkers<>() returns a bitboard of all pinned (against the
345 /// king) pieces for the given color. Or, when template parameter FindPinned is
346 /// false, the function return the pieces of the given color candidate for a
347 /// discovery check against the enemy king.
348
349 template<bool FindPinned>
350 Bitboard Position::hidden_checkers() const {
351
352   // Pinned pieces protect our king, dicovery checks attack the enemy king
353   Bitboard b, result = EmptyBoardBB;
354   Bitboard pinners = pieces(FindPinned ? flip(sideToMove) : sideToMove);
355   Square ksq = king_square(FindPinned ? sideToMove : flip(sideToMove));
356
357   // Pinners are sliders, that give check when candidate pinned is removed
358   pinners &=  (pieces(ROOK, QUEEN) & RookPseudoAttacks[ksq])
359             | (pieces(BISHOP, QUEEN) & BishopPseudoAttacks[ksq]);
360
361   while (pinners)
362   {
363       b = squares_between(ksq, pop_1st_bit(&pinners)) & occupied_squares();
364
365       // Only one bit set and is an our piece?
366       if (b && !(b & (b - 1)) && (b & pieces(sideToMove)))
367           result |= b;
368   }
369   return result;
370 }
371
372
373 /// Position:pinned_pieces() returns a bitboard of all pinned (against the
374 /// king) pieces for the side to move.
375
376 Bitboard Position::pinned_pieces() const {
377
378   return hidden_checkers<true>();
379 }
380
381
382 /// Position:discovered_check_candidates() returns a bitboard containing all
383 /// pieces for the side to move which are candidates for giving a discovered
384 /// check.
385
386 Bitboard Position::discovered_check_candidates() const {
387
388   return hidden_checkers<false>();
389 }
390
391 /// Position::attackers_to() computes a bitboard of all pieces which attacks a
392 /// given square. Slider attacks use occ bitboard as occupancy.
393
394 Bitboard Position::attackers_to(Square s, Bitboard occ) const {
395
396   return  (attacks_from<PAWN>(s, BLACK) & pieces(PAWN, WHITE))
397         | (attacks_from<PAWN>(s, WHITE) & pieces(PAWN, BLACK))
398         | (attacks_from<KNIGHT>(s)      & pieces(KNIGHT))
399         | (rook_attacks_bb(s, occ)      & pieces(ROOK, QUEEN))
400         | (bishop_attacks_bb(s, occ)    & pieces(BISHOP, QUEEN))
401         | (attacks_from<KING>(s)        & pieces(KING));
402 }
403
404 /// Position::attacks_from() computes a bitboard of all attacks of a given piece
405 /// put in a given square. Slider attacks use occ bitboard as occupancy.
406
407 Bitboard Position::attacks_from(Piece p, Square s, Bitboard occ) {
408
409   assert(square_is_ok(s));
410
411   switch (p)
412   {
413   case WB: case BB: return bishop_attacks_bb(s, occ);
414   case WR: case BR: return rook_attacks_bb(s, occ);
415   case WQ: case BQ: return bishop_attacks_bb(s, occ) | rook_attacks_bb(s, occ);
416   default: return StepAttacksBB[p][s];
417   }
418 }
419
420
421 /// Position::move_attacks_square() tests whether a move from the current
422 /// position attacks a given square.
423
424 bool Position::move_attacks_square(Move m, Square s) const {
425
426   assert(is_ok(m));
427   assert(square_is_ok(s));
428
429   Bitboard occ, xray;
430   Square f = move_from(m), t = move_to(m);
431
432   assert(!square_is_empty(f));
433
434   if (bit_is_set(attacks_from(piece_on(f), t), s))
435       return true;
436
437   // Move the piece and scan for X-ray attacks behind it
438   occ = occupied_squares();
439   do_move_bb(&occ, make_move_bb(f, t));
440   xray = ( (rook_attacks_bb(s, occ)   & pieces(ROOK, QUEEN))
441           |(bishop_attacks_bb(s, occ) & pieces(BISHOP, QUEEN)))
442          & pieces(color_of(piece_on(f)));
443
444   // If we have attacks we need to verify that are caused by our move
445   // and are not already existent ones.
446   return xray && (xray ^ (xray & attacks_from<QUEEN>(s)));
447 }
448
449
450 /// Position::pl_move_is_legal() tests whether a pseudo-legal move is legal
451
452 bool Position::pl_move_is_legal(Move m, Bitboard pinned) const {
453
454   assert(is_ok(m));
455   assert(pinned == pinned_pieces());
456
457   Color us = side_to_move();
458   Square from = move_from(m);
459
460   assert(color_of(piece_on(from)) == us);
461   assert(piece_on(king_square(us)) == make_piece(us, KING));
462
463   // En passant captures are a tricky special case. Because they are rather
464   // uncommon, we do it simply by testing whether the king is attacked after
465   // the move is made.
466   if (is_enpassant(m))
467   {
468       Color them = flip(us);
469       Square to = move_to(m);
470       Square capsq = to + pawn_push(them);
471       Square ksq = king_square(us);
472       Bitboard b = occupied_squares();
473
474       assert(to == ep_square());
475       assert(piece_on(from) == make_piece(us, PAWN));
476       assert(piece_on(capsq) == make_piece(them, PAWN));
477       assert(piece_on(to) == PIECE_NONE);
478
479       clear_bit(&b, from);
480       clear_bit(&b, capsq);
481       set_bit(&b, to);
482
483       return   !(rook_attacks_bb(ksq, b) & pieces(ROOK, QUEEN, them))
484             && !(bishop_attacks_bb(ksq, b) & pieces(BISHOP, QUEEN, them));
485   }
486
487   // If the moving piece is a king, check whether the destination
488   // square is attacked by the opponent. Castling moves are checked
489   // for legality during move generation.
490   if (type_of(piece_on(from)) == KING)
491       return is_castle(m) || !(attackers_to(move_to(m)) & pieces(flip(us)));
492
493   // A non-king move is legal if and only if it is not pinned or it
494   // is moving along the ray towards or away from the king.
495   return   !pinned
496         || !bit_is_set(pinned, from)
497         ||  squares_aligned(from, move_to(m), king_square(us));
498 }
499
500
501 /// Position::move_is_legal() takes a random move and tests whether the move
502 /// is legal. This version is not very fast and should be used only
503 /// in non time-critical paths.
504
505 bool Position::move_is_legal(const Move m) const {
506
507   for (MoveList<MV_LEGAL> ml(*this); !ml.end(); ++ml)
508       if (ml.move() == m)
509           return true;
510
511   return false;
512 }
513
514
515 /// Position::is_pseudo_legal() takes a random move and tests whether the move
516 /// is pseudo legal. It is used to validate moves from TT that can be corrupted
517 /// due to SMP concurrent access or hash position key aliasing.
518
519 bool Position::is_pseudo_legal(const Move m) const {
520
521   Color us = sideToMove;
522   Color them = flip(sideToMove);
523   Square from = move_from(m);
524   Square to = move_to(m);
525   Piece pc = piece_on(from);
526
527   // Use a slower but simpler function for uncommon cases
528   if (is_special(m))
529       return move_is_legal(m);
530
531   // Is not a promotion, so promotion piece must be empty
532   if (promotion_piece_type(m) - 2 != PIECE_TYPE_NONE)
533       return false;
534
535   // If the from square is not occupied by a piece belonging to the side to
536   // move, the move is obviously not legal.
537   if (pc == PIECE_NONE || color_of(pc) != us)
538       return false;
539
540   // The destination square cannot be occupied by a friendly piece
541   if (color_of(piece_on(to)) == us)
542       return false;
543
544   // Handle the special case of a pawn move
545   if (type_of(pc) == PAWN)
546   {
547       // Move direction must be compatible with pawn color
548       int direction = to - from;
549       if ((us == WHITE) != (direction > 0))
550           return false;
551
552       // We have already handled promotion moves, so destination
553       // cannot be on the 8/1th rank.
554       if (rank_of(to) == RANK_8 || rank_of(to) == RANK_1)
555           return false;
556
557       // Proceed according to the square delta between the origin and
558       // destination squares.
559       switch (direction)
560       {
561       case DELTA_NW:
562       case DELTA_NE:
563       case DELTA_SW:
564       case DELTA_SE:
565       // Capture. The destination square must be occupied by an enemy
566       // piece (en passant captures was handled earlier).
567       if (color_of(piece_on(to)) != them)
568           return false;
569
570       // From and to files must be one file apart, avoids a7h5
571       if (abs(file_of(from) - file_of(to)) != 1)
572           return false;
573       break;
574
575       case DELTA_N:
576       case DELTA_S:
577       // Pawn push. The destination square must be empty.
578       if (!square_is_empty(to))
579           return false;
580       break;
581
582       case DELTA_NN:
583       // Double white pawn push. The destination square must be on the fourth
584       // rank, and both the destination square and the square between the
585       // source and destination squares must be empty.
586       if (   rank_of(to) != RANK_4
587           || !square_is_empty(to)
588           || !square_is_empty(from + DELTA_N))
589           return false;
590       break;
591
592       case DELTA_SS:
593       // Double black pawn push. The destination square must be on the fifth
594       // rank, and both the destination square and the square between the
595       // source and destination squares must be empty.
596       if (   rank_of(to) != RANK_5
597           || !square_is_empty(to)
598           || !square_is_empty(from + DELTA_S))
599           return false;
600       break;
601
602       default:
603           return false;
604       }
605   }
606   else if (!bit_is_set(attacks_from(pc, from), to))
607       return false;
608
609   // Evasions generator already takes care to avoid some kind of illegal moves
610   // and pl_move_is_legal() relies on this. So we have to take care that the
611   // same kind of moves are filtered out here.
612   if (in_check())
613   {
614       // In case of king moves under check we have to remove king so to catch
615       // as invalid moves like b1a1 when opposite queen is on c1.
616       if (type_of(piece_on(from)) == KING)
617       {
618           Bitboard b = occupied_squares();
619           clear_bit(&b, from);
620           if (attackers_to(move_to(m), b) & pieces(flip(us)))
621               return false;
622       }
623       else
624       {
625           Bitboard target = checkers();
626           Square checksq = pop_1st_bit(&target);
627
628           if (target) // double check ? In this case a king move is required
629               return false;
630
631           // Our move must be a blocking evasion or a capture of the checking piece
632           target = squares_between(checksq, king_square(us)) | checkers();
633           if (!bit_is_set(target, move_to(m)))
634               return false;
635       }
636   }
637
638   return true;
639 }
640
641
642 /// Position::move_gives_check() tests whether a pseudo-legal move gives a check
643
644 bool Position::move_gives_check(Move m, const CheckInfo& ci) const {
645
646   assert(is_ok(m));
647   assert(ci.dcCandidates == discovered_check_candidates());
648   assert(color_of(piece_on(move_from(m))) == side_to_move());
649
650   Square from = move_from(m);
651   Square to = move_to(m);
652   PieceType pt = type_of(piece_on(from));
653
654   // Direct check ?
655   if (bit_is_set(ci.checkSq[pt], to))
656       return true;
657
658   // Discovery check ?
659   if (ci.dcCandidates && bit_is_set(ci.dcCandidates, from))
660   {
661       // For pawn and king moves we need to verify also direction
662       if (  (pt != PAWN && pt != KING)
663           || !squares_aligned(from, to, king_square(flip(side_to_move()))))
664           return true;
665   }
666
667   // Can we skip the ugly special cases ?
668   if (!is_special(m))
669       return false;
670
671   Color us = side_to_move();
672   Bitboard b = occupied_squares();
673   Square ksq = king_square(flip(us));
674
675   // Promotion with check ?
676   if (is_promotion(m))
677   {
678       clear_bit(&b, from);
679       return bit_is_set(attacks_from(Piece(promotion_piece_type(m)), to, b), ksq);
680   }
681
682   // En passant capture with check ? We have already handled the case
683   // of direct checks and ordinary discovered check, the only case we
684   // need to handle is the unusual case of a discovered check through
685   // the captured pawn.
686   if (is_enpassant(m))
687   {
688       Square capsq = make_square(file_of(to), rank_of(from));
689       clear_bit(&b, from);
690       clear_bit(&b, capsq);
691       set_bit(&b, to);
692       return  (rook_attacks_bb(ksq, b) & pieces(ROOK, QUEEN, us))
693             ||(bishop_attacks_bb(ksq, b) & pieces(BISHOP, QUEEN, us));
694   }
695
696   // Castling with check ?
697   if (is_castle(m))
698   {
699       Square kfrom, kto, rfrom, rto;
700       kfrom = from;
701       rfrom = to;
702
703       if (rfrom > kfrom)
704       {
705           kto = relative_square(us, SQ_G1);
706           rto = relative_square(us, SQ_F1);
707       } else {
708           kto = relative_square(us, SQ_C1);
709           rto = relative_square(us, SQ_D1);
710       }
711       clear_bit(&b, kfrom);
712       clear_bit(&b, rfrom);
713       set_bit(&b, rto);
714       set_bit(&b, kto);
715       return bit_is_set(rook_attacks_bb(rto, b), ksq);
716   }
717
718   return false;
719 }
720
721
722 /// Position::do_move() makes a move, and saves all information necessary
723 /// to a StateInfo object. The move is assumed to be legal. Pseudo-legal
724 /// moves should be filtered out before this function is called.
725
726 void Position::do_move(Move m, StateInfo& newSt) {
727
728   CheckInfo ci(*this);
729   do_move(m, newSt, ci, move_gives_check(m, ci));
730 }
731
732 void Position::do_move(Move m, StateInfo& newSt, const CheckInfo& ci, bool moveIsCheck) {
733
734   assert(is_ok(m));
735   assert(&newSt != st);
736
737   nodes++;
738   Key key = st->key;
739
740   // Copy some fields of old state to our new StateInfo object except the ones
741   // which are recalculated from scratch anyway, then switch our state pointer
742   // to point to the new, ready to be updated, state.
743   struct ReducedStateInfo {
744     Key pawnKey, materialKey;
745     Value npMaterial[2];
746     int castleRights, rule50, pliesFromNull;
747     Score value;
748     Square epSquare;
749   };
750
751   memcpy(&newSt, st, sizeof(ReducedStateInfo));
752
753   newSt.previous = st;
754   st = &newSt;
755
756   // Update side to move
757   key ^= zobSideToMove;
758
759   // Increment the 50 moves rule draw counter. Resetting it to zero in the
760   // case of non-reversible moves is taken care of later.
761   st->rule50++;
762   st->pliesFromNull++;
763
764   if (is_castle(m))
765   {
766       st->key = key;
767       do_castle_move<true>(m);
768       return;
769   }
770
771   Color us = side_to_move();
772   Color them = flip(us);
773   Square from = move_from(m);
774   Square to = move_to(m);
775   Piece piece = piece_on(from);
776   PieceType pt = type_of(piece);
777   PieceType capture = is_enpassant(m) ? PAWN : type_of(piece_on(to));
778
779   assert(color_of(piece) == us);
780   assert(color_of(piece_on(to)) != us);
781   assert(capture != KING);
782
783   if (capture)
784   {
785       Square capsq = to;
786
787       // If the captured piece is a pawn, update pawn hash key, otherwise
788       // update non-pawn material.
789       if (capture == PAWN)
790       {
791           if (is_enpassant(m))
792           {
793               capsq += pawn_push(them);
794
795               assert(pt == PAWN);
796               assert(to == st->epSquare);
797               assert(relative_rank(us, to) == RANK_6);
798               assert(piece_on(to) == PIECE_NONE);
799               assert(piece_on(capsq) == make_piece(them, PAWN));
800
801               board[capsq] = PIECE_NONE;
802           }
803
804           st->pawnKey ^= zobrist[them][PAWN][capsq];
805       }
806       else
807           st->npMaterial[them] -= PieceValueMidgame[capture];
808
809       // Remove the captured piece
810       clear_bit(&byColorBB[them], capsq);
811       clear_bit(&byTypeBB[capture], capsq);
812       clear_bit(&occupied, capsq);
813
814       // Update piece list, move the last piece at index[capsq] position and
815       // shrink the list.
816       //
817       // WARNING: This is a not revresible operation. When we will reinsert the
818       // captured piece in undo_move() we will put it at the end of the list and
819       // not in its original place, it means index[] and pieceList[] are not
820       // guaranteed to be invariant to a do_move() + undo_move() sequence.
821       Square lastSquare = pieceList[them][capture][--pieceCount[them][capture]];
822       index[lastSquare] = index[capsq];
823       pieceList[them][capture][index[lastSquare]] = lastSquare;
824       pieceList[them][capture][pieceCount[them][capture]] = SQ_NONE;
825
826       // Update hash keys
827       key ^= zobrist[them][capture][capsq];
828       st->materialKey ^= zobrist[them][capture][pieceCount[them][capture]];
829
830       // Update incremental scores
831       st->value -= pst(make_piece(them, capture), capsq);
832
833       // Reset rule 50 counter
834       st->rule50 = 0;
835   }
836
837   // Update hash key
838   key ^= zobrist[us][pt][from] ^ zobrist[us][pt][to];
839
840   // Reset en passant square
841   if (st->epSquare != SQ_NONE)
842   {
843       key ^= zobEp[st->epSquare];
844       st->epSquare = SQ_NONE;
845   }
846
847   // Update castle rights if needed
848   if (    st->castleRights != CASTLES_NONE
849       && (castleRightsMask[from] & castleRightsMask[to]) != ALL_CASTLES)
850   {
851       key ^= zobCastle[st->castleRights];
852       st->castleRights &= castleRightsMask[from] & castleRightsMask[to];
853       key ^= zobCastle[st->castleRights];
854   }
855
856   // Prefetch TT access as soon as we know key is updated
857   prefetch((char*)TT.first_entry(key));
858
859   // Move the piece
860   Bitboard move_bb = make_move_bb(from, to);
861   do_move_bb(&byColorBB[us], move_bb);
862   do_move_bb(&byTypeBB[pt], move_bb);
863   do_move_bb(&occupied, move_bb);
864
865   board[to] = board[from];
866   board[from] = PIECE_NONE;
867
868   // Update piece lists, index[from] is not updated and becomes stale. This
869   // works as long as index[] is accessed just by known occupied squares.
870   index[to] = index[from];
871   pieceList[us][pt][index[to]] = to;
872
873   // If the moving piece is a pawn do some special extra work
874   if (pt == PAWN)
875   {
876       // Set en-passant square, only if moved pawn can be captured
877       if (   (to ^ from) == 16
878           && (attacks_from<PAWN>(from + pawn_push(us), us) & pieces(PAWN, them)))
879       {
880           st->epSquare = Square((from + to) / 2);
881           key ^= zobEp[st->epSquare];
882       }
883
884       if (is_promotion(m))
885       {
886           PieceType promotion = promotion_piece_type(m);
887
888           assert(relative_rank(us, to) == RANK_8);
889           assert(promotion >= KNIGHT && promotion <= QUEEN);
890
891           // Replace the pawn with the promoted piece
892           clear_bit(&byTypeBB[PAWN], to);
893           set_bit(&byTypeBB[promotion], to);
894           board[to] = make_piece(us, promotion);
895
896           // Update piece lists, move the last pawn at index[to] position
897           // and shrink the list. Add a new promotion piece to the list.
898           Square lastSquare = pieceList[us][PAWN][--pieceCount[us][PAWN]];
899           index[lastSquare] = index[to];
900           pieceList[us][PAWN][index[lastSquare]] = lastSquare;
901           pieceList[us][PAWN][pieceCount[us][PAWN]] = SQ_NONE;
902           index[to] = pieceCount[us][promotion];
903           pieceList[us][promotion][index[to]] = to;
904
905           // Update hash keys
906           key ^= zobrist[us][PAWN][to] ^ zobrist[us][promotion][to];
907           st->pawnKey ^= zobrist[us][PAWN][to];
908           st->materialKey ^=  zobrist[us][promotion][pieceCount[us][promotion]++]
909                             ^ zobrist[us][PAWN][pieceCount[us][PAWN]];
910
911           // Update incremental score
912           st->value +=  pst(make_piece(us, promotion), to)
913                       - pst(make_piece(us, PAWN), to);
914
915           // Update material
916           st->npMaterial[us] += PieceValueMidgame[promotion];
917       }
918
919       // Update pawn hash key
920       st->pawnKey ^= zobrist[us][PAWN][from] ^ zobrist[us][PAWN][to];
921
922       // Reset rule 50 draw counter
923       st->rule50 = 0;
924   }
925
926   // Prefetch pawn and material hash tables
927   Threads[threadID].pawnTable.prefetch(st->pawnKey);
928   Threads[threadID].materialTable.prefetch(st->materialKey);
929
930   // Update incremental scores
931   st->value += pst_delta(piece, from, to);
932
933   // Set capture piece
934   st->capturedType = capture;
935
936   // Update the key with the final value
937   st->key = key;
938
939   // Update checkers bitboard, piece must be already moved
940   st->checkersBB = EmptyBoardBB;
941
942   if (moveIsCheck)
943   {
944       if (is_special(m))
945           st->checkersBB = attackers_to(king_square(them)) & pieces(us);
946       else
947       {
948           // Direct checks
949           if (bit_is_set(ci.checkSq[pt], to))
950               st->checkersBB = SetMaskBB[to];
951
952           // Discovery checks
953           if (ci.dcCandidates && bit_is_set(ci.dcCandidates, from))
954           {
955               if (pt != ROOK)
956                   st->checkersBB |= attacks_from<ROOK>(king_square(them)) & pieces(ROOK, QUEEN, us);
957
958               if (pt != BISHOP)
959                   st->checkersBB |= attacks_from<BISHOP>(king_square(them)) & pieces(BISHOP, QUEEN, us);
960           }
961       }
962   }
963
964   // Finish
965   sideToMove = flip(sideToMove);
966   st->value += (sideToMove == WHITE ?  TempoValue : -TempoValue);
967
968   assert(pos_is_ok());
969 }
970
971
972 /// Position::undo_move() unmakes a move. When it returns, the position should
973 /// be restored to exactly the same state as before the move was made.
974
975 void Position::undo_move(Move m) {
976
977   assert(is_ok(m));
978
979   sideToMove = flip(sideToMove);
980
981   if (is_castle(m))
982   {
983       do_castle_move<false>(m);
984       return;
985   }
986
987   Color us = side_to_move();
988   Color them = flip(us);
989   Square from = move_from(m);
990   Square to = move_to(m);
991   Piece piece = piece_on(to);
992   PieceType pt = type_of(piece);
993   PieceType capture = st->capturedType;
994
995   assert(square_is_empty(from));
996   assert(color_of(piece) == us);
997   assert(capture != KING);
998
999   if (is_promotion(m))
1000   {
1001       PieceType promotion = promotion_piece_type(m);
1002
1003       assert(promotion == pt);
1004       assert(relative_rank(us, to) == RANK_8);
1005       assert(promotion >= KNIGHT && promotion <= QUEEN);
1006
1007       // Replace the promoted piece with the pawn
1008       clear_bit(&byTypeBB[promotion], to);
1009       set_bit(&byTypeBB[PAWN], to);
1010       board[to] = make_piece(us, PAWN);
1011
1012       // Update piece lists, move the last promoted piece at index[to] position
1013       // and shrink the list. Add a new pawn to the list.
1014       Square lastSquare = pieceList[us][promotion][--pieceCount[us][promotion]];
1015       index[lastSquare] = index[to];
1016       pieceList[us][promotion][index[lastSquare]] = lastSquare;
1017       pieceList[us][promotion][pieceCount[us][promotion]] = SQ_NONE;
1018       index[to] = pieceCount[us][PAWN]++;
1019       pieceList[us][PAWN][index[to]] = to;
1020
1021       pt = PAWN;
1022   }
1023
1024   // Put the piece back at the source square
1025   Bitboard move_bb = make_move_bb(to, from);
1026   do_move_bb(&byColorBB[us], move_bb);
1027   do_move_bb(&byTypeBB[pt], move_bb);
1028   do_move_bb(&occupied, move_bb);
1029
1030   board[from] = board[to];
1031   board[to] = PIECE_NONE;
1032
1033   // Update piece lists, index[to] is not updated and becomes stale. This
1034   // works as long as index[] is accessed just by known occupied squares.
1035   index[from] = index[to];
1036   pieceList[us][pt][index[from]] = from;
1037
1038   if (capture)
1039   {
1040       Square capsq = to;
1041
1042       if (is_enpassant(m))
1043       {
1044           capsq -= pawn_push(us);
1045
1046           assert(pt == PAWN);
1047           assert(to == st->previous->epSquare);
1048           assert(relative_rank(us, to) == RANK_6);
1049           assert(piece_on(capsq) == PIECE_NONE);
1050       }
1051
1052       // Restore the captured piece
1053       set_bit(&byColorBB[them], capsq);
1054       set_bit(&byTypeBB[capture], capsq);
1055       set_bit(&occupied, capsq);
1056
1057       board[capsq] = make_piece(them, capture);
1058
1059       // Update piece list, add a new captured piece in capsq square
1060       index[capsq] = pieceCount[them][capture]++;
1061       pieceList[them][capture][index[capsq]] = capsq;
1062   }
1063
1064   // Finally point our state pointer back to the previous state
1065   st = st->previous;
1066
1067   assert(pos_is_ok());
1068 }
1069
1070
1071 /// Position::do_castle_move() is a private method used to do/undo a castling
1072 /// move. Note that castling moves are encoded as "king captures friendly rook"
1073 /// moves, for instance white short castling in a non-Chess960 game is encoded
1074 /// as e1h1.
1075 template<bool Do>
1076 void Position::do_castle_move(Move m) {
1077
1078   assert(is_ok(m));
1079   assert(is_castle(m));
1080
1081   Square kto, kfrom, rfrom, rto, kAfter, rAfter;
1082
1083   Color us = side_to_move();
1084   Square kBefore = move_from(m);
1085   Square rBefore = move_to(m);
1086
1087   // Find after-castle squares for king and rook
1088   if (rBefore > kBefore) // O-O
1089   {
1090       kAfter = relative_square(us, SQ_G1);
1091       rAfter = relative_square(us, SQ_F1);
1092   }
1093   else // O-O-O
1094   {
1095       kAfter = relative_square(us, SQ_C1);
1096       rAfter = relative_square(us, SQ_D1);
1097   }
1098
1099   kfrom = Do ? kBefore : kAfter;
1100   rfrom = Do ? rBefore : rAfter;
1101
1102   kto = Do ? kAfter : kBefore;
1103   rto = Do ? rAfter : rBefore;
1104
1105   assert(piece_on(kfrom) == make_piece(us, KING));
1106   assert(piece_on(rfrom) == make_piece(us, ROOK));
1107
1108   // Remove pieces from source squares
1109   clear_bit(&byColorBB[us], kfrom);
1110   clear_bit(&byTypeBB[KING], kfrom);
1111   clear_bit(&occupied, kfrom);
1112   clear_bit(&byColorBB[us], rfrom);
1113   clear_bit(&byTypeBB[ROOK], rfrom);
1114   clear_bit(&occupied, rfrom);
1115
1116   // Put pieces on destination squares
1117   set_bit(&byColorBB[us], kto);
1118   set_bit(&byTypeBB[KING], kto);
1119   set_bit(&occupied, kto);
1120   set_bit(&byColorBB[us], rto);
1121   set_bit(&byTypeBB[ROOK], rto);
1122   set_bit(&occupied, rto);
1123
1124   // Update board
1125   Piece king = make_piece(us, KING);
1126   Piece rook = make_piece(us, ROOK);
1127   board[kfrom] = board[rfrom] = PIECE_NONE;
1128   board[kto] = king;
1129   board[rto] = rook;
1130
1131   // Update piece lists
1132   pieceList[us][KING][index[kfrom]] = kto;
1133   pieceList[us][ROOK][index[rfrom]] = rto;
1134   int tmp = index[rfrom]; // In Chess960 could be kto == rfrom
1135   index[kto] = index[kfrom];
1136   index[rto] = tmp;
1137
1138   if (Do)
1139   {
1140       // Reset capture field
1141       st->capturedType = PIECE_TYPE_NONE;
1142
1143       // Update incremental scores
1144       st->value += pst_delta(king, kfrom, kto);
1145       st->value += pst_delta(rook, rfrom, rto);
1146
1147       // Update hash key
1148       st->key ^= zobrist[us][KING][kfrom] ^ zobrist[us][KING][kto];
1149       st->key ^= zobrist[us][ROOK][rfrom] ^ zobrist[us][ROOK][rto];
1150
1151       // Clear en passant square
1152       if (st->epSquare != SQ_NONE)
1153       {
1154           st->key ^= zobEp[st->epSquare];
1155           st->epSquare = SQ_NONE;
1156       }
1157
1158       // Update castling rights
1159       st->key ^= zobCastle[st->castleRights];
1160       st->castleRights &= castleRightsMask[kfrom];
1161       st->key ^= zobCastle[st->castleRights];
1162
1163       // Reset rule 50 counter
1164       st->rule50 = 0;
1165
1166       // Update checkers BB
1167       st->checkersBB = attackers_to(king_square(flip(us))) & pieces(us);
1168
1169       // Finish
1170       sideToMove = flip(sideToMove);
1171       st->value += (sideToMove == WHITE ?  TempoValue : -TempoValue);
1172   }
1173   else
1174       // Undo: point our state pointer back to the previous state
1175       st = st->previous;
1176
1177   assert(pos_is_ok());
1178 }
1179
1180
1181 /// Position::do_null_move() is used to do/undo a "null move": It flips the side
1182 /// to move and updates the hash key without executing any move on the board.
1183 template<bool Do>
1184 void Position::do_null_move(StateInfo& backupSt) {
1185
1186   assert(!in_check());
1187
1188   // Back up the information necessary to undo the null move to the supplied
1189   // StateInfo object. Note that differently from normal case here backupSt
1190   // is actually used as a backup storage not as the new state. This reduces
1191   // the number of fields to be copied.
1192   StateInfo* src = Do ? st : &backupSt;
1193   StateInfo* dst = Do ? &backupSt : st;
1194
1195   dst->key      = src->key;
1196   dst->epSquare = src->epSquare;
1197   dst->value    = src->value;
1198   dst->rule50   = src->rule50;
1199   dst->pliesFromNull = src->pliesFromNull;
1200
1201   sideToMove = flip(sideToMove);
1202
1203   if (Do)
1204   {
1205       if (st->epSquare != SQ_NONE)
1206           st->key ^= zobEp[st->epSquare];
1207
1208       st->key ^= zobSideToMove;
1209       prefetch((char*)TT.first_entry(st->key));
1210
1211       st->epSquare = SQ_NONE;
1212       st->rule50++;
1213       st->pliesFromNull = 0;
1214       st->value += (sideToMove == WHITE) ?  TempoValue : -TempoValue;
1215   }
1216
1217   assert(pos_is_ok());
1218 }
1219
1220 // Explicit template instantiations
1221 template void Position::do_null_move<false>(StateInfo& backupSt);
1222 template void Position::do_null_move<true>(StateInfo& backupSt);
1223
1224
1225 /// Position::see() is a static exchange evaluator: It tries to estimate the
1226 /// material gain or loss resulting from a move. There are three versions of
1227 /// this function: One which takes a destination square as input, one takes a
1228 /// move, and one which takes a 'from' and a 'to' square. The function does
1229 /// not yet understand promotions captures.
1230
1231 int Position::see_sign(Move m) const {
1232
1233   assert(is_ok(m));
1234
1235   Square from = move_from(m);
1236   Square to = move_to(m);
1237
1238   // Early return if SEE cannot be negative because captured piece value
1239   // is not less then capturing one. Note that king moves always return
1240   // here because king midgame value is set to 0.
1241   if (PieceValueMidgame[piece_on(to)] >= PieceValueMidgame[piece_on(from)])
1242       return 1;
1243
1244   return see(m);
1245 }
1246
1247 int Position::see(Move m) const {
1248
1249   Square from, to;
1250   Bitboard occ, attackers, stmAttackers, b;
1251   int swapList[32], slIndex = 1;
1252   PieceType capturedType, pt;
1253   Color stm;
1254
1255   assert(is_ok(m));
1256
1257   // As castle moves are implemented as capturing the rook, they have
1258   // SEE == RookValueMidgame most of the times (unless the rook is under
1259   // attack).
1260   if (is_castle(m))
1261       return 0;
1262
1263   from = move_from(m);
1264   to = move_to(m);
1265   capturedType = type_of(piece_on(to));
1266   occ = occupied_squares();
1267
1268   // Handle en passant moves
1269   if (is_enpassant(m))
1270   {
1271       Square capQq = to - pawn_push(side_to_move());
1272
1273       assert(capturedType == PIECE_TYPE_NONE);
1274       assert(type_of(piece_on(capQq)) == PAWN);
1275
1276       // Remove the captured pawn
1277       clear_bit(&occ, capQq);
1278       capturedType = PAWN;
1279   }
1280
1281   // Find all attackers to the destination square, with the moving piece
1282   // removed, but possibly an X-ray attacker added behind it.
1283   clear_bit(&occ, from);
1284   attackers = attackers_to(to, occ);
1285
1286   // If the opponent has no attackers we are finished
1287   stm = flip(color_of(piece_on(from)));
1288   stmAttackers = attackers & pieces(stm);
1289   if (!stmAttackers)
1290       return PieceValueMidgame[capturedType];
1291
1292   // The destination square is defended, which makes things rather more
1293   // difficult to compute. We proceed by building up a "swap list" containing
1294   // the material gain or loss at each stop in a sequence of captures to the
1295   // destination square, where the sides alternately capture, and always
1296   // capture with the least valuable piece. After each capture, we look for
1297   // new X-ray attacks from behind the capturing piece.
1298   swapList[0] = PieceValueMidgame[capturedType];
1299   capturedType = type_of(piece_on(from));
1300
1301   do {
1302       // Locate the least valuable attacker for the side to move. The loop
1303       // below looks like it is potentially infinite, but it isn't. We know
1304       // that the side to move still has at least one attacker left.
1305       for (pt = PAWN; !(stmAttackers & pieces(pt)); pt++)
1306           assert(pt < KING);
1307
1308       // Remove the attacker we just found from the 'occupied' bitboard,
1309       // and scan for new X-ray attacks behind the attacker.
1310       b = stmAttackers & pieces(pt);
1311       occ ^= (b & (~b + 1));
1312       attackers |=  (rook_attacks_bb(to, occ)   & pieces(ROOK, QUEEN))
1313                   | (bishop_attacks_bb(to, occ) & pieces(BISHOP, QUEEN));
1314
1315       attackers &= occ; // Cut out pieces we've already done
1316
1317       // Add the new entry to the swap list
1318       assert(slIndex < 32);
1319       swapList[slIndex] = -swapList[slIndex - 1] + PieceValueMidgame[capturedType];
1320       slIndex++;
1321
1322       // Remember the value of the capturing piece, and change the side to
1323       // move before beginning the next iteration.
1324       capturedType = pt;
1325       stm = flip(stm);
1326       stmAttackers = attackers & pieces(stm);
1327
1328       // Stop before processing a king capture
1329       if (capturedType == KING && stmAttackers)
1330       {
1331           assert(slIndex < 32);
1332           swapList[slIndex++] = QueenValueMidgame*10;
1333           break;
1334       }
1335   } while (stmAttackers);
1336
1337   // Having built the swap list, we negamax through it to find the best
1338   // achievable score from the point of view of the side to move.
1339   while (--slIndex)
1340       swapList[slIndex-1] = Min(-swapList[slIndex], swapList[slIndex-1]);
1341
1342   return swapList[0];
1343 }
1344
1345
1346 /// Position::clear() erases the position object to a pristine state, with an
1347 /// empty board, white to move, and no castling rights.
1348
1349 void Position::clear() {
1350
1351   st = &startState;
1352   memset(st, 0, sizeof(StateInfo));
1353   st->epSquare = SQ_NONE;
1354
1355   memset(byColorBB,  0, sizeof(Bitboard) * 2);
1356   memset(byTypeBB,   0, sizeof(Bitboard) * 8);
1357   memset(pieceCount, 0, sizeof(int) * 2 * 8);
1358   memset(index,      0, sizeof(int) * 64);
1359
1360   for (int i = 0; i < 8; i++)
1361       for (int j = 0; j < 16; j++)
1362           pieceList[0][i][j] = pieceList[1][i][j] = SQ_NONE;
1363
1364   for (Square sq = SQ_A1; sq <= SQ_H8; sq++)
1365   {
1366       board[sq] = PIECE_NONE;
1367       castleRightsMask[sq] = ALL_CASTLES;
1368   }
1369   sideToMove = WHITE;
1370   nodes = 0;
1371   occupied = 0;
1372 }
1373
1374
1375 /// Position::put_piece() puts a piece on the given square of the board,
1376 /// updating the board array, pieces list, bitboards, and piece counts.
1377
1378 void Position::put_piece(Piece p, Square s) {
1379
1380   Color c = color_of(p);
1381   PieceType pt = type_of(p);
1382
1383   board[s] = p;
1384   index[s] = pieceCount[c][pt]++;
1385   pieceList[c][pt][index[s]] = s;
1386
1387   set_bit(&byTypeBB[pt], s);
1388   set_bit(&byColorBB[c], s);
1389   set_bit(&occupied, s);
1390 }
1391
1392
1393 /// Position::compute_key() computes the hash key of the position. The hash
1394 /// key is usually updated incrementally as moves are made and unmade, the
1395 /// compute_key() function is only used when a new position is set up, and
1396 /// to verify the correctness of the hash key when running in debug mode.
1397
1398 Key Position::compute_key() const {
1399
1400   Key result = zobCastle[st->castleRights];
1401
1402   for (Square s = SQ_A1; s <= SQ_H8; s++)
1403       if (!square_is_empty(s))
1404           result ^= zobrist[color_of(piece_on(s))][type_of(piece_on(s))][s];
1405
1406   if (ep_square() != SQ_NONE)
1407       result ^= zobEp[ep_square()];
1408
1409   if (side_to_move() == BLACK)
1410       result ^= zobSideToMove;
1411
1412   return result;
1413 }
1414
1415
1416 /// Position::compute_pawn_key() computes the hash key of the position. The
1417 /// hash key is usually updated incrementally as moves are made and unmade,
1418 /// the compute_pawn_key() function is only used when a new position is set
1419 /// up, and to verify the correctness of the pawn hash key when running in
1420 /// debug mode.
1421
1422 Key Position::compute_pawn_key() const {
1423
1424   Bitboard b;
1425   Key result = 0;
1426
1427   for (Color c = WHITE; c <= BLACK; c++)
1428   {
1429       b = pieces(PAWN, c);
1430       while (b)
1431           result ^= zobrist[c][PAWN][pop_1st_bit(&b)];
1432   }
1433   return result;
1434 }
1435
1436
1437 /// Position::compute_material_key() computes the hash key of the position.
1438 /// The hash key is usually updated incrementally as moves are made and unmade,
1439 /// the compute_material_key() function is only used when a new position is set
1440 /// up, and to verify the correctness of the material hash key when running in
1441 /// debug mode.
1442
1443 Key Position::compute_material_key() const {
1444
1445   Key result = 0;
1446
1447   for (Color c = WHITE; c <= BLACK; c++)
1448       for (PieceType pt = PAWN; pt <= QUEEN; pt++)
1449           for (int i = 0; i < piece_count(c, pt); i++)
1450               result ^= zobrist[c][pt][i];
1451
1452   return result;
1453 }
1454
1455
1456 /// Position::compute_value() compute the incremental scores for the middle
1457 /// game and the endgame. These functions are used to initialize the incremental
1458 /// scores when a new position is set up, and to verify that the scores are correctly
1459 /// updated by do_move and undo_move when the program is running in debug mode.
1460 Score Position::compute_value() const {
1461
1462   Bitboard b;
1463   Score result = SCORE_ZERO;
1464
1465   for (Color c = WHITE; c <= BLACK; c++)
1466       for (PieceType pt = PAWN; pt <= KING; pt++)
1467       {
1468           b = pieces(pt, c);
1469           while (b)
1470               result += pst(make_piece(c, pt), pop_1st_bit(&b));
1471       }
1472
1473   result += (side_to_move() == WHITE ? TempoValue / 2 : -TempoValue / 2);
1474   return result;
1475 }
1476
1477
1478 /// Position::compute_non_pawn_material() computes the total non-pawn middle
1479 /// game material value for the given side. Material values are updated
1480 /// incrementally during the search, this function is only used while
1481 /// initializing a new Position object.
1482
1483 Value Position::compute_non_pawn_material(Color c) const {
1484
1485   Value result = VALUE_ZERO;
1486
1487   for (PieceType pt = KNIGHT; pt <= QUEEN; pt++)
1488       result += piece_count(c, pt) * PieceValueMidgame[pt];
1489
1490   return result;
1491 }
1492
1493
1494 /// Position::is_draw() tests whether the position is drawn by material,
1495 /// repetition, or the 50 moves rule. It does not detect stalemates, this
1496 /// must be done by the search.
1497 template<bool SkipRepetition>
1498 bool Position::is_draw() const {
1499
1500   // Draw by material?
1501   if (   !pieces(PAWN)
1502       && (non_pawn_material(WHITE) + non_pawn_material(BLACK) <= BishopValueMidgame))
1503       return true;
1504
1505   // Draw by the 50 moves rule?
1506   if (st->rule50 > 99 && !is_mate())
1507       return true;
1508
1509   // Draw by repetition?
1510   if (!SkipRepetition)
1511   {
1512       int i = 4, e = Min(st->rule50, st->pliesFromNull);
1513
1514       if (i <= e)
1515       {
1516           StateInfo* stp = st->previous->previous;
1517
1518           do {
1519               stp = stp->previous->previous;
1520
1521               if (stp->key == st->key)
1522                   return true;
1523
1524               i +=2;
1525
1526           } while (i <= e);
1527       }
1528   }
1529
1530   return false;
1531 }
1532
1533 // Explicit template instantiations
1534 template bool Position::is_draw<false>() const;
1535 template bool Position::is_draw<true>() const;
1536
1537
1538 /// Position::is_mate() returns true or false depending on whether the
1539 /// side to move is checkmated.
1540
1541 bool Position::is_mate() const {
1542
1543   return in_check() && !MoveList<MV_LEGAL>(*this).size();
1544 }
1545
1546
1547 /// Position::init() is a static member function which initializes at startup
1548 /// the various arrays used to compute hash keys and the piece square tables.
1549 /// The latter is a two-step operation: First, the white halves of the tables
1550 /// are copied from PSQT[] tables. Second, the black halves of the tables are
1551 /// initialized by flipping and changing the sign of the white scores.
1552
1553 void Position::init() {
1554
1555   RKISS rk;
1556
1557   for (Color c = WHITE; c <= BLACK; c++)
1558       for (PieceType pt = PAWN; pt <= KING; pt++)
1559           for (Square s = SQ_A1; s <= SQ_H8; s++)
1560               zobrist[c][pt][s] = rk.rand<Key>();
1561
1562   for (Square s = SQ_A1; s <= SQ_H8; s++)
1563       zobEp[s] = rk.rand<Key>();
1564
1565   for (int i = 0; i < 16; i++)
1566       zobCastle[i] = rk.rand<Key>();
1567
1568   zobSideToMove = rk.rand<Key>();
1569   zobExclusion  = rk.rand<Key>();
1570
1571   for (Piece p = WP; p <= WK; p++)
1572   {
1573       Score ps = make_score(PieceValueMidgame[p], PieceValueEndgame[p]);
1574
1575       for (Square s = SQ_A1; s <= SQ_H8; s++)
1576       {
1577           pieceSquareTable[p][s] = ps + PSQT[p][s];
1578           pieceSquareTable[p+8][flip(s)] = -pieceSquareTable[p][s];
1579       }
1580   }
1581 }
1582
1583
1584 /// Position::flip_me() flips position with the white and black sides reversed. This
1585 /// is only useful for debugging especially for finding evaluation symmetry bugs.
1586
1587 void Position::flip_me() {
1588
1589   // Make a copy of current position before to start changing
1590   const Position pos(*this, threadID);
1591
1592   clear();
1593   threadID = pos.thread();
1594
1595   // Board
1596   for (Square s = SQ_A1; s <= SQ_H8; s++)
1597       if (!pos.square_is_empty(s))
1598           put_piece(Piece(pos.piece_on(s) ^ 8), flip(s));
1599
1600   // Side to move
1601   sideToMove = flip(pos.side_to_move());
1602
1603   // Castling rights
1604   if (pos.can_castle(WHITE_OO))
1605       set_castle_right(king_square(BLACK), flip(pos.castle_rook_square(WHITE_OO)));
1606   if (pos.can_castle(WHITE_OOO))
1607       set_castle_right(king_square(BLACK), flip(pos.castle_rook_square(WHITE_OOO)));
1608   if (pos.can_castle(BLACK_OO))
1609       set_castle_right(king_square(WHITE), flip(pos.castle_rook_square(BLACK_OO)));
1610   if (pos.can_castle(BLACK_OOO))
1611       set_castle_right(king_square(WHITE), flip(pos.castle_rook_square(BLACK_OOO)));
1612
1613   // En passant square
1614   if (pos.st->epSquare != SQ_NONE)
1615       st->epSquare = flip(pos.st->epSquare);
1616
1617   // Checkers
1618   st->checkersBB = attackers_to(king_square(sideToMove)) & pieces(flip(sideToMove));
1619
1620   // Hash keys
1621   st->key = compute_key();
1622   st->pawnKey = compute_pawn_key();
1623   st->materialKey = compute_material_key();
1624
1625   // Incremental scores
1626   st->value = compute_value();
1627
1628   // Material
1629   st->npMaterial[WHITE] = compute_non_pawn_material(WHITE);
1630   st->npMaterial[BLACK] = compute_non_pawn_material(BLACK);
1631
1632   assert(pos_is_ok());
1633 }
1634
1635
1636 /// Position::pos_is_ok() performs some consitency checks for the position object.
1637 /// This is meant to be helpful when debugging.
1638
1639 bool Position::pos_is_ok(int* failedStep) const {
1640
1641   // What features of the position should be verified?
1642   const bool debugAll = false;
1643
1644   const bool debugBitboards       = debugAll || false;
1645   const bool debugKingCount       = debugAll || false;
1646   const bool debugKingCapture     = debugAll || false;
1647   const bool debugCheckerCount    = debugAll || false;
1648   const bool debugKey             = debugAll || false;
1649   const bool debugMaterialKey     = debugAll || false;
1650   const bool debugPawnKey         = debugAll || false;
1651   const bool debugIncrementalEval = debugAll || false;
1652   const bool debugNonPawnMaterial = debugAll || false;
1653   const bool debugPieceCounts     = debugAll || false;
1654   const bool debugPieceList       = debugAll || false;
1655   const bool debugCastleSquares   = debugAll || false;
1656
1657   if (failedStep) *failedStep = 1;
1658
1659   // Side to move OK?
1660   if (side_to_move() != WHITE && side_to_move() != BLACK)
1661       return false;
1662
1663   // Are the king squares in the position correct?
1664   if (failedStep) (*failedStep)++;
1665   if (piece_on(king_square(WHITE)) != WK)
1666       return false;
1667
1668   if (failedStep) (*failedStep)++;
1669   if (piece_on(king_square(BLACK)) != BK)
1670       return false;
1671
1672   // Do both sides have exactly one king?
1673   if (failedStep) (*failedStep)++;
1674   if (debugKingCount)
1675   {
1676       int kingCount[2] = {0, 0};
1677       for (Square s = SQ_A1; s <= SQ_H8; s++)
1678           if (type_of(piece_on(s)) == KING)
1679               kingCount[color_of(piece_on(s))]++;
1680
1681       if (kingCount[0] != 1 || kingCount[1] != 1)
1682           return false;
1683   }
1684
1685   // Can the side to move capture the opponent's king?
1686   if (failedStep) (*failedStep)++;
1687   if (debugKingCapture)
1688   {
1689       Color us = side_to_move();
1690       Color them = flip(us);
1691       Square ksq = king_square(them);
1692       if (attackers_to(ksq) & pieces(us))
1693           return false;
1694   }
1695
1696   // Is there more than 2 checkers?
1697   if (failedStep) (*failedStep)++;
1698   if (debugCheckerCount && count_1s<CNT32>(st->checkersBB) > 2)
1699       return false;
1700
1701   // Bitboards OK?
1702   if (failedStep) (*failedStep)++;
1703   if (debugBitboards)
1704   {
1705       // The intersection of the white and black pieces must be empty
1706       if ((pieces(WHITE) & pieces(BLACK)) != EmptyBoardBB)
1707           return false;
1708
1709       // The union of the white and black pieces must be equal to all
1710       // occupied squares
1711       if ((pieces(WHITE) | pieces(BLACK)) != occupied_squares())
1712           return false;
1713
1714       // Separate piece type bitboards must have empty intersections
1715       for (PieceType p1 = PAWN; p1 <= KING; p1++)
1716           for (PieceType p2 = PAWN; p2 <= KING; p2++)
1717               if (p1 != p2 && (pieces(p1) & pieces(p2)))
1718                   return false;
1719   }
1720
1721   // En passant square OK?
1722   if (failedStep) (*failedStep)++;
1723   if (ep_square() != SQ_NONE)
1724   {
1725       // The en passant square must be on rank 6, from the point of view of the
1726       // side to move.
1727       if (relative_rank(side_to_move(), ep_square()) != RANK_6)
1728           return false;
1729   }
1730
1731   // Hash key OK?
1732   if (failedStep) (*failedStep)++;
1733   if (debugKey && st->key != compute_key())
1734       return false;
1735
1736   // Pawn hash key OK?
1737   if (failedStep) (*failedStep)++;
1738   if (debugPawnKey && st->pawnKey != compute_pawn_key())
1739       return false;
1740
1741   // Material hash key OK?
1742   if (failedStep) (*failedStep)++;
1743   if (debugMaterialKey && st->materialKey != compute_material_key())
1744       return false;
1745
1746   // Incremental eval OK?
1747   if (failedStep) (*failedStep)++;
1748   if (debugIncrementalEval && st->value != compute_value())
1749       return false;
1750
1751   // Non-pawn material OK?
1752   if (failedStep) (*failedStep)++;
1753   if (debugNonPawnMaterial)
1754   {
1755       if (st->npMaterial[WHITE] != compute_non_pawn_material(WHITE))
1756           return false;
1757
1758       if (st->npMaterial[BLACK] != compute_non_pawn_material(BLACK))
1759           return false;
1760   }
1761
1762   // Piece counts OK?
1763   if (failedStep) (*failedStep)++;
1764   if (debugPieceCounts)
1765       for (Color c = WHITE; c <= BLACK; c++)
1766           for (PieceType pt = PAWN; pt <= KING; pt++)
1767               if (pieceCount[c][pt] != count_1s<CNT32>(pieces(pt, c)))
1768                   return false;
1769
1770   if (failedStep) (*failedStep)++;
1771   if (debugPieceList)
1772       for (Color c = WHITE; c <= BLACK; c++)
1773           for (PieceType pt = PAWN; pt <= KING; pt++)
1774               for (int i = 0; i < pieceCount[c][pt]; i++)
1775               {
1776                   if (piece_on(piece_list(c, pt)[i]) != make_piece(c, pt))
1777                       return false;
1778
1779                   if (index[piece_list(c, pt)[i]] != i)
1780                       return false;
1781               }
1782
1783   if (failedStep) (*failedStep)++;
1784   if (debugCastleSquares)
1785       for (CastleRight f = WHITE_OO; f <= BLACK_OOO; f = CastleRight(f << 1))
1786       {
1787           if (!can_castle(f))
1788               continue;
1789
1790           Piece rook = (f & (WHITE_OO | WHITE_OOO) ? WR : BR);
1791
1792           if (   castleRightsMask[castleRookSquare[f]] != (ALL_CASTLES ^ f)
1793               || piece_on(castleRookSquare[f]) != rook)
1794               return false;
1795       }
1796
1797   if (failedStep) *failedStep = 0;
1798   return true;
1799 }