Method step

SPlisHSPlasH/AnimationField.cpp:42–143

}

void AnimationField::step()
{
	Simulation *sim = Simulation::getCurrent();
	TimeManager *tm = TimeManager::getCurrent();
	const Real t = tm->getTime();
	const Real dt = tm->getTimeStepSize();

	if (t >= m_startTime && t <= m_endTime)
	{
		// animate particles		
		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();

			// find animated field
			const FieldDescription *particleField = nullptr;
			for (unsigned int j = 0; j < fm->numberOfFields(); j++)
			{
				const FieldDescription &field = fm->getField(j);
				if (field.name == m_particleFieldName)
				{
					particleField = &field;
					break;
				}
			}

			if (particleField == nullptr)
				continue;

			#pragma omp parallel for schedule(static) default(shared)
			for (int i = 0; i < (int)numParticles; i++)
			{
				const Vector3r &xi = fm->getPosition(i);
				const Vector3r &vi = fm->getVelocity(i);

				const Eigen::Vector3d xi_double = xi.cast<double>();
				const Eigen::Vector3d vi_double = vi.cast<double>();

				if (inShape(m_type, xi, m_x, m_rotation, m_scale))
				{
					Eigen::Map<Vector3r> value((Real*) particleField->getFct(i));
					const Eigen::Vector3d value_double = Vector3r(value).cast<double>();

					const double t_double = static_cast<double>(t);
					const double dt_double = static_cast<double>(dt);

 					te_variable vars[] = { {"t", &t_double}, {"dt", &dt_double}, 
										   {"x", &xi_double[0]}, {"y", &xi_double[1]}, {"z", &xi_double[2]},
										   {"vx", &vi_double[0]}, {"vy", &vi_double[1]}, {"vz", &vi_double[2]},
										   {"valuex", &value_double[0]}, {"valuey", &value_double[1]}, {"valuez", &value_double[2]},
										 };
 					const int numVars = 11;
 					int err;

					//////////////////////////////////////////////////////////////////////////
					// v_x