2

我正在尝试为网格编写一个小解包器。这使用有限元方法来求解展平表面和原始表面之间的最小线性应力。目前有一些顶点被固定以获得结果。没有这个,三角形会随机旋转和平移......

但由于这个问题不需要固定,更好的解决方案是直接求解奇异矩阵。Petsc 通过提供有关零空间的一些信息,确实提供了一些解决奇异系统的方法。http://www.mcs.anl.gov/petsc/petsc-current/docs/manual.pdf#section.4.6我想知道在 Eigen 中是否有任何替代方案。如果没有,是否有任何其他可能性可以在不修复/固定顶点的情况下解决此问题。

谢谢,问候

另请参阅此链接以获取更多信息: 开发历史

4

1 回答 1

1

Eigen 提供了一种 SVD 分解算法:Jacobi SVD

SVD 分解给出零空间。按照维基百科文章的符号,让M = U D V是 M 的 SVD 分解,其中D是奇异值的对角矩阵。然后,从Range、null 空间和 rank

对应于 M 的消失奇异值的右奇异向量 [V] 跨越 M 的零空间

于 2016-11-15T07:56:20.120 回答