Function deformMesh

Tools/MeshSkinning/main.cpp:728–764

void deformMesh()
{
	const int numVertices = (int)mesh.numVertices();
	auto& meshX0 = mesh.getVertices();

	#pragma omp parallel default(shared)
	{
		#pragma omp for schedule(static)  
		for (int i = 0; i < numVertices; i++)
		{
			Vector3r &mesh_xi = meshX[i];
			const Vector3r& mesh_xi0 = meshX0[i];
			const Vector3f8 xi0_avx(mesh_xi0);

			const unsigned int numNeighbors = (unsigned int)initialMeshNeighbors[i].size();

			Vector3f8 xNew;
			xNew.setZero();
			for (unsigned int j = 0; j < numNeighbors; j += 8)
			{
				const unsigned int count = std::min(numNeighbors - j, 8u);

				const Vector3f8 xj_avx = convertVec_zero(&initialMeshNeighbors[i][j], &x[0], count);
				const Vector3f8 xj0_avx = convertVec_zero(&initialMeshNeighbors[i][j], &x0[0], count);
				const Matrix3f8 Fj_avx = convertMat_zero(&initialMeshNeighbors[i][j], &F[0], count);

				const Vector3f8 xi_xj_0 = xi0_avx - xj0_avx;
				const Scalarf8 restVolume_j_avx = convert_zero(&initialMeshNeighbors[i][j], &restVolumes[0], count);
				const Scalarf8 W = CubicKernel_AVX::W(xi_xj_0);

				xNew += (Fj_avx * xi_xj_0 + xj_avx) * restVolume_j_avx * W;
			}
			mesh_xi = xNew.reduce();
			mesh_xi *= shepard[i];
		}
	}
}

#else