Hessenberg reduction for solving families of shifted linear systems (σI + τJ)x = b at O(n²) per shift - View it on GitHub