We present a program for solving exactly the general pairing Hamiltonian based on diagonalization. The generates seniority-zero shell-model-like basis vectors via ‘01’ inversion algorithm. matrix is constructed in this space. evaluates all non-zero elements of “on fly” using scattering operator and search diagonalized by iterative Lanczos OpenMP parallel thus developed, PairDiag, can efficientl...