广义逆

在FETI方法中，需要一个稀疏矩阵K的广义逆。
在我的例子中，K是由四面体网格上的有限元离散产生的矩阵

线性弹性方程。我知道如何计算LDL^T的广义逆但我不确定如何有效地计算LDL^T分解。该算法同时使用矩阵的行和列。
我应该使用向左看还是向右看的算法?
如果有人能给我提供任何提示，我将非常感激。
谢谢你！