C++ TriangleMesh::nodes方法代码示例

本文整理汇总了C++中TriangleMesh::nodes方法的典型用法代码示例。如果您正苦于以下问题：C++ TriangleMesh::nodes方法的具体用法？C++ TriangleMesh::nodes怎么用？C++ TriangleMesh::nodes使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类TriangleMesh的用法示例。

在下文中一共展示了TriangleMesh::nodes方法的2个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。

示例1: TriangleMesh

TriangleMesh * CorrespondTemplates::computeSurface(char * exportFilename, bool useMapping)
{
	int i, j, counter, numTetNodes;
	TriangleMesh * resultMesh = new TriangleMesh();
	vector<int> surfaceIds;

	//Find the surface:
	findSurfaceNodes();

	numTetNodes = volumetricMeshNodes.numberOfItems();
	surfaceIds.resize(numTetNodes);

	g_NodeContainer nodesSurf;
	if ( useMapping ) 
	   nodesSurf = rigidlyAlignedSurfaceTemplate->nodes();

	//Add the nodes:
	counter = 0;
	for(i = 0; i < numTetNodes; i++)
	{
		if(surfaceNodes[i] == -1) 
			surfaceIds[i] = -1;
		else 
		{
			surfaceIds[i] = counter;
			counter++;

			g_Vector coord;
			if ( useMapping ) {
			  coord =  mappingWeightsTet[3*i] * nodesSurf[mappingIndicesTet[3*i]]->coordinate(); 
			  coord += mappingWeightsTet[3*i+1] * nodesSurf[mappingIndicesTet[3*i+1]]->coordinate(); 
			  coord += mappingWeightsTet[3*i+2] * nodesSurf[mappingIndicesTet[3*i+2]]->coordinate();
			} else {
			  coord = volumetricMeshNodes[i]->coordinate();
			}
			g_Node * newNode = new g_Node(coord); 
			resultMesh->node(newNode);
		}
	}

	//Add the triangles (find the correct orientation):
	for(i = 0; i < trianglesOnSurface.size(); i++)
	{
		g_Element * elem = new g_Element();
	
		for(j = 0; j < allTetrahedra.size(); j++)
		{
			set<int> triangle; triangle.insert(allTetrahedra[j][0]); triangle.insert(allTetrahedra[j][1]); triangle.insert(allTetrahedra[j][2]);
			if(trianglesOnSurface[i] == triangle)
			{
				elem->node(resultMesh->nodes()[surfaceIds[allTetrahedra[j][0]]]);
				elem->node(resultMesh->nodes()[surfaceIds[allTetrahedra[j][2]]]);
				elem->node(resultMesh->nodes()[surfaceIds[allTetrahedra[j][1]]]);
				break;
			}
			triangle.clear(); triangle.insert(allTetrahedra[j][0]); triangle.insert(allTetrahedra[j][1]); triangle.insert(allTetrahedra[j][3]);
			if(trianglesOnSurface[i] == triangle)
			{
				elem->node(resultMesh->nodes()[surfaceIds[allTetrahedra[j][0]]]);
				elem->node(resultMesh->nodes()[surfaceIds[allTetrahedra[j][1]]]);
				elem->node(resultMesh->nodes()[surfaceIds[allTetrahedra[j][3]]]);
				break;
			}
			triangle.clear(); triangle.insert(allTetrahedra[j][0]); triangle.insert(allTetrahedra[j][2]); triangle.insert(allTetrahedra[j][3]);
			if(trianglesOnSurface[i] == triangle)
			{
				elem->node(resultMesh->nodes()[surfaceIds[allTetrahedra[j][0]]]);
				elem->node(resultMesh->nodes()[surfaceIds[allTetrahedra[j][3]]]);
				elem->node(resultMesh->nodes()[surfaceIds[allTetrahedra[j][2]]]);
				break;
			}
			triangle.clear(); triangle.insert(allTetrahedra[j][1]); triangle.insert(allTetrahedra[j][2]); triangle.insert(allTetrahedra[j][3]);
			if(trianglesOnSurface[i] == triangle)
			{
				elem->node(resultMesh->nodes()[surfaceIds[allTetrahedra[j][1]]]);
				elem->node(resultMesh->nodes()[surfaceIds[allTetrahedra[j][2]]]);
				elem->node(resultMesh->nodes()[surfaceIds[allTetrahedra[j][3]]]);
				break;
			}
		}

		resultMesh->element(elem);
	}

	//Export:
	if(exportFilename != NULL)
	{
		exportMeshWrapper(exportFilename, resultMesh);
	}

	return resultMesh;
}

开发者ID:slyandys，项目名称:test_triangleMesh，代码行数:92，代码来源:CorrespondSurfaceAndVolumetricTemplate.cpp

示例2: computeMappings

void CorrespondTemplates::computeMappings(char * outfileMappingTet, char * outfileMappingSurf, g_Node * contactPoint)
{
	//Compute the mappings using barycentric coordinates
	findSurfaceNodes();
	maptet2surfID = surfaceNodes;
	// output tetmesh surface
	computeSurface( "tet_mesh_surf.wrl" );
	deformVolumetricMeshRigidly();

	int i, j, k, numNodesTet, numElemsTet, numNodesSurf, numElemsSurf;
	double dist, minDist, d, mind, distContactPoint, minDistContactPoint;
	numNodesTet = volumetricMeshNodes.numberOfItems();
	numElemsTet = trianglesOnSurface.size();
	numNodesSurf = rigidlyAlignedSurfaceTemplate->getNumberOfNodes();
	numElemsSurf = rigidlyAlignedSurfaceTemplate->getNumberOfTriangles();

	//First compute the mapping from the tetrahedral mesh to the surface mesh
	mappingIndicesTet.resize(3*numNodesTet);
	mappingWeightsTet.resize(3*numNodesTet);
	offsetsTet.resize(numNodesTet);

	for(i = 0; i < numNodesTet; i++)
	{
		if(surfaceNodes[i] == -1)
		{
			mappingIndicesTet[3*i] = mappingIndicesTet[3*i+1] = mappingIndicesTet[3*i+2] = -1;
			mappingWeightsTet[3*i] = mappingWeightsTet[3*i+1] = mappingWeightsTet[3*i+2] = -1;
		}
		else
		{
			for(j = 0; j < numElemsSurf; j++)
			{
			  // g_Vector closestPoint = computeClosestPoint(volumetricMeshNodes[i], rigidlyAlignedSurfaceTemplate->elements()[j]);
			  // dist = closestPoint.DistanceTo(volumetricMeshNodes[i]->coordinate());
		          g_Vector closestPoint;
			  dist = computeClosestPoint(*volumetricMeshNodes[i], *rigidlyAlignedSurfaceTemplate->elements()[j], closestPoint );
			  dist = sqrt(dist);

				if((j == 0) || (dist < minDist))
				{
					minDist = dist;
					
					mappingIndicesTet[3*i] = rigidlyAlignedSurfaceTemplate->elements()[j]->nodes()[0]->id()-1;
					mappingIndicesTet[3*i+1] = rigidlyAlignedSurfaceTemplate->elements()[j]->nodes()[1]->id()-1;
					mappingIndicesTet[3*i+2] = rigidlyAlignedSurfaceTemplate->elements()[j]->nodes()[2]->id()-1;

					//----------------------
					for(int n = 0; n < 3; n++){
						d = volumetricMeshNodes[i]->coordinate().DistanceTo(rigidlyAlignedSurfaceTemplate->elements()[j]->nodes()[n]->coordinate());
						if((n == 0) || (d < mind)){
							mind = d;
							maptet2surfID[i] = rigidlyAlignedSurfaceTemplate->elements()[j]->nodes()[n]->id()-1;
						}
					}
					//------------------------

					vector<double> barycentrics = computeBarycentrics(volumetricMeshNodes[i]->coordinate(), closestPoint, rigidlyAlignedSurfaceTemplate->elements()[j]);
					mappingWeightsTet[3*i] = barycentrics[0];
					mappingWeightsTet[3*i+1] = barycentrics[1];
					mappingWeightsTet[3*i+2] = barycentrics[2];
					offsetsTet[i] = barycentrics[3];
				}
			}
		}
	}

	//Export
	if(outfileMappingTet != NULL)
	{
#if 0
		FILE * fp = fopen(outfileMappingTet, "w");
		if(fp == NULL)
		{
			cout<<"Problem writing "<<outfileMappingTet<<endl;
			return;
		}
		for(i = 0; i < numNodesTet; i++)
			fprintf(fp, "%d %d %d %f %f %f %f\n", mappingIndicesTet[3*i], mappingIndicesTet[3*i+1], mappingIndicesTet[3*i+2], 
				mappingWeightsTet[3*i], mappingWeightsTet[3*i+1], mappingWeightsTet[3*i+2], offsetsTet[i]);
		fclose(fp);
#else
		// Make a surface mesh using the tetmesh surface plus the coordinates from the surface mesh
		computeSurface( outfileMappingTet, true );
#endif
	}

	//Second compute the mapping from the surface mesh to the tetrahedral mesh
	mappingIndicesSurf.resize(3*numNodesSurf);
	mappingWeightsSurf.resize(3*numNodesSurf);
	offsetsSurf.resize(numNodesSurf);
	contactTriangleId.resize(3);
	contactWeights.resize(3);
	contactoffset.resize(1);

	for(i = 0; i < numNodesSurf; i++)
	{
		for(j = 0; j < numElemsTet; j++)
		{
			g_Element elem;
			set<int>::iterator triangleIt;
//.........这里部分代码省略.........

开发者ID:slyandys，项目名称:test_triangleMesh，代码行数:101，代码来源:CorrespondSurfaceAndVolumetricTemplate.cpp

注：本文中的TriangleMesh::nodes方法示例由纯净天空整理自Github/MSDocs等开源代码及文档管理平台，相关代码片段筛选自各路编程大神贡献的开源项目，源码版权归原作者所有，传播和使用请参考对应项目的License；未经允许，请勿转载。