下一节:示例 8.6.3. 上一级:数值示例 上一节:示例 8.6.1.

示例 8.6.2.

即使初始向量经过H^{-1} B预乘,Ritz向量也可以通过实施如(8.27)所示的后处理步骤进一步改进。以下示例代表与示例8.6.1相同的轴,但网格更粗糙,使得\{A,B\}的阶数为160。使用了\sigma =1.e5i的移位,并保持了完全正交性。下表显示,后处理步骤将直接残差提高了多达四个数量级。
Ritz值 \Vert A\tilde{x} - \tilde{\lambda}B\tilde{x}\Vert \Vert A\tilde{x} - \tilde{\lambda}B\tilde{x}\Vert
\tilde{\lambda} (有后处理) (无后处理)
-1.96e\!-\!04 + 1.03e\!+\!05i 5.04e\!-\!08 5.04e\!-\!08
-1.30e\!-\!02 + 1.06e\!+\!05i 4.11e\!-\!08 4.11e\!-\!08
-7.16e\!-\!03 + 9.28e\!+\!04i 4.15e\!-\!08 4.18e\!-\!08
-7.60e\!-\!02 + 1.11e\!+\!05i 8.41e\!-\!08 8.51e\!-\!05
-1.01e\!-\!01 + 8.66e\!+\!04i 1.09e\!-\!07 5.06e\!-\!03




下一节:示例 8.6.3. 上一级:数值示例 上一节:示例 8.6.1.
Susan Blackford 2000-11-20