Function precomputeValues

Tools/MeshSkinning/main.cpp:841–921

void precomputeValues()
{
	const int numParticles = (int)activeParticles.size();
	unsigned int numVertices = mesh.numVertices();
	auto& meshX0 = mesh.getVertices();

#ifdef USE_AVX
	precomputed_indices8.clear();
	precomp_V_gradW8.clear();
	precomputed_indices8.resize(numParticles);
#else
	precomputed_indices.clear();
	precomp_V_gradW.clear();
	precomputed_indices.resize(numParticles);
#endif


	unsigned int sumNeighborParticles = 0;
	unsigned int sumNeighborParticles8 = 0;
	for (int i = 0; i < numParticles; i++)
	{
		const unsigned int index = activeParticles[i];
		const size_t numNeighbors = initialNeighbors[index].size();

#ifdef USE_AVX
		precomputed_indices8[i] = sumNeighborParticles8;

		// steps of 8 values due to avx
		sumNeighborParticles8 += (unsigned int) numNeighbors / 8u;
		if (numNeighbors % 8 != 0)
			sumNeighborParticles8++;
#else
		precomputed_indices[i] = sumNeighborParticles;
		sumNeighborParticles += (int) numNeighbors;
#endif
	}

#ifdef USE_AVX
	precomp_V_gradW8.resize(sumNeighborParticles8);
#else
	precomp_V_gradW.resize(sumNeighborParticles);
#endif

	#pragma omp parallel default(shared)
	{
		#pragma omp for schedule(static) 
		for (int i = 0; i < (int)numParticles; i++)
		{
			const unsigned int index = activeParticles[i];
			const Vector3r& xi0 = x0[index];
			const unsigned int numNeighbors = (unsigned int)initialNeighbors[index].size();

#ifdef USE_AVX
			const Vector3f8 xi0_avx(xi0);
			unsigned int base8 = precomputed_indices8[i];
			unsigned int idx = 0;
			Matrix3f8 L_avx(L[index]);
			for (unsigned int j = 0; j < numNeighbors; j += 8)