efficient k=0, pass y instead of Wy to Gaussian solver(Ryan)
cleaning a bit