对于已经计算出的L的数值稀疏结构,并且supernode也已经确定
计算前next、Lpos数组未定义,head数组均为-1,三者均在分解计算中更新
其中包含两层循环,外层为超节点遍历,内层为对该超节点有贡献的所有descendant的遍历
for s=0:nsuper-1 使用left-looking依次更新supernode k
dnext=head[s]
if dnext不为-1(即存在有贡献子孙)
d=dnext
使用Lpos定位子孙 d 相应Ls与Lx的offset
找到贡献矩阵C,包含C1与C2
C1至少包含1行(因为在对称正定问题中,对主元有贡献的子孙才加进next数组中)
C1 rank