B = torch.cat((B, y.new_zeros(Q.size(0))), 0)
M = torch.cat((M, math.sqrt(rho) * Q.t()), 0)
- return torch.lstsq(B, M).solution.view(-1)[:D+1]
+ return torch.lstsq(B, M).solution[:D+1, 0]
######################################################################