Two optimal orthogonalization processes are devised toorthogonalize,possibly approximately,the columns of a very large and possiblysparse matrix A∈C^(n×k).Algorithmically the aim is,at each step,to optimallydecr...Two optimal orthogonalization processes are devised toorthogonalize,possibly approximately,the columns of a very large and possiblysparse matrix A∈C^(n×k).Algorithmically the aim is,at each step,to optimallydecrease nonorthogonality of all the columns of A.One process relies on using translated small rank corrections.Another is a polynomial orthogonalization process forperforming the Löwdin orthogonalization.The steps rely on using iterative methods combined,preferably,with preconditioning which can have a dramatic effect on how fast thenonorthogonality decreases.The speed of orthogonalization depends on howbunched the singular values of A are,modulo the number of steps taken.These methods put the steps of the Gram-Schmidt orthogonalizationprocess into perspective regardingtheir(lack of)optimality.The constructions are entirely operatortheoretic and can be extended to infinite dimensional Hilbert spaces.展开更多
基金supported by the Academy of Finland(Grant No.288641)。
文摘Two optimal orthogonalization processes are devised toorthogonalize,possibly approximately,the columns of a very large and possiblysparse matrix A∈C^(n×k).Algorithmically the aim is,at each step,to optimallydecrease nonorthogonality of all the columns of A.One process relies on using translated small rank corrections.Another is a polynomial orthogonalization process forperforming the Löwdin orthogonalization.The steps rely on using iterative methods combined,preferably,with preconditioning which can have a dramatic effect on how fast thenonorthogonality decreases.The speed of orthogonalization depends on howbunched the singular values of A are,modulo the number of steps taken.These methods put the steps of the Gram-Schmidt orthogonalizationprocess into perspective regardingtheir(lack of)optimality.The constructions are entirely operatortheoretic and can be extended to infinite dimensional Hilbert spaces.
