{"id":"arxiv/2605.01202","url":"https://arxiv.org/abs/2605.01202","tags":["math.NA"],"score":0,"title":"Sampling Pfaffian point processes and the symplectic Arnoldi method","domain":"math.NA","source":"arxiv","authors":["Alan Edelman","Sungwoo Jeong","Simeon Schaub"],"reviews":0,"subject":"math.NA","summary":"We present an exact sampling algorithm for Pfaffian point processes based on a skew-symmetric analogue of the Cholesky factorization. This algorithm enables efficient sampling of a wide range of statistics arising in random matrix theory and combinatorics. For instance, we can sample eigenvalues of the orthogonal and symplectic ensembles ($β= 1,4$). In addition, we introduce a symplectic Arnoldi method for computing skew-orthogonal polynomials associated with a general weight function. This method can be used to efficiently construct the $2 \\times 2$ matrix valued skew-symmetric kernels that…","category":"research","comments":0,"disputed":false,"division":"math","mentions":0,"discussed":0,"published":"2026-05-02","reactions":{}}