C++ (Cpp) ParticleEventData Examples

Example #1

File: thermalCondSpeciesSpeciesE.cpp Project: pviswanathan/DynamO

 void 
 OPThermalConductivitySpeciesSpeciesE::updateConstDelG(const ParticleEventData& PDat)
 {
   double p1E = Sim->dynamics.getLiouvillean().getParticleKineticEnergy(PDat.getParticle());
 
   constDelG[PDat.getSpecies().getID()] += PDat.getParticle().getVelocity() * p1E 
     - PDat.getOldVel() * (p1E - PDat.getDeltaKE());
 }

Example #2

File: MFL.cpp Project: pviswanathan/DynamO

  void 
  OPMFL::A1ParticleChange(const ParticleEventData& PDat)
  {
    //We ignore stuff that hasn't had an event yet
    if (lastTime[PDat.getParticle().getID()] != 0.0)
      { 
	data[PDat.getSpecies().getID()]
	  .addVal(PDat.getParticle().getVelocity().nrm() 
		  * (Sim->dSysTime 
		     - lastTime[PDat.getParticle().getID()]));
      }
  
    lastTime[PDat.getParticle().getID()] = Sim->dSysTime;
  }

Example #3

File: newtonian.cpp Project: vladutzik/DynamO

    BOOST_FOREACH(const size_t& ID, range2)
    {
        ParticleEventData tmpval
        (Sim->particles[ID],
         *Sim->species[Sim->particles[ID]],
         eType);

        Sim->particles[ID].getVelocity() += dP / structmass2;

        tmpval.setDeltaKE(0.5 * Sim->species[tmpval.getSpeciesID()]->getMass(ID)
                          * (Sim->particles[ID].getVelocity().nrm2()
                             - tmpval.getOldVel().nrm2()));

        retVal.L1partChanges.push_back(tmpval);
    }

Example #4

File: MFT.cpp Project: pviswanathan/DynamO

  void 
  OPMFT::A1ParticleChange(const ParticleEventData& PDat)
  {
    //We ignore stuff that hasn't had an event yet

    for (size_t collN = 0; collN < collisionHistoryLength; ++collN)
      if (lastTime[PDat.getParticle().getID()][collN] != 0.0)
	{
	  data[PDat.getSpecies().getID()][collN]
	    .addVal(Sim->dSysTime 
		    - lastTime[PDat.getParticle().getID()][collN]);
	}
  
    lastTime[PDat.getParticle().getID()].push_front(Sim->dSysTime);
  }

Example #5

File: vacf.cpp Project: armando-2011/DynamO

  void 
  OPVACF::newG(const ParticleEventData& PDat)
  {
    if (Sim->dynamics.liouvilleanTypeTest<LSLLOD>())
      Sim->dynamics.getLiouvillean().updateAllParticles();

    for (size_t i = 0; i < Sim->N; ++i)
      G[i].push_front(Sim->particleList[i].getVelocity());	      
  
    //Now correct the fact that the wrong velocity has been pushed
    G[PDat.getParticle().getID()].front() = PDat.getOldVel();

    //This ensures the list gets to accumilator size
    if (notReady)
      {
	if (++currCorrLen != CorrelatorLength)
	  return;
      
	notReady = false;
      }

    accPass();
  }

Example #6

File: thermaldiffE.cpp Project: armando-2011/DynamO

  void 
  OPThermalDiffusionE::updateConstDelG(const ParticleEventData& PDat) 
  {
    double p1E = Sim->dynamics.getLiouvillean().getParticleKineticEnergy(PDat.getParticle());
  
    constDelG += PDat.getParticle().getVelocity() * p1E 
      - PDat.getOldVel() * (p1E - PDat.getDeltaKE());

    sysMom += PDat.getDeltaP();
  
    if (PDat.getSpecies().getID() == species1)
      constDelGsp1 += PDat.getDeltaP();
  }

Example #7

File: viscosityE.cpp Project: armando-2011/DynamO

  void 
  OPViscosityE::updateConstDelG(const ParticleEventData& PDat)
  {
    for (size_t iDim = 0; iDim < NDIM; ++iDim)
      for (size_t jDim = 0; jDim < NDIM; ++jDim)
	constDelG[iDim][jDim] += 
	  (PDat.getParticle().getVelocity()[iDim]
	   * PDat.getParticle().getVelocity()[jDim]
	   - PDat.getOldVel()[iDim] * PDat.getOldVel()[jDim]
	   ) * PDat.getSpecies().getMass(PDat.getParticle().getID());
  }

Example #8

File: uenergy.cpp Project: armando-2011/DynamO

 void 
 OPUEnergy::A1ParticleChange(const ParticleEventData& PDat)
 {
   intECurrent += PDat.getDeltaU();
 }

Example #9

File: collenergychange.cpp Project: pviswanathan/DynamO

 void 
 OPCollEnergyChange::A1ParticleChange(const ParticleEventData& PDat)
 {
   data[PDat.getSpecies().getID()]
     .addVal(PDat.getDeltaKE());
 }

Example #10

File: kenergy.cpp Project: armando-2011/DynamO

 void
 OPKEnergy::A1ParticleChange(const ParticleEventData& PDat)
 {
   KECurrent += PDat.getDeltaKE();
 }

Example #11

File: momentum.cpp Project: herbsolo/DynamO

void 
OPMomentum::A1ParticleChange(const ParticleEventData& PDat)
{
  sysMom += PDat.getDeltaP();
}

Example #12

File: 2particleEventData.hpp Project: MarkRunWu/DynamO

 EEventType getType() const
 { return particle1_.getType(); }

Example #13

File: 2particleEventData.hpp Project: MarkRunWu/DynamO

 void setType(EEventType nType)
 { particle1_.setType(nType); particle2_.setType(nType); }