如何通过Eigen的CholmodSupport实现稀疏矩阵Cholesky秩1更新？

Eigen自带的CholmodSupport模块并未直接封装cholmod_updown函数，该模块目前仅对外暴露稀疏Cholesky分解、线性求解、行列式计算等基础能力，CHOLMOD原生支持的秩1/秩k更新、因子拼接等高级接口没有做上层API映射。
但你仍然可以借助该模块调用cholmod_updown，不需要重复做矩阵格式转换：

• 使用Eigen::CholmodSupernodalLLT或Eigen::CholmodSimplicialLLT完成初始分解后，类内部会维护CHOLMOD原生的cholmod_factor*类型因子对象，你可以通过继承分解类的方式访问保护成员m_cholmodFactor，将这个原生指针直接传入cholmod_updown即可完成秩更新操作。
• 调用前需要保证项目引入的CHOLMOD头文件、编译配置和Eigen编译时链接的CHOLMOD版本完全一致，尤其要匹配整数长度、浮点精度相关的宏定义，否则会出现内存越界、计算结果错误的问题。更新用的向量可以通过Eigen的Map接口直接映射为CHOLMOD要求的cholmod_dense/cholmod_sparse结构，不需要额外拷贝数据。

• 直接调用CHOLMOD原生接口：如果项目已经链接了CHOLMOD库，可以绕开Eigen的封装层，直接用CHOLMOD的C接口完成稀疏矩阵读入、Cholesky分解、因子拼接、秩更新全流程，不需要跨API访问内部成员，逻辑更可控，也没有额外性能损耗。
• 使用CHOLMOD官方C封装：新版SuiteSparse提供了CHOLMOD的面向对象C接口，原生封装了秩更新、因子修改方法，不需要手动操作C风格结构体，能大幅降低内存管理、格式匹配的出错概率。
• 自定义稀疏秩1更新逻辑：如果求解序列中的矩阵在拼接、秩更新后稀疏结构保持固定，可以基于Simplicial形式的稀疏Cholesky分解（即L为逐列存储的稀疏下三角矩阵）自行实现秩1更新，核心逻辑是对L的对应列做Givens旋转修正，代码量很小，且可以完全适配自定义的稀疏存储格式，没有额外第三方依赖。
• 避开稠密版本接口的坑：不要使用Eigen自带的稠密版cholupdate，该接口仅针对稠密矩阵设计，传入稀疏矩阵会触发全量稠密转换，对于大规模稀疏问题会产生远超实际需求的内存和时间开销，完全不适用。

内容的提问来源于stack exchange，提问作者Etienne Côme