ACCURATE COMPUTATION OF THE PRODUCT-INDUCED SINGULAR-VALUE DECOMPOSITION WITH APPLICATIONS

We present a new algorithm for floating-point computation of the singu lar value decomposition (SVD) of the product B-tau C, where B and C ar e full row rank matrices. The algorithm replaces the pair (B, C) with an equivalent pair (B', C') and then it uses the Jacobi SVD algorithm to compute the SVD of the explicitly computed matrix B'(tau) C'. In th is way, each nonzero singular value sigma is approximated with some si gma + delta sigma, where the relative error \delta sigma\/sigma is, up to a factor of the dimensions, of order epsilon{min(Delta is an eleme nt of D) kappa(2) (Delta B) + min(Delta is an element of D) kappa(2)(D elta C)}, where D denotes the set of diagonal nonsingular matrices, ka ppa(2) (.) denotes the spectral condition number, and epsilon is the r oundoff unit of floating-point arithmetic. The new algorithm is applie d to the eigenvalue problem HMx = lambda x with symmetric positive def inite H and M. It is shown that each eigenvalue lambda is computed wit h high relative accuracy and that the relative error \delta lambda\/la mbda of the computed approximation lambda + delta lambda is, up to a f actor of the dimension, of order epsilon{min(Delta is an element of D) kappa(2)(Delta H Delta) + min(Delta is an element of D) kappa(2)(Delt a M Delta)}. The new algorithm can also be used for accurate SVD compu tation of a single matrix G that admits an accurate factorization G = B-tau C.