git.sesse.net Git - stockfish/blob - src/tt.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-2015 Marco Costalba, Joona Kiiski, Tord Romstad
   5   Copyright (C) 2015-2019 Marco Costalba, Joona Kiiski, Gary Linscott, Tord Romstad
   6
   7   Stockfish is free software: you can redistribute it and/or modify
   8   it under the terms of the GNU General Public License as published by
   9   the Free Software Foundation, either version 3 of the License, or
  10   (at your option) any later version.
  11
  12   Stockfish is distributed in the hope that it will be useful,
  13   but WITHOUT ANY WARRANTY; without even the implied warranty of
  14   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15   GNU General Public License for more details.
  16
  17   You should have received a copy of the GNU General Public License
  18   along with this program.  If not, see <http://www.gnu.org/licenses/>.
  19 */
  20
  21 #include <cstring>   // For std::memset
  22 #include <iostream>
  23 #include <thread>
  24
  25 #include "bitboard.h"
  26 #include "misc.h"
  27 #include "thread.h"
  28 #include "tt.h"
  29 #include "uci.h"
  30
  31 TranspositionTable TT; // Our global transposition table
  32
  33 /// TTEntry::save saves a TTEntry
  34 void TTEntry::save(Key k, Value v, bool PvNode, Bound b, Depth d, Move m, Value ev) {
  35
  36   assert(d / ONE_PLY * ONE_PLY == d);
  37
  38   // Preserve any existing move for the same position
  39   if (m || (k >> 48) != key16)
  40       move16 = (uint16_t)m;
  41
  42   // Overwrite less valuable entries
  43   if (  (k >> 48) != key16
  44       || d / ONE_PLY > depth8 - 4
  45       || b == BOUND_EXACT)
  46   {
  47       key16     = (uint16_t)(k >> 48);
  48       value16   = (int16_t)v;
  49       eval16    = (int16_t)ev;
  50       genBound8 = (uint8_t)(TT.generation8 | PvNode << 2 | b);
  51       depth8    = (int8_t)(d / ONE_PLY);
  52   }
  53 }
  54
  55
  56 /// TranspositionTable::resize() sets the size of the transposition table,
  57 /// measured in megabytes. Transposition table consists of a power of 2 number
  58 /// of clusters and each cluster consists of ClusterSize number of TTEntry.
  59
  60 void TranspositionTable::resize(size_t mbSize) {
  61
  62   Threads.main()->wait_for_search_finished();
  63
  64   clusterCount = mbSize * 1024 * 1024 / sizeof(Cluster);
  65
  66   free(mem);
  67   mem = malloc(clusterCount * sizeof(Cluster) + CacheLineSize - 1);
  68
  69   if (!mem)
  70   {
  71       std::cerr << "Failed to allocate " << mbSize
  72                 << "MB for transposition table." << std::endl;
  73       exit(EXIT_FAILURE);
  74   }
  75
  76   table = (Cluster*)((uintptr_t(mem) + CacheLineSize - 1) & ~(CacheLineSize - 1));
  77   clear();
  78 }
  79
  80
  81 /// TranspositionTable::clear() initializes the entire transposition table to zero,
  82 //  in a multi-threaded way.
  83
  84 void TranspositionTable::clear() {
  85
  86   std::vector<std::thread> threads;
  87
  88   for (size_t idx = 0; idx < Options["Threads"]; ++idx)
  89   {
  90       threads.emplace_back([this, idx]() {
  91
  92           // Thread binding gives faster search on systems with a first-touch policy
  93           if (Options["Threads"] > 8)
  94               WinProcGroup::bindThisThread(idx);
  95
  96           // Each thread will zero its part of the hash table
  97           const size_t stride = clusterCount / Options["Threads"],
  98                        start  = stride * idx,
  99                        len    = idx != Options["Threads"] - 1 ?
 100                                 stride : clusterCount - start;
 101
 102           std::memset(&table[start], 0, len * sizeof(Cluster));
 103       });
 104   }
 105
 106   for (std::thread& th: threads)
 107       th.join();
 108 }
 109
 110 /// TranspositionTable::probe() looks up the current position in the transposition
 111 /// table. It returns true and a pointer to the TTEntry if the position is found.
 112 /// Otherwise, it returns false and a pointer to an empty or least valuable TTEntry
 113 /// to be replaced later. The replace value of an entry is calculated as its depth
 114 /// minus 8 times its relative age. TTEntry t1 is considered more valuable than
 115 /// TTEntry t2 if its replace value is greater than that of t2.
 116
 117 TTEntry* TranspositionTable::probe(const Key key, bool& found) const {
 118
 119   TTEntry* const tte = first_entry(key);
 120   const uint16_t key16 = key >> 48;  // Use the high 16 bits as key inside the cluster
 121
 122   for (int i = 0; i < ClusterSize; ++i)
 123       if (!tte[i].key16 || tte[i].key16 == key16)
 124       {
 125           tte[i].genBound8 = uint8_t(generation8 | tte[i].pv_hit() << 2 | tte[i].bound()); // Refresh
 126
 127           return found = (bool)tte[i].key16, &tte[i];
 128       }
 129
 130   // Find an entry to be replaced according to the replacement strategy
 131   TTEntry* replace = tte;
 132   for (int i = 1; i < ClusterSize; ++i)
 133       // Due to our packed storage format for generation and its cyclic
 134       // nature we add 263 (263 is the modulus plus 7 to keep the lowest
 135       // two bound bits from affecting the result) to calculate the entry
 136       // age correctly even after generation8 overflows into the next cycle.
 137       if (  replace->depth8 - ((263 + generation8 - replace->genBound8) & 0xF8)
 138           >   tte[i].depth8 - ((263 + generation8 -   tte[i].genBound8) & 0xF8))
 139           replace = &tte[i];
 140
 141   return found = false, replace;
 142 }
 143
 144
 145 /// TranspositionTable::hashfull() returns an approximation of the hashtable
 146 /// occupation during a search. The hash is x permill full, as per UCI protocol.
 147
 148 int TranspositionTable::hashfull() const {
 149
 150   int cnt = 0;
 151   for (int i = 0; i < 1000 / ClusterSize; ++i)
 152       for (int j = 0; j < ClusterSize; ++j)
 153           cnt += (table[i].entry[j].genBound8 & 0xFC) == generation8;
 154
 155   return cnt * 1000 / (ClusterSize * (1000 / ClusterSize));
 156 }