anki-3d-engine/src/scene/ParticleEmitter.cpp

#include "anki/scene/ParticleEmitter.h"
#include "anki/scene/SceneGraph.h"
#include "anki/scene/InstanceNode.h"
#include "anki/scene/Misc.h"
#include "anki/resource/Model.h"
#include "anki/util/Functions.h"
#include "anki/physics/PhysicsWorld.h"
#include "anki/Gl.h"

namespace anki {

//==============================================================================
// Misc                                                                        =
//==============================================================================

const U COMPONENTS = 3 + 1 + 1; // 3 position, 1 size, 1 alpha
const PtrSize VERT_SIZE = COMPONENTS * sizeof(F32);

//==============================================================================
static F32 getRandom(F32 initial, F32 deviation)
{
	return (deviation == 0.0) 
		? initial
		: initial + randFloat(deviation) * 2.0 - deviation;
}

//==============================================================================
static Vec3 getRandom(const Vec3& initial, const Vec3& deviation)
{
	if(deviation == Vec3(0.0))
	{
		return initial;
	}
	else
	{
		Vec3 out;
		for(U i = 0; i < 3; i++)
		{
			out[i] = getRandom(initial[i], deviation[i]);
		}
		return out;
	}
}

//==============================================================================
// ParticleBase                                                                =
//==============================================================================

//==============================================================================
void ParticleBase::revive(const ParticleEmitter& pe,
	F32 /*prevUpdateTime*/, F32 crntTime)
{
	ANKI_ASSERT(isDead());
	const ParticleEmitterProperties& props = pe;

	// life
	m_timeOfDeath = getRandom(crntTime + props.m_particle.m_life,
		props.m_particle.m_lifeDeviation);
	m_timeOfBirth = crntTime;
}

//==============================================================================
// ParticleSimple                                                              =
//==============================================================================

//==============================================================================
void ParticleSimple::simulate(const ParticleEmitter& pe,
	F32 prevUpdateTime, F32 crntTime)
{
	F32 dt = crntTime - prevUpdateTime;

	ANKI_ASSERT(
		static_cast<const ParticleEmitterProperties&>(pe).
		m_particle.m_gravity.getLength() > 0.0);

	Vec3 xp = m_position;
	Vec3 xc = m_acceleration * (dt * dt) + m_velocity * dt + xp;

	m_position = xc;

	m_velocity += m_acceleration * dt;
}

//==============================================================================
void ParticleSimple::revive(const ParticleEmitter& pe,
	F32 prevUpdateTime, F32 crntTime)
{
	ParticleBase::revive(pe, prevUpdateTime, crntTime);
	m_velocity = Vec3(0.0);

	const ParticleEmitterProperties& props = pe;

	m_acceleration = getRandom(props.m_particle.m_gravity,
			props.m_particle.m_gravityDeviation);

	// Set the initial position
	m_position = getRandom(props.m_particle.m_startingPos,
		props.m_particle.m_startingPosDeviation);

	m_position += pe.getWorldTransform().getOrigin();
}

//==============================================================================
// Particle                                                                    =
//==============================================================================

#if 0

//==============================================================================
Particle::Particle(
	const char* name, SceneGraph* scene, // SceneNode
	// RigidBody
	PhysicsWorld* masterContainer, const RigidBody::Initializer& init_)
	:	ParticleBase(name, scene, PT_PHYSICS)
{
	RigidBody::Initializer init = init_;

	getSceneGraph().getPhysics().newPhysicsObject<RigidBody>(body, init);

	sceneNodeProtected.rigidBodyC = body;
}

//==============================================================================
Particle::~Particle()
{
	getSceneGraph().getPhysics().deletePhysicsObject(body);
}

//==============================================================================
void Particle::revive(const ParticleEmitter& pe,
	F32 prevUpdateTime, F32 crntTime)
{
	ParticleBase::revive(pe, prevUpdateTime, crntTime);

	const ParticleEmitterProperties& props = pe;

	// pre calculate
	Bool forceFlag = props.forceEnabled;
	Bool worldGravFlag = props.wordGravityEnabled;

	// activate it (Bullet stuff)
	body->forceActivationState(ACTIVE_TAG);
	body->activate();
	body->clearForces();
	body->setLinearVelocity(btVector3(0.0, 0.0, 0.0));
	body->setAngularVelocity(btVector3(0.0, 0.0, 0.0));

	// force
	if(forceFlag)
	{
		Vec3 forceDir = getRandom(props.particle.forceDirection,
			props.particle.forceDirectionDeviation);
		forceDir.normalize();

		if(!pe.identityRotation)
		{
			// the forceDir depends on the particle emitter rotation
			forceDir = pe.getWorldTransform().getRotation() * forceDir;
		}

		F32 forceMag = getRandom(props.particle.forceMagnitude,
			props.particle.forceMagnitudeDeviation);

		body->applyCentralForce(toBt(forceDir * forceMag));
	}

	// gravity
	if(!worldGravFlag)
	{
		body->setGravity(toBt(getRandom(props.particle.gravity,
			props.particle.gravityDeviation)));
	}

	// Starting pos. In local space
	Vec3 pos = getRandom(props.particle.startingPos,
		props.particle.startingPosDeviation);

	if(pe.identityRotation)
	{
		pos += pe.getWorldTransform().getOrigin();
	}
	else
	{
		pos.transform(pe.getWorldTransform());
	}

	btTransform trf(
		toBt(Transform(pos, pe.getWorldTransform().getRotation(), 1.0)));
	body->setWorldTransform(trf);
}
#endif

//==============================================================================
// ParticleEmitter                                                             =
//==============================================================================

//==============================================================================
ParticleEmitter::ParticleEmitter(
	const char* name, SceneGraph* scene,
	const char* filename)
	:	SceneNode(name, scene),
		SpatialComponent(this),
		MoveComponent(this),
		RenderComponent(this),
		m_particles(getSceneAllocator()),
		m_transforms(getSceneAllocator())
{
	addComponent(static_cast<MoveComponent*>(this));
	addComponent(static_cast<SpatialComponent*>(this));
	addComponent(static_cast<RenderComponent*>(this));

	m_obb.setCenter(Vec3(0.0));
	m_obb.setExtend(Vec3(1.0));
	m_obb.setRotation(Mat3::getIdentity());

	// Load resource
	m_particleEmitterResource.load(filename);

	// copy the resource to me
	ParticleEmitterProperties& me = *this;
	const ParticleEmitterProperties& other =
		m_particleEmitterResource->getProperties();
	me = other;

	if(m_usePhysicsEngine)
	{
		createParticlesSimulation(scene);
		m_simulationType = SimulationType::PHYSICS_ENGINE;
	}
	else
	{
		createParticlesSimpleSimulation(scene);
		m_simulationType = SimulationType::SIMPLE;
	}

	m_timeLeftForNextEmission = 0.0;
	RenderComponent::init();

	// Create the vertex buffer and object
	//
	PtrSize buffSize = m_maxNumOfParticles * VERT_SIZE * 3;

	GlManager& gl = GlManagerSingleton::get();
	GlJobChainHandle jobs(&gl);

	m_vertBuff = GlBufferHandle(jobs, GL_ARRAY_BUFFER, buffSize, 
		GL_MAP_WRITE_BIT | GL_MAP_PERSISTENT_BIT | GL_MAP_COHERENT_BIT);

	// TODO Optimize that to avoir serialization
	jobs.finish();

	m_vertBuffMapping = (U8*)m_vertBuff.getPersistentMappingAddress();
}

//==============================================================================
ParticleEmitter::~ParticleEmitter()
{
	// Delete simple particles
	if(m_simulationType == SimulationType::SIMPLE)
	{
		for(ParticleBase* part : m_particles)
		{
			getSceneAllocator().deleteInstance(part);
		}
	}
}

//==============================================================================
void ParticleEmitter::buildRendering(RenderingBuildData& data)
{
	ANKI_ASSERT(data.m_subMeshIndicesCount <= m_transforms.size() + 1);

	if(m_aliveParticlesCount == 0)
	{
		return;
	}

	RenderingKey key = data.m_key;
	key.m_lod = 0;

	GlProgramPipelineHandle ppline = 
		m_particleEmitterResource->getMaterial().getProgramPipeline(key);

	ppline.bind(data.m_jobs);

	PtrSize offset = (getGlobTimestamp() % 3) * (m_vertBuff.getSize() / 3);

	// Position
	m_vertBuff.bindVertexBuffer(data.m_jobs, 
		3, GL_FLOAT, false, VERT_SIZE, offset + 0, 0);

	// Scale
	m_vertBuff.bindVertexBuffer(data.m_jobs, 
		1, GL_FLOAT, false, VERT_SIZE, offset + sizeof(F32) * 3, 6);

	// Alpha
	m_vertBuff.bindVertexBuffer(data.m_jobs, 
		1, GL_FLOAT, false, VERT_SIZE, offset + sizeof(F32) * 4, 7);

	GlDrawcallArrays dc = GlDrawcallArrays(
		GL_POINTS, 
		m_aliveParticlesCount,
		data.m_subMeshIndicesCount);

	dc.draw(data.m_jobs);
}

//==============================================================================
const Material& ParticleEmitter::getMaterial()
{
	return m_particleEmitterResource->getMaterial();
}

//==============================================================================
void ParticleEmitter::componentUpdated(SceneComponent& comp, 
	SceneComponent::UpdateType)
{
	if(comp.getType() == MoveComponent::getClassType())
	{
		m_identityRotation =
			getWorldTransform().getRotation() == Mat3::getIdentity();
	}
}

//==============================================================================
void ParticleEmitter::createParticlesSimulation(SceneGraph* scene)
{
#if 0
	collShape = getSceneAllocator().newInstance<btSphereShape>(particle.size);

	RigidBody::Initializer binit;
	binit.shape = collShape;
	binit.group = PhysicsWorld::CG_PARTICLE;
	binit.mask = PhysicsWorld::CG_MAP;

	particles.reserve(maxNumOfParticles);

	for(U i = 0; i < maxNumOfParticles; i++)
	{
		binit.mass = getRandom(particle.mass, particle.massDeviation);

		Particle* part = getSceneGraph().newSceneNode<Particle>(
			nullptr, &scene->getPhysics(), binit);

		part->size = getRandom(particle.size, particle.sizeDeviation);
		part->alpha = getRandom(particle.alpha, particle.alphaDeviation);
		part->getRigidBody()->forceActivationState(DISABLE_SIMULATION);

		particles.push_back(part);
	}
#endif
}

//==============================================================================
void ParticleEmitter::createParticlesSimpleSimulation(SceneGraph* scene)
{
	m_particles.resize(m_maxNumOfParticles);

	for(U i = 0; i < m_maxNumOfParticles; i++)
	{
		ParticleSimple* part = 
			getSceneAllocator().newInstance<ParticleSimple>();

		part->m_size = getRandom(m_particle.m_size, m_particle.m_sizeDeviation);
		part->m_alpha = 
			getRandom(m_particle.m_alpha, m_particle.m_alphaDeviation);

		m_particles[i] = part;
	}
}

//==============================================================================
void ParticleEmitter::frameUpdate(F32 prevUpdateTime, F32 crntTime, 
	SceneNode::UpdateType uptype)
{
	if(uptype == SceneNode::SYNC_UPDATE)
	{
		return;
	}

	// - Deactivate the dead particles
	// - Calc the AABB
	// - Calc the instancing stuff
	//
	Vec3 aabbmin(MAX_F32);
	Vec3 aabbmax(MIN_F32);
	m_aliveParticlesCount = 0;

	F32* verts = (F32*)(m_vertBuffMapping 
		+ (getGlobTimestamp() % 3) * (m_vertBuff.getSize() / 3));
	F32* verts_base = verts;
	(void)verts_base;

	for(ParticleBase* p : m_particles)
	{
		if(p->isDead())
		{
			// if its already dead so dont deactivate it again
			continue;
		}

		if(p->getTimeOfDeath() < crntTime)
		{
			// Just died
			p->kill();
		}
		else
		{
			// It's alive

			// Do checks
			ANKI_ASSERT(((PtrSize)verts + VERT_SIZE 
				- (PtrSize)m_vertBuffMapping) <= m_vertBuff.getSize());

			// This will calculate a new world transformation
			p->simulate(*this, prevUpdateTime, crntTime);

			const Vec3& origin = p->getPosition();

			for(U i = 0; i < 3; i++)
			{
				aabbmin[i] = std::min(aabbmin[i], origin[i]);
				aabbmax[i] = std::max(aabbmax[i], origin[i]);
			}

			F32 lifePercent = (crntTime - p->getTimeOfBirth())
				/ (p->getTimeOfDeath() - p->getTimeOfBirth());

			verts[0] = origin.x();
			verts[1] = origin.y();
			verts[2] = origin.z();

			// XXX set a flag for scale
			verts[3] = p->m_size + (lifePercent * m_particle.m_sizeAnimation);

			// Set alpha
			if(m_particle.m_alphaAnimation)
			{
				verts[4] = sin((lifePercent) * getPi<F32>()) * p->m_alpha;
			}
			else
			{
				verts[4] = p->m_alpha;
			}

			++m_aliveParticlesCount;
			verts += 5;
		}
	}

	if(m_aliveParticlesCount != 0)
	{
		Vec3 min = aabbmin - m_particle.m_size;
		Vec3 max = aabbmax + m_particle.m_size;
		Vec3 center = (min + max) / 2.0;

		m_obb = Obb(center, Mat3::getIdentity(), max - center);
	}
	else
	{
		m_obb = Obb(Vec3(0.0), Mat3::getIdentity(), Vec3(0.001));
	}
	SpatialComponent::markForUpdate();

	//
	// Emit new particles
	//
	if(m_timeLeftForNextEmission <= 0.0)
	{
		U particlesCount = 0; // How many particles I am allowed to emmit
		for(ParticleBase* pp : m_particles)
		{
			ParticleBase& p = *pp;
			if(!p.isDead())
			{
				// its alive so skip it
				continue;
			}

			p.revive(*this, prevUpdateTime, crntTime);

			// do the rest
			++particlesCount;
			if(particlesCount >= m_particlesPerEmittion)
			{
				break;
			}
		} // end for all particles

		m_timeLeftForNextEmission = m_emissionPeriod;
	} // end if can emit
	else
	{
		m_timeLeftForNextEmission -= crntTime - prevUpdateTime;
	}

	// Do something more
	doInstancingCalcs();
}

//==============================================================================
void ParticleEmitter::doInstancingCalcs()
{
	// Gather the move components of the instances
	SceneFrameVector<MoveComponent*> instanceMoves(getSceneFrameAllocator());
	Timestamp instancesTimestamp = 0;
	SceneObject::visitThisAndChildren<SceneNode>([&](SceneNode& sn)
	{	
		if(sn.tryGetComponent<InstanceComponent>())
		{
			MoveComponent& move = sn.getComponent<MoveComponent>();

			instanceMoves.push_back(&move);

			instancesTimestamp = 
				std::max(instancesTimestamp, move.getTimestamp());
		}
	});

	// If instancing
	if(instanceMoves.size() != 0)
	{
		Bool transformsNeedUpdate = false;

		// Check if an instance was added or removed and reset the spatials and
		// the transforms
		if(instanceMoves.size() != m_transforms.size())
		{
			transformsNeedUpdate = true;

			// Check if instances added or removed
			if(m_transforms.size() < instanceMoves.size())
			{
				// Instances added

				U diff = instanceMoves.size() - m_transforms.size();

				while(diff-- != 0)
				{
					ObbSpatialComponent* newSpatial = getSceneAllocator().
						newInstance<ObbSpatialComponent>(this);

					addComponent(newSpatial);
				}
			}
			else
			{
				// Instances removed

				// TODO
				ANKI_ASSERT(0 && "TODO");
			}

			m_transforms.resize(instanceMoves.size());
		}

		if(transformsNeedUpdate || m_transformsTimestamp < instancesTimestamp)
		{
			m_transformsTimestamp = instancesTimestamp;

			// Update the transforms
			for(U i = 0; i < instanceMoves.size(); i++)
			{
				m_transforms[i] = instanceMoves[i]->getWorldTransform();
			}
		}

		// Update the spatials anyway
		U count = 0;
		iterateComponentsOfType<SpatialComponent>([&](SpatialComponent& sp)
		{
			// Skip the first
			if(count != 0)	
			{
				ObbSpatialComponent* msp = 
					staticCastPtr<ObbSpatialComponent*>(&sp);
		
				Obb aobb = m_obb;
				aobb.setCenter(Vec3(0.0));
				msp->obb = aobb.getTransformed(m_transforms[count - 1]);

				msp->markForUpdate();
			}

			++count;
		});

		ANKI_ASSERT(count - 1 == m_transforms.size());
	} // end if instancing
}

//==============================================================================
Bool ParticleEmitter::getHasWorldTransforms()
{
	if(m_transforms.size() == 0)
	{
		return false;
	}
	else
	{
		return true;
	}
}

//==============================================================================
void ParticleEmitter::getRenderWorldTransform(U index, Transform& trf)
{
	ANKI_ASSERT(m_transforms.size() > 0);

	if(index == 0)
	{
		// Don't transform the particle positions. They are already in world 
		// space
		trf = Transform::getIdentity();
	}
	else
	{
		--index;
		ANKI_ASSERT(index < m_transforms.size());

		// The particle positions are already in word space. Move them back to
		// local space
		Transform invTrf = getWorldTransform().getInverse();
		trf = Transform::combineTransformations(m_transforms[index], invTrf);
	}
}

} // end namespace anki