[论文解读] Cohomologous Harmonic Cochains

Anil N. Hirani, Kaushik Kalyanaraman+2 more|arXiv (Cornell University)|Dec 13, 2010

Topological and Geometric Data Analysis参考文献 41被引用 4

一句话总结

本文提出了一种加权最小二乘法，用于在单纯复形网格上计算与给定上循环同调的调和共形，该方法利用离散Hodge-deRham同构。该方法通过生成更稀疏的线性系统，避免了先前方法的数值缺陷，显著提高了求解效率——尤其在使用Whitney Hodge算子时效果更显著，同时保持对Hodge算子选择的独立性，并适用于有限元和离散外微分计算框架。

ABSTRACT

We describe algorithms for finding harmonic cochains, an essential ingredient for solving elliptic partial differential equations in exterior calculus. Harmonic cochains are also useful in computational topology and computer graphics. We focus on finding harmonic cochains cohomologous to a given cocycle. Amongst other things, this allows localization near topological features of interest. We derive a weighted least squares method by proving a discrete Hodge-deRham theorem on the isomorphism between the space of harmonic cochains and cohomology. The solution obtained either satisfies the Whitney form finite element exterior calculus equations or the discrete exterior calculus equations for harmonic cochains, depending on the discrete Hodge star used.

研究动机与目标

• 开发一种在单纯复形网格上计算与给定上循环同调的调和共形的数值鲁棒且高效的方法。
• 克服现有基于Poisson的方法的数值低效问题，特别是使用Whitney Hodge算子时的问题。
• 提供一种无需计算完整调和基或同调基的方法，从而实现在拓扑特征附近的局部化求解。
• 建立离散Hodge-deRham同构，以证明加权最小二乘公式的合理性。
• 提升计算拓扑和椭圆PDE有限元方法中的可扩展性和性能。

提出的方法

• 该方法基于离散Hodge-deRham同构推导出加权最小二乘公式，将上同调类与调和共形联系起来。
• 求解一个等价于Demlow和Hirani中方程的混合有限元外微分计算系统（引理3.10），使用Hodge算子（Whitney或对偶原生型）。
• 线性系统被构造为在满足调和条件的前提下，最小化共形对给定上循环的偏差的加权范数。
• 该方法避免了解特征向量问题或计算完整调和基的需求，从而支持局部化计算。
• 该方法与Hodge算子的选择无关，而不同于先前方法在使用Whitney算子时出现病态或稠密系统的问题。
• 由于消除了多余项，所得到的线性系统相比Desbrun等人和Gu-Yau的方法更稀疏，尤其在使用Whitney Hodge算子时效果显著。

实验结果

研究问题

• RQ1 能否开发一种数值高效的算法，用于在不需完整调和基的前提下计算与给定上循环同调的调和共形？
• RQ2 线性系统的稀疏性和条件数如何影响求解时间，特别是在使用Whitney Hodge算子时？
• RQ3 能否直接从离散Hodge-deRham同构推导出最小二乘公式，以确保正确性和稳定性？
• RQ4 与Desbrun等人或Gu-Yau的现有Poisson求解器相比，该最小二乘法在性能上的定量差异是什么？
• RQ5 该方法在2D和3D单纯复形网格中，对Whitney与对偶原生型Hodge算子的选择是否均保持鲁棒性和高效性？

主要发现

• 最小二乘法生成的线性系统相比Desbrun等人方法显著更稀疏，尤其在使用Whitney Hodge算子时：在实心环形四面体网格中，非零元数量分别为56,000和1430万。
• 最小二乘法的求解时间分别为0.1181秒（DEC Hodge算子）和0.1355秒（Whitney Hodge算子），而Desbrun等人使用SuperLU求解的时间为13.05秒和1746秒。
• 该方法避免了Desbrun等人系统中的额外项，该额外项导致数值不稳定性和可扩展性差，尤其在使用Whitney Hodge算子时更为明显。
• 最小二乘公式与Hodge算子的选择无关，无论使用何种离散Hodge算子，均保持稳定和高效。
• 该方法无需特征向量计算或后处理，性能优于基于特征向量的替代方法。
• 以简化形式证明了离散Hodge-deRham同构，为最小二乘方法提供了理论依据，并实现了对同调调和共形问题的直接求解。