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