대칭 Lanczos 사각법을 이용한 확률적 트레이스 추정

초록

본 논문은 대칭 Lanczos 사각법이 실제로 구현 가능한지에 대한 이론적·실험적 근거를 제시한다. 저자는 대칭 사각법 존재에 필요한 필요조건과 충분조건을 정리하고, Jordan‑Wielandt 형태의 행렬에 대해 특수한 초기 벡터 분포를 선택하면 충분조건을 만족함을 증명한다. 이를 바탕으로 이분 그래프와 방향성 그래프의 Estrada 지수를 계산하는 무편향 추정기를 설계하고, 기존 비대칭 Lanczos 기반 추정기에 비해 분산이 감소함을 실험적으로 확인한다.

상세 분석

논문은 먼저 Lanczos 과정에서 얻어지는 Jacobi 행렬의 고유값(즉, 사각법 노드)과 첫 번째 성분의 제곱이 가중치가 되는 전통적인 Gaussian 사각법을 재검토한다. 여기서 “대칭 사각법”이란 모든 단계 m 에 대해 노드가 0을 중심으로 대칭(θₖ = –θⱼ)하고, 대칭되는 노드 쌍이 동일한 가중치를 갖는 경우를 의미한다. 기존 연구에서는 대칭성이 계산 효율성을 높일 수 있다는 기대만 있었고, 실제로 언제 성립하는지는 명확히 밝혀지지 않았다.

저자는 먼저 **필요조건(Theorem 3.1)**을 제시한다. 대칭 사각법이 모든 가능한 Lanczos 단계에서 유지되려면, 원 행렬 A 가 스펙트럼이 0을 중심으로 대칭이어야 하고, 초기 벡터 v^(1) 의 좌표 μ = Qᵀv^(1) 가 “r‑부분 절대 팰린드롬” 형태를 가져야 한다. 즉, 비제로 고유값의 절반에 해당하는 좌표들의 절댓값이 서로 대칭적으로 일치한다는 의미이다. 이 조건은 정확 연산 환경에서만 성립한다는 점을 강조하며, 실제 수치 계산에서는 부동소수점 오차와 Lanczos 붕괴 현상이 추가적인 제약을 만든다.

다음으로 **충분조건(Theorem 3.3, 3.4)**을 제시한다. 핵심 아이디어는 Jordan‑Wielandt 형태의 블록 행렬
A =