Method emitParticlesCircle

SPlisHSPlasH/Emitter.cpp:231–380

void Emitter::emitParticlesCircle(std::vector <unsigned int> &reusedParticles, unsigned int &indexReuse, unsigned int &numEmittedParticles)
{
	TimeManager *tm = TimeManager::getCurrent();
	const Real t = tm->getTime();
	const Real timeStepSize = tm->getTimeStepSize();
	const Vector3r & emitDir = m_rotation.col(0);
	Simulation *sim = Simulation::getCurrent();
	const Real particleRadius = sim->getParticleRadius();
	const Real diam = static_cast<Real>(2.0)*particleRadius;
	Vector3r emitVel = m_velocity * emitDir;

	// shortly before the emitter starts, cleanup the emitter from particles
	if (t < m_emitStartTime || t > m_emitEndTime)
		emitVel = emitDir * particleRadius * 10 / 0.25;

	if (t >= m_emitStartTime - 0.25 && t <= m_emitEndTime)
	{
		// animate emitted particles		
		const Vector3r & x0 = m_x;

		const Real animationMarginAhead = sim->getSupportRadius();
		const Vector3r size = getSize(static_cast<Real>(m_width), static_cast<Real>(m_height), m_type);
		const Real h = size[0];
 		const Real r = 0.5f * size[1];
		const Real r2 = r * r;
		const Vector3r pos = x0 + 0.5f * animationMarginAhead * emitDir;


		const unsigned int nModels = sim->numberOfFluidModels();
		for (unsigned int m = 0; m < nModels; m++)
		{
			FluidModel *fm = sim->getFluidModel(m);
			const unsigned int numParticles = fm->numActiveParticles();
 			#pragma omp parallel for schedule(static) default(shared)
 			for (int i = 0; i < (int)numParticles; i++)
 			{
 				Vector3r &xi = fm->getPosition(i);
 				if (inCylinder(xi, pos, m_rotation, h, r2))
 				{
 					fm->getVelocity(i) = emitVel;
 					fm->getPosition(i) += timeStepSize * emitVel;
 					fm->setParticleState(i, ParticleState::AnimatedByEmitter);
					fm->setObjectId(i, m_objectId);
 				}
 			}
		}
 	}
	if (t < m_nextEmitTime || t > m_emitEndTime)
	{
		return;
	}

	Vector3r axisHeight = m_rotation.col(1);
	Vector3r axisWidth = m_rotation.col(2);

	const Real radius = (static_cast<Real>(0.5) * (Real)m_width * diam);
	const Real radius2 = radius*radius;