I have a (non-symmetric) system that I would like to extract a kernel from, probably using QR factorisation or similar. What is the best way to do this at scale?