下一节:算法
上一级:收缩与重启
上一节:收缩
假设广义舒尔形式(参见8.12)按照\tau 进行排序,使得
\vert T^A_{1,1}/T^B_{1,1}-\tau\vert\leq\vert T^A_{2,2}/T^B_{2,2} 0 \tau\vert\leq\cdots \leq\vert T^A_{{m},{m}}/T^B_{{m},{m}}-\tau\vert,
其中{m}是\mathrm{span}(V)的维度。那么,对于i<{m},由 V S^R的前i列张成的空间
\mathrm{span} ( Vs^R_1, \ldots, Vs^R_i )
包含了最有可能的i个Petrov向量。对应的测试子空间由
\mathrm{span}(Ws^L,\ldots,Ws^L_i)给出。
因此,为了将子空间的维度(“隐式重启”)减少到{m}_{\min},{m}_{\min}<{m},可以简单地丢弃列
v_{{m}_{\min}+1}
到v_{m}以及
w_{{m}_{\min}+1}
到w_{m},然后继续使用Jacobi-Davidson算法,其中
V=[Vs^R_1,\ldots,Vs^R_{{m}_{\min}}]\quad \mathrm{and}\quad W=[Ws^L_1,\ldots,Ws^L_{{m}_{\min}}].
下一节:算法
上一级:收缩与重启
上一节:收缩
Susan Blackford
2000-11-20