当前位置: 首页>>代码示例>>C++>>正文


C++ TriMesh::area方法代码示例

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


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

示例1: area

Field<T,VertexHandle> mean_curvatures(const TriMesh& mesh) {
  Field<T,VertexHandle> area(mesh.n_vertices()); // Actually 8*area
  Field<TV,VertexHandle> normal(mesh.n_vertices());
  Field<TV,VertexHandle> Hn(mesh.n_vertices()); // Actually 4*Hn
  for (const auto f : mesh.face_handles()) {
    const auto v = mesh.vertex_handles(f);
    const TV x0 = mesh.point(v.x),
             x1 = mesh.point(v.y),
             x2 = mesh.point(v.z),
             x01 = x1-x0,
             x12 = x2-x1,
             x20 = x0-x2;
    const T cot0 = cot(x01,x20),
            cot1 = cot(x01,x12),
            cot2 = cot(x12,x20);
    // Compute A_mixed as in Meyer et al.
    const TV n = mesh.calc_face_normal(f);
    if (cot0<=0 && cot1<=0 && cot2<=0) { // Voronoi case
      const T area0 = cot0*sqr_magnitude(x12),
              area1 = cot1*sqr_magnitude(x20),
              area2 = cot2*sqr_magnitude(x01);
      area[v.x] -= (area1+area2);
      area[v.y] -= (area2+area0);
      area[v.z] -= (area0+area1);
      normal[v.x] -= (area1+area2)*n;
      normal[v.y] -= (area2+area0)*n;
      normal[v.z] -= (area0+area1)*n;
    } else { // One of the triangles is obtuse
      const T a = 2*mesh.area(f);
      area[v.x] += (1+(cot0>0))*a;
      area[v.y] += (1+(cot1>0))*a;
      area[v.z] += (1+(cot2>0))*a;
      normal[v.x] += (1+(cot0>0))*a*n;
      normal[v.y] += (1+(cot1>0))*a*n;
      normal[v.z] += (1+(cot2>0))*a*n;
    }
    // Accumulate mean curvature normal
    Hn[v.x] += cot2*x01-cot1*x20;
    Hn[v.y] += cot0*x12-cot2*x01;
    Hn[v.z] += cot1*x20-cot0*x12;
  }
  Field<T,VertexHandle> H(mesh.n_vertices(),uninit);
  for (const auto v : mesh.vertex_handles()) {
    const TV n = normal[v];
    H[v] = 2/area[v] * (mesh.is_boundary(v) ? dot(Hn[v],normalized(n))
                                            : copysign(magnitude(Hn[v]),dot(Hn[v],n)));
  }
  return H;
}
开发者ID:Haider-BA,项目名称:geode,代码行数:49,代码来源:curvature.cpp


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