Function init

Tools/MeshSkinning/main.cpp:276–432

void init()
{
	diameter = static_cast<Real>(2.0) * radius;
	mass = static_cast<Real>(1.0) * diameter * diameter * diameter * density0;

	restVolumes.resize(x0.size());
	shepard.resize(mesh.numVertices());
	densities.resize(x0.size());
	L.resize(x0.size());
	F.resize(x0.size());

	supportRadius = supportRadiusFactor * radius;

	CubicKernel::setRadius(supportRadius);
#ifdef USE_AVX
	CubicKernel_AVX::setRadius(supportRadius);
#endif
	W_zero = CubicKernel::W_zero();
	kernelFct = CubicKernel::W;
	gradKernelFct = CubicKernel::gradW;

	// Init neighborhood search
	neighborhoodSearch = new NeighborhoodSearch(supportRadius, false);
	neighborhoodSearch->set_radius(supportRadius);
	neighborhoodSearch->add_point_set(&x0[0][0], x0.size(), true, true);

	// find initial neighbors
	START_TIMING("neighborhoodSearch");
	neighborhoodSearch->find_neighbors();
	STOP_TIMING_AVG;

	// store initial neighbors and init rest volumes
	int numParticles = (int) x0.size();
	initialNeighbors.resize(numParticles);
	#pragma omp parallel default(shared)
	{
		#pragma omp for schedule(static) 
		for (int i = 0; i < numParticles; i++)
		{
			// only neighbors in same phase will influence elasticity
			const unsigned int numNeighbors = numberOfNeighbors(0, i);
			initialNeighbors[i].resize(numNeighbors);
			for (unsigned int j = 0; j < numNeighbors; j++)
				initialNeighbors[i][j] = getNeighbor(0, i, j);

			std::sort(initialNeighbors[i].begin(), initialNeighbors[i].end(), Comparator(x0[i], &x0));

			if (initialNeighbors[i].size() > maxNeighbors)
				initialNeighbors[i].resize(maxNeighbors);

			// compute volume
			Real d = W_zero;
			const Vector3r& xi = x0[i];
			for (unsigned int j = 0; j < numberOfNeighbors(0, i); j++)
			{
				const unsigned int neighborIndex = getNeighbor(0, i, j);
				const Vector3r& xj = x0[neighborIndex];
				d += kernelFct(xi - xj);