cur_rtol = np.nan
// Step A: QR decomposition of B page 1132 KNV
u, z = np.linalg.qr(B, mode="complete")
u0 = u[:, :B.shape[1]]
u1 = u[:, B.shape[1]:]
z = z[:B.shape[1], :]
After Change
//to debug with numpy qr uncomment the line below
//Q, _ = np.linalg.qr(pole_space_j, mode="complete")
Q, _ = s_qr(pole_space_j, mode="full")
ker_pole_j = Q[:, pole_space_j.shape[1]:]
// We want to select one vector in ker_pole_j to build the transfer
