main page modules namespaces classes files
Gecode home

search for 

support.icc

Go to the documentation of this file.

/*
 *  Main authors:
 *     Christian Schulte <schulte@gecode.org>
 *
 *  Copyright:
 *     Christian Schulte, 2004
 *
 *  Last modified:
 *     $Date: 2006-10-25 13:51:24 +0200 (Wed, 25 Oct 2006) $ by $Author: schulte $
 *     $Revision: 3787 $
 *
 *  This file is part of Gecode, the generic constraint
 *  development environment:
 *     http://www.gecode.org
 *
 *  See the file "LICENSE" for information on usage and
 *  redistribution of this file, and for a
 *     DISCLAIMER OF ALL WARRANTIES.
 *
 */

inline
Options::Options(const char* n)
  : icl(ICL_DEF),
    c_d(Search::Config::c_d),
    a_d(Search::Config::a_d),
    mode(EM_SOLUTION),
    quiet(false),
    samples(1),
    iterations(1),
    solutions(1),
    fails(-1),
    time(-1),
    naive(false),
    size(0),
    name(n) {}

namespace {
  class FailTimeStop : public Search::Stop {
  private:
    Search::TimeStop *ts;
    Search::FailStop *fs;
    FailTimeStop(int fails, int time) {
      ts = new Search::TimeStop(time);
      fs = new Search::FailStop(fails);
    }
  public:
    bool stop(const Search::Statistics& s) {
      return fs->stop(s) || ts->stop(s);
    }
    static Search::Stop* create(int fails, int time) {
      if (fails < 0 && time < 0) return NULL;
      if (fails < 0) return new Search::TimeStop( time);
      if (time  < 0) return new Search::FailStop(fails);
      return new FailTimeStop(fails, time);
    }
  };
}


template <class Script, template<class> class Engine>
void
Example::run(const Options& o) {
  using namespace std;
  try {
    if (!o.quiet)
      cout << o.name << endl;
    switch (o.mode) {
    case EM_SOLUTION:
      {
        Timer t;
        int i = o.solutions;
        t.start();
        Script* s = new Script(o);
        unsigned int n_p = 0;
        unsigned int n_b = 0;
        if (s->status() != SS_FAILED) {
          n_p = s->propagators();
          n_b = s->branchings();
        }
        Search::Stop* stop = FailTimeStop::create(o.fails, o.time);
        Engine<Script> e(s,o.c_d,o.a_d,stop);
        delete s;
        do {
          Example* ex = e.next();
          if (ex == NULL)
            break;
          ex->print();
          delete ex;
        } while (--i != 0);
        Search::Statistics stat = e.statistics();
        cout << endl;
        cout << "Initial" << endl
             << "\tpropagators:   " << n_p << endl
             << "\tbranchings:    " << n_b << endl
             << endl
             << "Summary" << endl
             << "\truntime:       " << t.stop() << endl
             << "\tsolutions:     " << abs(static_cast<int>(o.solutions) - i) << endl
             << "\tpropagations:  " << stat.propagate << endl
             << "\tfailures:      " << stat.fail << endl
             << "\tclones:        " << stat.clone << endl
             << "\tcommits:       " << stat.commit << endl
             << "\tpeak memory:   "
             << static_cast<int>((stat.memory+1023) / 1024) << " KB"
             << endl;
      }
      break;
    case EM_STAT:
      {
        int i = o.solutions;
        Script* s = new Script(o);
        unsigned int n_p = 0;
        unsigned int n_b = 0;
        if (s->status() != SS_FAILED) {
          n_p = s->propagators();
          n_b = s->branchings();
        }
        Engine<Script> e(s,o.c_d,o.a_d);
        delete s;
        do {
          Example* ex = e.next();
          if (ex == NULL)
            break;
          delete ex;
        } while (--i != 0);
        Search::Statistics stat = e.statistics();
        cout << endl
             << "\tpropagators:   " << n_p << endl
             << "\tbranchings:    " << n_b << endl
             << "\tsolutions:     " << abs(static_cast<int>(o.solutions) - i) << endl
             << "\tpropagations:  " << stat.propagate << endl
             << "\tfailures:      " << stat.fail << endl
             << "\tclones:        " << stat.clone << endl
             << "\tcommits:       " << stat.commit << endl
             << "\tpeak memory:   "
             << static_cast<int>((stat.memory+1023) / 1024) << " KB"
             << endl;
      }
      break;
    case EM_TIME:
      {
        Timer t;
        GECODE_AUTOARRAY(double,ts,o.samples);
        for (int s = o.samples; s--; ) {
          t.start();
          for (int k = o.iterations; k--; ) {
            unsigned int i = o.solutions;
            Script* s = new Script(o);
            Engine<Script> e(s,o.c_d,o.a_d);
            delete s;
            do {
              Example* ex = e.next();
              if (ex == NULL)
                break;
              delete ex;
            } while (--i != 0);
          }
          ts[s] = t.stop() / o.iterations;
        }
        double m = am(ts,o.samples);
        double d = dev(ts,o.samples) * 100.0;
        if (o.quiet) {
          cout << showpoint << setprecision(6) << m << " "
               << setprecision(2) << d
               << endl;
        } else {
          cout << "\tRuntime: "
               << setw(20) << right
               << showpoint << fixed
               << setprecision(6) << m << "ms"
               << setprecision(2) << " (" << d << "% deviation)"
               << endl;
        }
      }
      break;
    }
  } catch (Exception e) {
    cout << "Exception: " << e.what() << "." << endl
         << "Stopping..." << endl;
  }
}

// STATISTICS: example-any