下一节:算法 上一级:收缩与重启 上一节:收缩

重启

假设广义舒尔形式(参见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