tdiff = np.dot(inv_hs, dxi / lambda_[np.newaxis, :])