[논문 리뷰] Answers to Some Questions about Explicit Sinkhorn Limits posed by Mel Nathanson

Shalosh B. Ekhad, Doron Zeilberger|arXiv (Cornell University)|2019. 02. 27.

Computability, Logic, AI Algorithms인용 수 4

한 줄 요약

이 논문은 양의 성분을 가진 3×3 대칭 행렬의 싱크본 극한에 대해 명시적인 대수적 공식을 제공하며, 그레브너 기저를 사용하여 대각 스케일링 행렬을 구한다. 이는 원래 행렬 성분에 대한 4차 방정식의 근으로 표현되는 이중 스케일링 극한 행렬과 그 스케일링 인자에 대한 닫힌 형식의 표현을 유도하며, 멜 나탄슨이 제기한 명시적 싱크본 극한과 알고리즘의 유한한 종료에 관한 열린 질문을 해결한다.

ABSTRACT

At the Jan. 2018 Joint Mathematics Meetings, Avi Wigderson gave a series of three fascinating lectures whose starting point was the Sinkhorn algorithm. One of the people in the audience was Mel Nathanson, and this lead him to ask some intriguing questions, some of which we answer in this note.

연구 동기 및 목표

• 멜 나탄슨이 요청한, 양의 성분을 가진 3×3 대칭 행렬의 싱크본 극한에 대해 명시적인 대수적 표현을 제공하는 것.
• 싱크본 알고리즘이 유한한 단계 수로 종료될 수 있는지에 대한 나탄슨의 질문을 해결하기 위해, 이러한 경우가 발생하는 명시적 예를 구성하는 것.
• 그레브너 기저 기법을 사용하여 싱크본 극한과 스케일링 행렬에 대한 정확한 대수적 표현을 도출하는 것.
• 기호 또는 수치적 매개변수를 가진 특수한 3×3 대칭 행렬의 경우에 대해 분석을 확장하여 그 싱크본 극한에 대한 닫힌 형식의 해를 제공하는 것.

제안 방법

• 저자들은 그레브너 기저를 위한 부처버거 알고리즘을 구현한 Maple 패키지 SINKHORN.txt를 사용하여 싱크본 극한을 정의하는 방정식의 체계를 기호적으로 해결한다.
• 일반적인 3×3 대칭 행렬에 대해, 변수 $ z $ 에 대한 4차 방정식을 유도하며, 이 방정식의 양의 근이 대각 행렬 $ X $ 의 스케일링 인자를 결정한다. 여기서 $ s_{11} = a_{11}z $ 이다.
• 싱크본 알고리즘의 교차적 행과 열 정규화 연산을 기호적으로 적용한 후, 이중 스케일링 조건을 강제로 적용하여 다항식 제약 조건을 도출한다.
• 특정 행렬에서 싱크본 알고리즘이 두 단계 내에 종료되는 명시적 예를 구성함으로써 해가 검증되며, 일부 행렬에 대해 유한 수렴이 확인된다.
• 매개변수화된 성분을 가진 특수한 경우(예: $ K, L, M $)에 대해, 싱크본 극한과 스케일링 행렬을 얻는 데 사용되는 축소된 4차 방정식을 도출한다.
• 이 방법은 더 큰 $ n \times n $ 대칭 행렬로 일반화될 수 있으며, 정의 방정식의 차수는 $ 2^{n-1} $ 으로 증가하므로 $ n \geq 4 $ 에서는 기호 계산이 불가능해진다.

실험 결과

연구 질문

• RQ1 양의 성분을 가진 3×3 대칭 행렬의 싱크본 극한은 행렬 성분에 대한 대수적 표현을 통해 닫힌 형식으로 표현될 수 있는가?
• RQ2 싱크본 알고리즘이 유한한 단계 수로 종료되는 행렬이 존재하는가? 만약 그렇다면, 이러한 행렬은 명시적으로 구성될 수 있는가?
• RQ3 그레브너 기저 기법을 사용하여 싱크본 극한과 스케일링 행렬에 대한 정확한 대수적 표현을 도출할 수 있는가?
• RQ4 싱크본 극한에서 스케일링 인자의 정의 다항식의 구조는 어떠한가? 그리고 그 차수는 행렬 크기와 어떻게 관련되는가?
• RQ5 행렬 $ A(r) $ 에 대해, 성분이 $ \frac{r(r+1)}{2}, 1, 1 $ 인 경우, 일반적인 $ r $ 에 대해 싱크본 알고리즘이 유한한 단계 수로 종료되는가?

주요 결과

• 일반적인 3×3 대칭 행렬 $ A $ 의 싱크본 극한 $ S $ 는 $ s_{11} = a_{11}z $ 로 명시적으로 주어지며, 여기서 $ z $ 는 $ A $ 의 성분에 대한 4차 방정식의 양의 근이며, 계수는 차수 8의 동차 다항식이다.
• 대각 스케일링 행렬 $ X $ 는 $ x_{11} = \sqrt{z} $ 를 만족하며, 나머지 $ X $ 와 $ S $ 의 성분은 공개 URL을 통해 제공되며, 논문에 포함하기에는 너무 길다.
• 행렬 $ A = \begin{pmatrix} K & 1 & 1 \\ 1 & L & 1 \\ 1 & 1 & 1 \end{pmatrix} $ 에 대해, 싱크본 극한은 $ K $ 와 $ L $ 에 따라 달라지는 계수를 가진 $ z $ 에 대한 4차 방정식에 의해 결정되며, 특정 $ K, L $ 의 선택에서 알고리즘이 두 단계 내에 종료된다.
• 행렬 $ A = \begin{pmatrix} K & 1 & 1 \\ 1 & L & 1 \\ 1 & 1 & M \end{pmatrix} $ 에 대해, 싱크본 극한은 다시 $ z $ 에 대한 4차 방정식으로 주어지며, 계수는 $ K, L, M $ 를 포함한다. 이 경우 해는 명시적으로 계산 가능하다.
• 논문은 싱크본 극한이 단 두 번의 반복만으로 도달하는 행렬 $ M = \begin{pmatrix} 2 & 2 & 6 \\ 2 & 1 & 2 \\ 9 & 3 & 3 \end{pmatrix} $ 를 구성하여, 이중 스케일링이 아닌 행렬에 대해서도 유한 종료가 가능함을 확인한다.
• 행렬 $ A(r) $ 에 대해 $ a_{11} = \frac{r(r+1)}{2} $ 이면, 싱크본 극한은 $ S(r) $ 로 명시적으로 계산되며, 기호적 반복에서 지속적인 비영인 행합 차이가 나타나므로 알고리즘이 $ r = 1 $ 이나 $ r = -2 $ 가 아닐 경우 유한한 단계 수로 종료되지 않는다.