Os sistemas lineares de equações são difundidos nas estatísticas computacionais. Um sistema especial que encontrei (por exemplo, na análise fatorial) é o sistema
onde Aqui é uma matriz diagonal com uma diagonal estritamente positiva, é uma matriz semi-definida positiva simétrica (com ) e é uma matriz arbitrária . Somos solicitados a resolver um sistema linear diagonal (fácil) que foi perturbado por uma matriz de baixo escalão. A maneira ingênua de resolver o problema acima é inverter usando a fórmula de Woodbury . No entanto, isso não parece certo, pois as fatorações de Cholesky e QR geralmente podem acelerar drasticamente a solução de sistemas lineares (e equações normais). Eu vim recentemente no D n × n Ω m × m m ≪ n B n × m A