/*
	 NL reader
	
	 NL is a format for representing optimization problems such as linear,
	 quadratic, nonlinear, complementarity and constraint programming problems
	 in discrete or continuous variables. It is described in the technical report
	 "Writing .nl Files" (http://www.cs.sandia.gov/~dmgay/nlwrite.pdf).
	
	 This is a complete reusable C++ implementation of an NL reader.
	
	 Usage:
	   // Read an NL file:
	   ReadNLFile(filename, handler);
	
	   // Read an NL string:
	   ReadNLString(string, handler);
	
	 where handler is an object that receives notifications of problem
	 components. See NLHandler for an example of a handler class.
	
	 Copyright (C) 2013 - 2016 AMPL Optimization Inc
	
	 Permission to use, copy, modify, and distribute this software and its
	 documentation for any purpose and without fee is hereby granted,
	 provided that the above copyright notice appear in all copies and that
	 both that the copyright notice and this permission notice and warranty
	 disclaimer appear in supporting documentation.
	
	 The author and AMPL Optimization Inc disclaim all warranties with
	 regard to this software, including all implied warranties of
	 merchantability and fitness.  In no event shall the author be liable
	 for any special, indirect or consequential damages or any damages
	 whatsoever resulting from loss of use, data or profits, whether in an
	 action of contract, negligence or other tortious action, arising out
	 of or in connection with the use or performance of this software.
	
	 Author: Victor Zverovich
	 */
	
	#ifndef MP_NL_READER_H_
	#define MP_NL_READER_H_
	
	#include "mp/common.h"
	#include "mp/error.h"
	#include "mp/nl.h"
	#include "mp/os.h"
	
	#include <algorithm>
	#include <cctype>
	#include <cstdlib>
	#include <limits>
	#include <string>
	
	namespace mp {
	
	using fmt::internal::MakeUnsigned;
	
	/** A read error with location information. */
	class ReadError : public Error {
	 private:
	  std::string filename_;
	  int line_;
	  int column_;
	
	  void init(fmt::CStringRef filename, int line, int column,
	            fmt::CStringRef format_str, fmt::ArgList args);
	
	 public:
	  /** Constructs the exception object. */
	  ReadError(fmt::CStringRef filename, int line, int column,
	            fmt::CStringRef format_str, fmt::ArgList args) {
	    init(filename, line, column, format_str, args);
	  }
	
	  /** Destructs the exception object. */
	  ~ReadError() throw() {}
	
	  /** Returns the name of the file where error occurred. */
	  const std::string &filename() const { return filename_; }
	
	  /** Returns the line number where error occurred, starting from 1. */
	  int line() const { return line_; }
	
	  /** Returns the column number where error occurred, starting from 1. */
	  int column() const { return column_; }
	
	  FMT_VARIADIC_(char, , ReadError, init,
	                fmt::CStringRef, int, int, fmt::CStringRef)
	};
	
	/** A read error with information about offset in a binary input. */
	class BinaryReadError : public Error {
	 private:
	  std::string filename_;
	  std::size_t offset_;
	
	 public:
	  /** Constructs the exception object. */
	  BinaryReadError(
	      const std::string &filename, std::size_t offset, fmt::CStringRef message)
	  : Error(message), filename_(filename), offset_(offset) {}
	
	  /** Destructs the exception object. */
	  ~BinaryReadError() throw() {}
	
	  /** Returns the name of the file where error occurred. */
	  const std::string &filename() const { return filename_; }
	
	  /** Returns the offset in chars to the error location. */
	  std::size_t offset() const { return offset_; }
	};
	
	enum {
	  /** Maximum number of options reserved for AMPL use in NL and SOL formats. */
	  MAX_AMPL_OPTIONS = 9
	};
	
	namespace arith {
	
	/** Floating-point arithmetic kind. */
	enum Kind {
	
	  /** Unknown floating-point arithmetic. */
	  UNKNOWN = 0,
	
	  /**
	    \rst
	    Standard `IEEE-754 floating point
	    <http://en.wikipedia.org/wiki/IEEE_floating_point>`_ - little endian.
	    \endrst
	   */
	  IEEE_LITTLE_ENDIAN = 1,
	
	  /** Standard IEEE-754 floating point - big endian. */
	  IEEE_BIG_ENDIAN = 2,
	
	  /**
	    \rst
	    `IBM floating point
	    <http://en.wikipedia.org/wiki/IBM_Floating_Point_Architecture>`_.
	    \endrst
	   */
	  IBM = 3,
	
	  /** VAX floating point (legacy). */
	  VAX = 4,
	
	  /** Cray floating point. */
	  CRAY = 5,
	
	  /** Last floating point. */
	  LAST = CRAY
	};
	
	// Returns floating-point arithmetic kind used on the current system.
	Kind GetKind();
	
	inline bool IsIEEE(arith::Kind k) {
	  return k == IEEE_LITTLE_ENDIAN || k == IEEE_BIG_ENDIAN;
	}
	}  // namespace arith
	
	/**
	  \rst
	  An NL `header <http://en.wikipedia.org/wiki/Header_(computing)>`_
	  which contains information about problem dimensions, such as the number of
	  variables and constraints, and the input format.
	
	  Base class: `mp::ProblemInfo`
	  \endrst
	 */
	struct NLHeader : ProblemInfo {
	  /** Input/output format */
	  enum Format {
	    /**
	      Text format. The text format is fully portable meaning that an .nl file
	      can be written on a machine of one architecture and then read on a
	      machine of a different architecture.
	     */
	    TEXT = 0,
	
	    /**
	      Binary format. The binary format is not generally portable and should
	      normally be used on a single machine.
	     */
	    BINARY = 1
	  };
	
	  /** Input/output format. */
	  Format format;
	
	  /** The number of options reserved for AMPL use. */
	  int num_ampl_options;
	
	  /**
	    Values of options reserved for AMPL use. Leave the default values if not
	    using AMPL.
	   */
	  int ampl_options[MAX_AMPL_OPTIONS];
	
	  /**
	    Extra info for writing a solution reserved for AMPL use. Leave the default
	    value if not using AMPL.
	   */
	  double ampl_vbtol;
	
	  /**
	    \rst
	    Floating-point arithmetic kind used with binary format to check
	    if an .nl file is written using a compatible representation of
	    floating-point numbers. It is not used with the text format and normally
	    set to `mp::arith::UNKNOWN` there.
	    \endrst
	   */
	  arith::Kind arith_kind;
	
	  /** Flags. */
	  enum {
	    /** Flag that specifies whether to write output suffixes to a .sol file. */
	    WANT_OUTPUT_SUFFIXES = 1
	  };
	
	  /**
	    \rst
	    Flags. Can be either 0 or `mp::NLHeader::WANT_OUTPUT_SUFFIXES`.
	    \endrst
	   */
	  int flags;
	
	  NLHeader()
	    : ProblemInfo(), format(TEXT), num_ampl_options(0), ampl_vbtol(0),
	      arith_kind(arith::UNKNOWN), flags(0) {
	    std::fill(ampl_options, ampl_options + MAX_AMPL_OPTIONS, 0);
	  }
	};
	
	/** Writes NLHeader in the NL format. */
	fmt::Writer &operator<<(fmt::Writer &w, const NLHeader &h);
	
	/**
	  \rst
	  An NL handler.
	
	  `~mp::NLHandler` can be used as a base class for other handlers. Subclasses
	  only need to redefine methods that handle constructs they are interested
	  in and, possibly, the types used by these methods.
	
	  *Impl* is a type derived from `~mp::NLHandler` that will receive notifications
	  of unhandled constructs via `~mp::NLHandler::OnUnhandled`.
	
	  *ExprType* is a return type of expression handler methods such as
	  `~mp::NLHandler::OnUnary` useful for building expression objects. If not used
	  it can be any default-constructible type.
	  \endrst
	 */
	template <typename Impl, typename ExprType>
	class NLHandler {
	 public:
	  /** Destroys the object. */
	  virtual ~NLHandler() {}
	
	  /** An expression type. */
	  typedef ExprType Expr;
	
	  /**
	    \rst
	    .. |Expr| replace:: `Expr <mp::NLHandler::Expr>`
	
	    A numeric expression type.
	    It is a typedef of |Expr| but subclasses may define it as a different type
	    convertible to |Expr|.
	    \endrst
	   */
	  typedef Expr NumericExpr;
	
	  /**
	    \rst
	    A logical expression type.
	    It is a typedef of |Expr| but subclasses may define it as a different type
	    convertible to |Expr|.
	    \endrst
	   */
	  typedef Expr LogicalExpr;
	
	  /**
	    \rst
	    A count expression type.
	    It is a typedef of |Expr| but subclasses may define it as a different type
	    convertible to `NumericExpr <mp::NLHandler::NumericExpr>`.
	    \endrst
	   */
	  typedef Expr CountExpr;
	
	  /**
	    \rst
	    A reference expression type.
	    It is a typedef of |Expr| but subclasses may define it as a different type
	    convertible to `NumericExpr <mp::NLHandler::NumericExpr>`.
	    \endrst
	   */
	  typedef Expr Reference;
	
	  /**
	    Receives notification of an unhandled construct of the given kind.
	    Throws `UnsupportedError`.
	   */
	  void OnUnhandled(const char *kind) {
	    throw MakeUnsupportedError(kind);
	  }
	
	  /** Receives notification of an NL header. */
	  void OnHeader(const NLHeader &h) {
	    internal::Unused(&h);
	    MP_DISPATCH(OnUnhandled("NL header"));
	  }
	
	  /**
	    \rst
	    Returns true if the objective with index *obj_index* should be handled.
	    \endrst
	   */
	  bool NeedObj(int obj_index) const {
	    internal::Unused(obj_index);
	    return true;
	  }
	
	  /**
	    \rst
	    Returns final objective index
	    (e.g., returns 0 if objno>0 and multiobj=0)
	    \endrst
	   */
	  int resulting_obj_index(int obj_index) const { return obj_index; }
	
	  /**
	    Receives notification of an objective type and the nonlinear part of
	    an objective expression.
	   */
	  void OnObj(int index, obj::Type type, NumericExpr expr) {
	    internal::Unused(index, type, &expr);
	    MP_DISPATCH(OnUnhandled("objective"));
	  }
	
	  /**
	    Receives notification of the nonlinear part of an algebraic constraint
	    expression.
	   */
	  void OnAlgebraicCon(int index, NumericExpr expr) {
	    internal::Unused(index, &expr);
	    MP_DISPATCH(OnUnhandled("nonlinear constraint"));
	  }
	
	  /** Receives notification of a logical constraint expression. */
	  void OnLogicalCon(int index, LogicalExpr expr) {
	    internal::Unused(index, &expr);
	    MP_DISPATCH(OnUnhandled("logical constraint"));
	  }
	
	  /**
	    A class (struct) that receives notifications of terms in the linear part
	    of a common expression.
	   */
	  struct LinearExprHandler {
	    /** Receives notification of a term in the linear expression. */
	    void AddTerm(int var_index, double coef) {
	      internal::Unused(var_index, coef);
	    }
	  };
	
	  /**
	    Receives notification of the beginning of a common expression
	    (defined variable).
	   */
	  LinearExprHandler BeginCommonExpr(int index, int num_linear_terms) {
	    internal::Unused(index, num_linear_terms);
	    MP_DISPATCH(OnUnhandled("common expression"));
	    return LinearExprHandler();
	  }
	
	  /** Receives notification of the end of a common expression. */
	  void EndCommonExpr(int index, NumericExpr expr, int position) {
	    internal::Unused(index, &expr, position);
	    MP_DISPATCH(OnUnhandled("common expression"));
	  }
	
	  /**
	    \rst
	    Receives notification of a complementarity relation
	    ``var_lb <= x <= var_ub complements con_lb <= body <= con_ub``, where ``x``
	    is the variable at index *var_index* and ``body`` is the constraint body.
	    *info* gives the constraint bounds.
	    \endrst
	   */
	  void OnComplementarity(int con_index, int var_index, ComplInfo info) {
	    internal::Unused(con_index, var_index, &info);
	    MP_DISPATCH(OnUnhandled("complementarity constraint"));
	  }
	
	  /**
	    A typedef of a class that receives notifications of terms in the linear
	    part of an objective expression.
	   */
	  typedef LinearExprHandler LinearObjHandler;
	
	  /** Receives notification of the linear part of an objective expression. */
	  LinearObjHandler OnLinearObjExpr(int obj_index, int num_linear_terms) {
	    internal::Unused(obj_index, num_linear_terms);
	    MP_DISPATCH(OnUnhandled("linear objective"));
	    return LinearObjHandler();
	  }
	
	  /**
	    A typedef of a class that receives notifications of terms in the linear
	    part of a constraint expression.
	   */
	  typedef LinearExprHandler LinearConHandler;
	
	  /** Receives notification of the linear part of a constraint expression. */
	  LinearConHandler OnLinearConExpr(int con_index, int num_linear_terms) {
	    internal::Unused(con_index, num_linear_terms);
	    MP_DISPATCH(OnUnhandled("linear constraint"));
	    return LinearConHandler();
	  }
	
	  /** Receives notification of the linear part of a common expression. */
	  LinearExprHandler OnLinearCommonExpr(int expr_index, int num_linear_terms) {
	    internal::Unused(expr_index, num_linear_terms);
	    MP_DISPATCH(OnUnhandled("linear common expression"));
	    return LinearExprHandler();
	  }
	
	  /** Receives notification of variable bounds. */
	  void OnVarBounds(int index, double lb, double ub) {
	    internal::Unused(index, lb, ub);
	    MP_DISPATCH(OnUnhandled("variable bounds"));
	  }
	
	  /** Receives notification of constraint bounds (ranges). */
	  void OnConBounds(int index, double lb, double ub) {
	    internal::Unused(index, lb, ub);
	    MP_DISPATCH(OnUnhandled("constraint bounds"));
	  }
	
	  /** Receives notification of the initial value for a variable. */
	  void OnInitialValue(int var_index, double value) {
	    internal::Unused(var_index, value);
	    MP_DISPATCH(OnUnhandled("initial value"));
	  }
	
	  /** Receives notification of the initial value for a dual variable. */
	  void OnInitialDualValue(int con_index, double value) {
	    internal::Unused(con_index, value);
	    MP_DISPATCH(OnUnhandled("initial dual value"));
	  }
	
	  /** A class (struct) that receives notifications of Jacobian column sizes. */
	  struct ColumnSizeHandler {
	    /** Receives notification of a Jacobian column size. */
	    void Add(int size) { internal::Unused(size); }
	  };
	
	  /** Receives notification of Jacobian column sizes. */
	  ColumnSizeHandler OnColumnSizes() {
	    MP_DISPATCH(OnUnhandled("column sizes"));
	    return ColumnSizeHandler();
	  }
	
	  /**
	    \rst
	    Receives notification of a function.
	    The *name* argument is a function name and it is not null-terminated.
	    \endrst
	   */
	  void OnFunction(int index, fmt::StringRef name,
	                  int num_args, func::Type type) {
	    internal::Unused(index, &name, num_args, type);
	    MP_DISPATCH(OnUnhandled("function"));
	  }
	
	  /** A class (struct) that receives notifications of integer suffix values. */
	  struct IntSuffixHandler {
	    /** Receives notification of a suffix value. */
	    void SetValue(int index, int value) { internal::Unused(index, value); }
	  };
	
	  /**
	    \rst
	    Receives notification of an integer suffix.
	    The *name* argument is a suffix name and it is not null-terminated.
	    *kind* specifies the suffix kind.
	    \endrst
	   */
	  IntSuffixHandler OnIntSuffix(fmt::StringRef name, suf::Kind kind,
	                               int num_values) {
	    internal::Unused(&name, kind, num_values);
	    MP_DISPATCH(OnUnhandled("integer suffix"));
	    return IntSuffixHandler();
	  }
	
	  /** A class (struct) that receives notifications of double suffix values. */
	  struct DblSuffixHandler {
	    /** Receives notification of a suffix value. */
	    void SetValue(int index, double value) { internal::Unused(index, value); }
	  };
	
	  /**
	    \rst
	    Receives notification of a double suffix.
	    The *name* argument is a suffix name and it is not null-terminated.
	    *kind* specifies the suffix kind.
	    \endrst
	   */
	  DblSuffixHandler OnDblSuffix(fmt::StringRef name, suf::Kind kind,
	                               int num_values) {
	    internal::Unused(&name, kind, num_values);
	    MP_DISPATCH(OnUnhandled("double suffix"));
	    return DblSuffixHandler();
	  }
	
	  /**
	    \rst
	    A class (struct) that receives notifications of expression arguments.
	    All argument handlers in `mp::NLHandler` are typedefs of this class, but
	    subclasses of `mp::NLHandler` may define them as different classes.
	    \endrst
	   */
	  struct ArgHandler {
	    /** Receives notification of an argument. */
	    void AddArg(Expr arg) { internal::Unused(&arg); }
	  };
	
	  /**
	    \rst
	    A typedef of a class that receives notifications of `numeric
	    <mp::expr::FIRST_NUMERIC>` arguments.
	    \endrst
	   */
	  typedef ArgHandler NumericArgHandler;
	
	  /**
	    \rst
	    A typedef of a class that receives notifications of `vararg expression
	    <mp::expr::FIRST_VARARG>` arguments.
	    \endrst
	   */
	  typedef ArgHandler VarArgHandler;
	
	  /**
	    \rst
	    A typedef of a class that receives notifications of `call expression
	    <mp::expr::CALL>` arguments.
	    \endrst
	   */
	  typedef ArgHandler CallArgHandler;
	
	  /**
	    \rst
	    A typedef of a class that receives notifications of `numberof expression
	    <mp::expr::NUMBEROF>` arguments.
	    \endrst
	   */
	  typedef ArgHandler NumberOfArgHandler;
	
	  /**
	    \rst
	    A typedef of a class that receives notifications of `count expression
	    <mp::expr::COUNT>` arguments.
	    \endrst
	   */
	  typedef ArgHandler CountArgHandler;
	
	  /**
	    \rst
	    A typedef of a class that receives notifications of `logical
	    <mp::expr::FIRST_LOGICAL>` arguments.
	    \endrst
	   */
	  typedef ArgHandler LogicalArgHandler;
	
	  /**
	    \rst
	    A typedef of a class that receives notifications of `pairwise expression
	    <mp::expr::FIRST_PAIRWISE>` arguments.
	    \endrst
	   */
	  typedef ArgHandler PairwiseArgHandler;
	
	  /**
	    \rst
	    A typedef of a class that receives notifications of symbolic
	    (`numeric <mp::expr::FIRST_NUMERIC>` or `string <mp::expr::STRING>`)
	    arguments.
	    \endrst
	   */
	  typedef ArgHandler SymbolicArgHandler;
	
	  /**
	    \rst
	    Receives notification of a `number <mp::expr::FIRST_NUMERIC>` in a
	    nonlinear expression.
	    \endrst
	   */
	  NumericExpr OnNumber(double value) {
	    internal::Unused(value);
	    MP_DISPATCH(OnUnhandled("number"));
	    return NumericExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of a `variable reference <mp::expr::FIRST_REFERENCE>`.
	    \endrst
	   */
	  Reference OnVariableRef(int var_index) {
	    internal::Unused(var_index);
	    MP_DISPATCH(OnUnhandled("variable reference"));
	    return Reference();
	  }
	
	  /**
	    \rst
	    Receives notification of a `common expression <mp::expr::COMMON_EXPR>`
	    (defined variable) reference.
	    \endrst
	   */
	  Reference OnCommonExprRef(int expr_index) {
	    internal::Unused(expr_index);
	    MP_DISPATCH(OnUnhandled("common expression reference"));
	    return Reference();
	  }
	
	  /**
	    \rst
	    Receives notification of a `unary expression <mp::expr::FIRST_UNARY>`.
	    \endrst
	   */
	  NumericExpr OnUnary(expr::Kind kind, NumericExpr arg) {
	    internal::Unused(kind, &arg);
	    MP_DISPATCH(OnUnhandled("unary expression"));
	    return NumericExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of a `binary expression <mp::expr::FIRST_BINARY>`.
	    \endrst
	   */
	  NumericExpr OnBinary(expr::Kind kind, NumericExpr lhs, NumericExpr rhs) {
	    internal::Unused(kind, &lhs, &rhs);
	    MP_DISPATCH(OnUnhandled("binary expression"));
	    return NumericExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of an `if expression <mp::expr::IF>`.
	    \endrst
	   */
	  NumericExpr OnIf(LogicalExpr condition,
	      NumericExpr then_expr, NumericExpr else_expr) {
	    internal::Unused(&condition, &then_expr, &else_expr);
	    MP_DISPATCH(OnUnhandled("if expression"));
	    return NumericExpr();
	  }
	
	  /**
	    \rst
	    A class (struct) that receives notifications of slopes and breakpoints in a
	    `piecewise-linear term <mp::expr::PLTERM>`.
	    \endrst
	   */
	  struct PLTermHandler {
	    /** Receives notification of a slope. */
	    void AddSlope(double slope) { internal::Unused(slope); }
	
	    /** Receives notification of a breakpoint. */
	    void AddBreakpoint(double breakpoint) { internal::Unused(breakpoint); }
	  };
	
	  /**
	    \rst
	    Receives notification of the beginning of a `piecewise-linear term
	    <mp::expr::PLTERM>`.
	    \endrst
	   */
	  PLTermHandler BeginPLTerm(int num_breakpoints) {
	    internal::Unused(num_breakpoints);
	    MP_DISPATCH(OnUnhandled("piecewise-linear term"));
	    return PLTermHandler();
	  }
	
	  /**
	    \rst
	    Receives notification of the end of a `piecewise-linear term
	    <mp::expr::PLTERM>`.
	
	    *arg*: argument that is a variable or a common expression reference.
	    \endrst
	   */
	  NumericExpr EndPLTerm(PLTermHandler handler, Reference arg) {
	    internal::Unused(&handler, &arg);
	    MP_DISPATCH(OnUnhandled("piecewise-linear term"));
	    return NumericExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of the beginning of a `call expression
	    <mp::expr::CALL>`.
	    \endrst
	   */
	  CallArgHandler BeginCall(int func_index, int num_args) {
	    internal::Unused(func_index, num_args);
	    MP_DISPATCH(OnUnhandled("call expression"));
	    return CallArgHandler();
	  }
	
	  /**
	    \rst
	    Receives notification of the end of a `call expression <mp::expr::CALL>`.
	    \endrst
	   */
	  NumericExpr EndCall(CallArgHandler handler) {
	    internal::Unused(handler);
	    MP_DISPATCH(OnUnhandled("call expression"));
	    return NumericExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of the beginning of a `vararg expression
	    <mp::expr::FIRST_VARARG>`.
	    \endrst
	   */
	  VarArgHandler BeginVarArg(expr::Kind kind, int num_args) {
	    internal::Unused(kind, num_args);
	    MP_DISPATCH(OnUnhandled("vararg expression"));
	    return NumericArgHandler();
	  }
	
	  /**
	    \rst
	    Receives notification of the end of a `vararg expression
	    <mp::expr::FIRST_VARARG>`.
	    \endrst
	   */
	  NumericExpr EndVarArg(VarArgHandler handler) {
	    internal::Unused(handler);
	    MP_DISPATCH(OnUnhandled("vararg expression"));
	    return NumericExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of the beginning of a `summation <mp::expr::SUM>`.
	    \endrst
	   */
	  NumericArgHandler BeginSum(int num_args) {
	    internal::Unused(num_args);
	    MP_DISPATCH(OnUnhandled("summation"));
	    return NumericArgHandler();
	  }
	
	  /**
	    \rst
	    Receives notification of the end of a `summation <mp::expr::SUM>`.
	    \endrst
	   */
	  NumericExpr EndSum(NumericArgHandler handler) {
	    internal::Unused(handler);
	    MP_DISPATCH(OnUnhandled("summation"));
	    return NumericExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of the beginning of a `count expression
	    <mp::expr::COUNT>`.
	    \endrst
	   */
	  CountArgHandler BeginCount(int num_args) {
	    internal::Unused(num_args);
	    MP_DISPATCH(OnUnhandled("count expression"));
	    return CountArgHandler();
	  }
	
	  /**
	    \rst
	    Receives notification of the end of a `count expression <mp::expr::COUNT>`.
	    \endrst
	   */
	  CountExpr EndCount(CountArgHandler handler) {
	    internal::Unused(handler);
	    MP_DISPATCH(OnUnhandled("count expression"));
	    return NumericExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of the beginning of a `numberof expression
	    <mp::expr::NUMBEROF>`.
	    \endrst
	   */
	  NumberOfArgHandler BeginNumberOf(int num_args, NumericExpr arg0) {
	    internal::Unused(num_args, &arg0);
	    MP_DISPATCH(OnUnhandled("numberof expression"));
	    return NumberOfArgHandler();
	  }
	
	  /**
	    \rst
	    Receives notification of the end of a `numberof expression
	    <mp::expr::NUMBEROF>`.
	    \endrst
	   */
	  NumericExpr EndNumberOf(NumberOfArgHandler handler) {
	    internal::Unused(handler);
	    MP_DISPATCH(OnUnhandled("numberof expression"));
	    return NumericExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of the beginning of a `symbolic numberof
	    expression <mp::expr::NUMBEROF_SYM>`.
	    \endrst
	   */
	  SymbolicArgHandler BeginSymbolicNumberOf(int num_args, Expr arg0) {
	    internal::Unused(num_args, &arg0);
	    MP_DISPATCH(OnUnhandled("symbolic numberof expression"));
	    return SymbolicArgHandler();
	  }
	
	  /**
	    \rst
	    Receives notification of the end of a `symbolic numberof expression
	    <mp::expr::NUMBEROF_SYM>`.
	    \endrst
	   */
	  NumericExpr EndSymbolicNumberOf(SymbolicArgHandler handler) {
	    internal::Unused(handler);
	    MP_DISPATCH(OnUnhandled("symbolic numberof expression"));
	    return NumericExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of a `Boolean value <mp::expr::BOOL>`.
	    \endrst
	   */
	  LogicalExpr OnBool(bool value) {
	    internal::Unused(value);
	    MP_DISPATCH(OnUnhandled("bool"));
	    return LogicalExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of a `logical not <mp::expr::NOT>`.
	    \endrst
	   */
	  LogicalExpr OnNot(LogicalExpr arg) {
	    internal::Unused(&arg);
	    MP_DISPATCH(OnUnhandled("logical not"));
	    return LogicalExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of a `binary logical expression
	    <mp::expr::FIRST_BINARY_LOGICAL>`.
	    \endrst
	   */
	  LogicalExpr OnBinaryLogical(
	      expr::Kind kind, LogicalExpr lhs, LogicalExpr rhs) {
	    internal::Unused(kind, &lhs, &rhs);
	    MP_DISPATCH(OnUnhandled("binary logical expression"));
	    return LogicalExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of a `relational expression
	    <mp::expr::FIRST_RELATIONAL>`.
	    \endrst
	   */
	  LogicalExpr OnRelational(expr::Kind kind, NumericExpr lhs, NumericExpr rhs) {
	    internal::Unused(kind, &lhs, &rhs);
	    MP_DISPATCH(OnUnhandled("relational expression"));
	    return LogicalExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of a `logical count expression
	    <mp::expr::FIRST_LOGICAL_COUNT>`.
	    \endrst
	   */
	  LogicalExpr OnLogicalCount(expr::Kind kind, NumericExpr lhs, CountExpr rhs) {
	    internal::Unused(kind, &lhs, &rhs);
	    MP_DISPATCH(OnUnhandled("logical count expression"));
	    return LogicalExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of an `implication expression
	    <mp::expr::IMPLICATION>`.
	    \endrst
	   */
	  LogicalExpr OnImplication(
	      LogicalExpr condition, LogicalExpr then_expr, LogicalExpr else_expr) {
	    internal::Unused(&condition, &then_expr, &else_expr);
	    MP_DISPATCH(OnUnhandled("implication expression"));
	    return LogicalExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of the beginning of an `iterated logical
	    expression <mp::expr::FIRST_ITERATED_LOGICAL>`.
	    \endrst
	   */
	  LogicalArgHandler BeginIteratedLogical(expr::Kind kind, int num_args) {
	    internal::Unused(kind, num_args);
	    MP_DISPATCH(OnUnhandled("iterated logical expression"));
	    return LogicalArgHandler();
	  }
	
	  /**
	    \rst
	    Receives notification of the end of an `iterated logical expression
	    <mp::expr::FIRST_ITERATED_LOGICAL>`.
	    \endrst
	   */
	  LogicalExpr EndIteratedLogical(LogicalArgHandler handler) {
	    internal::Unused(handler);
	    MP_DISPATCH(OnUnhandled("iterated logical expression"));
	    return LogicalExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of the beginning of a `pairwise expression
	    <mp::expr::FIRST_PAIRWISE>`.
	    \endrst
	   */
	  PairwiseArgHandler BeginPairwise(expr::Kind kind, int num_args) {
	    internal::Unused(kind, num_args);
	    MP_DISPATCH(OnUnhandled("pairwise expression"));
	    return PairwiseArgHandler();
	  }
	
	  /**
	    \rst
	    Receives notification of the end of a `pairwise expression
	    <mp::expr::FIRST_PAIRWISE>`.
	    \endrst
	   */
	  LogicalExpr EndPairwise(PairwiseArgHandler handler) {
	    internal::Unused(handler);
	    MP_DISPATCH(OnUnhandled("pairwise expression"));
	    return LogicalExpr();
	  }
	
	  /**
	    \rst
	    Receives notification of a `string <mp::expr::STRING>`.
	    The *value* argument is a string value and it is not null-terminated.
	    \endrst
	   */
	  Expr OnString(fmt::StringRef value) {
	    internal::Unused(&value);
	    MP_DISPATCH(OnUnhandled("string"));
	    return Expr();
	  }
	
	  /**
	    \rst
	    Receives notification of a `symbolic if expression <mp::expr::IFSYM>`.
	    \endrst
	   */
	  Expr OnSymbolicIf(LogicalExpr condition, Expr then_expr, Expr else_expr) {
	    internal::Unused(&condition, &then_expr, &else_expr);
	    MP_DISPATCH(OnUnhandled("symbolic if expression"));
	    return Expr();
	  }
	
	  /** Receives notification of the end of the input. */
	  void EndInput() {}
	};
	
	/**
	  \rst
	  An NL handler that ignores all input.
	
	  `~mp::NullNLHandler` can be used as a base class when only a subset of
	  constructs needs to be handled. Unhandled constructs will be ignored, not
	  reported.
	
	  *ExprType* is a return type of expression handler methods such as
	  `~mp::NLHandler::OnUnary` useful for building expression objects. If not used
	  it can be any default-constructible type.
	  \endrst
	 */
	template <typename ExprType>
	class NullNLHandler : public NLHandler<NullNLHandler<ExprType>, ExprType> {
	 public:
	  /** Receives notification of an unhandled construct and ignores it. */
	  void OnUnhandled(const char *) {}
	};
	
	namespace internal {
	
	class ReaderBase {
	 protected:
	  const char *ptr_, *start_, *end_;
	  const char *token_;  // start of the current token
	  std::string name_;
	
	  ~ReaderBase() {}
	
	 public:
	  ReaderBase(NLStringRef data, fmt::CStringRef name);
	
	  char ReadChar() {
	    token_ = ptr_;
	    return *ptr_++;
	  }
	

	  const char *ptr() const { return ptr_; }
	  void set_ptr(const char *ptr) { token_ = ptr_ = ptr; }
	
	  bool IsEOF(const char *ptr) const { return ptr == end_ + 1; }
	  bool IsEOF() const { return IsEOF(ptr_); }
	};
	
	#ifdef FMT_LOCALE
	typedef fmt::Locale Locale;
	#else
	# warning "Parsing is locale-dependent"
	struct Locale {
	  double strtod(const char *&str) const {
	    char *end = 0;
	    double result = std::strtod(str, &end);
	    str = end;
	    return result;
	  }
	};
	#endif
	
	template <typename Locale = Locale>
	class TextReader : public ReaderBase {
	 private:
	  const char *line_start_;
	  int line_;
	
	  class CopyableLocale : public Locale {
	   public:
	    CopyableLocale() {}
	    // All Locale objects refer to the "C" locale, so no need to copy.
	    CopyableLocale(const CopyableLocale &) : Locale() {}
	    CopyableLocale &operator=(const CopyableLocale &) { return *this; }
	  };
	  CopyableLocale locale_;
	
	  // Reads an integer without a sign.
	  // Int: signed or unsigned integer type.
	  template <typename Int>
	  bool ReadIntWithoutSign(Int& value) {
	    char c = *ptr_;
	    if (c < '0' || c > '9')
	      return false;
	    typedef typename MakeUnsigned<Int>::Type UInt;
	    UInt result = 0;
	    do {
	      UInt new_result = result * 10 + (c - '0');
	      if (new_result < result)
	        ReportError("number is too big");
	      result = new_result;
	      c = *++ptr_;
	    } while (c >= '0' && c <= '9');
	    UInt max = std::numeric_limits<Int>::max();
	    if (result > max)
	      ReportError("number is too big");
	    value = result;
	    return true;
	  }
	
	  template <typename Int>
	  bool DoReadOptionalInt(Int &value) {
	    SkipSpace();
	    char sign = *ptr_;
	    if (sign == '+' || sign == '-')
	      ++ptr_;
	    typedef typename MakeUnsigned<Int>::Type UInt;
	    UInt result = 0;
	    if (!ReadIntWithoutSign<UInt>(result))
	      return false;
	    UInt max = std::numeric_limits<Int>::max();
	    if (result > max && !(sign == '-' && result == max + 1))
	      ReportError("number is too big");
	    value = sign != '-' ? result : 0 - result;
	    return true;
	  }
	
	  // Reads a nonnegative integer and checks that adding it to accumulator
	  // doesn't overflow.
	  int ReadUInt(int &accumulator) {
	    int value = ReadUInt();
	    if (accumulator > std::numeric_limits<int>::max() - value)
	      ReportError("integer overflow");
	    accumulator += value;
	    return value;
	  }
	
	  template <typename Int>
	  Int ReadUInt() {
	    SkipSpace();
	    Int value = 0;
	    if (!ReadIntWithoutSign(value))
	      ReportError("expected unsigned integer");
	    return value;
	  }
	
	  bool ReadOptionalInt(int &value) { return DoReadOptionalInt(value); }
	
	  bool ReadOptionalUInt(int &value) {
	    SkipSpace();
	    return ReadIntWithoutSign(value);
	  }
	
	  bool ReadOptionalDouble(double &value);
	
	  void DoReportError(
	      const char *loc, fmt::CStringRef format_str,
	      const fmt::ArgList &args = fmt::ArgList());
	
	  void SkipSpace() {
	    while (std::isspace(*ptr_) && *ptr_ != '\n')
	      ++ptr_;
	    token_ = ptr_;
	  }
	
	 public:
	  TextReader(NLStringRef data, fmt::CStringRef name);
	
	  Locale &locale() { return locale_; }
	
	  void ReportError(fmt::CStringRef format_str, const fmt::ArgList &args) {
	    DoReportError(token_, format_str, args);
	  }
	  FMT_VARIADIC(void, ReportError, fmt::CStringRef)
	
	  void ReadTillEndOfLine() {
	    while (char c = *ptr_) {
	      ++ptr_;
	      if (c == '\n') {
	        line_start_ = ptr_;
	        ++line_;
	        return;
	      }
	    }
	    DoReportError(ptr_, "expected newline");
	  }
	
	  template <typename Int>
	  Int ReadInt() {
	    Int value = 0;
	    if (!DoReadOptionalInt(value))
	      ReportError("expected integer");
	    return value;
	  }
	
	  int ReadUInt() { return ReadUInt<int>(); }
	
	  double ReadDouble() {
	    SkipSpace();
	    const char *start = ptr_;
	    double value = 0;
	    if (*ptr_ != '\n')
	      value = locale_.strtod(ptr_);
	    if (ptr_ == start)
	      ReportError("expected double");
	    return value;
	  }
	
	  fmt::StringRef ReadString();
	
	  // Reads a function or suffix name.
	  fmt::StringRef ReadName();
	
	  // Reads an .nl file header. The header is always in text format, so this
	  // function doesn't have a counterpart in BinaryReader.
	  void ReadHeader(NLHeader &header);
	};
	
	// Converter that doesn't change the input.
	class IdentityConverter {
	 public:
	  template <typename T>
	  T Convert(T value) { return value; }
	};
	
	// Converter that changes the input endianness.
	class EndiannessConverter {
	 private:
	  void Convert(char *data, std::size_t size) {
	    std::reverse(data, data + size);
	  }
	
	 public:
	  template <typename T>
	  T Convert(T value) {
	    Convert(reinterpret_cast<char*>(&value), sizeof(T));
	    return value;
	  }
	};
	
	class BinaryReaderBase : public ReaderBase {
	 protected:
	  explicit BinaryReaderBase(const ReaderBase &base) : ReaderBase(base) {}
	
	  // Reads length chars.
	  const char *Read(int length) {
	    if (end_ - ptr_ < length) {
	      token_ = end_;
	      ReportError("unexpected end of file");
	    }
	    const char *start = ptr_;
	    ptr_ += length;
	    return start;
	  }
	
	 public:
	  void ReportError(fmt::CStringRef format_str, const fmt::ArgList &args);
	  FMT_VARIADIC(void, ReportError, fmt::CStringRef)
	
	  void ReadTillEndOfLine() {
	    // Do nothing.
	  }
	};
	
	// Binary reader.
	template <typename InputConverter = IdentityConverter>
	class BinaryReader : private InputConverter, public BinaryReaderBase {
	 public:
	  explicit BinaryReader(const ReaderBase &base) : BinaryReaderBase(base) {}
	
	  template <typename Int>
	  Int ReadInt() {
	    token_ = ptr_;
	    Int val;
	    memcpy(&val, Read(sizeof(Int)), sizeof(Int));
	    return this->Convert(val);
	  }
	
	  int ReadUInt() {
	    int value = ReadInt<int>();
	    if (value < 0)
	      ReportError("expected unsigned integer");
	    return value;
	  }
	
	  double ReadDouble() {
	    token_ = ptr_;
	    double val;
	    memcpy(&val, Read(sizeof(double)), sizeof(double));
	    return this->Convert(val);
	  }
	
	  fmt::StringRef ReadString() {
	    int length = ReadUInt();
	    return fmt::StringRef(length != 0 ? Read(length) : 0, length);
	  }
	
	  // Reads a function or suffix name.
	  fmt::StringRef ReadName() { return ReadString(); }
	};
	
	// An NLHandler that forwards notification of variable bounds to another
	// handler and ignores all other notifications.
	template <typename Handler>
	class VarBoundHandler : public NullNLHandler<typename Handler::Expr> {
	 private:
	  Handler &handler_;
	
	 public:
	  explicit VarBoundHandler(Handler &h) : handler_(h) {}
	
	  void OnVarBounds(int index, double lb, double ub) {
	    handler_.OnVarBounds(index, lb, ub);
	  }
	};
	
	// Linear expression handler that ignores input.
	struct NullLinearExprHandler {
	  void AddTerm(int, double) {}
	};
	
	// An NL reader.
	// Handler: a class implementing the NLHandler concept that receives
	//          notifications of NL constructs
	template <typename Reader, typename Handler>
	class NLReader {
	 private:
	  Reader &reader_;
	  const NLHeader &header_;
	  Handler &handler_;
	  int flags_;
	  int num_vars_and_exprs_;  // Number of variables and common expressions.
	
	  typedef typename Handler::Expr Expr;
	  typedef typename Handler::NumericExpr NumericExpr;
	  typedef typename Handler::LogicalExpr LogicalExpr;
	  typedef typename Handler::Reference Reference;
	
	  double ReadConstant(char code);
	  double ReadConstant() { return ReadConstant(reader_.ReadChar()); }
	
	  // Reads a nonnegative integer and checks that it is less than ub.
	  // ub is unsigned so that it can hold value INT_MAX + 1u.
	  int ReadUInt(unsigned ub) {
	    int value = reader_.ReadUInt();
	    unsigned unsigned_value = value;
	    if (unsigned_value >= ub)
	      reader_.ReportError("integer {} out of bounds", value);
	    return value;
	  }
	
	  // Reads a nonnegative integer and checks that it is in the range [lb, ub).
	  int ReadUInt(unsigned lb, unsigned ub) {
	    int value = reader_.ReadUInt();
	    unsigned unsigned_value = value;
	    if (unsigned_value < lb || unsigned_value >= ub)
	      reader_.ReportError("integer {} out of bounds", value);
	    return value;
	  }
	
	  // Minimum number of arguments for an iterated expression that has a
	  // binary counterpart. Examples: sum (+), forall (&&), exists (||).
	  enum {MIN_ITER_ARGS = 3};
	
	  int ReadNumArgs(int min_args = MIN_ITER_ARGS) {
	    int num_args = reader_.ReadUInt();
	    if (num_args < min_args)
	      reader_.ReportError("too few arguments");
	    return num_args;
	  }
	
	  // Reads a variable or a common expression reference.
	  Reference DoReadReference() {
	    int index = ReadUInt(num_vars_and_exprs_);
	    reader_.ReadTillEndOfLine();
	    return index < header_.num_vars ?
	          handler_.OnVariableRef(index) :
	          handler_.OnCommonExprRef(index - header_.num_vars);
	  }
	
	  // Reads a variable or a common expression reference.
	  Reference ReadReference() {
	    if (reader_.ReadChar() != 'v')
	      reader_.ReportError("expected reference");
	    return DoReadReference();
	  }
	
	  template <typename ExprReader, typename ArgHandler>
	  void DoReadArgs(int num_args, ArgHandler &arg_handler) {
	    ExprReader expr_reader;
	    for (int i = 0; i < num_args; ++i)
	      arg_handler.AddArg(expr_reader.Read(*this));
	  }
	
	  template <typename ExprReader, typename ArgHandler>
	  void ReadArgs(int num_args, ArgHandler &arg_handler) {
	    reader_.ReadTillEndOfLine();
	    DoReadArgs<ExprReader>(num_args, arg_handler);
	  }
	
	  int ReadOpCode() {
	    int opcode = reader_.ReadUInt();
	    if (opcode > internal::MAX_OPCODE)
	      reader_.ReportError("invalid opcode {}", opcode);
	    reader_.ReadTillEndOfLine();
	    return opcode;
	  }
	
	  typename Handler::CountExpr ReadCountExpr() {
	    int num_args = ReadNumArgs(1);
	    typename Handler::CountArgHandler args = handler_.BeginCount(num_args);
	    ReadArgs<LogicalExprReader>(num_args, args);
	    return handler_.EndCount(args);
	  }
	
	  // Helper structs to provide a uniform interface to Read{Numeric,Logical}Expr
	  // since it is not possible to overload on expression type as NumericExpr
	  // and LogicalExpr can be the same type.
	  struct NumericExprReader {
	    typedef NumericExpr Expr;
	    Expr Read(NLReader &r) const { return r.ReadNumericExpr(); }
	  };
	  struct LogicalExprReader {
	    typedef LogicalExpr Expr;
	    Expr Read(NLReader &r) const { return r.ReadLogicalExpr(); }
	  };
	  struct SymbolicExprReader {
	    typedef typename Handler::Expr Expr;
	    Expr Read(NLReader &r) const { return r.ReadSymbolicExpr(); }
	  };
	
	  // A helper struct used to make sure that the arguments to a binary
	  // expression are read in the correct order and avoid errors of the form:
	  //   MakeBinary(opcode, ReadNumericExpr(), ReadNumericExpr())
	  // The above code is incorrect as the order of evaluation of arguments is
	  // unspecified.
	  template <typename ExprReader = NumericExprReader>
	  struct BinaryArgReader {
	    typename ExprReader::Expr lhs;
	    typename ExprReader::Expr rhs;
	    BinaryArgReader(NLReader &r)
	      : lhs(ExprReader().Read(r)), rhs(ExprReader().Read(r)) {}
	  };
	
	  // Reads a numeric or string expression.
	  Expr ReadSymbolicExpr();
	
	  // Reads a numeric expression.
	  // ignore_zero: if true, zero constants are ignored
	  NumericExpr ReadNumericExpr(bool ignore_zero = false) {
	    return ReadNumericExpr(reader_.ReadChar(), ignore_zero);
	  }
	  NumericExpr ReadNumericExpr(char code, bool ignore_zero);
	  NumericExpr ReadNumericExpr(int opcode);
	
	  // Reads a logical expression.
	  LogicalExpr ReadLogicalExpr();
	  LogicalExpr ReadLogicalExpr(int opcode);
	
	  enum ItemType { VAR, OBJ, CON, PROB };
	
	  template <ItemType T>
	  class ItemHandler {
	   protected:
	    NLReader &reader_;
	
	   public:
	    static const ItemType TYPE = T;
	    explicit ItemHandler(NLReader &r) : reader_(r) {}
	  };
	
	  struct VarHandler : ItemHandler<VAR> {
	    explicit VarHandler(NLReader &r) : ItemHandler<VAR>(r) {}
	
	    int num_items() const { return this->reader_.header_.num_vars; }
	
	    void SetBounds(int index, double lb, double ub) {
	      this->reader_.handler_.OnVarBounds(index, lb, ub);
	    }
	    void SetInitialValue(int index, double value) {
	      this->reader_.handler_.OnInitialValue(index, value);
	    }
	  };
	
	  struct ObjHandler : ItemHandler<OBJ> {
	    explicit ObjHandler(NLReader &r) : ItemHandler<OBJ>(r) {}
	
	    int num_items() const { return this->reader_.header_.num_objs; }
	
	    // Returns true if objective expression should be skipped.
	    bool SkipExpr(int obj_index) const {
	      return !this->reader_.handler_.NeedObj(obj_index);
	    }
	
	    typename Handler::LinearObjHandler OnLinearExpr(int index, int num_terms) {
	      auto& h = this->reader_.handler_;
	      return h.OnLinearObjExpr( h.resulting_obj_index(index), num_terms);
	    }
	  };
	
	  struct ConHandler : ItemHandler<CON> {
	    explicit ConHandler(NLReader &r) : ItemHandler<CON>(r) {}
	
	    int num_items() const {
	      return this->reader_.header_.num_algebraic_cons +
	          this->reader_.header_.num_logical_cons;
	    }
	  };
	
	  struct ProblemHandler : ItemHandler<PROB> {
	    explicit ProblemHandler(NLReader &r) : ItemHandler<PROB>(r) {}
	
	    // An NL input always contains one problem.
	    int num_items() const { return 1; }
	  };
	
	  // Reads the linear part of an objective or constraint expression.
	  template <typename LinearHandler>
	  void ReadLinearExpr();
	
	  template <typename LinearHandler>
	  void ReadLinearExpr(int num_terms, LinearHandler linear_expr);
	
	  // Reads column sizes, numbers of nonzeros in the first num_var − 1
	  // columns of the Jacobian sparsity matrix.
	  template <bool CUMULATIVE>
	  void ReadColumnSizes();
	
	  // Reads initial values for primal or dual variables.
	  template <typename ValueHandler>
	  void ReadInitialValues();
	
	  struct IntReader {
	    int operator()(Reader &r) const { return r.template ReadInt<int>(); }
	  };
	
	  struct DoubleReader {
	    double operator()(Reader &r) const { return r.ReadDouble(); }
	  };
	
	  template <typename ValueReader, typename SuffixHandler>
	  void ReadSuffixValues(int num_values, int num_items, SuffixHandler &handler) {
	    ValueReader read;
	    for (int i = 0; i < num_values; ++i) {
	      int index = ReadUInt(num_items);
	      handler.SetValue(index, read(reader_));
	      reader_.ReadTillEndOfLine();
	    }
	  }
	
	  template <typename ItemInfo>
	  void ReadSuffix(int info);
	
	 public:
	  NLReader(Reader &reader, const NLHeader &header, Handler &handler, int flags)
	    : reader_(reader), header_(header), handler_(handler), flags_(flags),
	      num_vars_and_exprs_(0) {}
	
	  // Algebraic constraint handler.
	  struct AlgebraicConHandler : ItemHandler<CON> {
	    explicit AlgebraicConHandler(NLReader &r) : ItemHandler<CON>(r) {}
	
	    int num_items() const { return this->reader_.header_.num_algebraic_cons; }
	
	    // Returns false because constraint expressions are always read.
	    bool SkipExpr(int) const { return false; }
	
	    typename Handler::LinearConHandler OnLinearExpr(int index, int num_terms) {
	      return this->reader_.handler_.OnLinearConExpr(index, num_terms);
	    }
	
	    void SetBounds(int index, double lb, double ub) {
	      this->reader_.handler_.OnConBounds(index, lb, ub);
	    }
	    void SetInitialValue(int index, double value) {
	      this->reader_.handler_.OnInitialDualValue(index, value);
	    }
	  };
	
	  // Reads variable or constraint bounds.
	  template <typename BoundHandler>
	  void ReadBounds();
	
	  // bound_reader: a reader after variable bounds section input
	  void Read(Reader *bound_reader);
	
	  void Read();
	};
	
	template <typename Reader, typename Handler>
	double NLReader<Reader, Handler>::ReadConstant(char code) {
	  double value = 0;
	  switch (code) {
	  case 'n':
	    value = reader_.ReadDouble();
	    break;
	  case 's':
	    value = reader_.template ReadInt<short>();
	    break;
	  case 'l':
	    // The following check is necessary for compatibility with ASL.
	    if (sizeof(double) == 2 * sizeof(int))
	      value = reader_.template ReadInt<int>();
	    else
	      value = reader_.template ReadInt<long>();
	    break;
	  default:
	    reader_.ReportError("expected constant");
	  }
	  reader_.ReadTillEndOfLine();
	  return value;
	}
	
	template <typename Reader, typename Handler>
	typename Handler::Expr NLReader<Reader, Handler>::ReadSymbolicExpr() {
	  char c = reader_.ReadChar();
	  switch (c) {
	  case 'h':
	    return handler_.OnString(reader_.ReadString());
	  case 'o': {
	    int opcode = ReadOpCode();
	    if (opcode != expr::nl_opcode(expr::IFSYM))
	      return ReadNumericExpr(opcode);
	    // Read symbolic if expression.
	    LogicalExpr condition = ReadLogicalExpr();
	    Expr then_expr = ReadSymbolicExpr();
	    Expr else_expr = ReadSymbolicExpr();
	    return handler_.OnSymbolicIf(condition, then_expr, else_expr);
	  }
	  }
	  return ReadNumericExpr(c, false);
	}
	
	template <typename Reader, typename Handler>
	typename Handler::NumericExpr
	    NLReader<Reader, Handler>::ReadNumericExpr(char code, bool ignore_zero) {
	  switch (code) {
	  case 'f': {
	    // Read a function call.
	    int func_index = ReadUInt(header_.num_funcs);
	    int num_args = reader_.ReadUInt();
	    reader_.ReadTillEndOfLine();
	    typename Handler::CallArgHandler args =
	        handler_.BeginCall(func_index, num_args);
	    for (int i = 0; i < num_args; ++i)
	      args.AddArg(ReadSymbolicExpr());
	    return handler_.EndCall(args);
	  }
	  case 'n': case 'l': case 's': {
	    // Read a number.
	    double value = ReadConstant(code);
	    if (ignore_zero && value == 0)
	      break;  // Ignore zero constant.
	    return handler_.OnNumber(value);
	  }
	  case 'o':
	    return ReadNumericExpr(ReadOpCode());
	  case 'v':
	    return DoReadReference();
	  default:
	    reader_.ReportError("expected expression");
	  }
	  return NumericExpr();
	}
	
	template <typename Reader, typename Handler>
	typename Handler::NumericExpr
	    NLReader<Reader, Handler>::ReadNumericExpr(int opcode) {
	  const internal::OpCodeInfo &info = internal::GetOpCodeInfo(opcode);
	  expr::Kind kind = info.kind;
	  switch (info.first_kind) {
	  case expr::FIRST_UNARY:
	    return handler_.OnUnary(kind, ReadNumericExpr());
	  case expr::FIRST_BINARY: {
	    BinaryArgReader<> args(*this);
	    return handler_.OnBinary(kind, args.lhs, args.rhs);
	  }
	  case expr::IF: {
	    LogicalExpr condition = ReadLogicalExpr();
	    NumericExpr then_expr = ReadNumericExpr();
	    NumericExpr else_expr = ReadNumericExpr();
	    return handler_.OnIf(condition, then_expr, else_expr);
	  }
	  case expr::PLTERM: {
	    // Read a piecewise-linear term.
	    int num_slopes = reader_.ReadUInt();
	    if (num_slopes <= 1)
	      reader_.ReportError("too few slopes in piecewise-linear term");
	    reader_.ReadTillEndOfLine();
	    typename Handler::PLTermHandler pl_handler =
	        handler_.BeginPLTerm(num_slopes - 1);
	    for (int i = 0; i < num_slopes - 1; ++i) {
	      pl_handler.AddSlope(ReadConstant());
	      pl_handler.AddBreakpoint(ReadConstant());
	    }
	    pl_handler.AddSlope(ReadConstant());
	    return handler_.EndPLTerm(pl_handler, ReadReference());
	  }
	  case expr::FIRST_VARARG: {
	    // Read a vararg expression (min or max).
	    int num_args = ReadNumArgs(1);
	    typename Handler::VarArgHandler args = handler_.BeginVarArg(kind, num_args);
	    ReadArgs<NumericExprReader>(num_args, args);
	    return handler_.EndVarArg(args);
	  }
	  case expr::SUM: {
	    int num_args = ReadNumArgs();
	    typename Handler::NumericArgHandler args = handler_.BeginSum(num_args);
	    ReadArgs<NumericExprReader>(num_args, args);
	    return handler_.EndSum(args);
	  }
	  case expr::COUNT:
	    return ReadCountExpr();
	  case expr::NUMBEROF: {
	    // Read a numberof expression.
	    int num_args = ReadNumArgs(1);
	    reader_.ReadTillEndOfLine();
	    typename Handler::NumberOfArgHandler args =
	        handler_.BeginNumberOf(num_args, ReadNumericExpr());
	    DoReadArgs<NumericExprReader>(num_args - 1, args);
	    return handler_.EndNumberOf(args);
	  }
	  case expr::NUMBEROF_SYM: {
	    // Read a symbolic numberof expression.
	    int num_args = ReadNumArgs(1);
	    reader_.ReadTillEndOfLine();
	    typename Handler::SymbolicArgHandler args =
	        handler_.BeginSymbolicNumberOf(num_args, ReadSymbolicExpr());
	    DoReadArgs<SymbolicExprReader>(num_args - 1, args);
	    return handler_.EndSymbolicNumberOf(args);
	  }
	  default:
	    reader_.ReportError("expected numeric expression opcode");
	  }
	  return NumericExpr();
	}
	
	template <typename Reader, typename Handler>
	typename Handler::LogicalExpr NLReader<Reader, Handler>::ReadLogicalExpr() {
	  switch (char c = reader_.ReadChar()) {
	  case 'n': case 'l': case 's':
	    return handler_.OnBool(ReadConstant(c) != 0);
	  case 'o':
	    return ReadLogicalExpr(ReadOpCode());
	  }
	  reader_.ReportError("expected logical expression");
	  return LogicalExpr();
	}
	
	template <typename Reader, typename Handler>
	typename Handler::LogicalExpr
	    NLReader<Reader, Handler>::ReadLogicalExpr(int opcode) {
	  const internal::OpCodeInfo &info = internal::GetOpCodeInfo(opcode);
	  expr::Kind kind = info.kind;
	  switch (info.first_kind) {
	  case expr::NOT:
	    return handler_.OnNot(ReadLogicalExpr());
	  case expr::FIRST_BINARY_LOGICAL: {
	    BinaryArgReader<LogicalExprReader> args(*this);
	    return handler_.OnBinaryLogical(kind, args.lhs, args.rhs);
	  }
	  case expr::FIRST_RELATIONAL: {
	    BinaryArgReader<> args(*this);
	    return handler_.OnRelational(kind, args.lhs, args.rhs);
	  }
	  case expr::FIRST_LOGICAL_COUNT: {
	    NumericExpr lhs = ReadNumericExpr();
	    char c = reader_.ReadChar();
	    if (c != 'o' || internal::GetOpCodeInfo(ReadOpCode()).kind != expr::COUNT)
	      reader_.ReportError("expected count expression");
	    return handler_.OnLogicalCount(kind, lhs, ReadCountExpr());
	  }
	  case expr::IMPLICATION: {
	    // Read an implication (=>).
	    LogicalExpr condition = ReadLogicalExpr();
	    LogicalExpr then_expr = ReadLogicalExpr();
	    LogicalExpr else_expr = ReadLogicalExpr();
	    return handler_.OnImplication(condition, then_expr, else_expr);
	  }
	  case expr::FIRST_ITERATED_LOGICAL: {
	    // Read an iterated logical expression (exists or forall).
	    int num_args = ReadNumArgs();
	    typename Handler::LogicalArgHandler args =
	        handler_.BeginIteratedLogical(kind, num_args);
	    ReadArgs<LogicalExprReader>(num_args, args);
	    return handler_.EndIteratedLogical(args);
	  }
	  case expr::FIRST_PAIRWISE: {
	    // Read a pairwise expression (alldiff or !alldiff).
	    int num_args = ReadNumArgs(1);
	    typename Handler::PairwiseArgHandler args =
	        handler_.BeginPairwise(kind, num_args);
	    ReadArgs<NumericExprReader>(num_args, args);
	    return handler_.EndPairwise(args);
	  }
	  default:
	    reader_.ReportError("expected logical expression opcode");
	  }
	  return LogicalExpr();
	}
	
	template <typename Reader, typename Handler>
	template <typename LinearHandler>
	void NLReader<Reader, Handler>::ReadLinearExpr() {
	  LinearHandler lh(*this);
	  int index = ReadUInt(lh.num_items());
	  // The number of terms should be less than num_vars because common
	  // expressions are not allowed in a linear expression.
	  int num_terms = ReadUInt(1, header_.num_vars + 1u);
	  reader_.ReadTillEndOfLine();
	  if (lh.SkipExpr(index))
	    ReadLinearExpr(num_terms, NullLinearExprHandler());
	  else
	    ReadLinearExpr(num_terms, lh.OnLinearExpr(index, num_terms));
	}
	
	template <typename Reader, typename Handler>
	template <typename LinearHandler>
	void NLReader<Reader, Handler>::ReadLinearExpr(
	    int num_terms, LinearHandler linear_expr) {
	  for (int i = 0; i < num_terms; ++i) {
	    // Variable index should be less than num_vars because common
	    // expressions are not allowed in a linear expression.
	    int var_index = ReadUInt(header_.num_vars);
	    double coef = reader_.ReadDouble();
	    reader_.ReadTillEndOfLine();
	    linear_expr.AddTerm(var_index, coef);
	  }
	}
	
	template <typename Reader, typename Handler>
	template <typename BoundHandler>
	void NLReader<Reader, Handler>::ReadBounds() {
	  enum BoundType {
	    RANGE,     // Both lower and upper bounds: l <= body <= u.
	    UPPER,     // Only upper bound: body <= u.
	    LOWER,     // Only lower bound: l <= body.
	    FREE,      // No constraints on body (free variable or constraint).
	    CONSTANT,  // Equal to constant: body = c.
	    COMPL      // Body complements variable v[i - 1].
	  };
	  reader_.ReadTillEndOfLine();
	  double lb = 0, ub = 0;
	  BoundHandler bh(*this);
	  int num_bounds = bh.num_items();
	  double infinity = std::numeric_limits<double>::infinity();
	  for (int i = 0; i < num_bounds; ++i) {
	    switch (reader_.ReadChar() - '0') {
	    case RANGE:
	      lb = reader_.ReadDouble();
	      ub = reader_.ReadDouble();
	      break;
	    case UPPER:
	      lb = -infinity;
	      ub = reader_.ReadDouble();
	      break;
	    case LOWER:
	      lb = reader_.ReadDouble();
	      ub = infinity;
	      break;
	    case FREE:
	      lb = -infinity;
	      ub =  infinity;
	      break;
	    case CONSTANT:
	      lb = ub = reader_.ReadDouble();
	      break;
	    case COMPL:
	      if (BoundHandler::TYPE == CON) {
	        int flags = reader_.template ReadInt<int>();
	        int var_index = reader_.ReadUInt();
	        // Don't use NLReader::ReadUInt(int, int) as num_vars + 1 may overflow.
	        if (var_index == 0 || var_index > header_.num_vars)
	          reader_.ReportError("integer {} out of bounds", var_index);
	        --var_index;
	        int mask = ComplInfo::INF_LB | ComplInfo::INF_UB;
	        handler_.OnComplementarity(i, var_index, ComplInfo(flags & mask));
	        reader_.ReadTillEndOfLine();
	        continue;
	      }
	      reader_.ReportError("COMPL bound type is invalid for variables");
	      break;
	    default:
	      reader_.ReportError("expected bound");
	    }
	    reader_.ReadTillEndOfLine();
	    bh.SetBounds(i, lb, ub);
	  }
	}
	
	template <typename Reader, typename Handler>
	template <bool CUMULATIVE>
	void NLReader<Reader, Handler>::ReadColumnSizes() {
	  int num_sizes = header_.num_vars - 1;
	  if (reader_.ReadUInt() != num_sizes)
	    reader_.ReportError("expected {}", num_sizes);
	  reader_.ReadTillEndOfLine();
	  typename Handler::ColumnSizeHandler size_handler = handler_.OnColumnSizes();
	  int prev_size = 0;
	  for (int i = 0; i < num_sizes; ++i) {
	    int size = reader_.ReadUInt();
	    if (CUMULATIVE) {
	      if (size < prev_size)
	        reader_.ReportError("invalid column offset");
	      size -= prev_size;
	      prev_size += size;
	    }
	    size_handler.Add(size);
	    reader_.ReadTillEndOfLine();
	  }
	}
	
	template <typename Reader, typename Handler>
	template <typename ValueHandler>
	void NLReader<Reader, Handler>::ReadInitialValues() {
	  int num_values = reader_.ReadUInt();
	  ValueHandler vh(*this);
	  if (num_values > vh.num_items())
	    reader_.ReportError("too many initial values");
	  reader_.ReadTillEndOfLine();
	  for (int i = 0; i < num_values; ++i) {
	    int index = ReadUInt(vh.num_items());
	    vh.SetInitialValue(index, reader_.ReadDouble());
	    reader_.ReadTillEndOfLine();
	  }
	}
	
	template <typename Reader, typename Handler>
	template <typename ItemInfo>
	void NLReader<Reader, Handler>::ReadSuffix(int info) {
	  int num_items = ItemInfo(*this).num_items();
	  int num_values = ReadUInt(1, num_items + 1);
	  fmt::StringRef name = reader_.ReadName();
	  reader_.ReadTillEndOfLine();
	  suf::Kind kind = static_cast<suf::Kind>(info & internal::SUFFIX_KIND_MASK);
	  if ((info & suf::FLOAT) != 0) {
	    typename Handler::DblSuffixHandler
	        suffix_handler = handler_.OnDblSuffix(name, kind, num_values);
	    ReadSuffixValues<DoubleReader>(num_values, num_items, suffix_handler);
	  } else {
	    typename Handler::IntSuffixHandler
	        suffix_handler = handler_.OnIntSuffix(name, kind, num_values);
	    ReadSuffixValues<IntReader>(num_values, num_items, suffix_handler);
	  }
	}
	
	template <typename Reader, typename Handler>
	void NLReader<Reader, Handler>::Read(Reader *bound_reader) {
	  bool read_bounds = bound_reader == 0;
	  // TextReader::ReadHeader checks that this doesn't overflow.
	  num_vars_and_exprs_ = header_.num_vars +
	      header_.num_common_exprs_in_both +
	      header_.num_common_exprs_in_cons +
	      header_.num_common_exprs_in_objs +
	      header_.num_common_exprs_in_single_cons +
	      header_.num_common_exprs_in_single_objs;
	  for (;;) {
	    char c = reader_.ReadChar();
	    switch (c) {
	    case 'C': {
	      // Nonlinear part of an algebraic constraint body.
	      int index = ReadUInt(header_.num_algebraic_cons);
	      reader_.ReadTillEndOfLine();
	      handler_.OnAlgebraicCon(index, ReadNumericExpr(true));
	      break;
	    }
	    case 'L': {
	      // Logical constraint expression.
	      int index = ReadUInt(header_.num_logical_cons);
	      reader_.ReadTillEndOfLine();
	      handler_.OnLogicalCon(index, ReadLogicalExpr());
	      break;
	    }
	    case 'O': {
	      // Objective type and nonlinear part of an objective expression.
	      int index = ReadUInt(header_.num_objs);
	      int obj_type = reader_.ReadUInt();
	      reader_.ReadTillEndOfLine();
	      NumericExpr expr = ReadNumericExpr(true);
	      if (handler_.NeedObj(index))
	        handler_.OnObj(handler_.resulting_obj_index(index),
	                       obj_type != 0 ? obj::MAX : obj::MIN,
	                       expr);
	      break;
	    }
	    case 'V': {
	      // Defined variable definition (must precede V, C, L, O segments
	      // where used).
	      int expr_index = ReadUInt(header_.num_vars, num_vars_and_exprs_);
	      expr_index -= header_.num_vars;
	      int num_linear_terms = reader_.ReadUInt();
	      int position = reader_.ReadUInt();
	      reader_.ReadTillEndOfLine();
	      typename Handler::LinearExprHandler
	          expr_handler(handler_.BeginCommonExpr(expr_index, num_linear_terms));
	      if (num_linear_terms != 0)
	        ReadLinearExpr(num_linear_terms, expr_handler);
	      handler_.EndCommonExpr(expr_index, ReadNumericExpr(), position);
	      break;
	    }
	    case 'F': {
	      // Imported function description.
	      int index = ReadUInt(header_.num_funcs);
	      int type = reader_.ReadUInt();
	      if (type != func::NUMERIC && type != func::SYMBOLIC)
	        reader_.ReportError("invalid function type");
	      int num_args = reader_.template ReadInt<int>();
	      fmt::StringRef name = reader_.ReadName();
	      reader_.ReadTillEndOfLine();
	      handler_.OnFunction(index, name, num_args, static_cast<func::Type>(type));
	      break;
	    }
	    case 'G':
	      // Linear part of an objective expression & gradient sparsity.
	      ReadLinearExpr<ObjHandler>();
	      break;
	    case 'J':
	      // Jacobian sparsity & linear terms in constraints.
	      ReadLinearExpr<AlgebraicConHandler>();
	      break;
	    case 'S': {
	      // Suffix values.
	      int info = reader_.ReadUInt();
	      if (info > (internal::SUFFIX_KIND_MASK | suf::FLOAT))
	        reader_.ReportError("invalid suffix kind");
	      switch (info & internal::SUFFIX_KIND_MASK) {
	      case suf::VAR:
	        ReadSuffix<VarHandler>(info);
	        break;
	      case suf::CON:
	        ReadSuffix<ConHandler>(info);
	        break;
	      case suf::OBJ:
	        ReadSuffix<ObjHandler>(info);
	        break;
	      case suf::PROBLEM:
	        ReadSuffix<ProblemHandler>(info);
	        break;
	      }
	      break;
	    }
	    case 'b':
	      // Bounds on variables.
	      if (read_bounds) {
	        ReadBounds<VarHandler>();

	        reader_.ptr();
	        if ((flags_ & READ_BOUNDS_FIRST) != 0)
	          return;
	        read_bounds = false;
	        break;
	      }
	      if (!bound_reader)
	        reader_.ReportError("duplicate 'b' segment");
	      reader_ = *bound_reader;
	      bound_reader = 0;
	      break;
	    case 'r':
	      // Bounds on algebraic constraint bodies ("ranges").
	      ReadBounds<AlgebraicConHandler>();
	      break;
	    case 'K':
	      // Jacobian sparsity & linear constraint term matrix column sizes
	      // (must precede all J segments).
	      ReadColumnSizes<false>();
	      break;
	    case 'k':
	      // Jacobian sparsity & linear constraint term matrix cumulative column
	      // sizes (must precede all J segments).
	      ReadColumnSizes<true>();
	      break;
	    case 'x':
	      // Primal initial guess.
	      ReadInitialValues<VarHandler>();
	      break;
	    case 'd':
	      // Dual initial guess.
	      ReadInitialValues<AlgebraicConHandler>();
	      break;
	    case '\0':
	      if (reader_.IsEOF()) {
	        if (read_bounds)
	          reader_.ReportError("segment 'b' missing");
	        return;
	      }
	      // Fall through.
	    default:
	      reader_.ReportError("invalid segment type");
	    }
	  }
	}
	
	template <typename Reader, typename Handler>
	void NLReader<Reader, Handler>::Read() {
	  if ((flags_ & READ_BOUNDS_FIRST) != 0) {
	    // Read variable bounds first because this allows more efficient
	    // problem construction.
	    VarBoundHandler<Handler> bound_handler(handler_);
	    Reader bound_reader(reader_);
	    NLReader< Reader, VarBoundHandler<Handler> >
	        reader(bound_reader, header_, bound_handler, flags_);
	    reader.Read(0);
	    // Read everything else.
	    Read(&bound_reader);
	  } else {
	    Read(0);
	  }
	  handler_.EndInput();
	}
	
	// An .nl file reader.
	template <typename File = fmt::File>
	class NLFileReader {
	 private:
	  File file_;
	  std::size_t size_;
	  std::size_t rounded_size_;  // Size rounded up to a multiple of page size.
	
	  void Open(fmt::CStringRef filename);
	
	  // Reads the file into an array.
	  void Read(fmt::internal::MemoryBuffer<char, 1> &array);
	
	 public:
	  NLFileReader() : size_(0), rounded_size_(0) {}
	
	  const File &file() { return file_; }
	
	  // Opens and reads the file.
	  template <typename Handler>
	  void Read(fmt::CStringRef filename, Handler &handler, int flags) {
	    Open(filename);
	    if (size_ == rounded_size_) {
	      // Don't use mmap, because the file size is a multiple of the page size
	      // and therefore the mmap'ed buffer won't be null-terminated.
	      fmt::internal::MemoryBuffer<char, 1> array;
	      Read(array);
	      return ReadNLString(
	            NLStringRef(&array[0], size_), handler, filename, flags);
	    }
	    MemoryMappedFile<File> mapped_file(file_, rounded_size_);
	    ReadNLString(
	          NLStringRef(mapped_file.start(), size_), handler, filename, flags);
	  }
	};
	
	template <typename InputConverter, typename Handler>
	void ReadBinary(TextReader<> &reader, const NLHeader &header,
	                Handler &handler, int flags) {
	  BinaryReader<InputConverter> bin_reader(reader);
	  NLReader<BinaryReader<InputConverter>, Handler>(
	        bin_reader, header, handler, flags).Read();
	}
	
	template <typename NameHandler>
	void ReadNames(fmt::CStringRef filename, fmt::StringRef data,
	               NameHandler &handler) {
	  int line = 1;
	  const char *start = data.data();
	  const char *end = start + data.size();
	  for (const char *ptr = start; ptr != end; ++ptr) {
	    if (*ptr == '\n') {
	      handler.OnName(fmt::StringRef(start, ptr - start));
	      start = ptr + 1;
	      ++line;
	    }
	  }
	  if (start != end) {
	    int column = static_cast<int>(end - start + 1);
	    throw ReadError(filename, line, column, "missing newline");
	  }
	}
	
	// A name file reader.
	class NameReader {
	 private:
	  MemoryMappedFile<> mapped_file_;
	
	 public:
	  // Reads names from the file *filename* sending the names to the *handler*
	  // object by calling ``handler.OnName(name)``. The name argument to
	  // ``OnName`` is a ``fmt::StringRef`` object and the string it refers to
	  // is not null-terminated.
	  // Each name in the input file should be on a separate line ended with a
	  // newline character ('\n').
	  template <typename NameHandler>
	  void Read(fmt::CStringRef filename, NameHandler &handler) {
	    fmt::File file(filename, fmt::File::RDONLY);
	    mapped_file_.map(file, filename);
	    fmt::StringRef data(mapped_file_.start(), mapped_file_.size());
	    ReadNames(filename, data, handler);
	  }
	};
	
	// An NL handler that constructs an optimization problem using ProblemBuilder.
	template <typename ProblemBuilder>
	class NLProblemBuilder {
	 public:
	  typedef typename ProblemBuilder::Function Function;
	  typedef typename ProblemBuilder::Expr Expr;
	  typedef typename ProblemBuilder::NumericExpr NumericExpr;
	  typedef typename ProblemBuilder::LogicalExpr LogicalExpr;
	  typedef typename ProblemBuilder::CountExpr CountExpr;
	  typedef typename ProblemBuilder::Reference Reference;
	
	 private:
	  ProblemBuilder &builder_;
	
	  template <typename Obj>
	  void SetObj(const Obj &obj, obj::Type type, NumericExpr expr) {
	    obj.set_type(type);
	    obj.set_nonlinear_expr(expr);
	  }
	
	  // Sets bounds on a problem item (objective or constraint).
	  template <typename Item>
	  void SetBounds(const Item &item, double lb, double ub) {
	    item.set_lb(lb);
	    item.set_ub(ub);
	  }
	
	  template <typename CommonExpr>
	  void SetCommonExpr(const CommonExpr &common_expr,
	                     NumericExpr expr, int position) {
	    common_expr.set_nonlinear_expr(expr);
	    common_expr.set_position(position);
	  }
	
	 public:
	  explicit NLProblemBuilder(ProblemBuilder &builder): builder_(builder) {}
	
	  ProblemBuilder &builder() { return builder_; }
	
	  void OnHeader(const NLHeader &h) {
	    builder_.SetInfo(h);
	
	    // As nl-benchmark shows, adding problem components at once and then
	    // updating them is faster than adding them incrementally. The latter
	    // requires additional checks to make sure that prolbem components are
	    // in the correct order.
	    if (int n = h.num_continuous_vars())
	      builder_.AddVars(n, var::CONTINUOUS);
	    if (int n = h.num_integer_vars())
	      builder_.AddVars(n, var::INTEGER);
	    if (int n = h.num_common_exprs())
	      builder_.AddCommonExprs(n);
	    int n_objs = resulting_nobj( h.num_objs );
	    if (n_objs != 0)
	      builder_.AddObjs( n_objs );
	    if (h.num_algebraic_cons != 0)
	      builder_.AddAlgebraicCons(h.num_algebraic_cons);
	    if (h.num_logical_cons != 0)
	      builder_.AddLogicalCons(h.num_logical_cons);
	    if (h.num_funcs != 0)
	      builder_.AddFunctions(h.num_funcs);
	  }
	
	  virtual int objno() const { return 1; }
	  virtual bool multiobj() const { return true; }
	
	  int resulting_nobj(int nobj_header) const {
	    return multiobj() ? nobj_header :
	                        std::min( (objno()>0), (nobj_header>0) );
	  }
	  bool NeedObj(int obj_index) const {
	    return multiobj() || objno()-1==obj_index;
	  }
	  int resulting_obj_index(int index) const {
	    if (multiobj())
	      return index;
	    assert(objno()-1==index);
	    return 0;
	  }
	
	  void OnObj(int index, obj::Type type, NumericExpr expr) {
	    SetObj(builder_.obj(index), type, expr);
	  }
	
	  void OnAlgebraicCon(int index, NumericExpr expr) {
	    builder_.algebraic_con(index).set_nonlinear_expr(expr);
	  }
	
	  void OnLogicalCon(int index, LogicalExpr expr) {
	    builder_.logical_con(index).set_expr(expr);
	  }
	
	  void OnComplementarity(int con_index, int var_index, ComplInfo info) {
	    builder_.SetComplementarity(con_index, var_index, info);
	  }
	
	  typedef typename ProblemBuilder::LinearObjBuilder LinearObjHandler;
	
	  LinearObjHandler OnLinearObjExpr(int obj_index, int num_linear_terms) {
	    return builder_.obj(obj_index).set_linear_expr(num_linear_terms);
	  }
	
	  typedef typename ProblemBuilder::LinearConBuilder LinearConHandler;
	
	  LinearConHandler OnLinearConExpr(int con_index, int num_linear_terms) {
	    return builder_.algebraic_con(con_index).set_linear_expr(num_linear_terms);
	  }
	
	  typedef typename ProblemBuilder::LinearExprBuilder LinearExprHandler;
	
	  LinearExprHandler BeginCommonExpr(int index, int num_linear_terms) {
	    return builder_.common_expr(index).set_linear_expr(num_linear_terms);
	  }
	  void EndCommonExpr(int index, NumericExpr expr, int position) {
	    SetCommonExpr(builder_.common_expr(index), expr, position);
	  }
	
	  void OnVarBounds(int index, double lb, double ub) {
	    SetBounds(builder_.var(index), lb, ub);
	  }
	
	  void OnConBounds(int index, double lb, double ub) {
	    SetBounds(builder_.algebraic_con(index), lb, ub);
	  }
	
	  void OnInitialValue(int var_index, double value) {
	    builder_.var(var_index).set_value(value);
	  }
	
	  void OnInitialDualValue(int con_index, double value) {
	    builder_.algebraic_con(con_index).set_dual(value);
	  }
	
	  struct ColumnSizeHandler {
	    void Add(int) {
	      // Ignore column sizes as the constraints are stored row-wise by default.
	    }
	  };
	
	  ColumnSizeHandler OnColumnSizes() {
	    return ColumnSizeHandler();
	  }
	
	  void OnFunction(int index, fmt::StringRef name,
	                  int num_args, func::Type type) {
	    builder_.DefineFunction(index, name, num_args, type);
	  }
	
	  typedef typename ProblemBuilder::IntSuffixHandler IntSuffixHandler;
	
	  IntSuffixHandler OnIntSuffix(fmt::StringRef name, suf::Kind kind,
	                               int num_values) {
	    return builder_.AddIntSuffix(name, kind, num_values);
	  }
	
	  typedef typename ProblemBuilder::DblSuffixHandler DblSuffixHandler;
	
	  DblSuffixHandler OnDblSuffix(fmt::StringRef name, suf::Kind kind,
	                               int num_values) {
	    return builder_.AddDblSuffix(name, kind, num_values);
	  }
	
	  NumericExpr OnNumber(double value) {
	    return builder_.MakeNumericConstant(value);
	  }
	
	  Reference OnVariableRef(int var_index) {
	    return builder_.MakeVariable(var_index);
	  }
	
	  Reference OnCommonExprRef(int expr_index) {
	    return builder_.MakeCommonExpr(expr_index);
	  }
	
	  NumericExpr OnUnary(expr::Kind kind, NumericExpr arg) {
	    return builder_.MakeUnary(kind, arg);
	  }
	
	  NumericExpr OnBinary(expr::Kind kind, NumericExpr lhs, NumericExpr rhs) {
	    return builder_.MakeBinary(kind, lhs, rhs);
	  }
	
	  NumericExpr OnIf(LogicalExpr condition,
	      NumericExpr then_expr, NumericExpr else_expr) {
	    return builder_.MakeIf(condition, then_expr, else_expr);
	  }
	
	  typedef typename ProblemBuilder::PLTermBuilder PLTermHandler;
	
	  PLTermHandler BeginPLTerm(int num_breakpoints) {
	    return builder_.BeginPLTerm(num_breakpoints);
	  }
	  NumericExpr EndPLTerm(PLTermHandler handler, Reference arg) {
	    return builder_.EndPLTerm(handler, arg);
	  }
	
	  typedef typename ProblemBuilder::CallExprBuilder CallArgHandler;
	
	  CallArgHandler BeginCall(int func_index, int num_args) {
	    // Check if the function is defined.
	    if (Function func = builder_.function(func_index))
	      return builder_.BeginCall(func, num_args);
	    throw Error("function {} is not defined", func_index);
	  }
	  NumericExpr EndCall(CallArgHandler handler) {
	    return builder_.EndCall(handler);
	  }
	
	  typedef typename ProblemBuilder::IteratedExprBuilder VarArgHandler;
	
	  VarArgHandler BeginVarArg(expr::Kind kind, int num_args) {
	    return builder_.BeginIterated(kind, num_args);
	  }
	  NumericExpr EndVarArg(VarArgHandler handler) {
	    return builder_.EndIterated(handler);
	  }
	
	  typedef typename ProblemBuilder::IteratedExprBuilder NumericArgHandler;
	
	  NumericArgHandler BeginSum(int num_args) {
	    return builder_.BeginSum(num_args);
	  }
	  NumericExpr EndSum(NumericArgHandler handler) {
	    return builder_.EndSum(handler);
	  }
	
	  typedef typename ProblemBuilder::NumberOfExprBuilder NumberOfArgHandler;
	
	  NumberOfArgHandler BeginNumberOf(int num_args, NumericExpr arg0) {
	    return builder_.BeginNumberOf(num_args, arg0);
	  }
	  NumericExpr EndNumberOf(NumberOfArgHandler handler) {
	    return builder_.EndNumberOf(handler);
	  }
	
	  typedef typename ProblemBuilder::SymbolicNumberOfExprBuilder
	          SymbolicArgHandler;
	
	  SymbolicArgHandler BeginSymbolicNumberOf(int num_args, Expr arg0) {
	    return builder_.BeginSymbolicNumberOf(num_args, arg0);
	  }
	  NumericExpr EndSymbolicNumberOf(SymbolicArgHandler handler) {
	    return builder_.EndSymbolicNumberOf(handler);
	  }
	
	  typedef typename ProblemBuilder::CountExprBuilder CountArgHandler;
	
	  CountArgHandler BeginCount(int num_args) {
	    return builder_.BeginCount(num_args);
	  }
	  CountExpr EndCount(CountArgHandler handler) {
	    return builder_.EndCount(handler);
	  }
	
	  LogicalExpr OnBool(bool value) {
	    return builder_.MakeLogicalConstant(value);
	  }
	
	  LogicalExpr OnNot(LogicalExpr arg) {
	    return builder_.MakeNot(arg);
	  }
	
	  LogicalExpr OnBinaryLogical(
	      expr::Kind kind, LogicalExpr lhs, LogicalExpr rhs) {
	    return builder_.MakeBinaryLogical(kind, lhs, rhs);
	  }
	
	  LogicalExpr OnRelational(
	      expr::Kind kind, NumericExpr lhs, NumericExpr rhs) {
	    return builder_.MakeRelational(kind, lhs, rhs);
	  }
	
	  LogicalExpr OnLogicalCount(
	      expr::Kind kind, NumericExpr lhs, CountExpr rhs) {
	    return builder_.MakeLogicalCount(kind, lhs, rhs);
	  }
	
	  LogicalExpr OnImplication(
	      LogicalExpr condition, LogicalExpr then_expr, LogicalExpr else_expr) {
	    return builder_.MakeImplication(condition, then_expr, else_expr);
	  }
	
	  typedef typename ProblemBuilder::IteratedLogicalExprBuilder LogicalArgHandler;
	
	  LogicalArgHandler BeginIteratedLogical(expr::Kind kind, int num_args) {
	    return builder_.BeginIteratedLogical(kind, num_args);
	  }
	  LogicalExpr EndIteratedLogical(LogicalArgHandler handler) {
	    return builder_.EndIteratedLogical(handler);
	  }
	
	  typedef typename ProblemBuilder::PairwiseExprBuilder PairwiseArgHandler;
	
	  PairwiseArgHandler BeginPairwise(expr::Kind kind, int num_args) {
	    return builder_.BeginPairwise(kind, num_args);
	  }
	  LogicalExpr EndPairwise(PairwiseArgHandler handler) {
	    return builder_.EndPairwise(handler);
	  }
	
	  Expr OnString(fmt::StringRef value) {
	    return builder_.MakeStringLiteral(value);
	  }
	
	  Expr OnSymbolicIf(LogicalExpr condition, Expr then_expr, Expr else_expr) {
	    return builder_.MakeSymbolicIf(condition, then_expr, else_expr);
	  }
	
	  void EndInput() {}
	};
	
	template <typename Handler, bool>
	struct NLAdapter {
	  typedef Handler Type;
	  typedef Handler &RefType;
	};
	
	template <typename Handler>
	struct NLAdapter<Handler, true> {
	  typedef NLProblemBuilder<typename Handler::Builder> Type;
	  typedef Type RefType;
	};
	
	// SV added to suppress warning in following class
	// see also https://gcc.gnu.org/bugzilla/show_bug.cgi?id=71484#c5
	#ifdef __GNUC__
	#pragma GCC diagnostic ignored "-Wctor-dtor-privacy"
	#endif
	
	// Checks if T has a member type Builder.
	template <typename T>
	class HasBuilder {
	 private:
	  template <typename U> static fmt::internal::Yes &test(typename U::Builder *);
	  template <typename U> static fmt::internal::No &test(...);
	 public:
	  enum {value = sizeof(test<T>(0)) == sizeof(fmt::internal::Yes)};
	};
	}  // namespace internal
	
	template <typename Handler>
	void ReadNLString(NLStringRef str, Handler &handler,
	                  fmt::CStringRef name, int flags) {
	  // If handler is a Problem-like object (Problem::Builder type is defined)
	  // then use ProblemBuilder API to populate it.
	  typedef internal::NLAdapter<
	      Handler, internal::HasBuilder<Handler>::value> Adapter;
	  typename Adapter::RefType adapter(handler);
	  internal::TextReader<> reader(str, name);
	  NLHeader header = NLHeader();
	  reader.ReadHeader(header);
	  adapter.OnHeader(header);
	  switch (header.format) {
	  case NLHeader::TEXT:
	    internal::NLReader<internal::TextReader<>, typename Adapter::Type>(
	          reader, header, adapter, flags).Read();
	    break;
	  case NLHeader::BINARY: {
	      using internal::ReadBinary;
	    arith::Kind arith_kind = arith::GetKind();
	    if (arith_kind == header.arith_kind) {
	      ReadBinary<internal::IdentityConverter>(reader, header, adapter, flags);
	      break;
	    }
	    if (!IsIEEE(arith_kind) || !IsIEEE(header.arith_kind))
	      throw ReadError(name, 0, 0, "unsupported floating-point arithmetic");
	    ReadBinary<internal::EndiannessConverter>(reader, header, adapter, flags);
	    break;
	  }
	  }
	}
	
	template <typename Handler>
	inline void ReadNLFile(fmt::CStringRef filename, Handler &handler, int flags) {
	  internal::NLFileReader<>().Read(filename, handler, flags);
	}
	}  // namespace mp
	
	#endif  // MP_NL_READER_H_