Efficient time propagation of U-D covariance factors

Time propagation of the Kalman filter covariance matrix involves an operation of the form ?P?T where for many applications ? is a sparse transition matrix. When the filter implementation employs U-D covariance factorization (i.e., recursions for U and D are used, where P = UDUT with U unit upper triangular and D diagonal) the corresponding time propagation involves W = ?U. Both the ?P?T and ?U computations can exploit transition matrix sparseness. If, however, the structure of W is not exploited, the computation involved with transforming W to an equivalent triangular form can be prohibitively expensive. The contribution of this paper is a streamlined Gram-Schmidt orthogonalization algorithm that can dramatically reduce UD time propagation computation costs.