o815/sim.cpp

   1 #include "o815.h"
   2
   3 #include "latlib/hypercache.h"
   4 #include "latlib/progress.h"
   5
   6 o815::sim::sim(o815 *_O815, const int& confmemSize) {
   7   O815 = _O815;
   8   log = O815->out->log;
   9   hypercache::initC(O815->programid, O815->comargs.nequi, O815->comargs.nskip, O815->comargs.confcache.first, &confMem, confmemSize, O815->comargs.confcache.second, log);
  10   toEquilibrate = true;
  11   nequi = O815->comargs.nequi;
  12   nskip = O815->comargs.nskip;
  13 }
  14
  15 void o815::sim::nextConfig() {
  16   int nequileft = hypercache::readC();
  17
  18   if ( nequileft != -1 ) {
  19     if(toEquilibrate) {
  20       _reset();
  21       *log << "SIM: Starting equilibration." << endl << flush;
  22       if (nequileft != nequi)
  23         *log << "SIM: " << nequileft << " of " << nequi << " equilibration steps left after virtual equilibration." << endl << flush;
  24       progress equiProg(nequileft);
  25       for( int iequi=0; iequi<nequileft; iequi++ ) {
  26         _makeSweep(1);
  27         while( equiProg.madeStep(iequi) ) *log << "SIM: " << equiProg.done()*100 << "% of equilibration done." << endl << flush;
  28       }
  29     }
  30     _makeSweep(nskip);
  31     hypercache::writeC();
  32   }
  33
  34   toEquilibrate = false;
  35 }