The cross-product matrix-based CJ-FEAST SVDsolver proposed previously by the authors is shown to compute the left singular vector possibly much less accurately than the right singular vector and may be numerically backward unstable when a desired singular value is small. In this paper, an alternative augmented matrix-based CJ-FEAST SVDsolver is considered to compute the singular triplets of a large matrix $A$ with the singular values in an interval $[a,b]$ contained in the singular spectrum. The new CJ-FEAST SVDsolver is a subspace iteration applied to an approximate spectral projector of the augmented matrix $[0, A^T; A, 0]$ associated with the eigenvalues in $[a,b]$, and constructs approximate left and right singular subspaces with the desired singular values independently, onto which $A$ is projected to obtain the Ritz approximations to the desired singular triplets. Compact estimates are given for the accuracy of the approximate spectral projector, and a number of convergence results are established. The new solver is proved to be always numerically backward stable. A convergence comparison of the cross-product and augmented matrix-based CJ-FEAST SVDsolvers is made, and a general-purpose choice strategy between the two solvers is proposed for the robustness and overall efficiency. Numerical experiments confirm all the results.
翻译:作者先前提议的基于跨产品矩阵的 CJ-FTEST SVDSolver, 显示其计算左单向量的精确度可能远低于右单向量, 当想要的单值小时, 可能是数字向后的不稳定性。 在本文中, 另一种基于 CJ- FTEST SVDSolver 的增强矩阵基数的CJ- FTEST SVDSolver, 被视为独立计算一个大型矩阵的奇特三重值, 以单频谱中包含的单值为 $[a, b] 。 新的 CJ- FTET SVDsolver 是用于扩展矩阵光谱投影仪的近似频谱投影器 $0; A, 0, 0, 与 $[a, b] 相联的eigen值相关联的A, 且以所期望的单值为左右的单值为左右的单值, $A 。 用于光谱投影投影投影投影仪的精确度, 和一些趋同结果。 新的解结果结果被确定为新的解算为V 和S- fl- greal- decal- develeval- decal- decalentalental develvical degreal prisal pride sal pressal pressal develviolviolvical pressal pressal gal