Clean-up some shifting in space calculation (#1955)
[stockfish] / 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].genBound8 & 0x7)); // 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 (256 is the modulus plus 7 to keep the unrelated
135       // lowest three 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 & 0xF8) == generation8;
154
155   return cnt * 1000 / (ClusterSize * (1000 / ClusterSize));
156 }