Method initModel

SPlisHSPlasH/FluidModel.cpp:285–327

}

void FluidModel::initModel(const std::string &id, const unsigned int nFluidParticles, Vector3r* fluidParticles, Vector3r* fluidVelocities, unsigned int* fluidObjectIds, const unsigned int nMaxEmitterParticles)
{
	m_id = id;
	init();
	releaseFluidParticles();
	resizeFluidParticles(nFluidParticles + nMaxEmitterParticles);

	// copy fluid positions
	#pragma omp parallel default(shared)
	{
		#pragma omp for schedule(static)  
		for (int i = 0; i < (int)nFluidParticles; i++)
		{
			getPosition0(i) = fluidParticles[i];
			getPosition(i) = fluidParticles[i];
			getVelocity0(i) = fluidVelocities[i];
			getVelocity(i) = fluidVelocities[i];
			getAcceleration(i).setZero();
			m_density[i] = 0.0;
			m_particleId[i] = i;
			m_objectId[i] = fluidObjectIds[i];
			m_objectId0[i] = fluidObjectIds[i];
			if (m_particleState[i] != ParticleState::Fixed)
				m_particleState[i] = ParticleState::Active;
		}
	}
	// set IDs for emitted particles
	for (unsigned int i = nFluidParticles; i < (nFluidParticles + nMaxEmitterParticles); i++)
	{
		m_particleId[i] = i;
		m_objectId[i] = 0;
	}

	// initialize masses
	initMasses();

	// Fluids 
	NeighborhoodSearch *neighborhoodSearch = Simulation::getCurrent()->getNeighborhoodSearch();
	m_pointSetIndex = neighborhoodSearch->add_point_set(&getPosition(0)[0], nFluidParticles, true, true, true, this);

	m_numActiveParticles0 = nFluidParticles;
	m_numActiveParticles = m_numActiveParticles0;
}


void FluidModel::performNeighborhoodSearchSort()