denom = self.diag_y + self.eps_f * float(tau2)
Source: torchmatch/transport/samples/kernels/cg_dense.py:501