def _KNV0(B, ker_pole, transfer_matrix, j, poles):
 transfer_matrix_not_j = np.delete(transfer_matrix, j, axis=1)
 (Q, R) = s_qr(transfer_matrix_not_j, mode='full')
 mat_ker_pj = np.dot(ker_pole[j], ker_pole[j].T)
 yj = np.dot(mat_ker_pj, Q[:, (-1)])
 if (not np.allclose(yj, 0)):
  xj = (yj / np.linalg.norm(yj))
  transfer_matrix[:, j] = xj

def _KNV0 ( B , ker_pole , transfer_matrix , j , poles ) : transfer_matrix_not_j = np . delete ( transfer_matrix , j , axis = 1 ) ( Q , R ) = s_qr ( transfer_matrix_not_j , mode = 'full' ) mat_ker_pj = np . dot ( ker_pole [ j ] , ker_pole [ j ] . T ) yj = np . dot ( mat_ker_pj , Q [ : , ( - 1 ) ] ) if ( not np . allclose ( yj , 0 ) ) : xj = ( yj / np . linalg . norm ( yj ) ) transfer_matrix [ : , j ] = xj
algorithm "knv0" kautsky et al .