Fix profile build for syzygy
[stockfish] / src / timeman.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-2014 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 <algorithm>
21 #include <cfloat>
22 #include <cmath>
23
24 #include "search.h"
25 #include "timeman.h"
26 #include "uci.h"
27
28 namespace {
29
30   enum TimeType { OptimumTime, MaxTime };
31
32   const int MoveHorizon   = 50;   // Plan time management at most this many moves ahead
33   const double MaxRatio   = 7.0;  // When in trouble, we can step over reserved time with this ratio
34   const double StealRatio = 0.33; // However we must not steal time from remaining moves over this ratio
35
36   const double xscale     = 9.3;
37   const double xshift     = 59.8;
38   const double skewfactor = 0.172;
39
40
41   // move_importance() is a skew-logistic function based on naive statistical
42   // analysis of "how many games are still undecided after n half-moves". Game
43   // is considered "undecided" as long as neither side has >275cp advantage.
44   // Data was extracted from CCRL game database with some simple filtering criteria.
45
46   double move_importance(int ply) {
47
48     return pow((1 + exp((ply - xshift) / xscale)), -skewfactor) + DBL_MIN; // Ensure non-zero
49   }
50
51   template<TimeType T>
52   int remaining(int myTime, int movesToGo, int currentPly, int slowMover)
53   {
54     const double TMaxRatio   = (T == OptimumTime ? 1 : MaxRatio);
55     const double TStealRatio = (T == OptimumTime ? 0 : StealRatio);
56
57     double thisMoveImportance = (move_importance(currentPly) * slowMover) / 100;
58     double otherMovesImportance = 0;
59
60     for (int i = 1; i < movesToGo; ++i)
61         otherMovesImportance += move_importance(currentPly + 2 * i);
62
63     double ratio1 = (TMaxRatio * thisMoveImportance) / (TMaxRatio * thisMoveImportance + otherMovesImportance);
64     double ratio2 = (thisMoveImportance + TStealRatio * otherMovesImportance) / (thisMoveImportance + otherMovesImportance);
65
66     return int(myTime * std::min(ratio1, ratio2));
67   }
68
69 } // namespace
70
71
72 void TimeManager::init(const Search::LimitsType& limits, int currentPly, Color us)
73 {
74   /* We support four different kinds of time controls:
75
76       increment == 0 && movesToGo == 0 means: x basetime  [sudden death!]
77       increment == 0 && movesToGo != 0 means: x moves in y minutes
78       increment >  0 && movesToGo == 0 means: x basetime + z increment
79       increment >  0 && movesToGo != 0 means: x moves in y minutes + z increment
80
81     Time management is adjusted by following parameters:
82
83       emergencyMoveHorizon: Be prepared to always play at least this many moves
84       emergencyBaseTime   : Always attempt to keep at least this much time (in ms) at clock
85       emergencyMoveTime   : Plus attempt to keep at least this much time for each remaining emergency move
86       minThinkingTime     : No matter what, use at least this much thinking before doing the move
87   */
88
89   int hypMTG, hypMyTime, t1, t2;
90
91   // Read uci parameters
92   int moveOverhead    = Options["Move Overhead"];
93   int minThinkingTime = Options["Minimum Thinking Time"];
94   int slowMover       = Options["Slow Mover"];
95
96   // Initialize unstablePvFactor to 1 and search times to maximum values
97   unstablePvFactor = 1;
98   optimumSearchTime = maximumSearchTime = std::max(limits.time[us], minThinkingTime);
99
100   // We calculate optimum time usage for different hypothetical "moves to go"-values and choose the
101   // minimum of calculated search time values. Usually the greatest hypMTG gives the minimum values.
102   for (hypMTG = 1; hypMTG <= (limits.movestogo ? std::min(limits.movestogo, MoveHorizon) : MoveHorizon); ++hypMTG)
103   {
104       // Calculate thinking time for hypothetical "moves to go"-value
105       hypMyTime =  limits.time[us]
106                  + limits.inc[us] * (hypMTG - 1)
107                  - moveOverhead * (2 + std::min(hypMTG, 40));
108
109       hypMyTime = std::max(hypMyTime, 0);
110
111       t1 = minThinkingTime + remaining<OptimumTime>(hypMyTime, hypMTG, currentPly, slowMover);
112       t2 = minThinkingTime + remaining<MaxTime>(hypMyTime, hypMTG, currentPly, slowMover);
113
114       optimumSearchTime = std::min(optimumSearchTime, t1);
115       maximumSearchTime = std::min(maximumSearchTime, t2);
116   }
117
118   if (Options["Ponder"])
119       optimumSearchTime += optimumSearchTime / 4;
120
121   // Make sure that maxSearchTime is not over absoluteMaxSearchTime
122   optimumSearchTime = std::min(optimumSearchTime, maximumSearchTime);
123 }