在 PETSc 共轭梯度求解器中如何实现零空间校正？我搜索了源代码和文档，但找不到有关实际实现的参考资料。

作为背景，我正在解决由双线性形式产生的线性系统，该双线性形式在商空间中是椭圆的$\mathbb{R}^n\big/\ker{\mathbf{Z}}$， 在哪里$\mathbf{Z}$是一个$m \times n$具有非平凡零空间的矩阵。

我计算 SVD$\mathbf{Z}$并将对应于“大”奇异值的右奇异向量作为零空间向量提供给 PETSc 共轭梯度求解器。它就像一个魅力，除了当$\mathbf{Z}$不是学术规模。