Exemple #1
0
  void
  SSnapshot::runEvent() const
  {
    double locdt = dt;
  
#ifdef DYNAMO_DEBUG 
    if (boost::math::isnan(dt))
      M_throw() << "A NAN system event time has been found";
#endif
    

    Sim->dSysTime += locdt;

    Sim->ptrScheduler->stream(locdt);
  
    //dynamics must be updated first
    Sim->dynamics.stream(locdt);
  
    dt += _period;
  
    locdt += Sim->freestreamAcc;
    Sim->freestreamAcc = 0;

    //This is done here as most ticker properties require it
    Sim->dynamics.getLiouvillean().updateAllParticles();

    BOOST_FOREACH(shared_ptr<OutputPlugin>& Ptr, Sim->outputPlugins)
      Ptr->eventUpdate(*this, NEventData(), locdt);
  
    std::string filename = magnet::string::search_replace("Snapshot.%i.xml.bz2", "%i", boost::lexical_cast<std::string>(_saveCounter++));
    Sim->writeXMLfile(filename, _applyBC);
  }
Exemple #2
0
  void
  SysTicker::runEvent() const
  {
    double locdt = dt;
  
#ifdef DYNAMO_DEBUG 
    if (std::isnan(dt))
      M_throw() << "A NAN system event time has been found";
#endif
    

    Sim->systemTime += locdt;

    Sim->ptrScheduler->stream(locdt);
  
    //dynamics must be updated first
    Sim->stream(locdt);
  
    dt += period;
  
    //This is done here as most ticker properties require it
    Sim->dynamics->updateAllParticles();

    for (shared_ptr<OutputPlugin>& Ptr : Sim->outputPlugins)
      {
	shared_ptr<OPTicker> ptr = std::dynamic_pointer_cast<OPTicker>(Ptr);
	if (ptr) ptr->ticker();
      }

    for (shared_ptr<OutputPlugin>& Ptr : Sim->outputPlugins)
      Ptr->eventUpdate(*this, NEventData(), locdt);
  }
Exemple #3
0
  NEventData
  SysTicker::runEvent()
  {
    dt += period;  
    //This is done here as most ticker properties require it
    Sim->dynamics->updateAllParticles();
    for (shared_ptr<OutputPlugin>& Ptr : Sim->outputPlugins)
      {
	shared_ptr<OPTicker> ptr = std::dynamic_pointer_cast<OPTicker>(Ptr);
	if (ptr) ptr->ticker();
      }
    return NEventData();
  }
Exemple #4
0
  NEventData
  SysSnapshot::runEvent()
  {
    if (_eventPeriod) 
      dt = HUGE_VAL;
    else
      dt += _period;

    Sim->dynamics->updateAllParticles();
  
    std::string filename = magnet::string::search_replace("Snapshot."+_format+".xml.bz2", "%COUNT", boost::lexical_cast<std::string>(_saveCounter));
    filename = magnet::string::search_replace(filename, "%ID", boost::lexical_cast<std::string>(Sim->stateID));
    Sim->writeXMLfile(filename, _applyBC);
    
    dout << "Printing SNAPSHOT" << std::endl;
    
    filename = magnet::string::search_replace("Snapshot.output."+_format+".xml.bz2", "%COUNT", boost::lexical_cast<std::string>(_saveCounter++));
    filename = magnet::string::search_replace(filename, "%ID", boost::lexical_cast<std::string>(Sim->stateID));
    Sim->outputData(filename);
    return NEventData();
  }
Exemple #5
0
  void
  SysSnapshot::runEvent()
  {
    double locdt = dt;
    if (_eventPeriod) 
      {
	locdt = 0;
	dt = HUGE_VAL;
      }
    else
      dt += _period;

#ifdef DYNAMO_DEBUG 
    if (std::isnan(locdt))
      M_throw() << "A NAN system event time has been found";
#endif
    
    Sim->systemTime += locdt;

    Sim->ptrScheduler->stream(locdt);
  
    //dynamics must be updated first
    Sim->stream(locdt);
    
    //This is done here as most ticker properties require it
    Sim->dynamics->updateAllParticles();

    for (shared_ptr<OutputPlugin>& Ptr : Sim->outputPlugins)
      Ptr->eventUpdate(*this, NEventData(), locdt);
  
    std::string filename = magnet::string::search_replace("Snapshot."+_format+".xml.bz2", "%COUNT", boost::lexical_cast<std::string>(_saveCounter));
    filename = magnet::string::search_replace(filename, "%ID", boost::lexical_cast<std::string>(Sim->simID));
    Sim->writeXMLfile(filename, _applyBC);
    
    dout << "Printing SNAPSHOT" << std::endl;
    
    filename = magnet::string::search_replace("Snapshot.output."+_format+".xml.bz2", "%COUNT", boost::lexical_cast<std::string>(_saveCounter++));
    filename = magnet::string::search_replace(filename, "%ID", boost::lexical_cast<std::string>(Sim->simID));
    Sim->outputData(filename);
  }
Exemple #6
0
  NEventData
  SVisualizer::runEvent()
  {
    if (_window->dynamoParticleSync())
      Sim->dynamics->updateAllParticles();
  
    for (shared_ptr<System>& system : Sim->systems)
      {
	shared_ptr<SysRotateGravity> rgrav =  std::dynamic_pointer_cast<SysRotateGravity>(system);
	if (!rgrav) continue;

	std::shared_ptr<DynGravity> dynamics = std::dynamic_pointer_cast<DynGravity>(Sim->dynamics);
	if (dynamics)
	  _window->getGLContext()->queueTask(std::bind(&magnet::GL::Camera::setUp, &(_window->getCamera()), -dynamics->getGravityVector(), rgrav->getAxis()));
	break;
      }

    _window->simupdateTick(Sim->systemTime / Sim->units.unitTime());

    //Now that the update has been performed, set up the next "tick"
    dt = _window->getUpdateInterval() * Sim->units.unitTime();
    _lastUpdate = boost::posix_time::microsec_clock::local_time();
    return NEventData();
  }
Exemple #7
0
  void
  SVisualizer::runEvent()
  {
    //Dont rewind in time, the -HUGE_VAL time is only used to ensure
    //the event takes place before any event, including negative time events.
    if (dt == -HUGE_VAL) dt = 0;
  
    //Actually move forward the system time
    Sim->systemTime += dt;
    Sim->ptrScheduler->stream(dt);
    //dynamics must be updated first
    Sim->stream(dt);

    if (_window->dynamoParticleSync())
      Sim->dynamics->updateAllParticles();

    for (shared_ptr<OutputPlugin>& Ptr : Sim->outputPlugins)
      Ptr->eventUpdate(*this, NEventData(), dt);
  
    for (shared_ptr<System>& system : Sim->systems)
      {
	shared_ptr<SysRotateGravity> rgrav =  std::dynamic_pointer_cast<SysRotateGravity>(system);
	if (!rgrav) continue;

	std::shared_ptr<DynGravity> dynamics = std::dynamic_pointer_cast<DynGravity>(Sim->dynamics);
	if (dynamics)
	  _window->getGLContext()->queueTask(std::bind(&magnet::GL::Camera::setUp, &(_window->getCamera()), -dynamics->getGravityVector(), rgrav->getAxis()));
	break;
      }

    _window->simupdateTick(Sim->systemTime / Sim->units.unitTime());

    //Now that the update has been performed, set up the next "tick"
    dt = _window->getUpdateInterval() * Sim->units.unitTime();
    _lastUpdate = boost::posix_time::microsec_clock::local_time();
  }
Exemple #8
0
  void 
  CSRingDSMC::runEvent() const
  {
    double locdt = dt;
  
#ifdef DYNAMO_DEBUG 
    if (boost::math::isnan(locdt))
      M_throw() << "A NAN system event time has been found";
#endif
    
    Sim->dSysTime += locdt;
    
    Sim->ptrScheduler->stream(locdt);
  
    //dynamics must be updated first
    Sim->dynamics.stream(locdt);

    dt = tstep;

    locdt += Sim->freestreamAcc;
    Sim->freestreamAcc = 0;

    BOOST_FOREACH(std::tr1::shared_ptr<OutputPlugin>& Ptr, Sim->outputPlugins)
      Ptr->eventUpdate(*this, NEventData(), locdt);

    //////////////////// T(1,2) operator
    double intPart;
    double fracpart = std::modf(maxprob12 * range1->size(), &intPart);
 
    size_t nmax = static_cast<size_t>(intPart) + (Sim->uniform_sampler() < fracpart);
  
    {
      boost::variate_generator
	<dynamo::baseRNG&, boost::uniform_int<size_t> >
	id1sampler(Sim->ranGenerator, 
		   boost::uniform_int<size_t>(0, (range1->size()/2) - 1));
    
      for (size_t n = 0; n < nmax; ++n)
	{
	  size_t pairID(id1sampler());
	  const Particle& p1(Sim->particleList[*(range1->begin() + 2 * pairID)]);
	  const Particle& p2(Sim->particleList[*(range1->begin() + 2 * pairID + 1)]);
	
	  Sim->dynamics.getLiouvillean().updateParticlePair(p1, p2);
	
	  CPDData PDat;
	
	  for (size_t iDim(0); iDim < NDIM; ++iDim)
	    PDat.rij[iDim] = Sim->normal_sampler();
	
	  PDat.rij *= diameter / PDat.rij.nrm();
	
	  if (Sim->dynamics.getLiouvillean().DSMCSpheresTest
	      (p1, p2, maxprob12, factor12, PDat))
	    {
	      ++Sim->eventCount;
	      ++n12;

	      const PairEventData
		SDat(Sim->dynamics.getLiouvillean().DSMCSpheresRun(p1, p2, e, PDat));
	    
	      Sim->signalParticleUpdate(SDat);
	    
	      Sim->ptrScheduler->fullUpdate(p1, p2);
	    
	      BOOST_FOREACH(std::tr1::shared_ptr<OutputPlugin>& Ptr, Sim->outputPlugins)
		Ptr->eventUpdate(*this, SDat, 0.0);
	    }
	}
    }

    //////////////////// T(1,3) operator
    {
      fracpart = std::modf(maxprob13 * range1->size(), &intPart);
    
      nmax = static_cast<size_t>(intPart) + (Sim->uniform_sampler() < fracpart);
    
      boost::variate_generator
	<dynamo::baseRNG&, boost::uniform_int<size_t> >
	id1sampler(Sim->ranGenerator, 
		   boost::uniform_int<size_t>(0, range1->size() - 1));
    
      for (size_t n = 0; n < nmax; ++n)
	{
	  const Particle& p1(Sim->particleList[*(range1->begin() + id1sampler())]);
	
	  size_t secondID(id1sampler());
	
	  while ((secondID == p1.getID())
		 || ((secondID % 2) 
		     ? ((secondID-1) == p1.getID())
		     : ((secondID+1) == p1.getID())))
	    secondID = id1sampler();
	
	  const Particle& p2(Sim->particleList[*(range1->begin() + secondID)]);
	
	  Sim->dynamics.getLiouvillean().updateParticlePair(p1, p2);
	
	  CPDData PDat;
	
	  for (size_t iDim(0); iDim < NDIM; ++iDim)
	    PDat.rij[iDim] = Sim->normal_sampler();
	
	  PDat.rij *= diameter / PDat.rij.nrm();
	
	  if (Sim->dynamics.getLiouvillean().DSMCSpheresTest
	      (p1, p2, maxprob13, factor13, PDat))
	    {
	      ++Sim->eventCount;
	      ++n13;

	      const PairEventData
		SDat(Sim->dynamics.getLiouvillean().DSMCSpheresRun(p1, p2, e, PDat));

	      Sim->signalParticleUpdate(SDat);
	    
	      Sim->ptrScheduler->fullUpdate(p1, p2);
	    
	      BOOST_FOREACH(std::tr1::shared_ptr<OutputPlugin>& Ptr, Sim->outputPlugins)
		Ptr->eventUpdate(*this, SDat, 0.0);
	    }
	}
    }
  }