def METHOD_NAME(self, parents_params: np.ndarray, n: int) -> np.ndarray:
    # Compute orthogonal basis vectors for the subspace orthogonal to psl.
    e_12 = UNDXCrossover._normalized_x1_to_x2(
        parents_params
    )  # Normalized vector from x1 to x2.
    basis_matrix = np.identity(n)
    if np.count_nonzero(e_12) != 0:
        basis_matrix[0] = e_12
    basis_matrix_t = basis_matrix.T
    Q, _ = np.linalg.qr(basis_matrix_t)
    return Q.T[1:]