矩阵逆的模(dropping strategy)

2011-10-17

下面的文本主要用于算法设计和算法实现的需要。

主要内容参考了Golub的书籍。

适合任何三角矩阵（无论对角元是否为1，属于通用算法）

U (1)

1 估计1-

1.1 方法1： (1)

1.2 方法2： (2)

2 估计1-L (2)

2.1 方法1：estfullnrm_forward(T,D) (2)

2.2 方法2： (3)

1估计1-U

见Golub的书：

1.1 方法1：

1.2 方法2：

2 估计1

-L 参考Golub 的书，在估计1-U 的模算法基础上得到的对于下三角矩阵的算法：

2.1 方法1：

0):1(=n p

for n k :1=

Choose )(k d

),(/))()(()(k k T k p k d k y -=

),:1()():1():1(k n k T k y n k p n k p +++=+

end

2.2 方法2：

0):1(=n p

for n k :1=

),(/))(1()(k k T k p k y -=+

),(/))(1()(k k T k p k y --=-

),:1()():1()(k n k T k y n k p k p +++=++ ),:1()():1()(k n k T k y n k p k p +++=-- if 11||)(|||)(|||)(|||)(|--+++≥+k p k y k p k y

+=)()(k y k y

+=+)():1(k p n k p else

-=)()(k y k y

-=+)():1(k p n k p end

end