flatzinc.hh

Go to the documentation of this file.

/* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */
/*
 *  Main authors:
 *     Guido Tack <tack@gecode.org>
 *
 *  Contributing authors:
 *     Gabriel Hjort Blindell <gabriel.hjort.blindell@gmail.com>
 *
 *  Copyright:
 *     Guido Tack, 2007-2012
 *     Gabriel Hjort Blindell, 2012
 *
 *  Last modified:
 *     $Date: 2015-03-19 11:47:57 +0100 (Thu, 19 Mar 2015) $ by $Author: tack $
 *     $Revision: 14463 $
 *
 *  This file is part of Gecode, the generic constraint
 *  development environment:
 *     http://www.gecode.org
 *
 *  Permission is hereby granted, free of charge, to any person obtaining
 *  a copy of this software and associated documentation files (the
 *  "Software"), to deal in the Software without restriction, including
 *  without limitation the rights to use, copy, modify, merge, publish,
 *  distribute, sublicense, and/or sell copies of the Software, and to
 *  permit persons to whom the Software is furnished to do so, subject to
 *  the following conditions:
 *
 *  The above copyright notice and this permission notice shall be
 *  included in all copies or substantial portions of the Software.
 *
 *  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
 *  EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
 *  MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
 *  NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
 *  LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
 *  OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
 *  WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
 *
 */

#ifndef __GECODE_FLATZINC_HH__
#define __GECODE_FLATZINC_HH__

#include <iostream>

#include <gecode/kernel.hh>
#include <gecode/int.hh>
#ifdef GECODE_HAS_SET_VARS
#include <gecode/set.hh>
#endif
#ifdef GECODE_HAS_FLOAT_VARS
#include <gecode/float.hh>
#endif
#include <map>

/*
 * Support for DLLs under Windows
 *
 */

#if !defined(GECODE_STATIC_LIBS) && \
    (defined(__CYGWIN__) || defined(__MINGW32__) || defined(_MSC_VER))

#ifdef GECODE_BUILD_FLATZINC
#define GECODE_FLATZINC_EXPORT __declspec( dllexport )
#else
#define GECODE_FLATZINC_EXPORT __declspec( dllimport )
#endif

#else

#ifdef GECODE_GCC_HAS_CLASS_VISIBILITY

#define GECODE_FLATZINC_EXPORT __attribute__ ((visibility("default")))

#else

#define GECODE_FLATZINC_EXPORT

#endif
#endif

// Configure auto-linking
#ifndef GECODE_BUILD_FLATZINC
#define GECODE_LIBRARY_NAME "FlatZinc"
#include <gecode/support/auto-link.hpp>
#endif

#include <gecode/driver.hh>

#include <gecode/flatzinc/conexpr.hh>
#include <gecode/flatzinc/ast.hh>
#include <gecode/flatzinc/varspec.hh>

namespace Gecode { namespace FlatZinc {

  class GECODE_FLATZINC_EXPORT Printer {
  private:
    AST::Array* _output;
    void printElem(std::ostream& out,
                   AST::Node* ai,
                   const Gecode::IntVarArray& iv,
                   const Gecode::BoolVarArray& bv
#ifdef GECODE_HAS_SET_VARS
                   ,
                   const Gecode::SetVarArray& sv
#endif
#ifdef GECODE_HAS_FLOAT_VARS
                  ,
                  const Gecode::FloatVarArray& fv
#endif
                   ) const;
    void printElemDiff(std::ostream& out,
                       AST::Node* ai,
                       const Gecode::IntVarArray& iv1,
                       const Gecode::IntVarArray& iv2,
                       const Gecode::BoolVarArray& bv1,
                       const Gecode::BoolVarArray& bv2
#ifdef GECODE_HAS_SET_VARS
                       ,
                       const Gecode::SetVarArray& sv1,
                       const Gecode::SetVarArray& sv2
#endif
#ifdef GECODE_HAS_FLOAT_VARS
                       ,
                       const Gecode::FloatVarArray& fv1,
                       const Gecode::FloatVarArray& fv2
#endif
                       ) const;
  public:
    Printer(void) : _output(NULL) {}
    void init(AST::Array* output);

    void print(std::ostream& out,
               const Gecode::IntVarArray& iv,
               const Gecode::BoolVarArray& bv
#ifdef GECODE_HAS_SET_VARS
               ,
               const Gecode::SetVarArray& sv
#endif
#ifdef GECODE_HAS_FLOAT_VARS
               ,
               const Gecode::FloatVarArray& fv
#endif
               ) const;

    void printDiff(std::ostream& out,
               const Gecode::IntVarArray& iv1, const Gecode::IntVarArray& iv2,
               const Gecode::BoolVarArray& bv1, const Gecode::BoolVarArray& bv2
#ifdef GECODE_HAS_SET_VARS
               ,
               const Gecode::SetVarArray& sv1, const Gecode::SetVarArray& sv2
#endif
#ifdef GECODE_HAS_FLOAT_VARS
               ,
               const Gecode::FloatVarArray& fv1,
               const Gecode::FloatVarArray& fv2
#endif
               ) const;

  
    ~Printer(void);
    
    void shrinkElement(AST::Node* node,
                       std::map<int,int>& iv, std::map<int,int>& bv, 
                       std::map<int,int>& sv, std::map<int,int>& fv);

    void shrinkArrays(Space& home,
                      int& optVar, bool optVarIsInt,
                      Gecode::IntVarArray& iv,
                      Gecode::BoolVarArray& bv
#ifdef GECODE_HAS_SET_VARS
                      ,
                      Gecode::SetVarArray& sv
#endif
#ifdef GECODE_HAS_FLOAT_VARS
                      ,
                      Gecode::FloatVarArray& fv
#endif
                     );
    
  private:
    Printer(const Printer&);
    Printer& operator=(const Printer&);
  };

  class FlatZincOptions : public Gecode::BaseOptions {
  protected:

      Gecode::Driver::IntOption         _solutions; 
      Gecode::Driver::BoolOption        _allSolutions; 
      Gecode::Driver::DoubleOption      _threads;   
      Gecode::Driver::BoolOption        _free; 
      Gecode::Driver::DoubleOption      _decay;       
      Gecode::Driver::UnsignedIntOption _c_d;       
      Gecode::Driver::UnsignedIntOption _a_d;       
      Gecode::Driver::UnsignedIntOption _node;      
      Gecode::Driver::UnsignedIntOption _fail;      
      Gecode::Driver::UnsignedIntOption _time;      
      Gecode::Driver::IntOption         _seed;      
      Gecode::Driver::StringOption      _restart;   
      Gecode::Driver::DoubleOption      _r_base;    
      Gecode::Driver::UnsignedIntOption _r_scale;   
      Gecode::Driver::BoolOption        _nogoods;   
      Gecode::Driver::UnsignedIntOption _nogoods_limit; 
      Gecode::Driver::BoolOption        _interrupt; 
      Gecode::Driver::DoubleOption      _step;        

    

      Gecode::Driver::StringOption      _mode;       
      Gecode::Driver::BoolOption        _stat;       
      Gecode::Driver::StringValueOption _output;     

  public:
    FlatZincOptions(const char* s)
    : Gecode::BaseOptions(s),
      _solutions("-n","number of solutions (0 = all, -1 = one/best)",-1),
      _allSolutions("-a", "return all solutions (equal to -n 0)"),
      _threads("-p","number of threads (0 = #processing units)",
               Gecode::Search::Config::threads),
      _free("-f", "free search, no need to follow search-specification"),
      _decay("-decay","decay factor",0.99),
      _c_d("-c-d","recomputation commit distance",Gecode::Search::Config::c_d),
      _a_d("-a-d","recomputation adaption distance",Gecode::Search::Config::a_d),
      _node("-node","node cutoff (0 = none, solution mode)"),
      _fail("-fail","failure cutoff (0 = none, solution mode)"),
      _time("-time","time (in ms) cutoff (0 = none, solution mode)"),
      _seed("-r","random seed",0),
      _restart("-restart","restart sequence type",RM_NONE),
      _r_base("-restart-base","base for geometric restart sequence",1.5),
      _r_scale("-restart-scale","scale factor for restart sequence",250),
      _nogoods("-nogoods","whether to use no-goods from restarts",false),
      _nogoods_limit("-nogoods-limit","depth limit for no-good extraction",
                     Search::Config::nogoods_limit),
      _interrupt("-interrupt","whether to catch Ctrl-C (true) or not (false)",
                 true),
      _step("-step","step distance for float optimization",0.0),
      _mode("-mode","how to execute script",Gecode::SM_SOLUTION),
      _stat("-s","emit statistics"),
      _output("-o","file to send output to") {

      _mode.add(Gecode::SM_SOLUTION, "solution");
      _mode.add(Gecode::SM_STAT, "stat");
      _mode.add(Gecode::SM_GIST, "gist");
      _restart.add(RM_NONE,"none");
      _restart.add(RM_CONSTANT,"constant");
      _restart.add(RM_LINEAR,"linear");
      _restart.add(RM_LUBY,"luby");
      _restart.add(RM_GEOMETRIC,"geometric");

      add(_solutions); add(_threads); add(_c_d); add(_a_d);
      add(_allSolutions);
      add(_free);
      add(_decay);
      add(_node); add(_fail); add(_time); add(_interrupt);
      add(_seed);
      add(_step);
      add(_restart); add(_r_base); add(_r_scale); 
      add(_nogoods); add(_nogoods_limit);
      add(_mode); add(_stat);
      add(_output);
    }

    void parse(int& argc, char* argv[]) {
      Gecode::BaseOptions::parse(argc,argv);
      if (_allSolutions.value() && _solutions.value()==-1) {
        _solutions.value(0);
      }
      if (_stat.value())
        _mode.value(Gecode::SM_STAT);
    }
  
    virtual void help(void) {
      std::cerr << "Gecode FlatZinc interpreter" << std::endl
                << " - Supported FlatZinc version: " << GECODE_FLATZINC_VERSION
                << std::endl << std::endl;
      Gecode::BaseOptions::help();
    }
  
    int solutions(void) const { return _solutions.value(); }
    bool allSolutions(void) const { return _allSolutions.value(); }
    double threads(void) const { return _threads.value(); }
    bool free(void) const { return _free.value(); }
    unsigned int c_d(void) const { return _c_d.value(); }
    unsigned int a_d(void) const { return _a_d.value(); }
    unsigned int node(void) const { return _node.value(); }
    unsigned int fail(void) const { return _fail.value(); }
    unsigned int time(void) const { return _time.value(); }
    int seed(void) const { return _seed.value(); }
    double step(void) const { return _step.value(); }
    const char* output(void) const { return _output.value(); }
    Gecode::ScriptMode mode(void) const {
      return static_cast<Gecode::ScriptMode>(_mode.value());
    }

    double decay(void) const { return _decay.value(); }
    RestartMode restart(void) const {
      return static_cast<RestartMode>(_restart.value());
    }
    double restart_base(void) const { return _r_base.value(); }
    unsigned int restart_scale(void) const { return _r_scale.value(); }
    bool nogoods(void) const { return _nogoods.value(); }
    unsigned int nogoods_limit(void) const { return _nogoods_limit.value(); }
    bool interrupt(void) const { return _interrupt.value(); }

    void allSolutions(bool b) { _allSolutions.value(b); }
  };

  class BranchInformation : public SharedHandle {
  public:
    BranchInformation(void);
    BranchInformation(const BranchInformation& bi);
    void init(void);
    void add(const BrancherHandle& bh,
             const std::string& rel0,
             const std::string& rel1,
             const std::vector<std::string>& n);
    void print(const BrancherHandle& bh,
               int a, int i, int n, std::ostream& o) const;
#ifdef GECODE_HAS_FLOAT_VARS

    void print(const BrancherHandle& bh,
               int a, int i, const FloatNumBranch& nl, std::ostream& o) const;
#endif
  };

 class GECODE_FLATZINC_EXPORT FznRnd {
 protected:
   Gecode::Support::RandomGenerator random;
   Gecode::Support::Mutex mutex;
 public:
   FznRnd(unsigned int s=1);
   unsigned int operator ()(unsigned int n);
 };

  class GECODE_FLATZINC_EXPORT FlatZincSpace : public Space {
  public:
    enum Meth {
      SAT, //< Solve as satisfaction problem
      MIN, //< Solve as minimization problem
      MAX  //< Solve as maximization problem
    };
  protected:
    int intVarCount;
    int boolVarCount;
    int floatVarCount;
    int setVarCount;

    int _optVar;
    bool _optVarIsInt;
  
    Meth _method;

    unsigned int _lns;
    
    FznRnd* _random;

    AST::Array* _solveAnnotations;

    FlatZincSpace(bool share, FlatZincSpace&);
  private:
    template<template<class> class Engine>
    void
    runEngine(std::ostream& out, const Printer& p, 
              const FlatZincOptions& opt, Gecode::Support::Timer& t_total);
    template<template<class> class Engine,
             template<template<class> class,class> class Meta>
    void
    runMeta(std::ostream& out, const Printer& p, 
            const FlatZincOptions& opt, Gecode::Support::Timer& t_total);
    void
    branchWithPlugin(AST::Node* ann);
  public:
    Gecode::IntVarArray iv;
    Gecode::IntVarArray iv_aux;
    
    Gecode::IntVarArray iv_lns;

    std::vector<bool> iv_introduced;
    int* iv_boolalias;
    Gecode::BoolVarArray bv;
    Gecode::BoolVarArray bv_aux;
    std::vector<bool> bv_introduced;
#ifdef GECODE_HAS_SET_VARS

    Gecode::SetVarArray sv;
    Gecode::SetVarArray sv_aux;
    std::vector<bool> sv_introduced;
#endif
#ifdef GECODE_HAS_FLOAT_VARS

    Gecode::FloatVarArray fv;
    Gecode::FloatVarArray fv_aux;
    std::vector<bool> fv_introduced;
    Gecode::FloatNum step;
#endif

    bool needAuxVars;
    FlatZincSpace(FznRnd* random = NULL);
  
    ~FlatZincSpace(void);
  
    void init(int intVars, int boolVars, int setVars, int floatVars);

    void newIntVar(IntVarSpec* vs);
    void aliasBool2Int(int iv, int bv);
    int aliasBool2Int(int iv);
    void newBoolVar(BoolVarSpec* vs);
    void newSetVar(SetVarSpec* vs);
    void newFloatVar(FloatVarSpec* vs);
  
    void postConstraints(std::vector<ConExpr*>& ces);
  
    void solve(AST::Array* annotation);
    void minimize(int var, bool isInt, AST::Array* annotation);
    void maximize(int var, bool isInt, AST::Array* annotation);

    void run(std::ostream& out, const Printer& p, 
             const FlatZincOptions& opt, Gecode::Support::Timer& t_total);
  
    void print(std::ostream& out, const Printer& p) const;

    void compare(const Space& s, std::ostream& out) const;
    void compare(const FlatZincSpace& s, std::ostream& out,
                 const Printer& p) const;

    void shrinkArrays(Printer& p);

    Meth method(void) const;

    int optVar(void) const;
    bool optVarIsInt(void) const;

    void createBranchers(AST::Node* ann,
                         int seed, double decay,
                         bool ignoreUnknown,
                         std::ostream& err = std::cerr);

    AST::Array* solveAnnotations(void) const;

    BranchInformation branchInfo;

    virtual void constrain(const Space& s);
    virtual Gecode::Space* copy(bool share);

    virtual bool slave(const CRI& cri);
    


    IntArgs arg2intargs(AST::Node* arg, int offset = 0);
    IntArgs arg2boolargs(AST::Node* arg, int offset = 0);
    IntSet arg2intset(AST::Node* n);
    IntSetArgs arg2intsetargs(AST::Node* arg, int offset = 0);
    IntVarArgs arg2intvarargs(AST::Node* arg, int offset = 0);
    BoolVarArgs arg2boolvarargs(AST::Node* arg, int offset = 0, int siv=-1);
    BoolVar arg2BoolVar(AST::Node* n);
    IntVar arg2IntVar(AST::Node* n);
    bool isBoolArray(AST::Node* b, int& singleInt);
#ifdef GECODE_HAS_SET_VARS

    SetVar arg2SetVar(AST::Node* n);
    SetVarArgs arg2setvarargs(AST::Node* arg, int offset = 0, int doffset = 0,
                              const IntSet& od=IntSet::empty);
#endif
#ifdef GECODE_HAS_FLOAT_VARS

    FloatValArgs arg2floatargs(AST::Node* arg, int offset = 0);
    FloatVar arg2FloatVar(AST::Node* n);
    FloatVarArgs arg2floatvarargs(AST::Node* arg, int offset = 0);
#endif

    IntConLevel ann2icl(AST::Node* ann);
  };

  class GECODE_VTABLE_EXPORT Error {
  private:
    const std::string msg;
  public:
    Error(const std::string& where, const std::string& what)
    : msg(where+": "+what) {}
    const std::string& toString(void) const { return msg; }
  };

  GECODE_FLATZINC_EXPORT
  FlatZincSpace* parse(const std::string& fileName,
                       Printer& p, std::ostream& err = std::cerr,
                       FlatZincSpace* fzs=NULL, FznRnd* rnd=NULL);

  GECODE_FLATZINC_EXPORT
  FlatZincSpace* parse(std::istream& is,
                       Printer& p, std::ostream& err = std::cerr,
                       FlatZincSpace* fzs=NULL, FznRnd* rnd=NULL);

}}

#endif

// STATISTICS: flatzinc-any