]> git.sesse.net Git - stockfish/blobdiff - src/search.cpp
Reformat sending of PV information
[stockfish] / src / search.cpp
index ba31bd0221bc3a811a861c99c56d7c1df0025b2b..679ebbbd1e457a8c111bfdde2aa3b75ed1202af6 100644 (file)
@@ -17,6 +17,7 @@
   along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */
 
+#include <algorithm>
 #include <cassert>
 #include <cmath>
 #include <cstring>
@@ -24,7 +25,6 @@
 #include <iostream>
 #include <sstream>
 #include <vector>
-#include <algorithm>
 
 #include "book.h"
 #include "evaluate.h"
@@ -42,7 +42,7 @@ namespace Search {
 
   volatile SignalsType Signals;
   LimitsType Limits;
-  std::vector<Move> RootMoves;
+  std::vector<Move> SearchMoves;
   Position RootPosition;
 }
 
@@ -60,15 +60,21 @@ namespace {
   enum NodeType { Root, PV, NonPV, SplitPointRoot, SplitPointPV, SplitPointNonPV };
 
   // RootMove struct is used for moves at the root of the tree. For each root
-  // move, we store a score, a node count, and a PV (really a refutation
-  // in the case of moves which fail low). Score is normally set at
-  // -VALUE_INFINITE for all non-pv moves.
+  // move we store a score, a node count, and a PV (really a refutation in the
+  // case of moves which fail low). Score is normally set at -VALUE_INFINITE for
+  // all non-pv moves.
   struct RootMove {
 
-    // RootMove::operator<() is the comparison function used when
-    // sorting the moves. A move m1 is considered to be better
-    // than a move m2 if it has an higher score
+    RootMove(){}
+    RootMove(Move m) {
+      nodes = 0;
+      score = prevScore = -VALUE_INFINITE;
+      pv.push_back(m);
+      pv.push_back(MOVE_NONE);
+    }
+
     bool operator<(const RootMove& m) const { return score < m.score; }
+    bool operator==(const Move& m) const { return pv[0] == m; }
 
     void extract_pv_from_tt(Position& pos);
     void insert_pv_in_tt(Position& pos);
@@ -79,15 +85,6 @@ namespace {
     std::vector<Move> pv;
   };
 
-  // RootMoveList struct is mainly a std::vector of RootMove objects
-  struct RootMoveList : public std::vector<RootMove> {
-
-    void init(Position& pos, Move rootMoves[]);
-    RootMove* find(const Move& m, int startIndex = 0);
-
-    int bestMoveChanges;
-  };
-
 
   /// Constants
 
@@ -147,26 +144,18 @@ namespace {
 
   /// Namespace variables
 
-  // Root move list
-  RootMoveList Rml;
-
-  // MultiPV mode
+  std::vector<RootMove> RootMoves;
   size_t MultiPV, UCIMultiPV, MultiPVIdx;
-
-  // Time management variables
   TimeManager TimeMgr;
-
-  // Skill level adjustment
+  int BestMoveChanges;
   int SkillLevel;
   bool SkillLevelEnabled;
-
-  // History table
   History H;
 
 
   /// Local functions
 
-  Move id_loop(Position& pos, Move rootMoves[], Move* ponderMove);
+  Move id_loop(Position& pos, Move* ponderMove);
 
   template <NodeType NT>
   Value search(Position& pos, Stack* ss, Value alpha, Value beta, Depth depth);
@@ -183,17 +172,14 @@ namespace {
   Value refine_eval(const TTEntry* tte, Value defaultEval, int ply);
   void update_history(const Position& pos, Move move, Depth depth, Move movesSearched[], int moveCount);
   void do_skill_level(Move* best, Move* ponder);
-
   int elapsed_time(bool reset = false);
   string score_to_uci(Value v, Value alpha = -VALUE_INFINITE, Value beta = VALUE_INFINITE);
-  string speed_to_uci(int64_t nodes);
-  string pv_to_uci(const Move pv[], int pvNum, bool chess960);
-  string pretty_pv(Position& pos, int depth, Value score, int time, Move pv[]);
-  string depth_to_uci(Depth depth);
-
-  // MovePickerExt template class extends MovePicker and allows to choose at compile
-  // time the proper moves source according to the type of node. In the default case
-  // we simply create and use a standard MovePicker object.
+  void pv_info_to_log(Position& pos, int depth, Value score, int time, Move pv[]);
+  void pv_info_to_uci(const Position& pos, int depth, Value alpha, Value beta);
+
+  // MovePickerExt class template extends MovePicker and allows to choose at
+  // compile time the proper moves source according to the type of node. In the
+  // default case we simply create and use a standard MovePicker object.
   template<bool SpNode> struct MovePickerExt : public MovePicker {
 
     MovePickerExt(const Position& p, Move ttm, Depth d, const History& h, Stack* ss, Value b)
@@ -218,17 +204,17 @@ namespace {
     return os << move_to_uci(m, chess960);
   }
 
-  // When formatting a move for std::cout we must know if we are in Chess960
-  // or not. To keep using the handy operator<<() on the move the trick is to
-  // embed this flag in the stream itself. Function-like named enum set960 is
-  // used as a custom manipulator and the stream internal general-purpose array,
-  // accessed through ios_base::iword(), is used to pass the flag to the move's
-  // operator<<() that will read it to properly format castling moves.
+  // When formatting a move for std::cout we must know if we are in Chess960 or
+  // not. To keep using the handy operator<<() on the move the trick is to embed
+  // this flag in the stream itself. Function-like named enum set960 is used as
+  // a custom manipulator and the stream internal general-purpose array, accessed
+  // through ios_base::iword(), is used to pass the flag to the move's operator<<
+  // that will read it to properly format castling moves.
   enum set960 {};
 
-  std::ostream& operator<< (std::ostream& os, const set960& f) {
+  std::ostream& operator<<(std::ostream& os, const set960& f) {
 
-    os.iword(0) = int(f);
+    os.iword(0) = f;
     return os;
   }
 
@@ -294,15 +280,12 @@ int64_t Search::perft(Position& pos, Depth depth) {
   StateInfo st;
   int64_t sum = 0;
 
-  // Generate all legal moves
   MoveList<MV_LEGAL> ml(pos);
 
-  // If we are at the last ply we don't need to do and undo
-  // the moves, just to count them.
+  // At the last ply just return the number of moves (leaf nodes)
   if (depth <= ONE_PLY)
       return ml.size();
 
-  // Loop through all legal moves
   CheckInfo ci(pos);
   for ( ; !ml.end(); ++ml)
   {
@@ -323,14 +306,12 @@ void Search::think() {
   static Book book; // Defined static to initialize the PRNG only once
 
   Position& pos = RootPosition;
-
-  // Reset elapsed search time
   elapsed_time(true);
+  TimeMgr.init(Limits, pos.startpos_ply_counter());
 
   // Set output stream mode: normal or chess960. Castling notation is different
   cout << set960(pos.is_chess960());
 
-  // Look for a book move
   if (Options["OwnBook"].value<bool>())
   {
       if (Options["Book File"].value<string>() != book.name())
@@ -351,9 +332,7 @@ void Search::think() {
   read_evaluation_uci_options(pos.side_to_move());
   Threads.read_uci_options();
 
-  // Set a new TT size if changed
   TT.set_size(Options["Hash"].value<int>());
-
   if (Options["Clear Hash"].value<bool>())
   {
       Options["Clear Hash"].set_value("false");
@@ -361,14 +340,13 @@ void Search::think() {
   }
 
   UCIMultiPV = Options["MultiPV"].value<size_t>();
-  SkillLevel = Options["Skill Level"].value<size_t>();
+  SkillLevel = Options["Skill Level"].value<int>();
 
   // Do we have to play with skill handicap? In this case enable MultiPV that
   // we will use behind the scenes to retrieve a set of possible moves.
   SkillLevelEnabled = (SkillLevel < 20);
   MultiPV = (SkillLevelEnabled ? std::max(UCIMultiPV, (size_t)4) : UCIMultiPV);
 
-  // Write current search header to log file
   if (Options["Use Search Log"].value<bool>())
   {
       Log log(Options["Search Log Filename"].value<string>());
@@ -390,8 +368,6 @@ void Search::think() {
 
   // Set best timer interval to avoid lagging under time pressure. Timer is
   // used to check for remaining available thinking time.
-   TimeMgr.init(Limits, pos.startpos_ply_counter());
-
   if (TimeMgr.available_time())
       Threads.set_timer(std::min(100, std::max(TimeMgr.available_time() / 8, 20)));
   else
@@ -399,15 +375,12 @@ void Search::think() {
 
   // We're ready to start thinking. Call the iterative deepening loop function
   Move ponderMove = MOVE_NONE;
-  Move bestMove = id_loop(pos, &RootMoves[0], &ponderMove);
+  Move bestMove = id_loop(pos, &ponderMove);
 
-  // Stop timer, no need to check for available time any more
+  // Stop timer and send all the slaves to sleep, if not already sleeping
   Threads.set_timer(0);
-
-  // This makes all the slave threads to go to sleep, if not already sleeping
   Threads.set_size(1);
 
-  // Write current search final statistics to log file
   if (Options["Use Search Log"].value<bool>())
   {
       int e = elapsed_time();
@@ -447,32 +420,33 @@ namespace {
   // with increasing depth until the allocated thinking time has been consumed,
   // user stops the search, or the maximum search depth is reached.
 
-  Move id_loop(Position& pos, Move rootMoves[], Move* ponderMove) {
+  Move id_loop(Position& pos, Move* ponderMove) {
 
     Stack ss[PLY_MAX_PLUS_2];
-    Value bestValues[PLY_MAX_PLUS_2];
     int bestMoveChanges[PLY_MAX_PLUS_2];
-    int depth, aspirationDelta;
-    Value bestValue, alpha, beta;
+    int depth;
+    Value bestValue, alpha, beta, delta;
     Move bestMove, skillBest, skillPonder;
     bool bestMoveNeverChanged = true;
 
-    // Initialize stuff before a new search
     memset(ss, 0, 4 * sizeof(Stack));
     TT.new_search();
     H.clear();
+    RootMoves.clear();
     *ponderMove = bestMove = skillBest = skillPonder = MOVE_NONE;
-    depth = aspirationDelta = 0;
-    bestValue = alpha = -VALUE_INFINITE, beta = VALUE_INFINITE;
+    depth = 0;
+    bestValue = alpha = -VALUE_INFINITE, beta = delta = VALUE_INFINITE;
     ss->currentMove = MOVE_NULL; // Hack to skip update gains
 
-    // Moves to search are verified and copied
-    Rml.init(pos, rootMoves);
+    for (MoveList<MV_LEGAL> ml(pos); !ml.end(); ++ml)
+        if (   SearchMoves.empty()
+            || std::count(SearchMoves.begin(), SearchMoves.end(), ml.move()))
+            RootMoves.push_back(RootMove(ml.move()));
 
     // Handle special case of searching on a mate/stalemate position
-    if (Rml.empty())
+    if (RootMoves.empty())
     {
-        cout << "info" << depth_to_uci(DEPTH_ZERO)
+        cout << "info depth 0"
              << score_to_uci(pos.in_check() ? -VALUE_MATE : VALUE_DRAW, alpha, beta) << endl;
 
         return MOVE_NONE;
@@ -482,25 +456,20 @@ namespace {
     while (!Signals.stop && ++depth <= PLY_MAX && (!Limits.maxDepth || depth <= Limits.maxDepth))
     {
         // Save now last iteration's scores, before Rml moves are reordered
-        for (size_t i = 0; i < Rml.size(); i++)
-            Rml[i].prevScore = Rml[i].score;
+        for (size_t i = 0; i < RootMoves.size(); i++)
+            RootMoves[i].prevScore = RootMoves[i].score;
 
-        Rml.bestMoveChanges = 0;
+        BestMoveChanges = 0;
 
         // MultiPV loop. We perform a full root search for each PV line
-        for (MultiPVIdx = 0; MultiPVIdx < std::min(MultiPV, Rml.size()); MultiPVIdx++)
+        for (MultiPVIdx = 0; MultiPVIdx < std::min(MultiPV, RootMoves.size()); MultiPVIdx++)
         {
-            // Calculate dynamic aspiration window based on previous iterations
-            if (depth >= 5 && abs(Rml[MultiPVIdx].prevScore) < VALUE_KNOWN_WIN)
+            // Aspiration window
+            if (depth >= 5 && abs(RootMoves[MultiPVIdx].prevScore) < VALUE_KNOWN_WIN)
             {
-                int prevDelta1 = bestValues[depth - 1] - bestValues[depth - 2];
-                int prevDelta2 = bestValues[depth - 2] - bestValues[depth - 3];
-
-                aspirationDelta = std::min(std::max(abs(prevDelta1) + abs(prevDelta2) / 2, 16), 24);
-                aspirationDelta = (aspirationDelta + 7) / 8 * 8; // Round to match grainSize
-
-                alpha = std::max(Rml[MultiPVIdx].prevScore - aspirationDelta, -VALUE_INFINITE);
-                beta  = std::min(Rml[MultiPVIdx].prevScore + aspirationDelta,  VALUE_INFINITE);
+                delta = Value(16);
+                alpha = RootMoves[MultiPVIdx].prevScore - delta;
+                beta  = RootMoves[MultiPVIdx].prevScore + delta;
             }
             else
             {
@@ -521,19 +490,19 @@ namespace {
                 // we want to keep the same order for all the moves but the new
                 // PV that goes to the front. Note that in case of MultiPV search
                 // the already searched PV lines are preserved.
-                sort<RootMove>(Rml.begin() + MultiPVIdx, Rml.end());
+                sort<RootMove>(RootMoves.begin() + MultiPVIdx, RootMoves.end());
 
                 // In case we have found an exact score and we are going to leave
                 // the fail high/low loop then reorder the PV moves, otherwise
                 // leave the last PV move in its position so to be searched again.
                 // Of course this is needed only in MultiPV search.
                 if (MultiPVIdx && bestValue > alpha && bestValue < beta)
-                    sort<RootMove>(Rml.begin(), Rml.begin() + MultiPVIdx);
+                    sort<RootMove>(RootMoves.begin(), RootMoves.begin() + MultiPVIdx);
 
                 // Write PV back to transposition table in case the relevant entries
                 // have been overwritten during the search.
                 for (size_t i = 0; i <= MultiPVIdx; i++)
-                    Rml[i].insert_pv_in_tt(pos);
+                    RootMoves[i].insert_pv_in_tt(pos);
 
                 // If search has been stopped exit the aspiration window loop,
                 // note that sorting and writing PV back to TT is safe becuase
@@ -542,64 +511,43 @@ namespace {
                     break;
 
                 // Send full PV info to GUI if we are going to leave the loop or
-                // if we have a fail high/low and we are deep in the search. UCI
-                // protocol requires to send all the PV lines also if are still
-                // to be searched and so refer to the previous search's score.
+                // if we have a fail high/low and we are deep in the search.
                 if ((bestValue > alpha && bestValue < beta) || elapsed_time() > 2000)
-                    for (size_t i = 0; i < std::min(UCIMultiPV, Rml.size()); i++)
-                    {
-                        bool updated = (i <= MultiPVIdx);
-
-                        if (depth == 1 && !updated)
-                            continue;
-
-                        Depth d = (updated ? depth : depth - 1) * ONE_PLY;
-                        Value s = (updated ? Rml[i].score : Rml[i].prevScore);
-
-                        cout << "info"
-                             << depth_to_uci(d)
-                             << (i == MultiPVIdx ? score_to_uci(s, alpha, beta) : score_to_uci(s))
-                             << speed_to_uci(pos.nodes_searched())
-                             << pv_to_uci(&Rml[i].pv[0], i + 1, pos.is_chess960())
-                             << endl;
-                    }
+                    pv_info_to_uci(pos, depth, alpha, beta);
 
                 // In case of failing high/low increase aspiration window and
                 // research, otherwise exit the fail high/low loop.
                 if (bestValue >= beta)
                 {
-                    beta = std::min(beta + aspirationDelta, VALUE_INFINITE);
-                    aspirationDelta += aspirationDelta / 2;
+                    beta += delta;
+                    delta += delta / 2;
                 }
                 else if (bestValue <= alpha)
                 {
                     Signals.failedLowAtRoot = true;
                     Signals.stopOnPonderhit = false;
 
-                    alpha = std::max(alpha - aspirationDelta, -VALUE_INFINITE);
-                    aspirationDelta += aspirationDelta / 2;
+                    alpha -= delta;
+                    delta += delta / 2;
                 }
                 else
                     break;
 
+                assert(alpha >= -VALUE_INFINITE && beta <= VALUE_INFINITE);
+
             } while (abs(bestValue) < VALUE_KNOWN_WIN);
         }
 
-        // Collect info about search result
-        bestMove = Rml[0].pv[0];
-        *ponderMove = Rml[0].pv[1];
-        bestValues[depth] = bestValue;
-        bestMoveChanges[depth] = Rml.bestMoveChanges;
+        bestMove = RootMoves[0].pv[0];
+        *ponderMove = RootMoves[0].pv[1];
+        bestMoveChanges[depth] = BestMoveChanges;
 
         // Skills: Do we need to pick now the best and the ponder moves ?
         if (SkillLevelEnabled && depth == 1 + SkillLevel)
             do_skill_level(&skillBest, &skillPonder);
 
         if (Options["Use Search Log"].value<bool>())
-        {
-            Log log(Options["Search Log Filename"].value<string>());
-            log << pretty_pv(pos, depth, bestValue, elapsed_time(), &Rml[0].pv[0]) << endl;
-        }
+             pv_info_to_log(pos, depth, bestValue, elapsed_time(), &RootMoves[0].pv[0]);
 
         // Filter out startup noise when monitoring best move stability
         if (depth > 2 && bestMoveChanges[depth])
@@ -743,7 +691,7 @@ namespace {
     excludedMove = ss->excludedMove;
     posKey = excludedMove ? pos.get_exclusion_key() : pos.get_key();
     tte = TT.probe(posKey);
-    ttMove = RootNode ? Rml[MultiPVIdx].pv[0] : tte ? tte->move() : MOVE_NONE;
+    ttMove = RootNode ? RootMoves[MultiPVIdx].pv[0] : tte ? tte->move() : MOVE_NONE;
 
     // At PV nodes we check for exact scores, while at non-PV nodes we check for
     // a fail high/low. Biggest advantage at probing at PV nodes is to have a
@@ -932,7 +880,6 @@ namespace {
 
 split_point_start: // At split points actual search starts from here
 
-    // Initialize a MovePicker object for the current position
     MovePickerExt<SpNode> mp(pos, ttMove, depth, H, ss, PvNode ? -VALUE_INFINITE : beta);
     CheckInfo ci(pos);
     ss->bestMove = MOVE_NONE;
@@ -941,7 +888,7 @@ split_point_start: // At split points actual search starts from here
                            && !SpNode
                            && depth >= SingularExtensionDepth[PvNode]
                            && ttMove != MOVE_NONE
-                           && !excludedMove // Do not allow recursive singular extension search
+                           && !excludedMove // Recursive singular search is not allowed
                            && (tte->type() & VALUE_TYPE_LOWER)
                            && tte->depth() >= depth - 3 * ONE_PLY;
     if (SpNode)
@@ -967,7 +914,7 @@ split_point_start: // At split points actual search starts from here
       // At root obey the "searchmoves" option and skip moves not listed in Root
       // Move List, as a consequence any illegal move is also skipped. In MultiPV
       // mode we also skip PV moves which have been already searched.
-      if (RootNode && !Rml.find(move, MultiPVIdx))
+      if (RootNode && !std::count(RootMoves.begin() + MultiPVIdx, RootMoves.end(), move))
           continue;
 
       // At PV and SpNode nodes we want all moves to be legal since the beginning
@@ -987,12 +934,10 @@ split_point_start: // At split points actual search starts from here
           // This is used by time management
           Signals.firstRootMove = (moveCount == 1);
 
-          // Save the current node count before the move is searched
           nodes = pos.nodes_searched();
 
-          // For long searches send current move info to GUI
           if (pos.thread() == 0 && elapsed_time() > 2000)
-              cout << "info" << depth_to_uci(depth)
+              cout << "info depth " << depth / ONE_PLY
                    << " currmove " << move
                    << " currmovenumber " << moveCount + MultiPVIdx << endl;
       }
@@ -1155,30 +1100,28 @@ split_point_start: // At split points actual search starts from here
       // be trusted, and we don't update the best move and/or PV.
       if (RootNode && !Signals.stop)
       {
-          // Remember searched nodes counts for this move
-          RootMove* rm = Rml.find(move);
-          rm->nodes += pos.nodes_searched() - nodes;
+          RootMove& rm = *std::find(RootMoves.begin(), RootMoves.end(), move);
+          rm.nodes += pos.nodes_searched() - nodes;
 
           // PV move or new best move ?
           if (isPvMove || value > alpha)
           {
-              // Update PV
-              rm->score = value;
-              rm->extract_pv_from_tt(pos);
+              rm.score = value;
+              rm.extract_pv_from_tt(pos);
 
               // We record how often the best move has been changed in each
               // iteration. This information is used for time management: When
               // the best move changes frequently, we allocate some more time.
               if (!isPvMove && MultiPV == 1)
-                  Rml.bestMoveChanges++;
+                  BestMoveChanges++;
           }
           else
               // All other moves but the PV are set to the lowest value, this
               // is not a problem when sorting becuase sort is stable and move
               // position in the list is preserved, just the PV is pushed up.
-              rm->score = -VALUE_INFINITE;
+              rm.score = -VALUE_INFINITE;
 
-      } // RootNode
+      }
 
       if (value > bestValue)
       {
@@ -1219,7 +1162,7 @@ split_point_start: // At split points actual search starts from here
     if (!moveCount)
         return excludedMove ? oldAlpha : inCheck ? value_mated_in(ss->ply) : VALUE_DRAW;
 
-    // We have pruned all the moves, so return a fail-low score
+    // If we have pruned all the moves without searching return a fail-low score
     if (bestValue == -VALUE_INFINITE)
     {
         assert(!playedMoveCount);
@@ -1228,8 +1171,7 @@ split_point_start: // At split points actual search starts from here
     }
 
     // Step 21. Update tables
-    // If the search is not aborted, update the transposition table,
-    // history counters, and killer moves.
+    // Update transposition table entry, history and killers
     if (!SpNode && !Signals.stop && !thread.cutoff_occurred())
     {
         move = bestValue <= oldAlpha ? MOVE_NONE : ss->bestMove;
@@ -1264,6 +1206,7 @@ split_point_start: // At split points actual search starts from here
     return bestValue;
   }
 
+
   // qsearch() is the quiescence search function, which is called by the main
   // search function when the remaining depth is zero (or, to be more precise,
   // less than ONE_PLY).
@@ -1344,7 +1287,6 @@ split_point_start: // At split points actual search starts from here
         if (PvNode && bestValue > alpha)
             alpha = bestValue;
 
-        // Futility pruning parameters, not needed when in check
         futilityBase = ss->eval + evalMargin + FutilityMarginQS;
         enoughMaterial = pos.non_pawn_material(pos.side_to_move()) > RookValueMidgame;
     }
@@ -1426,7 +1368,6 @@ split_point_start: // At split points actual search starts from here
       if (!pos.pl_move_is_legal(move, ci.pinned))
           continue;
 
-      // Update current move
       ss->currentMove = move;
 
       // Make and search the move
@@ -1647,8 +1588,8 @@ split_point_start: // At split points actual search starts from here
   }
 
 
-  // can_return_tt() returns true if a transposition table score
-  // can be used to cut-off at a given point in search.
+  // can_return_tt() returns true if a transposition table score can be used to
+  // cut-off at a given point in search.
 
   bool can_return_tt(const TTEntry* tte, Depth depth, Value beta, int ply) {
 
@@ -1663,8 +1604,8 @@ split_point_start: // At split points actual search starts from here
   }
 
 
-  // refine_eval() returns the transposition table score if
-  // possible otherwise falls back on static position evaluation.
+  // refine_eval() returns the transposition table score if possible, otherwise
+  // falls back on static position evaluation.
 
   Value refine_eval(const TTEntry* tte, Value defaultEval, int ply) {
 
@@ -1680,8 +1621,8 @@ split_point_start: // At split points actual search starts from here
   }
 
 
-  // update_history() registers a good move that produced a beta-cutoff
-  // in history and marks as failures all the other moves of that ply.
+  // update_history() registers a good move that produced a beta-cutoff in
+  // history and marks as failures all the other moves of that ply.
 
   void update_history(const Position& pos, Move move, Depth depth,
                       Move movesSearched[], int moveCount) {
@@ -1737,55 +1678,49 @@ split_point_start: // At split points actual search starts from here
   }
 
 
-  // speed_to_uci() returns a string with time stats of current search suitable
-  // to be sent to UCI gui.
+  // pv_info_to_uci() sends search info to GUI. UCI protocol requires to send all
+  // the PV lines also if are still to be searched and so refer to the previous
+  // search score.
 
-  string speed_to_uci(int64_t nodes) {
+  void pv_info_to_uci(const Position& pos, int depth, Value alpha, Value beta) {
 
-    std::stringstream s;
     int t = elapsed_time();
+    int selDepth = 0;
 
-    s << " nodes " << nodes
-      << " nps " << (t > 0 ? int(nodes * 1000 / t) : 0)
-      << " time "  << t;
-
-    return s.str();
-  }
-
+    for (int i = 0; i < Threads.size(); i++)
+        if (Threads[i].maxPly > selDepth)
+            selDepth = Threads[i].maxPly;
 
-  // pv_to_uci() returns a string with information on the current PV line
-  // formatted according to UCI specification.
+    for (size_t i = 0; i < std::min(UCIMultiPV, RootMoves.size()); i++)
+    {
+        bool updated = (i <= MultiPVIdx);
 
-  string pv_to_uci(const Move pv[], int pvNum, bool chess960) {
+        if (depth == 1 && !updated)
+            continue;
 
-    std::stringstream s;
+        int d = (updated ? depth : depth - 1);
+        Value s = (updated ? RootMoves[i].score : RootMoves[i].prevScore);
 
-    s << " multipv " << pvNum << " pv " << set960(chess960);
+        cout << "info"
+             << " depth " << d
+             << " seldepth " << selDepth
+             << (i == MultiPVIdx ? score_to_uci(s, alpha, beta) : score_to_uci(s))
+             << " nodes " << pos.nodes_searched()
+             << " nps " << (t > 0 ? pos.nodes_searched() * 1000 / t : 0)
+             << " time " << t
+             << " multipv " << i + 1 << " pv";
 
-    for ( ; *pv != MOVE_NONE; pv++)
-        s << *pv << " ";
+        for (int j = 0; RootMoves[i].pv[j] != MOVE_NONE; j++)
+            cout <<  " " << RootMoves[i].pv[j];
 
-    return s.str();
+        cout << endl;
+    }
   }
 
 
-  // depth_to_uci() returns a string with information on the current depth and
-  // seldepth formatted according to UCI specification.
-
-  string depth_to_uci(Depth depth) {
-
-    std::stringstream s;
-
-    // Retrieve max searched depth among threads
-    int selDepth = 0;
-    for (int i = 0; i < Threads.size(); i++)
-        if (Threads[i].maxPly > selDepth)
-            selDepth = Threads[i].maxPly;
-
-     s << " depth " << depth / ONE_PLY << " seldepth " << selDepth;
-
-    return s.str();
-  }
+  // pv_info_to_log() writes human-readable search information to the log file
+  // (which is created when the UCI parameter "Use Search Log" is "true"). It
+  // uses the two below helpers to pretty format time and score respectively.
 
   string time_to_string(int millisecs) {
 
@@ -1801,7 +1736,8 @@ split_point_start: // At split points actual search starts from here
     if (hours)
         s << hours << ':';
 
-    s << std::setfill('0') << std::setw(2) << minutes << ':' << std::setw(2) << seconds;
+    s << std::setfill('0') << std::setw(2) << minutes << ':'
+                           << std::setw(2) << seconds;
     return s.str();
   }
 
@@ -1814,30 +1750,23 @@ split_point_start: // At split points actual search starts from here
     else if (v <= VALUE_MATED_IN_PLY_MAX)
         s << "-#" << (VALUE_MATE + v) / 2;
     else
-        s << std::setprecision(2) << std::fixed << std::showpos << float(v) / PawnValueMidgame;
+        s << std::setprecision(2) << std::fixed << std::showpos
+          << float(v) / PawnValueMidgame;
 
     return s.str();
   }
 
-
-  // pretty_pv() creates a human-readable string from a position and a PV.
-  // It is used to write search information to the log file (which is created
-  // when the UCI parameter "Use Search Log" is "true").
-
-  string pretty_pv(Position& pos, int depth, Value value, int time, Move pv[]) {
+  void pv_info_to_log(Position& pos, int depth, Value value, int time, Move pv[]) {
 
     const int64_t K = 1000;
     const int64_t M = 1000000;
-    const int startColumn = 28;
-    const size_t maxLength = 80 - startColumn;
 
     StateInfo state[PLY_MAX_PLUS_2], *st = state;
     Move* m = pv;
-    string san;
+    string san, padding;
+    size_t length;
     std::stringstream s;
-    size_t length = 0;
 
-    // First print depth, score, time and searched nodes...
     s << set960(pos.is_chess960())
       << std::setw(2) << depth
       << std::setw(8) << score_to_string(value)
@@ -1845,31 +1774,37 @@ split_point_start: // At split points actual search starts from here
 
     if (pos.nodes_searched() < M)
         s << std::setw(8) << pos.nodes_searched() / 1 << "  ";
+
     else if (pos.nodes_searched() < K * M)
         s << std::setw(7) << pos.nodes_searched() / K << "K  ";
+
     else
         s << std::setw(7) << pos.nodes_searched() / M << "M  ";
 
-    // ...then print the full PV line in short algebraic notation
+    padding = string(s.str().length(), ' ');
+    length = padding.length();
+
     while (*m != MOVE_NONE)
     {
         san = move_to_san(pos, *m);
-        length += san.length() + 1;
 
-        if (length > maxLength)
+        if (length + san.length() > 80)
         {
-            length = san.length() + 1;
-            s << "\n" + string(startColumn, ' ');
+            s << "\n" + padding;
+            length = padding.length();
         }
+
         s << san << ' ';
+        length += san.length() + 1;
 
         pos.do_move(*m++, *st++);
     }
 
-    // Restore original position before to leave
-    while (m != pv) pos.undo_move(*--m);
+    while (m != pv)
+        pos.undo_move(*--m);
 
-    return s.str();
+    Log l(Options["Search Log Filename"].value<string>());
+    l << s.str() << endl;
   }
 
 
@@ -1882,79 +1817,41 @@ split_point_start: // At split points actual search starts from here
 
     static RKISS rk;
 
-    // Rml list is already sorted by score in descending order
-    int s;
-    size_t size = std::min(MultiPV, Rml.size());
-    int max_s = -VALUE_INFINITE;
-    int max = Rml[0].score;
-    int var = std::min(max - Rml[size - 1].score, int(PawnValueMidgame));
-    int wk = 120 - 2 * SkillLevel;
-
-    // PRNG sequence should be non deterministic
+    // PRNG sequence should be not deterministic
     for (int i = abs(get_system_time() % 50); i > 0; i--)
         rk.rand<unsigned>();
 
-    // Choose best move. For each move's score we add two terms both dependent
-    // on wk, one deterministic and bigger for weaker moves, and one random,
+    // Rml list is already sorted by score in descending order
+    size_t size = std::min(MultiPV, RootMoves.size());
+    int variance = std::min(RootMoves[0].score - RootMoves[size - 1].score, PawnValueMidgame);
+    int weakness = 120 - 2 * SkillLevel;
+    int max_s = -VALUE_INFINITE;
+
+    // Choose best move. For each move score we add two terms both dependent on
+    // weakness, one deterministic and bigger for weaker moves, and one random,
     // then we choose the move with the resulting highest score.
     for (size_t i = 0; i < size; i++)
     {
-        s = Rml[i].score;
+        int s = RootMoves[i].score;
 
         // Don't allow crazy blunders even at very low skills
-        if (i > 0 && Rml[i-1].score > s + EasyMoveMargin)
+        if (i > 0 && RootMoves[i-1].score > s + EasyMoveMargin)
             break;
 
-        // This is our magical formula
-        s += ((max - s) * wk + var * (rk.rand<unsigned>() % wk)) / 128;
+        // This is our magic formula
+        s += (  weakness * int(RootMoves[0].score - s)
+              + variance * (rk.rand<unsigned>() % weakness)) / 128;
 
         if (s > max_s)
         {
             max_s = s;
-            *best = Rml[i].pv[0];
-            *ponder = Rml[i].pv[1];
+            *best = RootMoves[i].pv[0];
+            *ponder = RootMoves[i].pv[1];
         }
     }
   }
 
 
-  /// RootMove and RootMoveList method's definitions
-
-  void RootMoveList::init(Position& pos, Move rootMoves[]) {
-
-    Move* sm;
-    bestMoveChanges = 0;
-    clear();
-
-    // Generate all legal moves and add them to RootMoveList
-    for (MoveList<MV_LEGAL> ml(pos); !ml.end(); ++ml)
-    {
-        // If we have a rootMoves[] list then verify the move
-        // is in the list before to add it.
-        for (sm = rootMoves; *sm && *sm != ml.move(); sm++) {}
-
-        if (sm != rootMoves && *sm != ml.move())
-            continue;
-
-        RootMove rm;
-        rm.pv.push_back(ml.move());
-        rm.pv.push_back(MOVE_NONE);
-        rm.score = rm.prevScore = -VALUE_INFINITE;
-        rm.nodes = 0;
-        push_back(rm);
-    }
-  }
-
-  RootMove* RootMoveList::find(const Move& m, int startIndex) {
-
-    for (size_t i = startIndex; i < size(); i++)
-        if ((*this)[i].pv[0] == m)
-            return &(*this)[i];
-
-    return NULL;
-  }
-
-
   // extract_pv_from_tt() builds a PV by adding moves from the transposition table.
   // We consider also failing high nodes and not only VALUE_TYPE_EXACT nodes. This
   // allow to always have a ponder move even when we fail high at root and also a
@@ -2002,7 +1899,7 @@ split_point_start: // At split points actual search starts from here
     Value v, m = VALUE_NONE;
     int ply = 0;
 
-    assert(pv[0] != MOVE_NONE && pos.is_pseudo_legal(pv[0]));
+    assert(pv[ply] != MOVE_NONE && pos.is_pseudo_legal(pv[ply]));
 
     do {
         k = pos.get_key();
@@ -2024,9 +1921,9 @@ split_point_start: // At split points actual search starts from here
 } // namespace
 
 
-// Thread::idle_loop() is where the thread is parked when it has no work to do.
-// The parameter 'sp', if non-NULL, is a pointer to an active SplitPoint object
-// for which the thread is the master.
+/// Thread::idle_loop() is where the thread is parked when it has no work to do.
+/// The parameter 'sp', if non-NULL, is a pointer to an active SplitPoint object
+/// for which the thread is the master.
 
 void Thread::idle_loop(SplitPoint* sp) {
 
@@ -2040,7 +1937,6 @@ void Thread::idle_loop(SplitPoint* sp) {
       {
           assert((!sp && threadID) || Threads.use_sleeping_threads());
 
-          // Slave thread should exit as soon as do_terminate flag raises
           if (do_terminate)
           {
               assert(!sp);
@@ -2115,15 +2011,16 @@ void Thread::idle_loop(SplitPoint* sp) {
 }
 
 
-// do_timer_event() is called by the timer thread when the timer triggers
+/// do_timer_event() is called by the timer thread when the timer triggers. It
+/// is used to print debug info and, more important, to detect when we are out of
+/// available time and so stop the search.
 
 void do_timer_event() {
 
   static int lastInfoTime;
   int e = elapsed_time();
 
-  // Print debug information every one second
-  if (!lastInfoTime || get_system_time() - lastInfoTime >= 1000)
+  if (get_system_time() - lastInfoTime >= 1000 || !lastInfoTime)
   {
       lastInfoTime = get_system_time();
 
@@ -2131,7 +2028,6 @@ void do_timer_event() {
       dbg_print_hit_rate();
   }
 
-  // Should we stop the search?
   if (Limits.ponder)
       return;