given that the next part is rather large, it would be useful to have a code comment on these lines that would help in understanding what rvir, rvrb, xtvix, ... are.
although I might guess     `resid @ inv(V) @ resid`, and so on
(new dot operator @ :)