This function applies V inverse via the Woodbury matrix identity
Vinv(G, qrM, tau_over_sigma, sigma, RHS)
| G | Matrix of genotype markers (size n x L) |
|---|---|
| qrM | Pre-computation for LxL linear system solve |
| tau_over_sigma | Tau over sigma from precomputation |
| sigma | Variance component from model noise epsilon |
| RHS | Matrix or vector on right-hand side of V inverse |
Matrix or vector resulting from left multiplication of Vinv with input RHS