Method updateVisMesh

Simulation/TetModel.cpp:286–324

}
 
void TetModel::updateVisMesh(const ParticleData &pd)
{
	if (m_attachments.size() == 0)
		return;

	// The collision mesh is the boundary of the tet mesh
	unsigned int *faces = m_surfaceMesh.getFaces().data();
	const unsigned int nFaces = m_surfaceMesh.numFaces();

	const Vector3r *normals = m_surfaceMesh.getVertexNormals().data();

	#pragma omp parallel default(shared)
	{
		#pragma omp for schedule(static)  
		for (int i = 0; i < (int) m_attachments.size(); i++)
		{
			const unsigned int pindex = m_attachments[i].m_index;
			const unsigned int triindex = m_attachments[i].m_triIndex;
			const Real *bary = m_attachments[i].m_bary;

			const unsigned int indexA = faces[3 * triindex] + m_indexOffset;
			const unsigned int indexB = faces[3 * triindex + 1] + m_indexOffset;
			const unsigned int indexC = faces[3 * triindex + 2] + m_indexOffset;

			const Vector3r &a = pd.getPosition(indexA);
			const Vector3r &b = pd.getPosition(indexB);
			const Vector3r &c = pd.getPosition(indexC);
			Vector3r p2 = bary[0] * a + bary[1] * b + bary[2] * c;
			Vector3r n = bary[0] * normals[faces[3 * triindex]] + bary[1] * normals[faces[3 * triindex + 1]] + bary[2] * normals[faces[3 * triindex + 2]];
			n.normalize();

			Vector3r &p = m_visVertices.getPosition(pindex);
			p = p2 - n*m_attachments[i].m_dist;
		}
	}

	m_visMesh.updateNormals(m_visVertices, 0);
	m_visMesh.updateVertexNormals(m_visVertices);
}
 
 
 bool TetModel::pointInTriangle(const Vector3r& p0, const Vector3r& p1, const Vector3r& p2, const Vector3r& p,